[发明专利]基于时序遥感图像的油田污染水体识别方法及装置

权利要求书

1.一种基于时序遥感图像的油田污染水体识别方法，其特征在于，包括：

根据研究区的多个时相的原始卫星遥感图像数据建立时序遥感图像数据集；

对所述时序遥感图像数据集中的每幅原始卫星遥感图像数据进行辐射定标、大气校正以及正射校正和图像配准，得到各波段数据对应的反射率值；

根据各时序遥感图像中各波段数据对应的反射率值确定对应时序遥感图像的水体指数；

基于所述水体指数提取对应时序遥感图像中的水体目标区域；

根据所述水体目标区域对对应的时序遥感图像中的非水体区域进行掩膜处理；

根据掩膜处理后的时序遥感图像数据集得到时序聚合遥感图像数据集；

获取所述时序聚合遥感图像数据集中各聚合遥感数据中水体目标区域内各像素点的生产废液指数；

根据水体目标在时序遥感图像集中的出现次数，基于预设的阈值-次数对应关系，获取对应的阈值；

根据水体目标区域内各像素点的生产废液指数以及对应的阈值识别该水体目标是否为生产污染水体。

2.根据权利要求1所述的基于时序遥感图像的油田污染水体识别方法，其特征在于，所述对所述时序遥感图像数据集中的每幅原始卫星遥感图像数据进行正射校正、辐射定标、大气校正和图像配准，包括：

根据多个时相的原始卫星遥感数据的头文件信息得到定标参数和校正参数；

根据所述定标参数和所述校正参数对所述时序遥感图像数据集中的每幅原始卫星遥感图像数据进行辐射定标、大气校正以及正射校正；

对校正后的图像进行地理配准。

3.根据权利要求1所述的基于时序遥感图像的油田污染水体识别方法，其特征在于，按照下述公式确定时序遥感图像中各个区域对应的水体指数：

NDWI＝((p(Green)-p(NIR))/((p(Green)+p(NIR))

其中，p(Green)和p(NIR)分别表示绿光波段和近红外波段数据对应的反射率值。

4.根据权利要求1所述的基于时序遥感图像的油田污染水体识别方法，其特征在于，还包括：

设置阈值-次数对应关系。

5.根据权利要求4所述的基于时序遥感图像的油田污染水体识别方法，其特征在于，所述设置阈值-次数对应关系，包括：

对掩膜处理后的时序遥感图像数据集进行时相叠加处理，得到叠加遥感图像数据集；

获取各叠加遥感图像中典型水体目标区域的生产废液指数的平均值以及对应的出现次数；

对各叠加遥感图像中典型水体目标区域的生产废液指数平均值以及对应的出现次数进行拟合，得到阈值-次数对应关系。

6.根据权利要求1所述的基于时序遥感图像的油田污染水体识别方法，其特征在于，所述根据掩膜处理后的时序遥感图像数据集得到时序聚合遥感图像数据集，包括：

对掩膜处理后的时序遥感图像进行叠加得到所述聚合遥感图像。

7.根据权利要求1至6任一项所述的基于时序遥感图像的油田污染水体识别方法，其特征在于，按照下述公式确定水体目标区域内各像素点的生产废液指数：

PWI＝(p(Green)-p(Red))/p(NIR)

其中，p(Green)表示绿光波段数据对应的反射率值，p(Red)表示红光波段数据对应的反射率值，p(NIR)表示近红外波段数据对应的反射率值。