一种基于时序遥感图像的水体提取方法技术

本发明专利技术公开一种基于时序遥感图像的水体提取方法，其步骤：首先收集覆盖研究区域的短期时间序列遥感图像，其次对收集的时间序列遥感图像进行几何校正、图像配准和辐射定标，然后采用趋势分析法分析近红外波段辐亮度随太阳高度角变化的趋势特征，得到近红外波段辐亮度随太阳高度角的变化斜率图，最后对近红外波段辐亮度随太阳高度角的变化斜率设定阈值进行水体提取，得到水体提取结果图。方法易行，操作简便，能够较好的将建设用地、山体阴影与水体区分开来，且对细小支流水体能够很好的检测出来，水体提取精度较高，同时本发明专利技术仅利用了近红外波段数据，减少了对多光谱遥感图像特定波段的依赖性。

【技术实现步骤摘要】
一种基于时序遥感图像的水体提取方法
本专利技术涉及遥感图像分类识别
，更具体地涉及一种利用短期时间序列遥感图像数据的遥感图像水体提取方法，适用于可以获取同一区域短期(地物类型未发生变化的时间段内，如数日)多个不同太阳高度角的多光谱遥感图像序列的情况下的地表水资源调查与监测。
技术介绍
地表水体分布信息的准确获取对水资源调查、河流综合治理、水利规划、洪涝旱灾监测以及灾害评估等领域具有重要意义。从卫星遥感影像中快速准确地提取水体信息已成为地表水资源调查与监测的一种重要手段。现阶段，利用遥感技术自动提取水体信息的方法主要有单波段阈值法、谱间关系法、监督分类法和水体指数法等。(1)单波段阈值法。该方法是利用水体在近红外或中红外波段处强吸收以及其它非水体在这两个波段处强反射的特点，来设定相应的阈值进行水体提取。这种方法原理较为简单，陆家驹等(陆家驹,李士鸿.TM资料水体识别技术的改进[J].环境遥感,1992,(01):17-23.)提出了基于红外单波段阈值识别水体的简单方法，但该方法难以精确提取细小水体。(2)监督分类法。该方法通过人工选取一定数量的代表性样本，根据已知类别的样本信息，利用统计学的方法对整幅影像进行分类，主要包括最大似然方法、支持向量机(SVM)方法等。段秋亚等(段秋亚,孟令奎,樊志伟,胡卫国,谢文君.GF-1卫星影像水体信息提取方法的适用性研究[J].国土资源遥感,2015,27(04):79-84.)通过对GF-1卫星影像水体提取方法的适用性研究，认为SVM法提取精度最高。该类方法对人工选取的样本具有很强的依赖性，且需要经过反复迭代，耗时较长。(3)谱间关系法。该方法通过对水体和背景地物的光谱曲线进行分析，确定水体有别于其它地物的光谱特征，进而进行水体信息的提取。该方法最早由周成虎等(周成虎，骆剑承，杨晓.1999，遥感影像地学理解与分析.北京：科学出版社)提出，他们提出的谱间关系水体提取模型为：(TM2+TM3)>(TM4+TM5)。张明华(张明华.用改进的谱间关系模型提取极高山地区水体信息[J].地理与地理信息科学,2008,24(2):14-16.)改进了谱间关系法，并和阈值法进行结合，构建了一种多条件谱间关系模型。(4)水体指数法。该方法利用水体及其它地物在不同波段的光谱差异的特征，通过构建特征指数从而达到水体提取的目的。常规的水体信息提取指数包括归一化差异水体指数(NormalizedDifferenceWaterIndex，NDWI)、修正归一化差异水体指数(ModifiedNormalizedDifferenceWaterIndex，MNDWI)等。NDWI利用水体特有的在绿波段的反射率高于近红外波段的特性，将这两个波段进行归一化差值处理，以此来突出水体信息，把水体信息与其背景地物更好的区分开来，该方法简单易操作，但容易将建设用地与水体混淆。徐涵秋(徐涵秋.利用改进的归一化差异水体指数(MNDWI)提取水体信息的研究[J].遥感学报,2005,9(5):589-595.)通过分析NDWI，将NDWI模型公式中的近红外波段修改为中红外波段，修改后的模型即为MNDWI，其水体提取效果优于NDWI，能更有效的提取城镇水体，并消除部分阴影对水体的影响等。虽然目前已经形成了一系列基于遥感的水体提取模型方法，但水体提取精度问题在某些背景地物反射率较低的地区仍然比较突出，如：建设用地、山体阴影等这些非水体地物与水体有相似的反射特征，目前的水体提取方法容易将其与水体混淆，从而导致水体的提取精度下降，MNDWI虽然能够一定程度上消除这些影响，但其利用的中红外波段目前多数光学卫星传感器并不具备。另外，当前水体提取使用的数据基本为中高分辨率的极轨卫星图像，如Landsat系列卫星影像，这些卫星的重访周期较长，且通常在一天中的固定时间拍摄，难以获取短期(地物类型未发生变化的时间段内，如数日)不同太阳高度角的图像，因此无法利用短期时间序列数据进行水体信息提取。高分四号(GF-4)卫星是我国第一颗地球同步轨道遥感卫星，具备高时间分辨率、中等空间分辨率观测能力，配置地球同步轨道五十米分辨率的凝视相机，可以完成分钟级别的遥感监测任务，能够获得同一区域短时内不同太阳高度角的图像，其时间维特征为基于遥感图像的地物识别和提取提供了新的思路。针对目前基于遥感图像的水体提取中存在的问题，结合GF-4卫星的优势，本专利技术提出了一种利用短期时序遥感图像数据的遥感图像水体提取方法，提高基于遥感图像的水体信息提取精度。
技术实现思路
为解决如何实现基于遥感图像的水体信息准确提取的技术问题，本专利技术的目的在于提供了一种利用短期时序遥感图像数据的遥感图像水体提取方法，方法易行，操作简便，能够较好地将建设用地、山体阴影与水体区分开来，且对细小支流水体能够很好地检测出来，水体提取精度较高。同时本专利技术仅利用了近红外波段数据，减少了对多光谱遥感图像特定波段的依赖性。为了实现所述的目的，本专利技术的技术解决方案是：本专利技术技术构思是：首选是收集覆盖研究区域的不同太阳高度角的短期时间序列多光谱遥感图像；其次是对收集的时间序列遥感图像进行几何校正、图像配准和辐射定标；再是基于上述校正后的时间序列影像，采用趋势分析方法分析近红外波段辐亮度随太阳高度角变化的趋势特征，得到近红外波段辐亮度随太阳高度角的变化斜率图；最后是基于上述的近红外波段辐亮度随太阳高度角的变化斜率图，通过设定阈值进行水体提取，得到水体分布提取结果图。一种利用短期时序遥感图像数据的遥感图像水体提取方法，其步骤是：(A)收集覆盖研究区域的不同太阳高度角的短期时间序列多光谱遥感图像。所述的研究区域时间序列遥感图像收集的要求包括以下步骤：1)每景遥感图像均覆盖研究区域，且每景图像的云覆盖区域均不大于10％；2)每景遥感图像均包含近红外波段，最好是同一传感器拍摄的图像；3)收集不同时相遥感图像n景(n>＝2)，不同时相图像上的地物类型未发生变化，一般数据的时间跨度不超过10天；4)遥感图像拍摄时的太阳高度角在30－75度，且最大太阳高度角和最小太阳高度角之间的差异不小于10度；(B)对收集的时间序列遥感图像进行几何校正、图像配准和辐射定标。所述的时间序列遥感图像几何校正、图像配准和辐射定标包括以下步骤：1)利用卫星图像自带的有理多项式系数(RationalPolynomialCoefficients，RPC)，通过有理函数模型(RationalFunctionModel，RFM)对图像进行几何校正。有理函数模型是美国SpaceImaging公司提供的一种广义的新型遥感卫星传感器成像模型，利用有理函数逼近二维像平面与三维物平面的对应关系。其形式如下式：其中(P,L,H)和(X,Y)分别是地面坐标和影像坐标的正则化坐标，其取值范围为-1到1之间。其中，Longitude，Latitude，Height分别表示WGS-84坐标下的经度、纬度和高程；LONG_OFF，LAT_OFF，HEIGHT_OFF分别表示经度、纬度和高程的偏移量；LONG_SCALE，LAT_SCALE，HEIGHT_SCALE分别表示经度、纬度和高程的比例因子。其中，COL，ROW分别表示图本文档来自技高网...

