一种基于遥感的河口海岸水动力信息增强显示方法与流程

1.本发明涉及遥感技术领域，具体涉及一种基于遥感的河口海岸水动力信息增强显示方法。


背景技术：

2.河口海岸地区往往是人口密集、经济高度发达的地区，同时又受到径流、潮流、波浪、泥沙等多重因子共同作用，水流流速、流向多变，泥沙输运复杂。准确掌握河口海岸地区的潮流流向、流态以及泥沙输移特征等动力信息，是开展河口海岸相关科学研究、规划制定、工程论证等工作的前提。
3.相对于传统的站点式水文观测方法，利用遥感技术开展河口海岸地区的流态、泥沙等的水动力监测，具有监测范围广、成本低、效率高等优势。根据数据源的不同，利用遥感影像开展河口海岸区的水动力信息监测主要包括定性、半定量、定量等方法，其中定性、半定量方法基于光学卫星影像，定量方法主要基于雷达、微波数据源，但在河口海岸等近岸水体流速反演精度相对较差。
4.基于光学卫星的水动力定性分析方法，主要是在对遥感影像进行整体增强的基础上，以河口海岸区的泥沙为示踪剂，结合相关潮位、流速等实测资料，对研究区的宏观流场、流势及悬沙空间分布特征进行目视解译和分析。遥感影像波段组合和成像后拉伸增强效果会对解译结果产生直接影响，在进行影像处理时，一般以红光、绿光、蓝光的真彩色波段组合或近红外、红光、绿光的假彩色波段组合进行成像，同时对影像进行整体拉伸增强，最后制作流态分析的基础图件；若以红光、绿光、蓝光的真彩色波段组合进行成像，水体整体偏暗，层次感不强，难以突出显示河口流态信息；若以近红外、红光、绿光的假彩色波段组合进行成像，虽然河口流态信息可以得到突出显示，但由于植被对红光的强吸收作用，陆地区域一般显示为红色，不能真实反映陆域情况，同时基础图件不美观。


技术实现要素：

