一种面向遥感影像云影检测的云及云阴影距离估计方法技术

面向遥感影像云影检测的云及云阴影距离估计方法，包括：逐列提取影像灰度值；逐行计算累积灰度值，并除以像素数目，得到灰度值累积率；提取曲线上各段的斜率；将连续出现的低速累积与快速累积作为一个云阴影和云的区域，并将低速累积起始点和快速累积起始点之间的距离作为云和云阴影的纵向距离；逐行对相距Δr的两行ri和ri+Δr，计算灰度累积率，并识别其中的快速累积与低速累积区域；将标记为云和云阴影的像素进行种子生长，并根据连通区域进行云和云阴影标记，得到云和云阴影对象；在获得纵向距离Δr和横向Δc的基础上，设定阈值模板阈值r'和c'，在Δr±r'和Δc±c'范围移动云阴影模板，将云和云阴影相关系数达到最小时的距离作为该云和云阴影之间的最佳的距离估计

【技术实现步骤摘要】

本专利技术涉及一种面向遥感影像云影检测的云及云阴影距离估计方法。具体实现上，根据云和云阴影形成的内在关系，结合云和云阴影灰度值较高和较低的特点，将其转化为灰度值累积率，初步估计云和云阴影之间的距离；然后利用种子生长提取云和云阴影，再根据云阴影和云之间的灰度值呈现负相关的特点，通过极小化相关系数，获得云和云阴影的距离优化估计。本方法完全基于影像特征而无需太阳高度角、传感器位置与姿态等参数，准确估计云和阴影之间的距离，可用于云及云阴影的复合类别提取及其相互验证，提高云影检测精度。
技术介绍
被动遥感获取的影像上，云和云阴影区域(在不致混淆的情况下，以下简称云影)仅包含少量的地表信息，在遥感影像上形成信息的空白区域。一方面，云影的存在使得可用的无云数据匮乏，数据的匮乏成为面向地表应用的重要限制性因素。另一方面，由于云和云阴影引起影像统计特征发生畸变，制约影像色彩归一化、影像镶嵌、时间序列分析、变化检测等诸多应用，通过云影检测获得云影掩膜，从而将云影区域剔除，能够显著提高各种处理与分析的精度。因而，云影检测是遥感影像预处理的重要研究内容。云影检测的一种直观而简单方法是利用云亮度值较高而阴影亮度值较低的特点进行阈值分割，阈值选取是需要解决的关键问题。动态阈值(吴传庆、王桥、杨志峰.基于混合像元分解的水体遥感图像去云法[J].遥感学报.2006,(10)2:176-183.)和多阈值方法(王伟,宋卫国,刘士兴等.KMeans聚类与多光谱阈值相结合的MODIS云检测算法[J].光谱学与光谱分析,2011,31(4):1061-1064.)等先进阈值选取方法用于提高阈值选取的精度。然而，该类方法的一个无法解决的内在缺陷：难以有效区分云和高反射率地物以及云阴影和低反射率地物，阈值优化是在误检率和漏检率之间达到平衡。为了解决云和高反射率地物以及云阴影和低反射地物不可区分性问题，一方面，云顶温度较低等辅助信息用于提高云影检测精度。典型算法如针对LandsatETM+的ACCA(AutomaticCloud-CoverAssessment)算法通过云影在各个波段的反射特征构建决策树，，实现云影检测，从而提高云影检测精度(IrishRR,BarkerJL,GowardSN,etal.CharacterizationoftheLandsat-7ETM+AutomatedCloud-CoverAssessment(ACCA)Algorithm[J].PhotogrammetricEngineering&RemoteSensing,2006,72(10):1179–1188.)。多时相影像的灰度值比较方法是通过云和云阴影造成相近时间获取影像上灰度值突变，用于辅助云影检测(HagolleO,HucM,PascualDV,etal.AMulti-TemporalMethodforCloudDetection,AppliedtoFORMOSAT-2,VENμS,LANDSATandSENTINEL-2Images[J].RemoteSensingofEnvironment,2010,114(8):1747–1755.)。该类方法的缺点在于需要额外的辅助信息，如多期影像或特殊的光谱通道(如热红外)，在辅助信息不足的情况下，方法的实用性降低。云和云阴影形成上的内在联系，使得云和云阴影之间存在一定的空间关系，通过云和云阴影的相互验证，能够提高云和云阴影的检测的精度。云和云阴影在影像上的空间关系决定于影像获取时刻的太阳高度角、传感器位置和姿态、云的高度。影像元数据中一般提供了太阳高度角、传感器的位置和姿态等参数，但是云高度确定是一个难于获得参数，因而无法通过物理模型直接获得云和阴影位置关系的参数。对此一般的解决方法是在影像分割获得云和云阴影的基础上，通过太阳高度角确定云和云之间的对应关系的面向对象方法(Zhu,Z.,&Woodcock,C.E.(2012).Object-basedCloudandCloudShadowDetectioninLandsatImagery.RemoteSensingofEnvironment,118,83–94.；WatmoughGr,AtkinsonPm,HuttonCW.ACombinedSpectralandObject-BasedApproachtoTransparentCloudRemovalinanOperationalSettingforLandsatETM+[J].InternationalJournalOfAppliedEarthObservationAndGeoinformation,2011,13(2):220–227.)。该类方法取得了较好的效果的同时，也存在一些不足，如云和云阴影提取的精度决定了距离估计的精度，给算法的应用带来不确定性。对此本专利技术提出一种基于影像的云及云阴影的距离估计方法。
技术实现思路
针对传统的遥感影像云和云阴影距离估计方法的不足，本专利技术提出一种基于影像的云及云阴影距离估计方法。本专利技术的原理为：在可见光与近红外通道，云对遥感过程的影响主要体现在(如图1)：(I)云对可见光以及近红外波段的高反射，一方面将地表信息的地球长波辐射反射回地表，使得到达传感器的包含地表信息辐射量减少；另一方面将太阳辐射大量反射并被传感器接收，相当于在影像中加入强度较大的白噪声，在影像上形成云。(II)同时，云对太阳辐射的大量反射，使得到达地表的辐射量减少，从而使得地表反射的长波辐射相应减少，在影像上形成云阴影。可以看出云和云阴影之间存在共生关系且灰度值之间表现为负相关。搭载光学传感器的卫星平台一般采用太阳同步轨道，以获取北半球上午影像(为了论述的方便，这里仅考虑北半球上午影像，其他情形可以根据太阳高度角进行初始估计)。由于阳光从东南方向照射过来，导致遥感影像上的云阴影一般出现在云的西北方向。同时，由于太阳和云的距离远远大于传感器与云的距离，且传感器一般是垂直对地观测，使得南北方向的距离远远大于东西方向的距离。因而，典型北半球的云和云阴影之间的距离如图2所示。根据以上原理，为了准确估计云和云阴影之间的距离，本专利技术的一种面向遥感影像云影检测的云及云阴影距离估计方法采用如下技术方案处理一景影像：步骤1：逐列提取影像灰度值，将其看成是随行变化的函数。由于云的灰度值较高，云阴影灰度值较低，因而云在曲线上形成“高原”，而云阴影形成“洼地”；步骤2：从第一行开始，逐行计算累积灰度值，并除以像素数目，该数值可以视为灰度值累积率，对应的曲线称为灰度值累积率曲线；步骤3：采用道格拉斯方法对上述曲线进行简化，用不同长度的线段进行表示，提取曲线上各段的斜率。云阴影部分灰度值较低，累积缓慢，导致相应线段的斜率较小。而云的亮度值较高，快速累积，导致曲线相应线段的斜率较大；步骤4：将连续出现的低速累积与快速累积作为一个云阴影和云的区域，并将低速累积起始点和快速累积起始点之间的距离作为云和云阴影的纵向距离Δr。将快速累积区域的像素标记为云，低速累积区域的像素标记为云阴影，将所有提取的云和云阴影纵向距离的均值作为整景影像纵向距离的估计；步骤5：逐行对相距Δr的两行ri和ri+Δr，按照步骤1和步骤2计算灰度本文档来自技高网...