【技术保护点】
一种利用短期时序遥感图像数据的遥感图像水体提取方法，其步骤是：(A)收集覆盖研究区域的不同太阳高度角的短期时间序列多光谱遥感图像：1)每景遥感图像均覆盖研究区域，每景图像的云覆盖区域均不大于10％；2)每景遥感图像均包含近红外波段，是同一传感器拍摄的图像；3)收集不同时相遥感图像n景(n>＝2)，不同时相图像上的地物类型未发生变化，数据的时间跨度不超过10天；4)遥感图像拍摄时的太阳高度角在30－75度，最大太阳高度角和最小太阳高度角之间的差异不小于10度；(B)对收集的时间序列遥感图像进行几何校正、图像配准和辐射定标：1)利用卫星图像自带的有理多项式系数，通过有理函数模型实现图像几何校正；2)选定一个时相的遥感图像为基准图像，其它(n‑1)个时相的遥感图像作为待配准图像，通过选取基准图像与待配准图像上的同名点，然后利用这些同名点建立基准图像与待配准图像之间的几何畸变模型，再利用该几何畸变模型进行几何校正，实现不同时相遥感图像的图像配准；3)通过从卫星图像的头文件中读取定标系数，使用以下计算公式将各时相影像的数字量化值DN转换成辐亮度，实现图像的辐射定标：Lt＝Gain×DN+Bias式中，Lt为图像的辐亮度，Gain为增益，Bias为偏置，单位均为W﹒m...

