diff --git a/example/equirectangle/output.png b/example/equirectangle/output.png
deleted file mode 100644
index f73f813..0000000
Binary files a/example/equirectangle/output.png and /dev/null differ
diff --git a/lib/Perspec2Equirec.py b/lib/Perspec2Equirec.py
deleted file mode 100644
index c843cf1..0000000
--- a/lib/Perspec2Equirec.py
+++ /dev/null
@@ -1,72 +0,0 @@
-import os
-import sys
-import cv2
-import numpy as np
-
-class Perspective:
-    def __init__(self, img_name , FOV, THETA, PHI ):
-        self._img = cv2.imread(img_name, cv2.IMREAD_COLOR)
-        [self._height, self._width, _] = self._img.shape
-        self.wFOV = FOV
-        self.THETA = THETA
-        self.PHI = PHI
-        self.hFOV = float(self._height) / self._width * FOV
-
-        self.w_len = np.tan(np.radians(self.wFOV / 2.0))
-        self.h_len = np.tan(np.radians(self.hFOV / 2.0))
-
-    
-
-    def GetEquirec(self,height,width):
-        #
-        # THETA is left/right angle, PHI is up/down angle, both in degree
-        #
-
-        x,y = np.meshgrid(np.linspace(-180, 180,width),np.linspace(90,-90,height))
-        
-        x_map = np.cos(np.radians(x)) * np.cos(np.radians(y))
-        y_map = np.sin(np.radians(x)) * np.cos(np.radians(y))
-        z_map = np.sin(np.radians(y))
-
-        xyz = np.stack((x_map,y_map,z_map),axis=2)
-
-        y_axis = np.array([0.0, 1.0, 0.0], np.float32)
-        z_axis = np.array([0.0, 0.0, 1.0], np.float32)
-        [R1, _] = cv2.Rodrigues(z_axis * np.radians(self.THETA))
-        [R2, _] = cv2.Rodrigues(np.dot(R1, y_axis) * np.radians(-self.PHI))
-
-        R1 = np.linalg.inv(R1)
-        R2 = np.linalg.inv(R2)
-
-        xyz = xyz.reshape([height * width, 3]).T
-        xyz = np.dot(R2, xyz)
-        xyz = np.dot(R1, xyz).T
-
-        xyz = xyz.reshape([height , width, 3])
-        inverse_mask = np.where(xyz[:,:,0]>0,1,0)
-
-        xyz[:,:] = xyz[:,:]/np.repeat(xyz[:,:,0][:, :, np.newaxis], 3, axis=2)
-        
-        
-        lon_map = np.where((-self.w_len<xyz[:,:,1])&(xyz[:,:,1]<self.w_len)&(-self.h_len<xyz[:,:,2])
-                    &(xyz[:,:,2]<self.h_len),(xyz[:,:,1]+self.w_len)/2/self.w_len*self._width,0)
-        lat_map = np.where((-self.w_len<xyz[:,:,1])&(xyz[:,:,1]<self.w_len)&(-self.h_len<xyz[:,:,2])
-                    &(xyz[:,:,2]<self.h_len),(-xyz[:,:,2]+self.h_len)/2/self.h_len*self._height,0)
-        mask = np.where((-self.w_len<xyz[:,:,1])&(xyz[:,:,1]<self.w_len)&(-self.h_len<xyz[:,:,2])
-                    &(xyz[:,:,2]<self.h_len),1,0)
-
-        persp = cv2.remap(self._img, lon_map.astype(np.float32), lat_map.astype(np.float32), cv2.INTER_CUBIC, borderMode=cv2.BORDER_WRAP)
-        
-        mask = mask * inverse_mask
-        mask = np.repeat(mask[:, :, np.newaxis], 3, axis=2)
-        persp = persp * mask
-        
-        
-        return persp , mask
-        
-
-
-
-
-
-
diff --git a/lib/multi_Perspec2Equirec.py b/lib/multi_Perspec2Equirec.py
deleted file mode 100644
index 19effc9..0000000
--- a/lib/multi_Perspec2Equirec.py
+++ /dev/null
@@ -1,41 +0,0 @@
-import os
-import sys
-import cv2
-import numpy as np
-import lib.Perspec2Equirec as P2E
-
-
-class Perspective:
-    def __init__(self, img_array , F_T_P_array ):
-        
-        assert len(img_array)==len(F_T_P_array)
-        
-        self.img_array = img_array
-        self.F_T_P_array = F_T_P_array
-    
-
-    def GetEquirec(self,height,width):
-        #
-        # THETA is left/right angle, PHI is up/down angle, both in degree
-        #
-        merge_image = np.zeros((height,width,3))
-        merge_mask = np.zeros((height,width,3))
-
-        for img_dir,[F,T,P] in zip (self.img_array,self.F_T_P_array):
-            per = P2E.Perspective(img_dir,F,T,P)        # Load equirectangular image
-            img , mask = per.GetEquirec(height,width)   # Specify parameters(FOV, theta, phi, height, width)
-            merge_image += img
-            merge_mask +=mask
-
-        merge_mask = np.where(merge_mask==0,1,merge_mask)
-        merge_image = (np.divide(merge_image,merge_mask))
-
-        
-        return merge_image
-        
-
-
-
-
-
-
diff --git a/lib/old_Perspec2Equirec.py b/lib/old_Perspec2Equirec.py
deleted file mode 100644
index e1b024d..0000000
--- a/lib/old_Perspec2Equirec.py
+++ /dev/null
@@ -1,89 +0,0 @@
-import os
-import sys
-import cv2
-import numpy as np
-
-class Perspective:
-    def __init__(self, img_name , FOV, THETA, PHI ):
-        self._img = cv2.imread(img_name, cv2.IMREAD_COLOR)
-        [self._height, self._width, _] = self._img.shape
-        self.FOV = FOV
-        self.THETA = THETA
-        self.PHI = PHI
-    
-
-    def GetEquirec(self,height,width):
-        #
-        # THETA is left/right angle, PHI is up/down angle, both in degree
-        #
-
-        equ_h = height
-        equ_w = width
-        equ_cx = (equ_w - 1) / 2.0
-        equ_cy = (equ_h - 1) / 2.0
-
-        wFOV = self.FOV
-        hFOV = float(self._height) / self._width * wFOV
-
-        w_len = np.tan(np.radians(wFOV / 2.0))
-        h_len = np.tan(np.radians(hFOV / 2.0))
-
-
-        x_map = np.ones([self._height, self._width], np.float32)
-        y_map = np.tile(np.linspace(-w_len, w_len,self._width), [self._height,1])
-        z_map = -np.tile(np.linspace(-h_len, h_len,self._height), [self._width,1]).T
-
-        print(z_map[0])
-
-        D = np.sqrt(x_map**2 + y_map**2 + z_map**2)
-        xyz = np.stack((x_map,y_map,z_map),axis=2)/np.repeat(D[:, :, np.newaxis], 3, axis=2)
-        print(xyz[0,:,2])
-        
-        y_axis = np.array([0.0, 1.0, 0.0], np.float32)
-        z_axis = np.array([0.0, 0.0, 1.0], np.float32)
-        [R1, _] = cv2.Rodrigues(z_axis * np.radians(self.THETA))
-        [R2, _] = cv2.Rodrigues(np.dot(R1, y_axis) * np.radians(-self.PHI))
-
-        xyz = xyz.reshape([self._height * self._width, 3]).T
-        xyz = np.dot(R1, xyz)
-        xyz = np.dot(R2, xyz).T
-        lat = np.arcsin(xyz[:, 2])
-        lon = np.arctan2(xyz[:, 1] , xyz[:, 0])
-
-        lon = lon / np.pi * 180
-        lat = -lat / np.pi * 180
-
-        print(lat.reshape([self._height , self._width])[0])
-        print(lon.reshape([self._height , self._width])[0])
-        
-        lon = (lon / 180 * equ_cx + equ_cx).astype(np.int)
-        lat = (lat / 90  * equ_cy + equ_cy).astype(np.int)
-        coordinate = (lat,lon)
-
-        x_map = np.repeat(np.arange(self._height), self._width)
-        y_map = np.tile(np.arange(self._width), self._height)
-
-        blank_map_x = np.zeros((height,width))
-        blank_map_y = np.zeros((height,width))
-        mask = np.zeros((height,width,3))
-
-        blank_map_x[coordinate] = x_map
-        blank_map_y[coordinate] = y_map
-        mask[coordinate] = [1,1,1]
-
-        # print(lat.reshape([self._height, self._width]))
-        # print(lon.reshape([self._height, self._width])[-1,1910:1930])
-
-
-        persp = cv2.remap(self._img, blank_map_y.astype(np.float32), blank_map_x.astype(np.float32), cv2.INTER_CUBIC, borderMode=cv2.BORDER_WRAP)
-        
-        persp = persp * mask
-        
-        return persp , mask
-        
-
-
-
-
-
-
diff --git a/panorama2cube.py b/panorama2cube.py
deleted file mode 100644
index 96f86c3..0000000
--- a/panorama2cube.py
+++ /dev/null
@@ -1,99 +0,0 @@
-import os
-import cv2 
-import lib.Equirec2Perspec as E2P
-import lib.Perspec2Equirec as P2E
-import lib.multi_Perspec2Equirec as m_P2E
-import glob
-import argparse
-
-
-
-def panorama2cube(input_dir,output_dir):
-
-    cube_size = 640
-
-    if not os.path.exists(output_dir):
-        os.mkdir(output_dir)
-
-    all_image = sorted(glob.glob(input_dir + '/*.*'))
-
-    print(all_image)
-
-
-    for index in range(len(all_image)):
-        # image = '../Opensfm/source/library/test-1/frame{:d}.png'.format(i)
-        equ = E2P.Equirectangular(all_image[index])    # Load equirectangular image
-        #
-        # FOV unit is degree
-        # theta is z-axis angle(right direction is positive, left direction is negative)
-        # phi is y-axis angle(up direction positive, down direction negative)
-        # height and width is output image dimension
-        #
-
-        out_dir = output_dir + '/%02d/'%(index)
-        if not os.path.exists(out_dir):
-            os.mkdir(out_dir)
-
-        img = equ.GetPerspective(90, 0, 0, cube_size, cube_size)  # Specify parameters(FOV, theta, phi, height, width)
-        output1 = out_dir +  'front.png'
-        cv2.imwrite(output1, img)
-
-        img = equ.GetPerspective(90, 90, 0, cube_size, cube_size)  # Specify parameters(FOV, theta, phi, height, width)
-        output2 = out_dir + 'right.png' 
-        cv2.imwrite(output2, img)
-
-
-        img = equ.GetPerspective(90, 180, 0, cube_size, cube_size)  # Specify parameters(FOV, theta, phi, height, width)
-        output3 = out_dir + 'back.png' 
-        cv2.imwrite(output3, img)
-
-        img = equ.GetPerspective(90, 270, 0, cube_size, cube_size)  # Specify parameters(FOV, theta, phi, height, width)
-        output4 = out_dir + 'left.png' 
-        cv2.imwrite(output4, img)
-
-        img = equ.GetPerspective(90, 0, 90, cube_size, cube_size)  # Specify parameters(FOV, theta, phi, height, width)
-        output5 = out_dir + 'top.png' 
-        cv2.imwrite(output5, img)
-
-        img = equ.GetPerspective(90, 0, -90, cube_size, cube_size)  # Specify parameters(FOV, theta, phi, height, width)
-        output6 = out_dir + 'bottom.png' 
-        cv2.imwrite(output6, img)
-
-
-def cube2panorama(input_dir,output_dir):
-
-    width = 1920
-    height = 960
-
-    if not os.path.exists(output_dir):
-        os.mkdir(output_dir)
-    
-    front = input_dir + '/front.png'
-    right = input_dir + '/right.png'
-    back = input_dir + '/back.png'
-    left = input_dir + '/left.png'
-    top = input_dir + '/top.png'
-    bottom = input_dir + '/bottom.png'
-
-    # this can turn cube to panorama
-    per = m_P2E.Perspective([front,right,back,left,top,bottom],
-                            [[90, 0, 0],[90, 90, 0],[90, 180, 0],
-                            [90, 270, 0],[90, 0, 90],[90, 0, -90]])    
-    
-    
-    img = per.GetEquirec(height,width)  
-    cv2.imwrite(output_dir + '/output.png', img)
-
-if __name__ == '__main__':
-    parser = argparse.ArgumentParser()
-    
-    parser.add_argument('--mode', type=str, default='panorama', choices=['panorama', 'cube'])
-    parser.add_argument('--input', type=str, default='./panorama')
-    parser.add_argument('--output', type=str, default='./output')
-
-    config = parser.parse_args()
-
-    if config.mode == 'panorama':
-        panorama2cube(config.input,config.output)
-    else:
-        cube2panorama(config.input,config.output)
diff --git a/pers2equir.py b/pers2equir.py
deleted file mode 100644
index 973da52..0000000
--- a/pers2equir.py
+++ /dev/null
@@ -1,43 +0,0 @@
-import os
-import cv2 
-import lib.Equirec2Perspec as E2P
-import lib.Perspec2Equirec as P2E
-import lib.multi_Perspec2Equirec as m_P2E
-import glob
-import argparse
-
-
-
-def pers2equir():
-    #
-    # FOV unit is degree
-    # theta is z-axis angle(right direction is positive, left direction is negative)
-    # phi is y-axis angle(up direction positive, down direction negative)
-    # height and width is output image dimension
-    #
-    
-    input_dir = './example/perspective'
-    output_dir = './example/equirectangle'
-
-    width = 1920
-    height = 960
-
-    if not os.path.exists(output_dir):
-        os.mkdir(output_dir)
-    
-    input1 = input_dir + '/perspective_1.png'
-    input2 = input_dir + '/perspective_2.png'
-
-    # this can turn cube to panorama
-    equ = m_P2E.Perspective([input1,input2],
-                            [[120, 0, 0],[120, 0, 90]])    
-    
-    
-    img = equ.GetEquirec(height,width)  
-    print(img.shape)
-    cv2.imwrite(output_dir + '/output.png', img)
-
-
-
-if __name__ == '__main__':
-    pers2equir()
\ No newline at end of file
diff --git a/spherical/circle-cropped.jpg b/spherical/circle-cropped.jpg
new file mode 100644
index 0000000..2202d86
Binary files /dev/null and b/spherical/circle-cropped.jpg differ
diff --git a/unwrapping b/unwrapping
new file mode 100644
index 0000000..e8ddab9
--- /dev/null
+++ b/unwrapping
@@ -0,0 +1,73 @@
+import cv2
+import argparse
+import numpy as np
+import math
+
+img = cv2.imread("./spherical/image.jpg",1)
+height,width,planes = (np.shape(img))
+ratio = height/width
+dim = (int(720/ratio),720)
+img_resize = cv2.resize(img, dim, interpolation = cv2.INTER_AREA) 
+height,width,plane = np.shape(img_resize)
+width_dst = height*2
+dst_img = np.zeros((height, width_dst,3),dtype= int)
+
+print(height)
+
+for k in range(0,3):
+    for i in range(1,height):
+    
+        altitude = abs(i-int(height/2))
+                
+        lat = np.sqrt(abs((height*height)/4 - (altitude*altitude)))
+        lat = int(lat)
+    
+        if(lat==0):
+            continue   
+    
+        left = (width/2) - lat
+        left = int(left)
+                    
+        right = (width/2) + lat
+        right = int(right)
+        
+        if(width%2==0):
+            left=left-1
+            right=right-1
+        
+        jump = width_dst/(2*lat)
+        if(jump<1):
+            jump=1
+        #jump = int(jump)
+        
+           
+        for j in range(left,right):
+            position = int((j-left)*jump)
+            #print("position is ",position,"j is ",j)
+            if(position>=width_dst or j>=width):
+                break;
+                
+            dst_img[i,position,k] = img_resize[i,j,k]
+          
+           
+    #print(dst_img[0:123,0:149])            
+    dst_img = dst_img.astype(np.uint8)
+    #print(dst_img[300,:])
+
+for c in range(0,3):    
+    for i in range(1,height):
+        pixel = 100;
+        for j in range(width_dst):
+            
+            if(dst_img[i,j,c]!=0):
+                pixel = dst_img[i,j,c]
+               
+        #print(pixel)
+            if(dst_img[i,j,c]==0):
+                dst_img[i,j,c] = pixel
+
+           
+
+
+
+cv2.imwrite('./panorama/input.jpg',dst_img)