【技术保护点】
一种面向遥感影像云影检测的云及云阴影距离估计方法，包括如下步骤：步骤1：灰度值提取；逐列提取各景影像的灰度值，第c列(c＝1,2,...,C)的灰度值可以表示为：Xc＝＜x1,c,x2,c,....,xR,c＞T；步骤2：灰度累积率曲线提取；第c列的灰度值累积率的计算方法为：Yc=<x1,cR,x1,c+x2,cR,...,x1,c+x2,c+...+xR,cR>T]]>其中，R为影像的行数，可以看出最后一项是该列灰度的算术平均值，该曲线被视为灰度值随行的累积速率，最后达到该列的均值；步骤3：灰度累积率曲线拟合；采用道格拉斯方法，对上述曲线进行简化，并将曲线表示为由各个线段组成的折线。假设一共拟合为n段折线，其中第i条线段表示为：L(i)＝(i,ki,bi,rbi,rei)其中ki和bi分别为第i曲线的斜率和截距，rbi和rei分别为该段曲线的起始行数和结束行数；步骤4：纵向距离估计与云阴影和云提取；将连续出现的低速累积线段和快速累积线段作为一个云阴影和云的候选区域，假设第i+1条线段检测为云，而i条线段检测为云阴影，将云和云阴影之间的纵向距离Δrk可以表示为：Δrk＝rbi+1‑rbirbi+1和rbi分别第i+1和i条线段的起始行数；在此基础上，将快速累积区域内的像素标记为云，低速累积区域内的像素标记为云阴影假设影像上共提取到K组对应云和云阴影，并纵向距离的均值该景影像上云和云阴影纵向距离的估计；Δr=1KΣk=1KΔrk.]]>步骤5：横向距离估计；逐行统计第r行和第r+Δr行的灰度累积率，将相距Δr上连续出现的低速累积和快速累积分别作为对应的云阴影和云，并按照步骤4估计云和云阴影的纵向距离估计方法估计横向距离，标记云和云阴影，并将云和云阴影的横向距离的均值作为整景影像云和云阴影距离的初步估计Δc。通过上述过程，可以获得遥感影像上云和云阴影距离的初步估计Δr和Δc；步骤6：云和云阴影种子生长；由于云和云阴影的边缘的数据可能未能有效提取，因而，将上述标记好的云和云阴影统计联通区域，并进行种子生长，获得云和云阴影；步骤7：云影距离优化；在获得纵向距离Δr和横向Δc的基础上，设定云阴影模板阈值r'和c'，在Δr±r'和Δc±c'范围内移动模板，并计算相关系数C，将负相关系数达到最小值的偏移作为最佳估计...

