一种气象卫星洪涝监测的云阴影去除方法、装置及设备与流程

1.本发明涉及卫星监测技术领域，特别是指一种气象卫星洪涝监测的云阴影去除方法、装置及设备。


背景技术：

2.风云三号极轨气象卫星具有全球的观测能力，是大范围且快速获取洪涝灾害信息的重要数据源。洪涝灾害的发生通常是由于降水造成的，在洪涝灾害发生期间，通常会有大量的云的存在，而云在地面上形成的阴影，在光谱特征上与水体很相似，难以区分。造成在利用光学传感器提取洪涝水体的时候，会把大量的云影识别成水体，产生大量的误判，引起较大的误差。


技术实现要素：

3.本发明提供了一种气象卫星洪涝监测的云阴影去除方法、装置及设备。实现减少了气象卫星洪涝监测时提取水体的误判，进而提高了洪涝监测的精度。
4.为解决上述技术问题，本发明的实施例提供以下方案：
5.一种气象卫星洪涝监测的云阴影去除方法，包括：
6.获取具有地理定位信息的卫星通道数据；
7.根据所述卫星通道数据，得到水体信息和对应的水体二值图；
8.获取云像元信息和太阳角度信息，并根据所述云像元信息和太阳角度信息，得到云阴影像元在卫星图像上的坐标数据；
9.根据所述云像元信息和云阴影像元在卫星图像上的坐标数据，得到云层二值图和预测云阴影二值图；
10.对所述水体二值图、云层二值图和预测云阴影二值图进行云阴影去除处理，得到目标水体二值图。
11.可选的，获取具有地理定位信息的卫星通道数据，包括：
12.获取卫星的原始数据；
13.对所述卫星的原始数据进行预处理，得到具有地理定位信息的卫星通道数据。
14.可选的，根据所述卫星通道数据，得到水体信息和对应的水体二值图，包括：
15.根据所述卫星通道数据，提取水体信息；
16.根据所述水体信息，通过公式得到水体二值图；其中，r
gre
为绿光波段的反射率，r
nir
为近红外波段的反射率，r
g_nth
为水体指数的预设阈值。
17.可选的，获取云像元信息和太阳角度信息，并根据所述云像元信息和太阳角度信息，得到云阴影像元在卫星图像上的坐标数据，包括：
18.获取云像元的高度、云像元信息和太阳角度信息；所述云像元信息包括云像元在卫星图像上的坐标数据、卫星传感器高度角以及卫星传感器方位角；所述太阳角度信息包
括太阳高度角和太阳方位角；
19.根据所述云像元的高度和所述云像元信息中的云像元在卫星图像上的坐标数据、卫星传感器高度角以及卫星传感器方位角，通过公式计算得到云垂直投影在地面上的位置地物在卫星图像上的坐标数据；其中，(x
c_nadir_img
，y
c_nadir_img
)为云垂直投影在地面上的位置地物在卫星图像上的坐标数据，(x
c_img
，y
c_img
)为云像元在卫星图像上的坐标数据，hc云像元的高度，θv为卫星传感器高度角，φv为卫星传感器方位角；
20.根据所述云垂直投影在地面上的位置地物在卫星图像上的坐标数据和所述太阳角度信息中的太阳高度角和太阳方位角，通过公式计算得到云阴影像元在卫星图像上的坐标数据；其中，(x
shadow_img
，y
shadow_img
)为云阴影像元在卫星图像上的坐标数据，(x
c_nadir_img
，y
c_nadir_img
)为云垂直投影在地面上的位置地物在卫星图像上的坐标数据，hc云像元的高度，θs为太阳高度角，φs为太阳方位角。
21.可选的，根据所述云像元信息和云阴影像元在卫星图像上的坐标数据，得到云层二值图，包括：
22.根据所述云像元信息，当云像元对应的卫星通道数据满足条件：
23.得到云层二值图；
24.其中，bt
far-inf
为远红外亮温通道的亮温值，bt
far-inf_shadow
为云阴影像元的亮温值，δbt取值为2k，k为温度单位，r
red
为红通道反射率，r
nir
为近红外通道反射率，r
red_shadow
为云阴影像元的红通道反射率，r
nir_shadow
为云阴影像元的近红外通道反射率。
25.可选的，根据所述云像元信息和云阴影像元在卫星图像上的坐标数据，得到预测云阴影二值图，包括：
26.通过所述云像元信息中云像元在卫星图像上的坐标数据，获取云像元的纬度信息；
27.根据所述云像元的纬度信息和所述云阴影像元在卫星图像上的坐标数据，按照预设间隔高度，得到云像元在预设间隔高度的所有预测云阴影范围；
28.将所述所有预测云阴影范围进行并集处理，得到预测云阴影二值图。
29.可选的，对所述水体二值图、云层二值图和预测云阴影二值图进行云阴影去除处理，得到目标水体二值图，包括：
30.通过公式计算得到去除云阴影影响的水体的空间分布区；其中，c1为云层二值图，s1为预测云阴影二值图，s2为云阴影二值图，w1为水体二值图，s
w1
为预测云阴影二值图中被误判成水体的云阴影区，w2为去除云阴影影响的水体的空间分布区；
31.根据至少两个时次的去除云阴影影响的水体的空间分布区进行并集处理，得到目标水体二值图。
32.可选的，根据至少两个时次的去除云阴影影响的水体的空间分布区进行并集处理，得到目标水体二值图，包括：
33.根据至少两个时次的去除云阴影影响的水体的空间分布区，通过公式w3＝w2
t1
∪w2
t2
∪w2
t3
∪
……
∪w2
tn
，得到目标水体二值图，其中，w3为目标水体二值图，w2
tn
为第n时次的去除云阴影影响的水体的空间分布区，n为大于1的整数。
34.本发明的实施例还提供一种气象卫星洪涝监测的云阴影去除装置，所述装置包括：
35.获取模块，用于获取具有地理定位信息的卫星通道数据；
36.处理模块，用于根据所述卫星通道数据，得到水体信息和对应的水体二值图；
37.所述获取模块，还用于获取云像元信息和太阳角度信息，并根据所述云像元信息和太阳角度信息，得到云阴影像元在卫星图像上的坐标数据；
38.所述处理模块，还用于根据所述云像元信息和云阴影像元在卫星图像上的坐标数据，得到云层二值图和预测云阴影二值图；对所述水体二值图、云层二值图和预测云阴影二值图进行云阴影去除处理，得到目标水体二值图。
39.本发明的实施例还提供一种计算设备，包括：处理器、存储器、通信接口和通信总线，所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信；
40.所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行如上述的气象卫星洪涝监测的云阴影去除方法对应的操作。
41.本发明的上述方案至少包括以下有益效果：
42.本发明的上述方案，通过获取具有地理定位信息的卫星通道数据；根据所述卫星通道数据，得到水体信息和对应的水体二值图；获取云像元信息和太阳角度信息，并根据所述云像元信息和太阳角度信息，得到云阴影像元在卫星图像上的坐标数据；根据所述云像元信息和云阴影像元在卫星图像上的坐标数据，得到云层二值图和预测云阴影二值图；对所述水体二值图、云层二值图和预测云阴影二值图进行云阴影去除处理，得到目标水体二值图。减少了气象卫星洪涝监测时提取水体的误判，进而提高了洪涝监测的精度。
附图说明
43.图1为本发明实施例的气象卫星洪涝监测的云阴影去除方法的流程示意图；
44.图2为本发明提供的具体的实施例中气象卫星洪涝监测的云阴影去除方法的示意图；
45.图3为本发明实施例的气象卫星洪涝监测的云阴影去除装置的结构示意图。
具体实施方式
46.下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
47.如图1所示，本发明提供一种气象卫星洪涝监测的云阴影去除方法，包括：
48.步骤11，获取具有地理定位信息的卫星通道数据；
49.步骤12，根据所述卫星通道数据，得到水体信息和对应的水体二值图；
50.步骤13，获取云像元信息和太阳角度信息，并根据所述云像元信息和太阳角度信息，得到云阴影像元在卫星图像上的坐标数据；
51.步骤14，根据所述云像元信息和云阴影像元坐标数据，得到云层二值图和预测云阴影二值图；
52.步骤15，对所述水体二值图、云层二值图和预测云阴影二值图进行云阴影去除处理，得到目标水体二值图。
53.该实施例利用卫星获取的具有地理定位信息的卫星通道数据，结合水体信息、云像元信息和太阳角度信息，确定水体二值图、云层二值图和预测云阴影二值图，进而通过云阴影去除处理，得到目标水体二值图；实现了减少气象卫星洪涝监测时提取水体的误判，进而提高了洪涝监测的精度。
54.本发明一可选的实施例中，步骤11，包括：
55.步骤111，获取卫星的原始数据；
56.步骤112，对所述卫星的原始数据进行预处理，得到具有地理定位信息的卫星通道数据。
57.本实施例中，通过卫星上搭载的成像仪获取原始数据，并对所述原始数据进行预处理，得到具有地理定位信息的卫星通道数据；
58.其中，预处理包括定位处理、定标处理以及质量检查处理：
59.定位处理用于对原始数据中的偏移位置进行平移矫正；
60.定标处理用于通过定标系数对卫星的原始数据进行参数换算，该参数优选的包括卫星测得的dn值(digital number，遥感影像像元亮度值)；
61.质量检查处理用于筛除原始数据中偏差值较大的数据和/或有误的数据。
62.一个具体的实施例1-1中，卫星为风云三号极轨气象卫星(fy-3)，卫星上搭载的中分辨成像仪mersi(medium resolution spectral imager)或mersi-ii(medium resolution spectral imager-ii)获得原始数据，该中分辨成像仪的通道中最高空间分辨率为250m，卫星fy-3a、卫星fy-3b和卫星fy-3c搭载mersi，卫星fy-3d上搭载的mersi-ii，其中，卫星fy-3a和卫星fy-3c为上午星，卫星fy-3b和卫星fy-3d为下午星，上午星是指极轨气象卫星在上午地方时飞过目标位置，下午星则是极轨气象卫星在下午地方时飞过目标位置，通过卫星fy-3a、卫星fy-3b、卫星fy-3c以及卫星fy-3d可实现对目标位置上午和下午的组网观测，一天内可实现对同一目标位置的多次监测，获取目标位置的原始数据，并通过对原始数据进行预处理，得到具有目标位置的地理定位信息的卫星通道数据。mersi获取的卫星通道数据，如下表所示：
[0063][0064][0065]
表1
[0066]
mersi-ii获取的卫星通道数据，如下表所示：
[0067][0068][0069]
表2
[0070]
上述表1示出了fy-3的mersi有20个通道，mersi-ii有25个通道，通过将卫星fy-3的原始数据经过预处理，生成具有地理定位信息的卫星通道数据。
[0071]
本发明一可选的实施例中，步骤12包括：
[0072]
步骤121，根据所述卫星通道数据，提取水体信息；
[0073]
步骤122，根据所述水体信息，通过公式得到水体二值图；其中，r
gre
为绿光波段的反射率，r
nir
为近红外波段的反射率，r
g_nth
为水体指数的预设阈值。
[0074]
本实施例中，通过归一化差分水体指数法ndwi(normalized difference water index)从卫星通道数据中提取水体信息，通过公式得到水体二值图；其中，水体信息包括绿光波段的反射率r
gre
和近红外波段的反射率r
nir
，水体指数的预设阈值r
g_nth
是依据卫星监测的目标位置的地理条件、天气条件以及卫星监测的时间等实际情况的差异而定；
[0075]
另外，需要说明的是，卫星监测得到的水体二值图会在像元内存在云阴影的情况下受到干扰，即当卫星监测的像元内存在云阴影时，所得到的水体二值图存在误差。
[0076]
一个具体的实施例1-2中，基于上述的具体实施例1-1提供的卫星通道数据，预设水体指数的预设阈值为r
g_th
，通过mersi或者mersi-ii中的2通道的数据，得到绿光波段的反射率r
gre
，通过mersi或者mersi-ii中的4通道的数据，得到近红外波段的反射率r
nir
；然后通过公式：
[0077]
公式计算得到水体二值图，满足大于或等于r
g_th
的条件为水体部分，不满足大于或等于r
g_th
的条件为非水体部分，得到水体二值图。
[0078]
本发明一可选的实施例中，步骤13包括：
[0079]
步骤131，获取云像元的高度、云像元信息和太阳角度信息；所述云像元信息包括云像元在卫星图像上的坐标数据、卫星传感器高度角以及卫星传感器方位角；所述太阳角度信息包括太阳高度角和太阳方位角；
[0080]
步骤132，根据所述云像元的高度和所述云像元信息中的云像元在卫星图像的坐标数据、卫星传感器高度角以及卫星传感器方位角，通过公式计算得到云垂直投影在地面上的位置地物在卫星图像上坐标数据；其中，(x
c_img
，y
c_img
)为云像元在卫星图像上坐标数据，(x
c_nadir_img
，y
c_nadir_img
)为云垂直投影在地面上的位置地物在卫星图像上的坐标数据，hc云像元的高度，θv为卫星传感器高度角，φv为卫星传感器方位角；
[0081]
步骤133，根据所述云垂直投影在地面上的位置地物在卫星图像上的坐标数据和所述太阳角度信息中的太阳高度角和太阳方位角，通过公式
计算得到云阴影像元在卫星图像上坐标数据；其中，(x
shadow_img
，y
shadow_img
)为云阴影像元在卫星图像上的坐标数据，(x
c_nadir_img
，y
c_nadir_img
)为云垂直投影在地面上的位置地物在卫星图像上的坐标数据，hc云像元的高度，θs为太阳高度角，φs为太阳方位角。
[0082]
本实施例中，卫星传感器高度角和太阳高度角表示云层和云阴影之间的水平距离，太阳和卫星传感器的方位角表示云层和云阴影之间的相对方向；
[0083]
通过公式：
[0084][0085]
可计算得到云垂直投影在地面上的位置地物在卫星图像上的坐标数据(x
c_nadir_img
，y
c_nadir_img
)；
[0086]
然后通过公式：
[0087][0088]
可计算得到云阴影像元在卫星图像上的坐标数据(x
shadow_img
，y
shadow_img
)。
[0089]
本发明一可选的实施例中，步骤14中得到云层二值图包括：
[0090]
步骤141，根据所述云像元信息和云阴影像元在卫星图像上的坐标数据，当云像元对应的卫星通道数据满足条件：
[0091]
得到云层二值图；
[0092]
其中，bt
far-inf
为远红外亮温通道的亮温值，bt
far-inf_shadow
为云阴影像元的亮温值，δbt取值为2k，k为温度单位，r
red
为红通道反射率，r
nir
为近红外通道反射率，r
red_shadow
为云阴影像元的红通道反射率，r
nir_shadow
为云阴影像元的近红外通道反射率。
[0093]
本实施例中，通过步骤13得到的云像元信息和云阴影像元在卫星图像上的坐标数据，对云层的云像元信息中的卫星通道数据进行判断，当满足上述条件时，为云层部分，当不满足上述条件时，为非云层部分；通过所述方法获得云层二值图的时效性更高，可满足气象卫星在进行洪涝监测时的更高的时效性需求。
[0094]
一个具体的实施例1-3中，基于上述的具体实施例1-1提供的卫星通道数据，得到云像元信息和云阴影像元坐标数据，当像元中的云层对应的卫星通道数据满足条件：
[0095]
时，为云层部分，当不满足该条件时，为非云层部分，得到云层二值图；其中，通过mersi的5通道或者mersi-ii中的24通道的数据可获得云像元信息中的远红外亮温通道的亮温值bt
far-inf
。
[0096]
本发明一可选的实施例中，步骤14中得到预测云阴影二值图包括：
[0097]
步骤142，通过所述云像元信息中云像元在卫星图像上的坐标数据，获取云像元的纬度信息；
[0098]
步骤143，根据所述云像元的纬度信息和所述云阴影像元在卫星图像上的坐标数据，按照预设间隔高度，得到云像元在预设间隔高度的所有预测云阴影范围；
[0099]
步骤144，将所述所有预测云阴影范围进行并集处理，得到预测云阴影二值图。
[0100]
本实施例中，由于卫星通道数据无法直接获取云的高度hc，因而，通过获取云像元的纬度信息，云的高度一般在0.2km(千米)至16km；不同纬度信息，云的高度也不同；当纬度信息为北纬30度至南纬30度之间，云的高度为0.2km至8km之间；当纬度信息为南北纬30度至60度之间，云的高度为0.2km至12km之间；当纬度信息为南北纬60度至90度之间，云的高度为0.2km至16km之间；
[0101]
依据云的高度不同，按照预设间隔高度，得到相同纬度信息下，每一个预设间隔高度的云阴影范围，将每个像元内对应的纬度信息下的所有预设间隔高度的云阴影范围进行并集处理，得到预测云阴影二值图，该方式得到的预测云阴影二值图可以实现云阴影范围的最大程度的识别，减少预测云阴影二值图的误差。
[0102]
本发明一可选的实施例中，步骤15包括：
[0103]
步骤151，通过公式计算得到去除云阴影影响的水体的空间分布区；其中，c1为云层二值图，s1为预测云阴影二值图，s2为云阴影二值图，w1为水体二值图，s
w1
为预测云阴影二值图中被误判成水体的云阴影区，w2为去除云阴影影响的水体的空间分布区；
[0104]
步骤152，根据至少两个时次的去除云阴影影响的水体的空间分布区进行并集处理，得到目标水体二值图。
[0105]
本实施例中，基于步骤11至步骤14得到的水体二值图、云层二值图以及预测云阴影二值图，得到目标水体二值图；由于云层二值图和预测云阴影二值图之间存在重叠部分，将云层二值图和预测云阴影二值图之间的重叠部分设置为云层的覆盖范围，即在预测云阴影二值图中删除云层二值图和预测云阴影二值图之间的重叠部分得到云阴影二值图，具体的，通过公式s2＝s1-(c1∩s1)，得到云阴影二值图s2，其中，s1为预测云阴影二值图，c1为云层二值图；
[0106]
由于水体二值图中存在将云阴影二值图中的一部分识别为水体，通过公式s
w1
＝(w1∩s2)得到预测云阴影二值图中被误判成水体的云阴影区s
w1
；其中，s2为云阴影二值图，
w1为水体二值图；
[0107]
将预测云阴影二值图中的云阴影区去除，即可得到去除云阴影影响的水体的空间分布区，具体的，通过公式w2＝w1-s
w1
得到去除云阴影影响的水体的空间分布区w2，其中，w1为水体二值图，s
w1
为云阴影区；
[0108]
去除云阴影影响的水体的空间分布区w2中还存在部分被误识别为云阴影区，而未被识别为水体的空间分布区的部分，因此，需要通过至少两个时次的去除云阴影影响的水体的空间分布区进行并集处理，得到目标水体二值图。
[0109]
本发明一可选的实施例中，步骤152包括：
[0110]
根据至少两个时次的去除云阴影影响的水体的空间分布区，通过公式w3＝w2
t1
∪w2
t2
∪w2
t3
∪
……
∪w2
tn
，得到目标水体二值图，其中，w3为目标水体二值图，w2
tn
为第n时次的去除云阴影影响的水体的空间分布区，n为大于1的整数。
[0111]
本实施例中，气象卫星通常会提取多时次的去除云阴影影响的水体的空间分布区，通过公式w3＝w2
t1
∪w2
t2
∪w2
t3
∪
……
∪w2
tn
，可进一步的减少云层和云阴影对水体监测的干扰；
[0112]
需要说明的是，基于目标水体二值图w3，还可以结合与土地利用数据中获得的背景水体wb进行叠加分析，可以得到更加精确的水体的空间分布范围。
[0113]
如图2所示，一个具体的实施例2中，运用风云三号气象卫星进行洪涝水体监测，风云三号气象卫星上挂载有中分辨成像仪mersi(或mersi-ii)，对卫星上获取的原始数据进行预处理，得到卫星通道数据，进而对卫星通道数据进行处理，分别得到水体提取后的水体二值图w1，云层覆盖信息提取后的云层二值图c1以及预测云阴影二值图s1，将云层二值图c1和预测云阴影二值图s1进行交集运算，得到云阴影二值图s2，根据水体二值图w1和云阴影二值图s2进行交集运算，可得到误判成水体的云阴影区s
w1
，进而得到去除云阴影区s
w1
后的水体的空间分布区域w2，对同一目标位置下的多时次的水体的空间分布区域w2和云层二值图进行并集运算，结合土地利用数据水体信息作为背景水体，可得到基于风云三号气象卫星的洪涝水体监测图。
[0114]
本发明的实施例，可运用于测量洪涝水体信息，通过获取具有地理定位信息的卫星通道数据；根据所述卫星通道数据，得到水体信息和对应的水体二值图；获取云像元信息和太阳角度信息，并根据所述云像元信息和太阳角度信息，得到云阴影像元在卫星图像上的坐标数据；根据所述云像元信息和云阴影像元在卫星图像上的坐标数据，得到云层二值图和预测云阴影二值图；对所述水体二值图、云层二值图和预测云阴影二值图进行云阴影去除处理，得到目标水体二值图；实现了减少了气象卫星洪涝监测时提取水体的误判，进而提高了洪涝监测的精度。
[0115]
如图3所示，本发明还提供一种气象卫星洪涝监测的云阴影去除装置，所述装置30包括：
[0116]
获取模块31，用于获取具有地理定位信息的卫星通道数据；
[0117]
处理模块32，用于根据所述卫星通道数据，得到水体信息和对应的水体二值图；
[0118]
所述获取模块31，还用于获取云像元信息和太阳角度信息，并根据所述云像元信息和太阳角度信息，得到云阴影像元在卫星图像上的坐标数据；
[0119]
所述处理模块32，还用于根据所述云像元信息和云阴影像元在卫星图像上的坐标
数据，得到云层二值图和预测云阴影二值图；对所述水体二值图、云层二值图和预测云阴影二值图进行云阴影去除处理，得到目标水体二值图。
[0120]
可选的，获取具有地理定位信息的卫星通道数据，包括：
[0121]
获取卫星的原始数据；
[0122]
对所述卫星的原始数据进行预处理，得到具有地理定位信息的卫星通道数据。
[0123]
可选的，根据所述卫星通道数据，得到水体信息和对应的水体二值图，包括：
[0124]
根据所述卫星通道数据，提取水体信息；
[0125]
根据所述水体信息，通过公式得到水体二值图；其中，r
gre
为绿光波段的反射率，r
nir
为近红外波段的反射率，r
g_nth
为水体指数的预设阈值。
[0126]
可选的，获取云像元信息和太阳角度信息，并根据所述云像元信息和太阳角度信息，得到云阴影像元在卫星图像上的坐标数据，包括：
[0127]
获取云像元的高度、云像元信息和太阳角度信息；所述云像元信息包括云像元在卫星图像上的的坐标数据、卫星传感器高度角以及卫星传感器方位角；所述太阳角度信息包括太阳高度角和太阳方位角；
[0128]
根据所述云像元的高度和所述云像元信息中的云像元在卫星图像上的坐标数据、卫星传感器高度角以及卫星传感器方位角，通过公式计算得到云垂直投影在地面上的位置地物在卫星图像上的坐标数据；其中，(x
c_img
，y
c_img
)为云像元在卫星图像上的坐标数据，(x
c_nadir_img
，y
c_nadir_img
)为云垂直投影在地面上的位置地物在卫星图像上的的坐标数据，hc云像元的高度，θv为卫星传感器高度角，φv为卫星传感器方位角；
[0129]
根据所述云垂直投影在地面上的位置地物在卫星图像上的坐标数据和所述太阳角度信息中的太阳高度角和太阳方位角，通过公式计算得到云阴影像元在卫星图像上的坐标数据；其中，(x
shadow_img
，y
shadow_img
)为云阴影像元在卫星图像上的坐标数据，(x
c_nadir_img
，y
c_nadir_img
)为云垂直投影在地面上的位置地物在卫星图像上的的坐标数据，hc云像元的高度，θs为太阳高度角，φs为太阳方位角。
[0130]
可选的，根据所述云像元信息和云阴影像元在卫星图像上的坐标数据，得到云层二值图，包括：
[0131]
根据所述云像元信息，当云像元对应的卫星通道数据满足条件：
[0132]
得到云层二值图；
[0133]
其中，bt
far-inf
为远红外亮温通道的亮温值，bt
far-inf_shadow
为云阴影像元的亮温值，
δbt取值为2k，k为温度单位，r
red
为红通道反射率，r
nir
为近红外通道反射率，r
red_shadow
为云阴影像元的红通道反射率，r
nir_shadow
为云阴影像元的近红外通道反射率。
[0134]
可选的，根据所述云像元信息和云阴影像元在卫星图像上的坐标数据，得到预测云阴影二值图，包括：
[0135]
通过所述云像元信息中云像元在卫星图像上的坐标数据，获取云像元的纬度信息；
[0136]
根据所述云像元的纬度信息和所述云阴影像元在卫星图像上的坐标数据，按照预设间隔高度，得到云像元在预设间隔高度的所有预测云阴影范围；
[0137]
将所述所有预测云阴影范围进行并集处理，得到预测云阴影二值图。
[0138]
可选的，对所述水体二值图、云层二值图和预测云阴影二值图进行交集处理和/或并集处理，得到目标水体二值图，包括：
[0139]
通过公式计算得到去除云阴影影响的水体的空间分布区；其中，c1为云层二值图，s1为预测云阴影二值图，s2为云阴影二值图，w1为水体二值图，s
w1
为预测云阴影二值图中被误判成水体的云阴影区，w2为去除云阴影影响的水体的空间分布区；
[0140]
根据至少两个时次的去除云阴影影响的水体的空间分布区进行并集处理，得到目标水体二值图。
[0141]
可选的，根据至少两个时次的去除云阴影影响的水体的空间分布区进行并集处理，得到目标水体二值图，包括：
[0142]
根据至少两个时次的去除云阴影影响的水体的空间分布区，通过公式w3＝w2
t1
∪w2
t2
∪w2
t3
∪
……
∪w2
tn
，得到目标水体二值图，其中，w3为目标水体二值图，w2
tn
为第n时次的去除云阴影影响的水体的空间分布区，n为大于1的整数。
[0143]
需要说明的是，该装置30是与上述方法对应的装置，上述气象卫星洪涝监测的云阴影去除方法实施例中的所有实现方式均适用于该装置的实施例中，也能达到相同的技术效果。
[0144]
本发明的实施例还提供一种计算设备，包括：处理器、存储器、通信接口和通信总线，所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信；
[0145]
所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行如上述的气象卫星洪涝监测的云阴影去除方法对应的操作。
[0146]
需要说明的是，该计算设备是与上述方法对应的设备，上述气象卫星洪涝监测的云阴影去除方法实施例中的所有实现方式均适用于该计算设备的实施例中，也能达到相同的技术效果。
[0147]
本发明的实施例还提供一种计算机可读存储介质，存储有指令，所述指令在计算机上运行时，使得计算机执行如上所述的方法。
[0148]
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员
可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
[0149]
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0150]
在本发明所提供的实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
[0151]
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
[0152]
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
[0153]
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、rom、ram、磁碟或者光盘等各种可以存储程序代码的介质。
[0154]
此外，需要指出的是，在本发明的装置和方法中，显然，各部件或各步骤是可以分解和/或重新组合的。这些分解和/或重新组合应视为本发明的等效方案。并且，执行上述系列处理的步骤可以自然地按照说明的顺序按时间顺序执行，但是并不需要一定按照时间顺序执行，某些步骤可以并行或彼此独立地执行。对本领域的普通技术人员而言，能够理解本发明的方法和装置的全部或者任何步骤或者部件，可以在任何计算装置(包括处理器、存储介质等)或者计算装置的网络中，以硬件、固件、软件或者它们的组合加以实现，这是本领域普通技术人员在阅读了本发明的说明的情况下运用他们的基本编程技能就能实现的。
[0155]
因此，本发明的目的还可以通过在任何计算装置上运行一个程序或者一组程序来实现。所述计算装置可以是公知的通用装置。因此，本发明的目的也可以仅仅通过提供包含实现所述方法或者装置的程序代码的程序产品来实现。也就是说，这样的程序产品也构成本发明，并且存储有这样的程序产品的存储介质也构成本发明。显然，所述存储介质可以是任何公知的存储介质或者将来所开发出来的任何存储介质。还需要指出的是，在本发明的装置和方法中，显然，各部件或各步骤是可以分解和/或重新组合的。这些分解和/或重新组合应视为本发明的等效方案。并且，执行上述系列处理的步骤可以自然地按照说明的顺序按时间顺序执行，但是并不需要一定按照时间顺序执行。某些步骤可以并行或彼此独立地执行。
[0156]
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。