5.针对现有技术的不足，本发明旨在提供一种基于遥感的河口海岸水动力信息增强显示方法。
6.为了实现上述目的，本发明采用如下技术方案：
7.一种基于遥感的河口海岸水动力信息增强显示方法，包括如下步骤：
8.s1、对原始的遥感影像进行辐射校正；
9.s2、对经过辐射校正的影像进行大气校正；
10.s3、对于经过步骤s2处理的影像，分别以红光、绿光、蓝光波段组合生成真彩色影像和以近红外、红光、绿光波段组合生成假彩色影像，并采用交互式线性拉伸方法分别对获得的真彩色影像和假彩色影像进行增强；
11.s4、根据归一化水体指数ndwi提取步骤s3中经过增强的假彩色影像的水域部分；
12.s5、将步骤s4获得的假彩色影像的水域部分的影像对位到步骤s3中经过增强的真
彩色影像的水域部分的上方进行镶嵌。
13.进一步地，步骤s1的具体过程为：获取卫星数据公布的辐射定标方法和定标参数，对原始的遥感影像进行辐射定标，获得辐射亮度影像。
14.进一步地，步骤s2的具体过程为：
15.s2.1、打开envi flaash大气校正模块，导入经过步骤s1处理的影像文件；
16.s2.2、输入传感器类型、图像区域平均海拔、成像时间；
17.s2.3、设置大气模型、气溶胶模型的类型，完成大气校正。
18.进一步地，步骤s3的具体过程为：
19.将经过步骤s2大气校正后的影像以红光、绿光、蓝光波段组合方式生成真彩色影像，然后采用交互式线性拉伸方法对真彩色影像进行增强，将增强后的真彩色影像保存为新的影像文件；
20.将经过步骤s2大气校正后的影像以近红外、红光、绿光波段组合方式生成假彩色影像，然后采用交互式线性拉伸方法对假彩色影像进行增强，将增强后的假彩色影像保存为新的影像文件。
21.进一步地，步骤s4的具体过程为：
22.s4.1、利用经过步骤s2大气校正后的影像的反射率计算归一化水体指数：
[0023][0024]
其中，r
w
(nir)和r
w
(green)分别表示近红外波段和绿光波段的反射率；
[0025]
s4.2、确定归一化水体指数阈值，生成水陆二值图像，之后将水陆二值图像通过栅格转矢量转换为矢量文件；
[0026]
s4.3、利用步骤s4.2得到的矢量文件对步骤s3中经过增强的假彩色影像的水域部分进行提取并保存为新的影像文件。
[0027]
更进一步地，步骤s5的具体过程为：将步骤s4生成的假彩色影像的水域部分的影像置于步骤s3中经过增强的真彩色影像的上方并进行镶嵌，忽略假彩色影像的水域部分的影像的0值，同时不考虑假彩色影像的水域部分的影像和真彩色影像之间的颜色平衡。
[0028]
本发明的有益效果在于：本发明方法基于光学遥感影像，根据陆域植被和水体光谱特征，首先将遥感影像分别以不同波段组合方式成像并进行增强显示，提取假彩色合成影像的水域，之后将水域和真彩色陆域进行镶嵌，最后得到突出显示水体流态信息而又真实表达陆域信息的河口海岸区水动力信息分析图件。相较于现有技术，利用本发明方法处理后的水动力分析影像色彩鲜明、图像丰富、层次好，不仅能使得水体水动力信息得到增强显示，同时又能真实表达陆域植被信息。
附图说明
[0029]
图1为本发明实施例的方法流程示意图；
[0030]
图2为本发明实施例中landsat8卫星影像增强后的真彩色影像；
[0031]
图3为本发明实施例中landsat8卫星影像增强后的假彩色影像；
[0032]
图4为本发明实施例中landsat8卫星影像经过本发明方法处理后水动力信息增强显示效果图；
[0033]
图5为本发明实施例中gf1卫星影像增强后的真彩色影像；
[0034]
图6为本发明实施例中gf1卫星影像增强后的假彩色影像；
[0035]
图7为本发明实施例中gf1卫星影像经本发明方法处理后水动力信息增强显示效果图。
具体实施方式
[0036]
以下将结合附图对本发明作进一步的描述，需要说明的是，本实施例以本技术方案为前提，给出了详细的实施方式和具体的操作过程，但本发明的保护范围并不限于本实施例。
[0037]
本实施例以广西钦江河口湾2013年12月4日landsat8卫星影像、广东珠江河口2016年12月27日gf1卫星影像数据作为示例，详细说明一种基于遥感的河口海岸水动力信息增强显示的方法的具体过程为，如图1所示：
[0038]
1、获取卫星数据的辐射定标参数：对于国外卫星影像数据，定标参数一般从头文件中读取；国产卫星数据定标参数变化比较频繁，可从中国资源卫星中心官网(url：http：//www.cresda.com/cn/)查询。
[0039]
本实施例中，对于2013年12月4日landsat8卫星影像，定标参数从头文件中读取；2016年12月27日gf1影像定标参数从中国资源卫星中心官网(url：http：//www.cresda.com/cn/)查询。
[0040]
2、对原始的遥感图像数据进行辐射校正：根据卫星数据公布的辐射定标方法及定标参数，对原始的遥感影像数据进行辐射定标，获得辐射亮度影像l(λ)。辐射亮度值按下式计算：
[0041]
l(λ)＝gain*dn bias
ꢀꢀꢀ
(1)
[0042]
式中，l(λ)为转换后的辐亮度值，单位为w/(m2.μm.sr)，λ为卫星传感器各通道的中心波长，单位μm；dn为卫星载荷观测值；gain为定标斜率，bias为定标截距，gf1号卫星参数均可从相关网站查询获取。
[0043]
3、打开envi flaash大气校正模块，导入辐射亮度影像；
[0044]
4、输入卫星数据传感器类型、图像区域平均海拔、成像时间等传感器与图像信息；
[0045]
5、设置大气模型、气溶胶模型的类型，完成大气校正。
[0046]
6、将大气校正后的影像以红光、绿光、蓝光波段组合方式生成真彩色影像，采用交互式线性拉伸方法对真彩色影像进行增强。真彩色影像重点突出显示陆域植被信息，同时确保图像丰富、色彩鲜明、层次好；将增强后的真彩色影像保存为新的影像文件。
[0047]
7、将大气校正后的影像以近红外、红光、绿光波段组合方式生成假彩色影像，采用交互式线性拉伸方法对假彩色影像进行增强。假彩色影像重点突出显示水体流态、悬沙等信息，同时确保图像丰富、色彩鲜明、层次好；将增强后的假彩色影像保存为新的影像文件。
[0048]
8、利用大气校正过遥感影像反射率计算归一化水体指数：
[0049][0050]
其中，r
w
(nir)、r
w
(green)分别表示近红外波段、绿光波段的反射率。计算过程可通过arcgis或envi波段运算功能实现。
[0051]
确定归一化水体指数阈值，本实施例中，对于两景影像阈值均取0，对于ndwi＞0判定为水体，生成水陆二值图像，计算过程可通过arcgis或envi波段运算功能实现。利用arcgis中栅格转矢量工具将水陆二值图像转换为矢量文件。
[0052]
9、利用步骤8获得的水域矢量文件对步骤7得到的假彩色影像的水域部分进行提取，将提取出来的假彩色影像的水域部分的影像保存为新的影像文件。
[0053]
10、将步骤9生成的将假彩色影像的水域部分对位到步骤6的真彩色影像的水域部分的上方然后进行镶嵌，忽略假彩色影像的水域部分的影像的0值，同时忽略假彩色影像和真彩色影像之间的颜色平衡；将镶嵌后的影像保存为新的影像文件。
[0054]
需要说明的是，步骤6、7、8都是对大气校正后的遥感图像进行处理或计算，三者的先后顺序可以调整。
[0055]
图2为landsat8影像红光、绿光、蓝光波段组合增强后的钦江河口湾流态图，图3为同景影像近红外、红光、绿光波段组合后流态图，图4为利用上述方法将两者镶嵌后的流态图。与图2相比，图4中钦江河口湾落潮阶段高含沙水体流向、流路得到了增强显示，同时外海低含沙水体颜色偏蓝，色调偏暗，与河口湾落潮水体对比鲜明，层次好，图像水体信息更为丰富，更好地显示了钦江河口湾落潮阶段的径潮相互作用。与图3相比，图4的水域信息与图3一致，但陆域以真彩色显示，真实反映了陆域尤其是植被信息，图像色彩更鲜明，成像效果更为美观。
[0056]
图5为gf1卫星影像红光、绿光、蓝光波段组合增强后的珠江河口湾流态图，图6为同景影像近红外、红光、绿光波段组合后流态图，图7为采用本方法将两者镶嵌后的流态图。与图5相比，图7中珠江河口湾涨潮阶段高含沙水体流向、流路得到了增强显示，同时外海低含沙水体颜色偏蓝，色调偏暗，与河口湾落潮水体对比鲜明，层次好，图像水体信息更为丰富，能清晰显示涨潮流流路，更好地展示了珠江河口湾落潮阶段的径潮相互作用。与图6相比，图7中的水域信息与图6一致，但陆域以真彩色显示，真实反映了陆域尤其是植被信息，图像色彩更鲜明，成像效果更为美观。
[0057]
对于本领域的技术人员来说，可以根据以上的技术方案和构思，给出各种相应的改变和变形，而所有的这些改变和变形，都应该包括在本发明权利要求的保护范围之内。