【技术特征摘要】
1.一种面向遥感影像云影检测的云及云阴影距离估计方法，包括如下步骤：步骤1：灰度值提取；逐列提取各景影像的灰度值，第c列(c＝1,2,...,C)的灰度值可以表示为：Xc＝＜x1,c,x2,c,....,xR,c＞T；步骤2：灰度累积率曲线提取；第c列的灰度值累积率的计算方法为：Yc=<x1,cR,x1,c+x2,cR,...,x1,c+x2,c+...+xR,cR>T]]>其中，R为影像的行数，可以看出最后一项是该列灰度的算术平均值，该曲线被视为灰度值随行的累积速率，最后达到该列的均值；步骤3：灰度累积率曲线拟合；采用道格拉斯方法，对上述曲线进行简化，并将曲线表示为由各个线段组成的折线。假设一共拟合为n段折线，其中第i条线段表示为：L(i)＝(i,ki,bi,rbi,rei)其中ki和bi分别为第i曲线的斜率和截距，rbi和rei分别为该段曲线的起始行数和结束行数；步骤4：纵向距离估计与云阴影和云提取；将连续出现的低速累积线段和快速累积线段作为一个云阴影和云的候选区域，假设第i+1条线段检测为云，而i条线段检测为云阴影，将云和云阴影之间的纵向距离Δrk可以表示为：Δrk＝rbi+1-rbirbi+1和rbi分别第i+1和i条线段的起始行数；在此基础上，将快速累积区域内的像素标记为云，低速累积区域内的像素标记为...

【专利技术属性】
技术研发人员：吴炜，夏列钢，沈瑛，杨海平，王卫红，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：浙江;33