【技术特征摘要】
1.一种利用短期时序遥感图像数据的遥感图像水体提取方法，其步骤是：(A)收集覆盖研究区域的不同太阳高度角的短期时间序列多光谱遥感图像：1)每景遥感图像均覆盖研究区域，每景图像的云覆盖区域均不大于10％；2)每景遥感图像均包含近红外波段，是同一传感器拍摄的图像；3)收集不同时相遥感图像n景(n>＝2)，不同时相图像上的地物类型未发生变化，数据的时间跨度不超过10天；4)遥感图像拍摄时的太阳高度角在30－75度，最大太阳高度角和最小太阳高度角之间的差异不小于10度；(B)对收集的时间序列遥感图像进行几何校正、图像配准和辐射定标：1)利用卫星图像自带的有理多项式系数，通过有理函数模型实现图像几何校正；2)选定一个时相的遥感图像为基准图像，其它(n-1)个时相的遥感图像作为待配准图像，通过选取基准图像与待配准图像上的同名点，然后利用这些同名点建立基准图像与待配准图像之间的几何畸变模型，再利用该几何畸变模型进行几何校正，实现不同时相遥感图像的图像配准；3)通过从卫星图像的头文件中读取定标系数，使用以下计算公式将各时相影像的数字量化值DN转换成辐亮度，实现图像的辐射定标：Lt＝Gain×DN+Bias式中，Lt为图像的辐亮度，Gain为增益，Bias为偏置，单位均为W﹒m-2﹒sr-1﹒μm-1，DN为卫星载荷观测值，无量纲；(C)基于上述校正后的时间序列影像，采用趋势分析方法...

【专利技术属性】
技术研发人员：任玉环，刘亚岚，
申请(专利权)人：中国科学院遥感与数字地球研究所，
类型：发明
国别省市：北京,11