from IPython.display import Image
Image('http://skeptic.smugmug.com/Nature/Costa-Rica-2011/i-xHQVnkv/0/L/IMG_1800-L.jpg')

!wget http://skeptic.smugmug.com/Nature/Costa-Rica-2011/i-xHQVnkv/0/L/IMG_1800-L.jpg

!convert IMG_1800-L.jpg sun:image.ras

!gmt2rgb image.ras -Lr -Gimage%c.grd -I1/1

! man grdfft

!grdfft imager.grd -Ew > fft

pspectrum=!cat fft

type(pspectrum)

%pylab inline

ax = pylab.axes()

ax.plot(pspectrum.fields(0),pspectrum.fields(1), 'm', lw=1)

ax.set_title('Radial Spectrum of IMG_1800-L.jpg')
ax.set_xlabel('Wavelength in pixels')
ax.set_ylabel('Power Spectral Density (pixel/Hz)')


img=matplotlib.image.imread('IMG_1800-L.jpg')
ax = pylab.axes()

plt.imshow(img)
ax.plot((200, 300),(400, 400), 'k', lw=2)
ax.set_xlim(0, 600)
ax.set_ylim(0, 600)

Image('http://skeptic.smugmug.com/Nature/Sea-Rim-State-Park-Texas-2003/i-Kg2WnjB/1/XL/DSC03045-XL.jpg')

!wget http://skeptic.smugmug.com/Nature/Sea-Rim-State-Park-Texas-2003/i-Kg2WnjB/1/XL/DSC03045-XL.jpg

!convert DSC03045-XL.jpg sun:image.ras
!gmt2rgb image.ras -Lr -Gimage%c.grd -I1/1
!grdfft imager.grd -Ew > fft
pspectrum=!cat fft

ax = pylab.axes()
ax.plot(pspectrum.fields(0),pspectrum.fields(1), 'm', lw=1)
ax.set_title('Radial Spectrum of DSC03045-XL.jpg')
ax.set_xlabel('Wavelength in pixels')
ax.set_ylabel('Power Spectral Density (pixel/Hz)')

ax = pylab.axes()
ax.semilogx(pspectrum.fields(0),pspectrum.fields(1), 'm', lw=1)
ax.set_title('Radial Spectrum of DSC03045-XL.jpg')
ax.set_xlabel('Wavelength in pixels')
ax.set_ylabel('Power Spectral Density (pixel/Hz)')

img=matplotlib.image.imread('DSC03045-XL.jpg')
ax = pylab.axes()

plt.imshow(img)
ax.plot((400, 400),(280, 340), 'k', lw=2)
ax.set_xlim(200, 700)
ax.set_ylim(200, 700)

Image('http://skeptic.smugmug.com/Nature/Sea-Rim-State-Park-Texas-2003/i-c5cC5gQ/1/XL/DSC03057-XL.jpg')

!wget http://skeptic.smugmug.com/Nature/Sea-Rim-State-Park-Texas-2003/i-c5cC5gQ/1/XL/DSC03057-XL.jpg

!convert DSC03057-XL.jpg sun:image.ras
!gmt2rgb image.ras -Lr -Gimage%c.grd -I1/1
!grdfft imager.grd -Ew > fft
pspectrum=!cat fft

ax = pylab.axes()
ax.semilogx(pspectrum.fields(0),pspectrum.fields(1), 'm', lw=1)
ax.set_title('Radial Spectrum of DSC03057-XL.jpg')
ax.set_xlabel('Wavelength in pixels')
ax.set_ylabel('Power Spectral Density (pixel/Hz)')

img=matplotlib.image.imread('DSC03057-XL.jpg')
ax = pylab.axes()

plt.imshow(img)
ax.plot((400, 350),(280, 380), 'k', lw=2) #100 pixel line
ax.plot((500, 470),(320, 380), 'b', lw=2) #60 pixel line
ax.plot((600, 580),(380, 420), 'r', lw=2) #30 pixel line

ax.set_xlim(200, 700)
ax.set_ylim(200, 700)