叶绿素浓度反演

原理

1.SeaDAS大气校正公式

2.叶绿素浓度计算公式：

IDL编程实现：

proautumn

;443波段autumn

file=filepath('lw443.2002autumn.flat',root_dir='E:',subdir='课件\\定量遥感\\实验二\\20121212-10级学生上机')

openr,lun,file,/get_lun

data443au=fltarr(700,1100,1)

readu,lun,data443au

data443au1=1.2386*data443au+0.0008574

help,data443au1,/str

free_lun,lun

;490 autumn

file=filepath('lw490.2002autumn.flat',root_dir='E:',subdir='课件\\定量遥感\\实验二\\20121212-10级学生上机')

openr,lun,file,/get_lun

data490au=fltarr(700,1100,1)

readu,lun,data490au

data490au1=0.92887*data490au+0.0015606

help,data490au1,/str

free_lun,lun

;510 autumn

file=filepath('lw510.2002autumn.flat',root_dir='E:',subdir='课件\\定量遥感\\实验二\\20121212-10级学生上机')

openr,lun,file,/get_lun

data510au=fltarr(700,1100,1)

readu,lun,data510au

data510au1=1.0118*data510au+0.00039303

help,data510au1,/str

free_lun,lun

;555 autumn

file=filepath('lw555.2002autumn.flat',root_dir='E:',subdir='课件\\定量遥感\\实验二\\20121212-10级学生上机')

openr,lun,file,/get_lun

data555au=fltarr(700,1100,1)

readu,lun,data555au

data555au1=1.0674*data555au+1.2241e-5

help,data555au1,/str

free_lun,lun

r1=data443au1>data490au1>data510au1/data555au1

help,r1

r=alog10(r1)

chla_au=10^(-2.2402*r^4+1.4345*r^3+0.15474*r^2-0.90456*r+0.025477)

help,chla_au

;window,1,xsize=700/2,ysize=1100/2

;new=rebin(chla_au,700/5,1100/5)

;print,min(new),max(new)

;tvscl,new,order=1

;

;

;window,2,xsize=500,ysize=700;打开窗口大于图像大小

;map_set,/isotropic,/noborder,

/cylindrical,limit=[18,114,40,128],xmargin=[4,4],ymargin=[8,8]$

;,title='autumn East China Sea chlorophyll inversion!c';设置投影，上下左右留白

;images=map_image(new,startx,starty,lonmin=114,lonmax=128,latmin=18,latmax=40,/bilinear,compress=1)

;

;device,decomposed = 0

;loadct,4

;tvscl,images,order=1,startx,starty

;map_grid,latdel=2,londel=2,/box,/label;先画图后加网格

im = IMAGE(chla_au, RGB_TABLE=4, $

POSITION=[0.25,0.05,0.95,0.9], $

FONT_COLOR='Green', FONT_SIZE=16, $

TITLE='autumn East China Sea chlorophyll inversion',/order)

c = COLORBAR(TARGET=im, ORIENTATION=1, $

POSITION=[0.3,0.05,0.35,0.5], $

TITLE='Chla_au (g/ml)')

end2.