本发明专利技术公开了一种基于成对HJ-1B影像时空信息的地表温度反演方法。该方法是针对现有的单窗（或单通道）算法在应用中所存在的参数获取困难而提出的一种新策略。实现该方法的步骤为：S1：获取研究区两个成对时相的HJ-1B影像，并进行相应的预处理；S2：对预处理后的影像以对称的方式进行边界延拓处理；S3：识别相似像元及其数目，确定相应的反演方法；S4：依据S3的判定结果，构建辐射传输方程组，并求解地表辐射值；S5：循环进行S3~S4,得的整个研究区在两个时相的地表辐射值，进行异常值识别和处理；S6：由处理后的地表辐射值反演地表温度。本发明专利技术能在相关参数无法准确获取的情况下准确反演地表温度，将有效促进HJ-1B热红外数据的定量应用。
【技术实现步骤摘要】
【专利摘要】本专利技术公开了一种基于成对HJ-1B影像时空信息的地表温度反演方法。该方法是针对现有的单窗（或单通道）算法在应用中所存在的参数获取困难而提出的一种新策略。实现该方法的步骤为：S1：获取研究区两个成对时相的HJ-1B影像，并进行相应的预处理；S2：对预处理后的影像以对称的方式进行边界延拓处理；S3：识别相似像元及其数目，确定相应的反演方法；S4：依据S3的判定结果，构建辐射传输方程组，并求解地表辐射值；S5：循环进行S3~S4,得的整个研究区在两个时相的地表辐射值，进行异常值识别和处理；S6：由处理后的地表辐射值反演地表温度。本专利技术能在相关参数无法准确获取的情况下准确反演地表温度，将有效促进HJ-1B热红外数据的定量应用。【专利说明】基于成对HJ-1B影像时空信息的地表温度反演方法
本专利技术涉及热红外遥感
，更具体的是在地表比辐射率和大气参数均未知 的情况下，利用两个成对时相HJ-IB影像的时空信息，从单一热红外波段的辐射亮度直接 定量反演地表温度的方法。
技术介绍
地表温度是影响生态系统的关键因子之一，它与许多地表过程有关。热红外卫星 遥感技术被认为是现阶段获取区域和全球尺度地表温度信息的一个可行且有效的手段。为 了实现从热红外卫星遥感影像准确提取地表温度信息，国内外研究人员提出了一系列的反 演方法，如单窗(或单通道)算法、劈窗算法、温度/比辐射率分离算法、白天/昼夜法。与其 它方法相比，单窗(或单通道）算法的优势在于它能从单一热红外波段影像反演地表温度。 因此，在实际应用中，对于Landsat TM/ETM+，CBERS，HJ-IB等仅具有单一热红外波段的传 感器，单窗(或单通道)算法是唯一的选择。由辐射传输方程可知，为了从单一热红外波段影 像准确地反演地表温度，需要同时对大气影响和地表比辐射率进行校正。已有的研究通过 大量的模拟实验和简化假设提出了相应的经验算法模型及其改进形式（Qin和Karnieli， 2001 ; Jimenez-Munoz 和 Sobrino, 2003;张勇等，2006;段四波等，2008 ;Jimenez_Munoz 等，2009;Cristobal 等，2009)。 然而,这些经验的算法均以地表比辐射率和大气参数(如大气平均温度、大气水分 含量）已知为前提。目前，对地表比辐射率的估算主要采用地表覆盖分类赋值法（Synder 等，1998)和植被指数法（Valor 和 Caselles，1996 ;覃志豪等，2004 ;Sobrino 等，2008)。 基于地表覆盖分类的方法需要预先知道不同地类的比辐射率，并认为每个类别的比辐射率 是固定不变的，这显然是不合理的；尽管基于植被指数的方法考虑了像元的组分差异，但依 然需要预先确定代表性组分的比辐射率，如土壤，植被等。另外，受准确大气参数获取条件 的限制，很多研究多假定整个研究区域的大气状态是均一分布的，在地表温度反演过程中 不考虑大气的影响，如Artis和Carnahan (1982)在早期所提出的一个单窗算法就没有考 虑大气的影响，因此利用该算法不能得到真实的地表温度。由于大气状态往往存在一定程 度的空间异质性，假定大气状况均质，不进行大气校正，或以某一点代替整个区域必将给遥 感地表温度反演带来误差（陈峰等，2010)。理论上，可以从MODIS的大气产品（如MOD 05) 提取现有经验算法所需的大气参数，但是该产品自身的不确定性可能会在有些情况下给最 终的地表温度结果带来不利影响。可见，由于在准确获取相关参数方面存在的困难，将在一 定程度上限制现有的单窗(或单通道）算法实施，从而也将影响HJ-IB单一热红外波段数据 的定量应用。 本专利技术针对HJ-IB星的高回访频率的特点，利用两个成对时相HJ-IB影像的时空 信息，基于两个合理的假设，提出一种直接从单一热红外波段影像反演地表温度的方法，而 无需过度依赖准确的地表比辐射率和大气参数信息。因此，该新方法可以视为对现有的单 窗(或单通道）算法的有益补充。
技术实现思路
要解决的技术问题本专利技术主要解决的技术问题是：在无法准确获取研究区地表 比辐射率和大气参数(如大气水分含量，大气平均作用温度）的情况下，即无法有效实施现 有的单窗(或单通道）算法时，实现从HJ-IB影像准确反演地表温度。 技术方案为了解决上述技术问题，本专利技术提出了一种基于成对HJ-IB影像时空 信息的地表温度反演方法，包括以下步骤： Sl获取研究区两个时相的HJ-IB影像(含C⑶多光谱和IRS热红外影像)，结合官方公 布资料和影像头文件信息，相应地，计算C⑶影像的各波段的星上反射率值和IRS热红外 影像辐射亮度值；并对HJ-IB C⑶反射率影像和热红外波段辐射影像进行地理校正等预处 理； S2对预处理后的影像以对称的方式进行边界延拓处理； S3依据每个目标像元所在的滑动窗口范围内各像元的光谱相似度指数，识别相似像 元及其数目，确定相应的反演方法； S4依据S3的结果，相应地提取相似或准相似像元的辐射值，构建一个简化的辐射传 输方程组，通过一个带约束的非线性优化模型实现对该方程组的求解，得到目标像元在两 个成对时相的地表辐射值； S5对延拓后的影像循环进行S3~S4，得的整个研究区两个时相的地表辐射值，并对其 进行异常值识别和处理； S6利用普朗克反方程或其简化模型从处理后的地表辐射值反演地表温度，并同时保 存反演方法的标识信息。 有益效果本专利技术可以在地表比辐射率和大气参数均无法准确获取的情况下，仅 利用两个成对时相HJ-IB影像的时空信息实现对地表温度的准确反演，克服了现有的方法 对辅助信息过于依赖的问题，进而能在一定程度上促进单一热红外波段影像的定量应用。 【专利附图】【附图说明】 图1是一种基于成对时相HJ-IB影像时空信息的地表温度反演方法示意图。 图2是基于本专利技术的两个时相HJ-IB地表温度反演结果图。 【具体实施方式】 以下结合附图和具体实施例，对本专利技术提供的基于成对HJ-IB影像时空信息的地 表温度反演方法进行详细描述。以下实施例仅用于说明本专利技术，而不是用来限制本专利技术的 范围。如图1所示，本专利技术基于成对HJ-IB影像时空信息的地表温度反演方法，主要包括以 下步骤： Sl获取研究区影像，并进行辐射校正、地理校正、重采样和裁剪预处理 在本实施例中，获取时相相近的两景HJ-IB IRS影像(获取时相分别为2011年12月 18日和2011年12月22 (path/row :452/91)，后文分别记为时相T18和T22)，并以珠江三 角洲的主要区域为实验区。由于数据质量的影响，仅有与T18相同时相的一景HJ-IB CCD2 多光谱影像（path/row: 454/90)是可用的(均来自http://www.cresda.cn)。因此，本实 施例中假定两个时相的地表状态相同，则可用T18时相的CCD影像充当T22时相的影像，而 不对最终的结果产生影响。一景Landsat5 TM影像（获取时相为2009年1月2日，path/ row: 122/44)，作为地理校正的参考影像。将HJ-IB CXD影像和IRS热红外本文档来自技高网...
【技术保护点】
一种基于成对HJ‑1B 影像时空信息的地表温度反演方法，其特征在于包括以下步骤： S1: 获取研究区两个成对时相的HJ‑1B影像（含CCD多光谱和IRS热红外影像），结合官方公布资料和影像头文件信息，相应地，计算CCD影像的各波段的星上反射率值和IRS热红外影像辐射亮度值；并对HJ‑1B CCD反射率影像和热红外波段辐射影像进行地理校正等预处理； S2: 对预处理后的影像以对称的方式进行边界延拓处理； S3: 依据每个目标像元所在的滑动窗口范围内各像元的光谱相似度指数，识别相似像元及其数目，确定相应的反演方法； S4: 依据S3的结果，相应地提取相似或准相似像元的辐射值，构建一个简化的辐射传输方程组，通过一个带约束的非线性优化模型实现对方程组的求解，得到目标像元在两个时相的地表辐射估计值； S5: 对延拓后的影像循环进行S3~S4, 得的整个研究区两个时相的地表辐射值，并对其进行异常值识别和处理； S6: 利用普朗克反方程或其简化模型从处理后的地表辐射值反演地表温度，并同时保存与反演方法相关的标识信息。
【技术特征摘要】
【专利技术属性】
技术研发人员：陈峰，赵小锋，全元，
申请(专利权)人：中国科学院城市环境研究所，
类型：发明
国别省市：福建;35
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。