(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211105590.5 (22)申请日 2022.09.06 (71)申请人 中国科学院空天信息创新研究院 地址 100101 北京市朝阳区大屯路甲20号 北 (72)发明人 孟庆岩　高健峰　张琳琳　胡新礼　 (74)专利代理 机构 北京亿腾知识产权代理事务 所(普通合伙) 11309 专利代理师 陈霁 (51)Int.Cl. G01N 25/20(2006.01) G06Q 10/06(2012.01) G06Q 50/26(2012.01) (54)发明名称 一种用于工业热污染的景观斑块降温速率 提取方法 (57)摘要 本申请提供了一种用于工业热污染的景观 斑块降温速率提取方法， 包括： 确定工业区域、 工 业区域边界、 周围区域的景观斑块； 基于所述工 业区域边界， 确定缓冲区， 从所述周围区域的景 观斑块中， 提取所述缓冲区内的景观斑块； 基于 已有工具或者数据， 提取所述景观斑块的地表温 度； 计算所述景观斑块的相对斑块宽度； 基于景 观斑块的所述地表温度和所述相对斑块宽度， 计 算所述景观斑块的斑块降温速率。 本申请提供的 方法， 研究了景观斑块与工业区域之间的位置关 系对地表温度变化率的影 响， 分析了工业区域周 围景观环 境对工业热污染的降温情况， 并构建了 定量化指标评估降温速率。 权利要求书2页 说明书5页 附图3页 CN 115308261 A 2022.11.08 CN 115308261 A 1.一种用于 工业热污染的景观斑块降温速率 提取方法， 其特 征在于， 包括： 确定工业区域、 工业区域 边界、 周围区域的景观斑块； 基于所述工业区域边界， 确定缓冲区， 从所述周围区域的景观斑块中， 提取所述缓冲区 内的景观斑块； 基于已有的遥感图像处理平台， 提取所述工业区域和缓冲区的遥感图像， 进行地表温 度反演， 得到所述工业区域和缓冲区的地表温度， 提取 所述景观斑块的地表温度； 基于所述工业区域、 所述景观斑块， 计算所述景观斑块的相对斑块宽度,其中， 所述相 对斑块宽度为所述工厂区域的中心与所述景观斑块的中心的连线及其延长线切割所述景 观斑块的线段的长度； 基于景观斑块的所述地表温度和所述相对斑块宽度， 计算所述景观斑块的斑块降温速 率。 2.根据权利要求1所述的提取方法， 其特征在于， 所述确定工业区域、 工业区域边界、 周 围区域的景观斑块， 包括： 提取工业区域 边界和全球土地覆盖数据； 对所述工业区域边界和全球土地覆盖数据进行相交处理， 得到工业区域、 工业区域边 界及周围区域的景观斑块。 3.根据权利要求2所述的提取方法， 其特征在于， 所述提取工业 区域边界和全球土地覆 盖数据， 包括： 利用Google  Earth中“添加多边形 ”工具， 通过目视解译提取工业区域边界； 通过 Google earth engine(GEE)下载ESA 2020年全球土地覆盖数据； 所述对所述工业区域边界和全球土地覆盖数据进行相交处理， 得到工业区域、 工业区 域边界及周围区域的景观斑块， 包括： 基于ArcGIS的 “相交”工具， 对所述工业区域边界和全球土地覆盖数据进行相交处理， 获取工业区域及周围区域的景观斑块。 4.根据权利要求1所述的提取方法， 其特征在于， 所述基于所述工业区域边界， 确定缓 冲区， 从所述周围区域的景观斑块中， 提取 所述缓冲区内的景观斑块， 包括： 基于所述工业区域边界， 确定缓冲区， 将所述缓冲区分割成若干个子缓冲区， 从所述周 围区域的景观斑块中， 提取 所述子缓冲区内的景观斑块。 5.根据权利要求4所述的提取方法， 其特征在于， 所述基于所述工业区域边界， 确定缓 冲区， 将所述缓冲区分割成若干个子缓冲区， 从所述周围区域的景观斑块中， 提取所述子缓 冲区内的景观斑块， 包括： 基于所述工业区域边界， 利用ArcGIS中的 “缓冲区工具 ”， 确定3千米缓冲区， 通过6等 分， 将所述缓冲区分割成6个500米缓冲区， 从所述周围区域的景观斑块中， 提取所述500米 缓冲区内的景观斑块。 6.根据权利要求1所述的提取方法， 其特征在于， 所述基于已有的遥感图像处理平台， 提取所述工业区域和缓冲区的遥感图像， 进行地表温度反演， 得到所述工业区域和缓冲区 的地表温度， 提取 所述景观斑块的地表温度， 包括： 基于已有的遥感图像处理平台， 提取所述工业区域和缓冲区的遥感图像， 进行预处理 和地表温度反演， 得到不同季节的所述工业区域和缓冲区的地表温度， 提取所述景观斑块权　利　要　求　书 1/2 页 2 CN 115308261 A 2的地表温度。 7.根据权利要求6所述的提取方法， 其特征在于， 所述基于已有的遥感图像处理平台， 提取所述工业区域和缓冲区的遥感图像， 进行预处理和地表温度反演， 得到不同季节的所 述工业区域和缓冲区的地表温度， 提取 所述景观斑块的地表温度， 包括： 基于已有的遥感图像处理平台ENVI， 提取所述工业区域和缓冲 区的遥感图像， 使用辐 射校准工具和大气校正工具实现进行预处理， 使用波段运算工具实现地表温度反演， 采用 辐射传输方程方法反演得到不同季节的所述工业区域和缓冲区的地表温度， 提取所述景观 斑块的地表温度。 8.根据权利要求1所述的提取方法， 其特征在于， 所述基于已有的遥感图像处理平台， 提取所述工业区域和缓冲区的遥感图像， 进行地表温度反演， 得到所述工业区域和缓冲区 的地表温度， 包括： 基于已有的遥感图像处理平台， 提取所述工业区域和缓冲区的遥感图像， 进行地表温 度反演， 得到所述工业区域和缓冲区的地表温度； 或者基于已发布的地表温度数据， 提取 所述工业区域和缓冲区的地表温度。 9.根据权利要求1所述的提取方法， 其特征在于， 所述基于所述工业区域、 所述景观斑 块， 计算所述景观斑块的相对斑块宽度,其中， 所述相对斑块宽度为所述工厂区域的中心与 所述景观斑块的中心的连线及其延长线切割所述景观斑块的线段的长度， 包括： 基于所述工业区域， 确定所述工业区域中心； 基于所述景观斑块， 确定所述景观斑块中心； 所述工厂区域中心与所述景观斑块中心的连线及其延长线切割所述景观斑块的线段 的长度， 即为相对斑块宽度。 10.根据权利要求1所述的提取方法， 其特征在于， 所述基于景观斑块的所述地表温度 和所述相对斑块宽度， 计算所述景观斑块的斑块降温速率， 包括： 基于景观斑块的所述 地表温度， 提取 景观斑块的最高地表温度和最低地表温度； 计算所述景观斑块的最高地表温度和最低 地表温度的差值， 将所述差值与 所述相对斑 块宽度相比， 得到所述景观斑块的斑块降温速率。权　利　要　求　书 2/2 页 3 CN 115308261 A 3