利用多期遥感影像提取砂质岸线的方法
【技术领域】
[0001] 本发明涉及一种卫星遥感海岸线探测方法，属于空间遥感技术领域，尤其涉及利 用多期遥感影像提取砂质岸线的方法。
【背景技术】
[0002] 海岸线附近蕴藏着丰富的生物、矿产、能源、土地等自然资源，它不仅是国防的前 口肖，又是海、陆交通的连接地，是人类经济活动频繁的地带。测定海岸线位置是人们研究海 陆相互作用、海洋地质研究、气候变化和环境保护等必须进行的一项技术活动，也是进行海 洋地形图测绘、海岸带调查等工作的重要内容。因此，对海岸线的快速准确探测有着重要的 价值。但因为自然和人为的作用，实际的海岸线是时刻处在连续、动态的变化过程。现有的 海岸线提取方法主要包括传统的现场测量踏勘和基于遥感与地理信息系统（GIS)技术的 综合，传统的现场测量踏勘方法提取海岸线，虽保证了岸线的精度，但费时、费力、成本高、 易受到气候天气影响，而且难以实现大尺度、长时间序列的监测，基于GIS技术的海岸线提 取因具有成本低、范围广和获取数据准确及时的优势，尤其是能够解决复杂地形下的海岸 线测量而被广泛应用。
[0003] 遥感影像上的水体是相对较特殊的图像部分，主要原因是水体的反射率随电磁波 波长的增大而递减，至近红外波段的清澈水体时反射率降至零；另外水体浑浊度的增加会 使水体的反射率有所上升，可能到短波红外波段，水体的反射率才降至零，正是根据这一特 性，可以实现水陆分割，从而提取海岸线。基于此状况，国内外专家学者作了大量的相关工 作：Bellomo、朱长明、Bouchahma、刘善伟、AmrYousef、马晓峰等基于遥感影像，结合NDWI、 MNDWI、SVM、DEM、边缘检测算子、图像纹理及数学形态学等方法自动提取了海岸线；同时 White、翁宇斌、孙伟富、Mujabar等专家学者通过人机交互方法提取了海岸线，这些方法均 是利用单景影像提取海岸线，虽可较迅速地提取出海岸线，但提取的岸线一般都是卫星或 航空飞行器过境时的瞬时水边线，未考虑同一年份不同时间内潮汐对瞬时水边线位置的影 响，并不是严格意义的海岸线，部分学者虽考虑到潮汐影响，但都是基于单景单期影像进行 潮汐校正。
[0004] 中国专利（授权公告号CN102013015B)公开了 "一种面向对象的遥感影像海岸线 提取方法"，属于全自动遥感影像海岸线提取领域。该发明主要是利用利用种子生长的方法 提先取出海水区域，最后利用相应的判别准则提取海岸线，但同样采用单景影像提取海岸 线，因而提取的海岸线结构不够精确。
[0005] 中国专利（申请公布号CN103530635A)公开了"一种基于卫星微波遥感图像 的海岸线提取方法"，该基于卫星微波遥感图像的海岸线提取方法从卫星微波遥感图像上 提取海岸线，包括预处理、图像分割和后处理三大步骤，该发明也是采用单景影像提取海岸 线，因而提取的海岸线结构不够精确，同时在图像分割及后处理过程中，环节较多，提取海 岸线后需要再对照原始图像手动修改，降低了海岸线的提取效率。
【发明内容】
[0006] 本发明提供一种利用多期遥感影像提取砂质岸线的方法，用于解决现有技术中单 景影像提取的海岸线不准确的问题，采用多期遥感影像作为海岸线的提取基础，获得瞬时 水边线再进行位置校正，经多期校正水边线的融合从而获取较为准确的砂质岸线。
[0007] 利用多期遥感影像提取砂质岸线的方法，包括以下步骤：
[0008] 步骤一：对通过卫星获得的遥感影像进行预处理；
[0009] 所述预处理包括采用图像回归法进行影像相对辐射归一化处理，再利用几何校正 方法对遥感影像进行配准，得到预处理后的遥感影像，相对辐射归一化公式：
[0010] yk=akxk+bk (1)
[0011] m
[0012] ⑶
[0013] 其中，波段的归一化参数，xk、yk*别为校正图像k波段归一化前后的像 元灰度值，分别为校正影像和参考影像k波段的像元平均值，分别为校正 影像与参考影像间k波段的协方差及参考影像本身的协方差；
[0014] 步骤二：获取NDWI图像；
[0015] 将预处理后的遥感影像进行NDWI指数计算，以获得NDWI图像，NDWI指数公式：
[0016](4) ..广
' +广.JN-Jtt
[0017] 其中，PG、PNIR分别表示绿光及近红外波段的反射率；
[0018] 步骤三：获取水陆分割的二值化影像；
[0019] 根据NDWI图像的直方图分布特征确定图像分割阈值，并对NDWI图像进行水陆分 割以获的二值化影像；
[0020] 步骤四：获取瞬时水边线；
[0021] 通过栅格矢量转换方法将二值化影像转为矢量数据，再通过碎斑剔除、面转线处 理及假边界剔除，得到瞬时水边线；
[0022] 步骤五：对瞬时水边线进行潮汐校正；
[0023] 利用潮汐校正公式对瞬时水边线进行校正，得到校正水边线上界，潮汐校正公 式；
[0024]
(5)
[0025] 其中，hl、h代表遥感影像成像时的潮高，Η代表平均大潮高潮位的潮高，ΔL为图 像上两水边线的距离，L表示水边线的校正距离，经潮汐校正后获得单期遥感影像的校正水 边线；
[0026] 输入单期瞬时水边线和校正距离，输出为单期校正水边线；
[0027] 步骤六：进行多期校正水边线的融合获取砂质岸线；
[0028] 校正水边线的上界为平均大潮高潮线，输入多期校正水边线进行融合，输出为获 得的砂质海岸线。
[0029] 如上所述的利用多期遥感影像提取砂质岸线的方法，所述步骤一中的配准是利用 已经过正射校正处理的高精度影像为基准对需校正的遥感影像进行几何配准，每景影像选 取均匀分布的同名地物点15个以上，控制点分布须均匀，影像的边缘、海岬顶点处和海湾 的凹形深处至少有一个控制点分布，并保证海岸线附近多设控制点，采用二次多项式进行 几何校正，双线性内插法进行亮度值的重采样，使校正精度控制在〇. 5个像素之内。
[0030] 如上所述的利用多期遥感影像提取砂质岸线的方法，所述步骤三中的分割阀值位 于NDWI图像的双峰之间的过度区域内。
[0031] 本发明的有益效果：
[0032] 1、本发明利用多期遥感影像，引用NDWI指数，提取瞬时水边线，同时考虑到潮汐 对岸线位置的影响，将多条潮汐校正的水边线上边线经过融合为较准确的砂质海岸线，能 够有效避免因潮间带地形（或者滩面地形）起伏造成的海岸线位置低估，同时可为我国砂 质岸线提取与变化分析提供参考依据。
[0033] 2、本发明提取的海岸线经过验证，基本和传统方法探测的海岸线重合，其均值和 均方根误差均处在较小的范围内，同时，海岸线的提取过程快速、准确，实现全自动化操作。
【附图说明】
[0034] 图1是本发明的流程图；
[0035] 图2是遥感影像经预处理后得到的NDWI图像直方图；
[0036] 图3是本发明中所述潮汐校正的原理模型图；
[0037] 图4是本发明中多期影像岸线的局部图；
[0038] 图5a是本发明中所述多期校正水边线融合的原理图；
[0039] 图5b是本发明提取的岸线示意图；
[0040] 图6是本发明提取的岸线结果图；
[0041] 图7是本发明提取的岸线精度验证图。
【具体实施方式】
[0042] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是 本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员 在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0043] 如图1所示，本实施例一种利用多期遥感影像提取砂质岸线的方法，以中国黄海 某处的多期卫星遥感影像进行砂质岸线探测做详细描述，包括以下步骤：
[0044] 步骤一：对通过卫星获得的遥感影像进行预处理；
[0045] 所述预处理包括采用图像回归法进行影像相对辐射归一化处理，再利用几何校正 方法对遥感影像进行配准，得到预处理后的遥感影像，由于影像辐射差异的不同，造成同一 地区的提取标准也有所不同，为消除同一地区不同时相影像间的辐射差异需对遥感影像进 行相对辐射归一化处理，相对辐射归一化公式如下（1)、（2)、（3):
[0046]yk=akxk+bk (1)
[0047] (2)
[0048] ⑶
[0049] 其中，ak、bk是k波段的归一化参数，xk