[[attachment:ASAR_seaice_mixed_20080421_f32_1000x1000.dat]]


{{{#!python
from scipy import float32,int16,tan,sin,array,zeros,arange,pi,fromfile,reshape,ones,clip,log10
from pylab import rc,imshow,show,gray,figure,cm,hold,colorbar,title,savefig,subplot,hist,title,close,xlabel,ylabel
import scipy.ndimage as ndi
import scipy.misc as misc

rc('image',origin='lower',interpolation='nearest')

def ql(img,vmin,vmax):
    figure()
    imshow(img,vmin=vmin,vmax=vmax)
    gray()
    show()
    
img=reshape(fromfile('ASAR_seaice_mixed_20080421_f32_1000x1000.dat',dtype=float32),(1000,1000))

close('all')
ql(img,vmin=0,vmax=0.3)
colorbar()
title('Sea ice radar backscatter intensity')
dbimg=10*log10(img)

figure()
subplot(2,1,1)
hist(img,bins=500)
xlabel('Intensity')
ylabel('Occurrence')
subplot(2,1,2)
hist(dbimg,bins=500)
xlabel('Intensity [dB]')
ylabel('Occurrence')


ql(dbimg,-30,0)
colorbar()
title('Sea ice radar backscatter intensity [dB]')
}}}