建国以来

新中国成立后，1954年特大洪水有多可怕？听老人们说，是继上古传说“洪水淹天”之后，受灾面积超575万公顷，京广铁路超百天停摆。

淮河流域罕见地下来近3个月的降雨，造成全流域大洪水，淮河干流中渡站三十日最大洪量为531亿立方米，超过1931年。湖南、湖北、江西、安徽、江苏五省，有123个县市受灾，淹没耕地4755万亩，受灾人数达到3000万。

淮北平原呈现出一大片洪泛区，其中河南省淮北平原几乎出现大洪水。安徽省洪涝灾害最严重，全省2620万亩农田被淹，房屋倒塌168万间，因灾死亡达1098人，死畜1052头。江苏省农田被淹1063万亩，受灾人口832人，冲毁桥梁1071座，涵洞156个。洞庭湖区900堤垸三分之二溃决，疫病流行近洞庭湖区死亡3.3万人，这是多么可怕的数字。受影响的河南省83县2市，淹没农田1342万亩，冲毁农田水利设施33970处，倒塌房屋30万间。最小的山东受灾，淮河下游里下河地区各河段连续暴雨水位上涨，内涝严重。

2012年百年一遇特大暴雨袭击华北，引发城市内涝，造成重大人员伤亡。7月21日，中国大部分地区遭受暴雨袭击，其中北京市及其周边受灾严重，是61年以来最强暴雨造成的洪涝灾害。

自入汛以来，北京市降水量达190.3毫米，受强降水的影响，北京、天津及河北涞源、廊坊、涿州等地出现城市内涝，受害人数达到79人。此次，暴雨造成10660间房屋倒塌，160.2万人受灾，造成116.4万元的经济损失。

其中溺水导致的死亡人数最多，占73%。各医疗机构共接治因雨灾导致的伤病员1388人次，住院治疗20人。

7月24日到8月24日，连续6个台风登陆，8月2日到8月8日，“苏拉”、“达维”、“海葵”三个台风登陆我国，直接影响15个省市，这也是近10年以来影响华北地区最严重的台风。

众所周知，台风天会引起强风暴雨，从而引发的风雹、洪涝、泥石流等灾害造成较大人员伤亡和财产损失。加强城市市政建设、提高灾前预警能力和完善城市应急救援体系是这场大雨留给我们的深刻教训与反思。

2021年河南暴雨致56人遇难5人失联，洪水无情同胞有爱，万众一心战洪魔。截至7月23日18时，受灾农作物面积达44209.7公顷，直接损失达655亿元，已转移避险39万余人，已致56人遇难，5人失联。

洪水检测卫星

Sentinel-1（哨兵1号）雷达卫星应用之浙江临海市洪灾监测

雷达卫星（SAR）的全天候、全天时及能穿透一些地物的成像特点，显示出它与光学遥感器相比的优越性。因此在洪灾监测中具有比较强的优势，可以穿透浓厚的云层来监测地表的淹没情况。

201910日1时45分左右，今年第9号台风“利奇马”在浙江台州温岭市城南镇登陆。台风给浙江台州临海市带来了严重的内涝，临海市灵江流域发生流域性洪水。临海全市被淹，临海古城被洪水围困，积水最深处一度可达5米。

从雷达影像图上可以看出，临海市多处被淹，积水区呈现黑色。

Planet 遥感卫星–Planet labs 卫星星座介绍与应用

2019年8月10日，受台风“利奇马”影响，寿光市及周边县市普降特大暴雨。

由上图可见：河道两旁的洼地出现了水浸泡过的痕迹，正常的田地应该是鲜绿色，水浸泡过的略显暗红色，另外主河道的水流明显变宽

关于SAR（合成孔径雷达）图像的洪水变化检测

1. 对于洪水区域的提取很多文献采用的变化检测技术，因为SAR图像全天时、全天候、不受云雨限制的原因，目前对提取洪水区域更多采用的是SAR图像，可选用的SAR图像也比较多，比如Sentinel-1、ASRA、TerraSAR、ALOS等。

2. 由于SAR图像中，水体的后向散射系数比较低，在图像上通常呈现暗黑色，所以通常选用的是SAR图像作为检测水体的影像。

3. 在变化检测技术当中，选用两幅图像作为主影像和副影像，通过对两幅经过预处理的影像进行作差，便可获得差异图，在变化检测过程中，获取差异图是重要的一步，在很多文献当中，很多人会对如何获取差异图做明确的解释，当然获取差异图的方法也有很多，比如差值法、比值法、对数比值法和均值比值法等等，这些方法主要针对SAR图像，这里对于光学图像的变化检测不多叙述。

4. 得到差异图之后，下面的工作就是从差异图中提取出洪水变化的地方，差异图属于变化检测的第二步，那么下面的工作就是变化检测的第三步，对于洪水区域的提取，可以使用阈值法、聚类法等。

5. 阈值的方法有OTSU(大津法)、EM算法、KI算法等，后两者需要假设差异图服从何种分布，使用算法对参数进行最优估计，之后在根据贝叶斯最小错误准则求出阈值，将差异图的变化区和非变化区区分出来。

6. 到这里，变化检测的大致步骤已经完成，之后需要对图像进行后处理，因为最终变化图会有很多的散斑，这对于之后与参考图进行精度评定的时候，会出现较大的误差，这些散斑可能是因为原图像中的斑点噪声，在预处理步骤中，使用的滤波器并不能所有的散斑消除掉，所以在后处理的步骤当中，还可以再次进行滤波处理。滤波窗口通常使用的是5*5，因为过大的窗口会使得真正变化区失真。

7. 最后的精度评价部分就是获取的变化检测结果图与参考图进行精度评定，精度评定的指标包括很多，比如kappa系数，总体精度（OA）,当然图像的漏检数、误检数也是作为评价指标之一。

8. 图像数据是使用matlab进行处理的。

印度洪水预报的进展

中 国 水 灾 年 表(1950-1979年)

洪水预警级别及划分确定

卫星数据下载地址整理（包含Sentinel、Modis、Landsat等）

1. 地理空间数据云

2. 哨兵一号数据网站

3. 国家气象科学数据中心

4. MODIS数据国外下载网

5. Python数据分析实战：降雨量统计分析报告分析

6. USGD土木