叶绿素浓度反演
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
原理
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,co mpress=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.