一种基于海洋遥感图像的海水污染监测预警方法及系统转让专利
申请号 : CN202110722560.8
文献号 : CN113177183B
文献日 : 2021-09-14
发明人 : 孙成琪 , 贾小平 , 安连彤
申请人 : 广东海洋大学
摘要 :
权利要求 :
1.一种基于海洋遥感图像的海水污染监测预警方法,其特征在于,所述方法包括以下步骤:
步骤1,获取海洋遥感数据并提取污染待定河流集合;
步骤2,对污染待定河流集合筛选得到疑似污染区域集合;
步骤3,利用疑似污染区域集合输出污染区域并发出警报信息;
其中,获取海洋遥感数据并提取污染待定河流集合的子步骤为:步骤1.1,获取目标区域的海洋遥感数据;所述目标区域包括有至少一段陆地海岸线,海岸线上包括至少一条或多条河流入海口,目标区域中至少包括一条或多条河流;海洋遥感数据至少包括目标区域的光谱数据、遥感图像,光谱数据包括光谱图像;
步骤1.2,对目标区域的遥感图像进行边缘检测得到边缘线所划分的图像区域作为河流区域,对河流区域集合Yr里的每个河流区域进行分区获得河流子区域,由所有河流区域构成河流区域集合Yr,Yr={Y1,Y2,Y3,……,Ysn},集合Yr的大小为sn,sn为边缘检测得到的边缘线分割的区域作为的河流区域的数量;
步骤1.3,利用步骤1.1获得的光谱数据分别判断每个河流子区域是否为清澈区域;
步骤1.4,如果一个河流区域有一个或以上河流子区域是非清澈区域,则此河流区域定义为非清澈河流,如果一个河流区域全部子区域是清澈区域,则此河流区域定义为清澈河流,所有非清澈河流构成污染待定河流集合YDr,YDr={YD1,YD2,YD3,……,YDdn},污染待定河流集合YDr的大小为dn,dn≤sn。
2.根据权利要求1所述的一种基于海洋遥感图像的海水污染监测预警方法,其特征在于,步骤1.2中,对目标区域的遥感图像进行边缘检测得到边缘线所划分的图像区域作为河流区域,对河流区域集合Yr里的每个河流区域进行分区获得河流子区域,由所有河流区域构成河流区域集合的子步骤为:
步骤1.2.1,读取光谱数据,把光谱图像转化为灰度图像,对灰度图像进行高斯滤波和提取轮廓得到边缘线;
步骤1.2.2,依次对步骤1.2.1得到的灰度图像的每个像素进行标记,标记方式为像素的灰度值大于设定的阈值标记为1,像素的灰度值小于设定的阈值标记为0;
步骤1.2.3,对步骤1.2.2得到的经过标记的灰度图像作为标记图像进行平滑边缘得到平滑图像;
步骤1.2.4,以平滑图像的图像矩阵的几何中心点作为锚点,对标记为1的像素利用锚点进行膨胀运算,获得待定河流区域分区图像;以待定河流区域分区图像中各个闭合区间的集合作为分区集合;筛选出分区集合中面积最大的闭合区间作为第一参照区域并将第一参照区域从分区集合中剔除;所述闭合区间为由边缘检测得到的边缘线分割待定河流区域分割图像得到的封闭区域;
步骤1.2.5,计算第一参照区域在遥感图像中相应区域中各个像素值的算术平均值作为参照像素值,计算分区集合中各个闭合区间内部所有像素的平均像素值;筛选出平均像素值小于参照像素值的闭合区间作为参照集合;选出参照集合中面积最大的闭合区间作为第二参照区域;
步骤1.2.6,通过Harris角点检测对第二参照区域进行检测得到多个角点,得到角点集合;将各个角点两两连接从而得到各个角点连线划分得到的新增区域,从而由各个新增区域得到新增区域集合;
步骤1.2.7,扫描新增区域集合中各个新增区域,如果新增区域属于第一参照区域的位置则依次拼接新增区域集合中与当前新增区域位置相邻的新增区域构成待拼接河流区域,并且将当前新增区域标记为河流子区域,待拼接河流区域拼接新增区域集合中与待拼接河流区域位置相邻的新增区域形成新的待拼接河流区域直到待拼接河流区域不存在位置相邻的新增区域,最终将待拼接河流区域作为河流区域;
步骤1.2.8,从而得到标记的河流子区域和由所有河流区域构成河流区域集合。
3.根据权利要求2所述的一种基于海洋遥感图像的海水污染监测预警方法,其特征在于,步骤1.3中,判断每个河流子区域是否为清澈区域的子步骤为:计算当前河流子区域的清澈系数:
式中,fx为第一波段波长,所述fx的范围为[743nm,753nm],fy为第二波段波长,所述fy的范围为[2105nm,2155nm],fz为第三波段波长,所述fz的范围为[1230nm,1250nm],Agc(fx)为河流子区域在波长fx的在大气顶部经过气体吸收校正的反射比,Ar(fx)为波长fx的瑞利散射反射率,Agc(fz)为河流子区域在波长fz的在大气顶部经过气体吸收校正的反射比,Ar(fz)为波长fz的瑞利散射反射率,D(z,y)为反射系数,exp为以自然常数e为底的指数函数;
其中D(z,y)的算法为:
Agc(fy)为河流子区域在波长fy的在大气顶部经过气体吸收校正的反射比,Ar(fy)为波长fy的瑞利散射反射率;
根据清澈系数判断当前河流子区域是否为清澈区域:如果Ic(fx,fz)>N则判定当前河流子区域为非清澈区域,如果Ic(fx,fz)≤N则判定当前河流子区域是清澈区域,其中,N为所有河流子区域的清澈系数Ic(fx,fz)的算术平均值或者将N设置为1.2。
4.根据权利要求3所述的一种基于海洋遥感图像的海水污染监测预警方法,其特征在于,步骤2中,对污染待定河流集合筛选得到疑似污染区域集合的子步骤为:步骤2.1,在污染待定河流集合YDr中,选取每个河流区域清澈系数最小的子区域作为污染待定河流基准子区域,组成污染待定河流基准子区域集合YEd,YEd={YE1,YE2,……,YEdn},污染待定河流基准子区域集合YEd的大小为dn,污染待定河流基准子区域集合YEd里的每个子区域分别对应污染待定河流集合YDr的每个非清澈河流,设置记录开始获得海洋遥感数据的时刻为T1;
步骤2.2,初始化变量m为1,以采样间隔Ti获得污染待定河流集合YDr中各个子区域的水污染系数VC;
步骤2.3,设置第一条件:在时刻T1+n×Ti和时刻T1之间第m个污染待定河流中子区域的水污染系数VC增量大于其所在污染待定河流中的污染待定河流基准子区域的水污染系数VC增量,n为设定的监测间隔值;
把符合第一条件的子区域的水污染系数VC增量放入集合DRm,DRm={DR1,DR2,DR3,……,DRm0},集合DRm表示第m个污染待定河流的符合上述条件的子区域的水污染系数VC增量集合,m0为集合DRm的大小,m0表示有m0个子区域在集合DRm中,对集合DRm按水污染系数VC增量值由大到小排序;
步骤2.4,计算集合DRm的污染增量平均值式中,m0为集合DRm的大小,DRu为集合DRm中的第u个数值;
步骤2.5,如果步骤2.4得到的污染增量平均值 大于 第m个污染待定河流的污染待定河流基准子区域,并且m0的值大于第m个污染待定河流中所有子区域数量数的30%,则判定第m个污染待定河流有污染倾向,把当前污染待定河流的污染倾向值加1;
步骤2.6,把第m个污染待定河流的污染待定河流基准子区域设置为集合DRm中第一个数值所在的子区域,把T1的值设置为T1+n×Ti,经过设定的时间间隔后,执行步骤1.1到步骤1.4以重新选出污染待定河流集合,并输出所有污染待定河流的污染倾向值构成疑似污染区域集合,否则重新执行步骤2.2到步骤2.6。
5.根据权利要求4所述的一种基于海洋遥感图像的海水污染监测预警方法,其特征在于,步骤2.2中,获得各个子区域的水污染系数VC的子步骤为:步骤2.2.1,获取各个子区域的遥感图像,对遥感图像进行几何校正和大气校正,融合遥感图像获得高空间分辨率数据;
步骤2.2.2,根据高空间分辨率数据计算每个子区域的水污染系数VC:式中,hj为子区域的高空间分辨率数据的第j波段的反射率,hi为子区域的高空间分辨率数据的第i波段的反射率,hk为子区域的高空间分辨率数据的第k波段的反射率,Δfij为子区域的高空间分辨率数据的第i波段的中心波长与第j波段的中心波长的差值,Δfjk为子区域的高空间分辨率数据的第j波段的中心波长与高空间分辨率数据的第k波段的中心波长的差值。
6.根据权利要求1所述的一种基于海洋遥感图像的海水污染监测预警方法,其特征在于,步骤3中,利用疑似污染区域集合输出污染区域并发出警报信息的子步骤为步骤3.1,分别选取疑似污染区域集合里每个污染待定河流的污染倾向值,如果在设定的时间区间里污染待定河流的污染倾向值超过设定的阈值,则把当前污染待定河流定义为污染区域,把当前污染待定河流的污染倾向值设置为0;
步骤3.2,输出步骤3.1得到的所有污染区域并发出警报信息,所述警报信息包括污染区域的坐标和污染区域连接的海洋区域。
7.一种基于海洋遥感图像的海水污染监测预警系统,其特征在于海洋遥感图像模块:用于获取海洋遥感数据并提取污染待定河流集合,海洋遥感图像模块包括海洋水色卫星,所述海洋水色卫星带有中分辨率成像光谱仪MODIS传感器和SeaWiFS传感器;
数据处理模块:用于对海洋遥感图像模块获得的海洋遥感数据进行处理;
污染预警模块:用于对所有污染区域发出警报;
其中,获取海洋遥感数据并提取污染待定河流集合的子步骤为:步骤1.1,获取目标区域的海洋遥感数据;所述目标区域包括有至少一段陆地海岸线,海岸线上包括至少一条或多条河流入海口,目标区域中至少包括一条或多条河流;海洋遥感数据至少包括目标区域的光谱数据、遥感图像,光谱数据包括光谱图像;
步骤1.2,对目标区域的遥感图像进行边缘检测得到边缘线所划分的图像区域作为河流区域,对河流区域集合Yr里的每个河流区域进行分区获得河流子区域,由所有河流区域构成河流区域集合Yr,Yr={Y1,Y2,Y3,……,Ysn},集合Yr的大小为sn,sn为边缘检测得到的边缘线分割的区域作为的河流区域的数量;
步骤1.3,利用步骤1.1获得的光谱数据分别判断每个河流子区域是否为清澈区域;
步骤1.4,如果一个河流区域有一个或以上河流子区域是非清澈区域,则此河流区域定义为非清澈河流,如果一个河流区域全部子区域是清澈区域,则此河流区域定义为清澈河流,所有非清澈河流构成污染待定河流集合YDr,YDr={YD1,YD2,YD3,……,YDdn},污染待定河流集合YDr的大小为dn,dn≤sn。
说明书 :
一种基于海洋遥感图像的海水污染监测预警方法及系统
技术领域
背景技术
自身的生存环境。
岸依旧在一些突出的水环境问题。例如,环境基础设施老化,污水管网不完善、污水溢流或
直排等,部分海洋水生态破坏严重,河湖生态流量保障不足,缺乏应有的水生植被和岸边缓
冲带,部分海洋水域蓝藻水华频发; 部分地区海洋沿岸高环境风险工业企业密集分布,存
在较高的环境风险,功能和效益不断下降, 成为制约社会经济可持续发展的瓶颈, 越来越
成为公众媒体关注的焦点。因此保护和监测海洋水域成为迫切需要解决的问题。
靠的多种观测资料。
FY‑1 系列卫星、FY‑2 系列卫星、美国 NOAA 系列卫星、SeaWiFS 卫星、EOS/MODIS 卫星,以
及日本 MTSAT 系列卫星等,能获取高精度数据,这些数据经过处理后得到海温、海冰(包括
极地海冰)、海洋水色、海洋污染、海上台风、海雾,以及海岸带动态等信息,可以用于环境监
测、灾害监测、环境数值预报、科学研究等不同领域。
发明内容
择或创造条件。
区域构成河流区域集合Yr,Yr={Y1,Y2,Y3,……,Ysn},集合Yr的大小为sn,sn为边缘检测得
到的边缘线分割的区域作为的河流区域的数量;
流,所有非清澈河流构成污染待定河流集合YDr,YDr={YD1,YD2,YD3,……,YDdn},污染待定
河流集合YDr的大小为dn,dn≤sn。
子步骤为:
区间的集合作为分区集合;筛选出分区集合中面积最大的闭合区间作为第一参照区域并将
第一参照区域从分区集合中剔除;所述闭合区间为由边缘检测得到的边缘线分割待定河流
区域分割图像得到的封闭区域;
区间的平均像素值;筛选出平均像素值小于参照像素值的闭合区间作为参照集合;选出参
照集合中面积最大的闭合区间作为第二参照区域;
如果是河流区域或者海湾则会分割第一参照区域,如果不是河流或者海湾则角点连线是分
割第二参照区域),从而由各个新增区域得到新增区域集合;
区域,并且将该新增区域标记为河流子区域,待拼接河流区域拼接新增区域集合中与待拼
接河流区域位置相邻的新增区域形成新的待拼接河流区域直到待拼接河流区域不存在位
置相邻的新增区域,最终将待拼接河流区域作为河流区域;
(fx)为海洋遥感卫星检测到的在波长fx的在大气顶部经过气体吸收校正的反射比(例如反
射比为0.62),Ar(fx)为波长fx的瑞利散射反射率(例如反射率为0.18),Agc(fz)为海洋遥感
卫星检测到的在波长fz的在大气顶部经过气体吸收校正的反射比(例如反射比为0.45),Ar
(fz)为波长fz的瑞利散射反射率(例如反射率为0.25),D(z,y)为反射系数,exp为以自然常
数e为底的指数函数;
域,其中,N为所有河流子区域的清澈系数Ic(fx,fz)的算术平均值或者将N设置为1.2。
YE2,……,YEdn},污染待定河流基准子区域集合YEd的大小为dn,污染待定河流基准子区域
集合YEd里的每个子区域分别对应污染待定河流集合YDr的每个非清澈河流,设置记录开始
获得海洋遥感数据的时刻为T1;
系数VC增量,n为设定的监测间隔值;
VC增量集合,m0为集合DRm的大小,m0表示有m0个子区域在集合DRm中,对集合DRm按大到小
排序;
则判定第m个污染待定河流有污染倾向,把当前污染待定河流的污染倾向值加1;
1.1到步骤1.3以重新选出污染待定河流集合,并输出所有污染待定河流的污染倾向值构成
疑似污染区域集合,否则重新执行步骤2.2到步骤2.6。
波段的中心波长与第j波段的中心波长的差值,Δfjk为高空间分辨率数据的第j波段的中心
波长与高空间分辨率数据的第k波段的中心波长的差值。
义为污染区域,把当前污染待定河流的污染倾向值设置为0;设定的时间区间为时间间隔
[0.5,4]小时;
Spectroradiometer)和Sea‑viewing Wide Field‑of‑view Sensor(SeaWiFS)传感器;
遥感数据至少包括目标区域的光谱数据、遥感图像,光谱数据包括光谱图像;
区域构成河流区域集合Yr,Yr={Y1,Y2,Y3,……,Ysn},集合Yr的大小为sn,sn为边缘检测得
到的边缘线分割的区域作为的河流区域的数量;
流,所有非清澈河流构成污染待定河流集合YDr,YDr={YD1,YD2,YD3,……,YDdn},污染待定
河流集合YDr的大小为dn,dn≤sn。
附图说明
同的参考标号表示相同或相似的元素,显而易见地,下面描述中的附图仅仅是本发明的一
些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些
附图获得其他的附图,在附图中:
具体实施方式
明。
均属于本发明的保护范围。下述示例具体的工艺参数等也仅是合适范围内的一个示例,即
本领域技术人员可以通过本文的说明做合适的范围内选择,而并非要限定于下文示例的具
体数值。
方法及系统,所述方法包括以下步骤:
区域构成河流区域集合Yr,Yr={Y1,Y2,Y3,……,Ysn},集合Yr的大小为sn,sn为边缘检测得
到的边缘线分割的区域作为的河流区域的数量;
流,所有非清澈河流构成污染待定河流集合YDr,YDr={YD1,YD2,YD3,……,YDdn},污染待定
河流集合YDr的大小为dn,dn≤sn。
子步骤为:
值为[128,225];
区间的集合作为分区集合;筛选出分区集合中面积最大的闭合区间作为第一参照区域并将
第一参照区域从分区集合中剔除;所述闭合区间为由边缘检测得到的边缘线分割待定河流
区域分割图像得到的封闭区域;
区间的平均像素值;筛选出平均像素值小于参照像素值的闭合区间作为参照集合;选出参
照集合中面积最大的闭合区间作为第二参照区域;
如果是河流区域或者海湾则会分割第一参照区域,如果不是河流或者海湾则角点连线是分
割第二参照区域),从而由各个新增区域得到新增区域集合;
区域,将待拼接河流区域作为河流区域;
拼接河流区域,并且将该新增区域标记为河流子区域,待拼接河流区域拼接新增区域集合
中与待拼接河流区域位置相邻的新增区域形成新的待拼接河流区域直到待拼接河流区域
不存在位置相邻的新增区域,最终将待拼接河流区域作为河流区域;
(fx)为海洋遥感卫星检测到的在波长fx的在大气顶部经过气体吸收校正的反射比(例如反
射比为0.62),Ar(fx)为波长fx的瑞利散射反射率(例如反射率为0.18),Agc(fz)为海洋遥感
卫星检测到的在波长fz的在大气顶部经过气体吸收校正的反射比(例如反射比为0.45),Ar
(fz)为波长fz的瑞利散射反射率(例如反射率为0.25),D(z,y)为反射系数,exp为以自然常
数e为底的指数函数;
域,其中,N为所有河流子区域的清澈系数Ic(fx,fz)的算术平均值或者将N设置为1.2。
YE2,……,YEdn},污染待定河流基准子区域集合YEd的大小为dn,污染待定河流基准子区域
集合YEd里的每个子区域分别对应污染待定河流集合YDr的每个非清澈河流,设置记录开始
获得海洋遥感数据的时刻为T1;
系数VC增量,n为设定的监测间隔值(n一般取3000毫秒);
VC增量集合,m0为集合DRm的大小,m0表示有m0个子区域在集合DRm中,对集合DRm按大到小
排序;
则判定第m个污染待定河流有污染倾向,把当前污染待定河流的污染倾向值加1;
隔一般取3000毫秒)后,执行步骤1.1到步骤1.3以重新选出污染待定河流集合,并输出所有
污染待定河流的污染倾向值构成疑似污染区域集合,否则重新执行步骤2.2到步骤2.6。
波段的中心波长与第j波段的中心波长的差值,Δfjk为高空间分辨率数据的第j波段的中心
波长与高空间分辨率数据的第k波段的中心波长的差值。
义为污染区域,把当前污染待定河流的污染倾向值设置为0;
Spectroradiometer)和Sea‑viewing Wide Field‑of‑view Sensor(SeaWiFS)传感器;
构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的
示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特
点可以在任何的一个或多个实施例或示例中以合适的方式结合。
测预警系统,可运行的系统可包括,但不仅限于,处理器、存储器。本领域技术人员可以理
解,所述例子仅仅是一种基于海洋遥感图像的海水污染监测预警系统的示例,并不构成对
一种基于海洋遥感图像的海水污染监测预警系统的限定,可以包括比例子更多或更少的部
件,或者组合某些部件,或者不同的部件,例如所述一种基于海洋遥感图像的海水污染监测
预警系统还可以包括输入输出设备、网络接入设备、总线等。
(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列 (Field‑
Programmable Gate Array,FPGA) 或者其他可编程逻辑器件、分立门或者晶体管逻辑器
件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理
器等,所述处理器是所述一种基于海洋遥感图像的海水污染监测预警系统运行系统的控制
中心,利用各种接口和线路连接整个一种基于海洋遥感图像的海水污染监测预警系统可运
行系统的各个部分。
一种基于海洋遥感图像的海水污染监测预警系统的各种功能。所述存储器可主要包括存储
程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序
(比如声音播放功能、图像播放功能等)等;存储数据区可存储根据手机的使用所创建的数
据(比如音频数据、电话本等)等。此外,存储器可以包括高速随机存取存储器,还可以包括
非易失性存储器,例如硬盘、内存、插接式硬盘,智能存储卡(Smart Media Card, SMC),安
全数字(Secure Digital, SD)卡,闪存卡(Flash Card)、至少一个磁盘存储器件、闪存器
件、或其他易失性固态存储器件。
发明的范围由权利要求及其等同物限定。