跪求高手用IDL编程将modis产品介绍的LST8天数据合成30天的,即四合一

MODIS+1B数据处理及反演地表温度和土壤湿度IDL实现_图文_百度文库
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
评价文档:
MODIS+1B数据处理及反演地表温度和土壤湿度IDL实现
阅读已结束,如果下载本文需要使用
想免费下载本文?
你可能喜欢32MODIS+1B数据处理及反演地表温度和土壤湿度IDL实现-第8页
上亿文档资料,等你来发现
32MODIS+1B数据处理及反演地表温度和土壤湿度IDL实现-8
4MODIS影像检测;物由于其结构、组成及理化性质的差异,导致其对电磁;CHI+CH4>CH2+CH6;其中,CHl、CH2、CH4和CH6分别为MOD;NDWI=(CH4.CH6)/(CH4+CH6);式中CH代表该通道的反射率值;当云层较薄时,虽然云的干扰依然很严重,但也会有部;R=K幸CH2/CHl’;其中,K为放大倍数;MODIS1B数据处理及反演地
4MODIS影像检测物由于其结构、组成及理化性质的差异,导致其对电磁波的反射信息和地物本身存在的差异。云和雪在可见光波段的反射率远远高于其它地物,所以很容易将云、雪从其它地物中分辨出来;土壤、植被、水体在可见光和近红外波长范围的光谱特征有着较大的差异。植被和土壤在红光(0.6"--0.7um)的反射率较近红外(0.7~1.1pm)的反射率低;水体的反射主要在蓝绿光波段,其它波段吸收率很强,特别在近红外、中红外有很强的吸收带,反射率几乎为零,因此在遥感中常用近红外波段确定水体的位置和轮廓。对于MoDIS资料,丁莉东根据周成虎等人针对TM数据建立用于识别水体的谱间关系模型建立了基于MODIS资料的的水体光谱特征模型【24】:CHI+CH4>CH2+CH6其中,CHl、CH2、CH4和CH6分别为MODIS数据第1、2、4和6通道的反射率。根据水体的光谱特征,通过蓝光和近红外波段的组合,可以建立增大水与植被和土壤反差的归一化差分水指数(NDwI)采用第4和6通道的反射率计算归一化差分水指NDWI:?NDWI=(CH4.CH6)/(CH4+CH6)式中CH代表该通道的反射率值。可见光波段0.645lam,云雪的反射率一般大于O.3,因此当影像第1波段反射率<20%,可以隔离此波段较高反射率的云雾雪和沙漠等在。通过NDWI>=O,以及NDVI<=0,隔离水与植被和土壤等地物,对水体进行检测。4.3.2薄云雾下的水体的识别当云层较薄时,虽然云的干扰依然很严重,但也会有部分水体信息透过云层被传感器接收,而且水体和陆地的光谱特征相差较大,通过多光谱处理,增强这种差异,有可能识别薄云雾下的水体。由于NDWI是在利用第1通道消除云雪的干挠前提下对水体的识别,所以无法识别出薄云正气的水体。范伟提出利用MODIS波段2与1的比值增强模式来识别水体【25】,其公式为:R=K幸CH2/CHl’其中,K为放大倍数。该方法可以增强水体和陆地的光谱差异,但该模式还需要估测K和R两个阈值,对自动识别水体有一定的复杂。本文通过对影像进行云检测,确定有云的区域;然后对云区进行归一化植被指数NDVI光谱分析,发现当NDVI<=0时,能够很好的识别薄云下的水体,甚至是一些水面宽度较窄的河流等小水体混合像元。31MODIS1B数据处理及反演地表温度和土壤湿度IDL实现4.3.3消除阴影对水体的误判山体阴影和云阴影在影像中的特征与水体极为相此,在可见光、近红外波道的反射率都比较低,因此在进行水体识别时,常常会将阴影误判为水体。然而MODIS有20"25,27"36共16个热红外通道,主要反映地物自身的热辐射特性。其中通道27,28是太阳反射光和地物红外热辐射的交叉区。在白天,既接收来自地物的反射辐射,又接收来自地物的热辐射,因而此波段对地表温度特别敏感。通过对MODIS影像中20和31热红外通道中山体阴影与水体的亮温值比较,得到如下结论:11阴影区域比一般水体温度低。虽然水体在白天温度较低,反射率小,故呈现浅色调;但阴影部分由于热辐射受到山体、云层、高大建筑的阻挡,其亮温值要比水体小的多。2)MODIS的热红外通道20"-25,29~36,他们分别处在热红外的两个大气窗口区,适合于探测地物表面温度。水体在这些波段的发射率近似为常量,所以从理论上讲,水体在这些波段的差异应当很小。经过检测,同一水体像元在20,31波段下的亮温值不会超过5K,而阴影的差异相比水体而言就大的多。根据上述对与水体相关的光谱分析,得出水体识别流程图(见图4.4)。Chart4-4:Waterbodyidentificationflowchart4.4植被覆盖率4.4.1方法研究植被覆盖率作为反映地表信息的重要参数,一直是植被遥感领域的重要研究课题。324MODIS影像检测由于植被覆盖率具有显著的时空分异特性,因而,遥感已成为估算植被覆盖率的主要技术手段。植被覆盖率遥感研究的主要方法有回归模型法、混合光谱模型方法和植被指数方法。回归模型法主要基于对卫星同步实地观测数据的统计分析,利用遥感数据计算出的植被指数与实测植被覆盖率进行回归,建立线性回归或非线性回归的经验模型,由于回归模型区域性的限制,模型对所用遥感数据的空间分辨率要求较高,一般需要用TIVI或SPOT卫星资料。该模型只适用于特定的区域和特定的植被类型,并受观测时大气状况、土壤状况的影响,具有较大的局限性,但对小范围局部区域的植被覆盖率测量具有较高的精度【261。张良培【271对线性混合理论做了很好地论述,他认为,由不同物质形成的混合光谱可以表达成每种单物质(Endmember)光谱的线性组合,即:R=aRa+bRb+cRc+…式中,a、b、c等是加权系数;Ra、Rb、Rc等是单物质的反射比。在实际应用中,加权系数被认为是各种单成分物质在混合区域内(一个像元内)所占的面积,与其覆盖率有很大关系。其物理含义是,各个单成分物质的光谱是单独存在的,混合光谱是在单成分光谱进入传感器后形成的。基于这样的理论,遥感学家提出了许多用线性光谱模型计算植被覆盖率,这些模型都需要纯土壤和纯植被的反射率来计算覆盖率,而且方法要考虑的参数比较多,这些参数的获取又相对比较复杂,因此应用有一定的难度。在遥感应用领域,植被指数作为反映地表植被信息的最重要信息源,已广泛用来定性和定量评价植被覆盖及其生长活力。本文将重点阐述利用植被指数方法进行植被覆盖率检测。4.4.2植被指数植被生长状况的好坏可通过植被指数(Vl,Vegetationlndex)来反映,它是对地表植被活动的简单有效的经验度量,是对地表植被分布和变化情况的简单、有效和经验的衡量标准。利用遥感技术识别地物通常是基于其光谱差异的,将两个或多个光谱观测通道组合可得到植被指数,这一指数在一定程度上反映着植被的演化信息【14】。指数是通过地表覆盖物在可见波谱段的吸收和在近红外波谱的反射特性,建立的用于描述植被数量和质量的参数。许多研究表明红色通道和近红外通道反射能量与地表植被量有关。由于色素吸收在蓝色和红色波段最敏感。无植被或少植被区反差最小,中等植被区反差是红色和近红外波段变化结果,而高植被区则只有近红外波段对反差有贡献,红色波段趋于33MODIS1B数据处理及反演地表温度和土壤湿度IDL实现饱和。最早的植被指数是Jordan于1969年提出的比值植被指数(RvI),即RVI=knk/Lred(其中,),air表示近红外波段的反射率,kred表示红光波段的反射率)[2s】。RVI对大气的影响敏感,而且当植被覆盖不够浓密时(小于50%),它的分辨能力很弱,只有在植被覆盖浓密的情况下效果最好。Derring于1978年提出的归一化差分植被指数(NDVI,NormalizeDifferenceVI),将比值限定在【-1,1】范围内,可以消除大部分与太阳角、地形、云/暗影和大气条件有关的辐照度条件的变化。NDVI=(knir-钆red)|(9Lnir+Lred)虽然NDVI是植被生长状态及植被覆盖度提最常用的指示因子,但也存在不足:增强了NDVI低值部分,抑制了高值部分,导致NDVI数据容易饱和,对高植被密度区敏感性降低;另一方面NDVI容易受到植冠背景的影响,没有消除大气散射、吸心对光谱产生的影响。‘增强植被指数EVI,考虑到大气对光谱影响的波段相关性,用红、蓝波段的组合,代替NDVI中的红波段,减小大气影响的效果129瑚】:EVI=G*(Lnir-kred)/(Lnir+C1?Lred.C2*九blue+L)其中,G为调节因子,一般取G=2.5;C1,C2为抵抗大气调节系数,C1=6,C2=7.5,主要是通过Red与Blue波段的合理组合,调节大气气溶对植被指数的影响;L为土壤调节因子,用来调节土壤背景对植被指数的影响,取值一般为1.7-1。植被指数NDVI和EVI在MODIS数据中得到了广泛的应用,其MODIS.NDVI和MODIS.EVI表达式为:NDVI=(CH2.CHl)/(CH2+CHl)EVI=2.5?(CH2.CHl)/(CH2+6宰CHl.7.5"CH3+1)其中CHl、CH2、CH3分别为MODIS数据中第1、2、3波段的反射率。4.4.3植被指数合成在利用MODIS传感器进行观测时,太阳光照角度和观测角度以及云的条件的变化都非常大,因此得到的是来自地表的双向反射率信息。要构造大面积区域的晴空下的植被指数(vD,需要把给定时间段内的多幅VI图像合成一张图像,并且要使大气效应和角度效应的影响最小。合成期可根据研究需要进行按天合成:也即将每天多条轨道拼接4MODIS影像检测在一起,日合成产品也是MODIS植被指数产品系列中的第一个产品。也可以按旬或更长时间合成植被指数。MODIS植被指数的合成是在像元水平上进行的,采用最广泛使用的最大值合成(MVC),即通过云检测、质量检查等步骤后,逐像元地比较几张NDVI图像并选取最大的NDVI值合成后的NDVI值。MVC在近朗伯体表面的假定条件下效果很好,但是不适用于各向异性表面。归一化差分植被指数不能消除对表面各向异性,所选择的NDVI最大值与地表和大气的双向特性有关,前向散射产生较高NDVI值,后向散射则相反。因此当朗伯体为各向异性表面所取代时,MVC选择结果可信度降低。尽管MVC受到这些因素的限制,但是该方法对于未作大气订正的资料还是最好的方法,也是目前业务上使用的主要方法。4.4.4基于NDVI的植被覆盖率检测用植被指数方法估算植被覆盖率广泛研究于90年代初。DuncanJ,StowD[311等曾研究了墨西哥荒漠地区灌木林覆盖率与NDVI的关系,得到了较好的关系模型;LarssonH[32】分别从TM、MSS和SPOT卫星图像数据估算植被指数,并建立了阿拉伯森林地区植被指数与覆盖率的关系模型:中科院植物所的池宏康【331等通过分析沙地反射机理,建立了沙地油蒿群落盖度与修正后的土壤调节植被指数(MSAVI)之间的关系模型。众所周知,关系模型只适用于特定地区和特定的时间,因此应用起来有很大的局限性。与关系模型相比,下面中科院地理所张仁华提出的基于植被指数法计算植被覆盖率的公式模型I蚓应用更加方便也更加普遍为:pv=(NDVI―NDVIs)/(NDVIv?NDVIs)式中,NDVI是植被指数,NDVIv和NDVIs分别是茂密植被覆盖和完全裸土像元的NDVI值,通常取NDVIv=0.9,NDVIs=0.15。4.5对MODIS影像的检测实例所研究的云南区域MODIS影像,经过经纬度投影变换,切割范围为东经97031’36.65”---,106011’40.19”,北纬21008’33.87”"-'29013’26.15”。选择6、2、4波段进行合成,得到便于人工识别伪彩色遥感图像(图4.5)。根据色彩可以大体看出,图像右边及上、下角白色团状部分为云层,最右上角浅灰色团状为薄雾层,左上角蓝色树枝状35包含各类专业文献、外语学习资料、高等教育、生活休闲娱乐、各类资格考试、行业资料、应用写作文书、32MODIS+1B数据处理及反演地表温度和土壤湿度IDL实现等内容。 
 地表温度反演IDL程序_理学_高等教育_教育专区。PRO LST ; 从文件夹中读取输入数据...因本计算数据中亮温为实际值 的 10 倍,所以多乘一个 0.1 dt_31[index_...  4、3、2 波段组合影像 6 波段 数据处理:地表温度反演公式化简后的单窗体算法模型公式如下: Ts={a(1-C-D)+[b(1-C+D)+C+D]T6-DTa}/C 式中: C =...  接收的 MODIS HDF 格式原始影像数据,进行去条带、 几何校正和大气校正等预处理...3.2.1 地表温度反演――分裂窗算法 3.2.2 IDL 编程算法实现流程 20-21 3...  tp compile_opt idl2 envi,/restore_base_save_files envi_batch_init,log_...(title='打开红光波段影像') irimg_file='F:\HJ-1B 温度反演\...  数据进行操作反演, 实现襄樊市地区的地表真实温度的...真实温度 TS: TS = K2/ln(K1/ B(TS)+ 1) ...如城市热岛监测、计算土壤湿度指数(NDVI/T) 等。 ...  数据进行操作反演, 实现襄樊市地区的地表真实温度的...(TS)为: B(TS) = [Lλ - L↑- τ?(1-...如城市热岛监测、计算土壤湿度指数(NDVI/T) 等。 ...  地表温度反演单窗算法_建筑/土木_工程科技_专业资料。一:单窗算法 1、TM/ETM...(1+K1/L6) 其中:K1=607.76,K2=1260.56(覃志豪,用陆地卫星 TM6 数据演算...  利用envi反演地表温度_天文/地理_自然科学_专业资料。1. 遥感数据预处理 影像数据...2. 地表温度反演 Landsat TM5影像第6波段为热红外辐射波段,接收的热红外强度与...  的数据和算法介绍: 1、MODIS 数据是采用的 MODIS ...4、气溶胶反演的查找表是利用 IDL 调用 6S 辐射...处理 4.1 云检测 云检测工具是扩展工具,该工具实现...MODIS 1B数据处理及反演地表温度和土壤湿度IDL实现(可编辑)
扫扫二维码,随身浏览文档
手机或平板扫扫即可继续访问
MODIS 1B数据处理及反演地表温度和土壤湿度IDL实现(可编辑)
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由:
将文档分享至:
分享完整地址
文档地址:
粘贴到BBS或博客
flash地址:
支持嵌入FLASH地址的网站使用
html代码:
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布,请您等待!
3秒自动关闭窗口 上传我的文档
 下载
 收藏
该文档贡献者很忙,什么也没留下。
 下载此文档
正在努力加载中...
modis数据地表温度反演分裂窗算法的idl实现
下载积分:1000
内容提示:modis数据地表温度反演分裂窗算法的idl实现
文档格式:PDF|
浏览次数:2|
上传日期: 07:31:50|
文档星级:
该用户还上传了这些文档
modis数据地表温度反演分裂窗算法的idl实现
官方公共微信MODIS+1B数据处理及反演地表温度和土壤湿度IDL实现_图文_百度文库
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
评价文档:
MODIS+1B数据处理及反演地表温度和土壤湿度IDL实现
阅读已结束,如果下载本文需要使用
想免费下载本文?
你可能喜欢}

我要回帖

更多关于 modis数据下载 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信