|
Size: 1143
Comment:
|
← Revision 4 as of 2008-04-27 16:55:42 ⇥
Size: 990
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| {{attachement:ASAR_seaice_20080421_f32_5000x5000.dat}} | [[attachment:ASAR_seaice_mixed_20080421_f32_1000x1000.dat]] |
| Line 17: | Line 18: |
| y0,x0,wy,wx=2300,2800,1000,1000 img=reshape(fromfile('ASAR_seaice_20080421_f32_5000x5000.dat',dtype=float32),(5000,5000))[y0:y0+wy,x0:x0+wx] |
img=reshape(fromfile('ASAR_seaice_mixed_20080421_f32_1000x1000.dat',dtype=float32),(1000,1000)) |
| Line 24: | Line 24: |
| savefig('seaice_intensity.png',dpi=100) | |
| Line 36: | Line 35: |
| savefig('seaice_histogram.png',dpi=100) | |
| Line 42: | Line 40: |
| savefig('seaice_dB.png',dpi=100)}}} | }}} |
ASAR_seaice_mixed_20080421_f32_1000x1000.dat
1 from scipy import float32,int16,tan,sin,array,zeros,arange,pi,fromfile,reshape,ones,clip,log10
2 from pylab import rc,imshow,show,gray,figure,cm,hold,colorbar,title,savefig,subplot,hist,title,close,xlabel,ylabel
3 import scipy.ndimage as ndi
4 import scipy.misc as misc
5
6 rc('image',origin='lower',interpolation='nearest')
7
8 def ql(img,vmin,vmax):
9 figure()
10 imshow(img,vmin=vmin,vmax=vmax)
11 gray()
12 show()
13
14 img=reshape(fromfile('ASAR_seaice_mixed_20080421_f32_1000x1000.dat',dtype=float32),(1000,1000))
15
16 close('all')
17 ql(img,vmin=0,vmax=0.3)
18 colorbar()
19 title('Sea ice radar backscatter intensity')
20 dbimg=10*log10(img)
21
22 figure()
23 subplot(2,1,1)
24 hist(img,bins=500)
25 xlabel('Intensity')
26 ylabel('Occurrence')
27 subplot(2,1,2)
28 hist(dbimg,bins=500)
29 xlabel('Intensity [dB]')
30 ylabel('Occurrence')
31
32
33 ql(dbimg,-30,0)
34 colorbar()
35 title('Sea ice radar backscatter intensity [dB]')