In [1]:

    

%run ../common.ipynb


    

    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

image = imread('../sample.tif')[1]


    

In [7]:

    

F = fftshift(fft2(image))


    

In [8]:

    

imshow(log(abs(F)))


    

    
Out[8]:





(<matplotlib.figure.Figure at 0x11a63d7d0>,
 <matplotlib.axes._subplots.AxesSubplot at 0x1251c9550>,
 <matplotlib.image.AxesImage at 0x1251ef3d0>)






    

In [9]:

    

rF = rot90(F)


    

In [14]:

    

rF = ifftshift(rF)
rimg = ifft2(rF)


    

In [16]:

    

rimg.shape


    

    
Out[16]:





(3875, 3846)

In [19]:

    

%%time
fft2


    

    




CPU times: user 18.3 s, sys: 315 ms, total: 18.6 s
Wall time: 18.7 s

In [20]:

    

fft2??