地表温度反演IDL程序.docxVIP

下载本文档

14
0
约4.25千字
约 5页
2017-08-29 发布于重庆
举报
版权申诉

地表温度反演IDL程序.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

地表温度反演IDL程序

PRO LST?; 从文件夹中读取输入数据的文件名，数据顺序为2、4、6.....?imagenames=findfile(E:\rs\ld2\热岛\*, count=count)?outpath = E:\rs\ld2\200802\sz_lst\for k =3,count-1,2 do begin?print,imagenames[k]; 得到图像的行列数目及投影信息?envi_open_file,imagenames[k],R_FID=image_fid, /no_realize?if (image_fid eq -1) then return?envi_file_query, image_fid, ns=ns, nl=nl?map_info = envi_get_map_info(fid=image_fid)?dims = [-1, 0, ns-1, 0, nl-1];1,2,19波段的反射率?data_1 = envi_get_data(fid=image_fid,dims=dims,pos=0)?data_2 = envi_get_data(fid=image_fid,dims=dims,pos=1)?data_19 = envi_get_data(fid=image_fid,dims=dims,pos=20);31,32波段的亮温?data_31 = envi_get_data(fid=image_fid,dims=dims,pos=32)?data_32 = envi_get_data(fid=image_fid,dims=dims,pos=33)index_bad1= where((data_1 eq 65534) or (data_2 eq 65534) or (data_19 eq 65534) or $?????????????????? (data_31 eq 65534) or (data_32 eq 65534))envi_file_mng, id=image_fid, /remove;计算大气水含量?;w=((alfa-ln(ref19/ref2))/beta)2 2次幂?;alfa=0.02 beta=0.651?a=data_19*1.0/data_2?b=ALOG(a)?w=((0.02-b)/0.651)^2;计算大气透过率;传感器视角为10度的星下大气透过率，在水汽含量为0.4-2.0，2.0-4.0，4.0-6.0的计算方程?index_w1 = where((w ge 0.4) and (w lt 2.0))?index_w2 = where((w ge 2.0) and (w lt 4.0))?index_w3 = where((w ge 4.0) and (w le 6.0))?index_bad2 = where((w lt 0.4) or (w gt 6.0))t10_31 = fltarr(ns,nl)?t10_32 = fltarr(ns,nl)if index_w1[0] ne -1 then begin??? t10_31[index_w1] = 0.99513-0.08082*w[index_w1]??? t10_32[index_w1] = 0.99376-0.11369*w[index_w1]?endifif index_w2[0] ne -1 then begin??? t10_31[index_w2] = 1.08692-0.12759*w[index_w2]??? t10_32[index_w2] = 1.07900-0.15925*w[index_w2]?endifif index_w3[0] ne -1 then begin?t10_31[index_w3] = 1.07268-0.12571*w[index_w3]?t10_32[index_w3] = 0.93821-0.12613*w[index_w3]?endif;大气透过率温度校正函数?;因数据中亮温为实际亮温的10倍，阈值都乘以10倍，如318k在计算中用3180?index_31_t1 = where(data_31 gt 3180)?index_31_t2 = where((data_31 le 3180) and (data_31 ge 2780))?index_31_t3 = where(data_31 lt 2780)dt_31 = fltarr(ns,nl)if index_31_t1[0] ne -1 then begin?dt_31[index_31_t1] = 0.08?endifif index_31_t2[0] ne -1 then begin?;下面的公式中的0.