CN107967696B - 一种水面船只雷达遥感检测方法、电子设备 - Google Patents
一种水面船只雷达遥感检测方法、电子设备 Download PDFInfo
- Publication number
- CN107967696B CN107967696B CN201711178445.9A CN201711178445A CN107967696B CN 107967696 B CN107967696 B CN 107967696B CN 201711178445 A CN201711178445 A CN 201711178445A CN 107967696 B CN107967696 B CN 107967696B
- Authority
- CN
- China
- Prior art keywords
- image
- gray level
- ship
- pixels
- remote sensing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 21
- 238000001514 detection method Methods 0.000 title claims abstract description 20
- 238000001914 filtration Methods 0.000 claims abstract description 15
- 238000012880 independent component analysis Methods 0.000 claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims abstract description 12
- 238000000034 method Methods 0.000 claims abstract description 12
- 238000000926 separation method Methods 0.000 claims abstract description 12
- 230000000877 morphologic effect Effects 0.000 claims abstract description 11
- 238000010606 normalization Methods 0.000 claims abstract description 11
- 230000001186 cumulative effect Effects 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 4
- 239000002131 composite material Substances 0.000 abstract description 2
- 238000012544 monitoring process Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
- G06T7/41—Analysis of texture based on statistical description of texture
- G06T7/44—Analysis of texture based on statistical description of texture using image operators, e.g. filters, edge density metrics or local histograms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Probability & Statistics with Applications (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种水面船只雷达遥感检测方法、电子设备,方法包括:输入原始灰度图像,对其直方图均衡化,以实现图像增强;对图像增强后的原始灰度图像进行3×3窗口大小的Gamma滤波;计算灰度共生矩阵,获得纹理特征方差统计量灰度图像;进行形态学滤波的闭运算并与原始灰度图像进行波段组合生成新的两波段合成图像;对两波段合成图像进行独立成分分析,然后将独立分量IC2求反,进行归一化处理,并统计其直方图,取直方图中累计百分比达到预设值时对应的像元值为船只信息分离阈值K;将像元值大于船只信息分离阈值K的像元判定为船只像元,统计得到船只的数量、大小及空间位置。本发明具有易实现,高效率和高精度的技术特点。
Description
技术领域
本发明涉及卫星遥感目标提取技术领域,特别是指一种水面船只雷达遥感检测方法、电子设备。
背景技术
水面船只目标的卫星遥感自动检测、识别和定位,对于船只搜救,渔业管理、航运监控、安全监测、非法采挖和非法走私监测、海域态势感知以及海防预警等多个领域具有十分广阔的应用前景,在民用和军事领域均有重要意义。合成孔径雷达(SAR)具有全天时、全天候、大范围及高重复观测频率等对地观测特点,是海洋舰船监测的重要手段之一。近年来,以我国高分三号为代表的一批新型的高空间分辨率、短重访周期、多极化的合成孔径雷达成像卫星的涌现,为水面船只的遥感监测提供了极为丰富的雷达遥感数据源。但是目前使用的基于雷达遥感影像的水面船只目标检测方法普遍存在工作效率低下、检测精度不佳的技术问题。
发明内容
有鉴于此,本发明的目的在于提出一种易于实施、且具有高工作效率和高检测精度的水面船只雷达遥感检测方法、电子设备。
基于上述目的本发明提供的一种水面船只雷达遥感检测方法,包括:
输入一幅单极化的合成孔径雷达2级产品原始灰度图像,对其进行直方图均衡化,以实现图像增强;
对图像增强后的所述原始灰度图像进行3×3窗口大小的Gamma滤波;
基于所述3×3窗口大小计算灰度共生矩阵,获得纹理特征方差统计量灰度图像;
统计所述纹理特征方差统计量灰度图像的最小值和最大值,进行归一化处理,并将其扩展至8比特的数值范围;
基于所述3×3窗口大小对扩展至8比特的所述纹理特征方差统计量灰度图像进行形态学滤波的闭运算;
将形态学滤波后的所述纹理特征方差统计量灰度图像与所述原始灰度图像进行波段组合生成新的两波段合成图像;
对所述两波段合成图像进行独立成分分析;对所述独立成分分析的结果中的独立分量IC2求反,将暗像元转换为亮像元,进行归一化处理,然后将其扩展至8比特的数值范围,并统计其直方图,取直方图中累计百分比达到预设值时对应的像元值为船只信息分离阈值K;
将像元值大于所述船只信息分离阈值K的像元判定为船只像元,基于所述船只像元统计得到船只的数量、大小及空间位置。
在一些实施方式中,所述灰度共生矩阵的计算公式为:
P(i,j)=#{[(r1,c1),(r2,c2)]∈M×M|I(r1,c1)=i,I(r2,c2)=j};
其中,#{x}表示集合x的元素个数;i和j分别为两个像素的灰度值;(r1,c1)与(r2,c2)表示图像空间上的一个点对的像素坐标。
在一些实施方式中,基于所述灰度共生矩阵,纹理特征均值统计量的计算公式为:
所述纹理特征方差统计量的计算公式为:
在一些实施方式中,所述预设值为98.5%。
另一方面,本发明实施例还提供了一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上任意一项所述的方法。
从上面所述可以看出,本发明提供的水面船只雷达遥感检测方法、电子设备,通过对水面船只雷达影像纹理特征的形态学滤波构建船只雷达影像异常体,并利用独立成分分析方法实现船只与背景水域的分离,具有易实现,高效率和高精度的技术特点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的水面船只雷达遥感检测方法流程图;
图2为本发明实施例中Gamma滤波后的灰度图像;
图3为本发明实施例中纹理特征方差统计量灰度图像;
图4为本发明实施例中形态学滤波闭运算后的灰度图像;
图5为本发明实施例中第二个独立分量IC2的灰度图像;
图6为本发明实施例中求反后8比特的第二个独立分量IC2的灰度图像;
图7为本发明实施例中提取的船只信息二值图像。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
本发明实施例提供了一种水面船只雷达遥感检测方法。参考图1,所述水面船只雷达遥感检测方法,包括以下步骤:
步骤101、输入一幅单极化的合成孔径雷达2级产品原始灰度图像,对其进行直方图均衡化,以实现图像增强。
本步骤中,所述的原始灰度图像来自于合成孔径雷达(SAR),即SAR灰度图像。
步骤102、对图像增强后的所述原始灰度图像进行3×3窗口大小的Gamma滤波。
图像增强后的所述原始灰度图像进行Gamma滤波后的结果如图2所示。
步骤103、基于所述3×3窗口大小计算灰度共生矩阵(Gray Level Co-occurrenceMatrix,GLCM),获得纹理特征方差统计量灰度图像。
本步骤中,所述灰度共生矩阵的计算公式为:
P(i,j)=#{[(r1,c1),(r2,c2)]∈M×M|I(r1,c1)=i,I(r2,c2)=j};
其中,#{x}表示集合x的元素个数;i和j分别为两个像素的灰度值;(r1,c1)与(r2,c2)表示图像空间上的一个点对的像素坐标。取像元对间距为d,两像元连线的方向为α,取不同间距d,不同方向α的灰度共生矩阵为P(i,j,d,α)。
基于所述灰度共生矩阵,纹理特征均值统计量的计算公式为:
所述纹理特征方差统计量的计算公式为:
通过计算灰度共生矩阵获得方差统计量的灰度图像如图3所示。
步骤104、统计所述纹理特征方差(Variance)统计量灰度图像的最小值和最大值,进行归一化处理,并将其扩展至8比特(0~255)的数值范围。
步骤105、基于所述3×3窗口大小对扩展至8比特的所述纹理特征方差统计量灰度图像进行形态学滤波的闭运算。
本步骤中,经过形态学滤波的闭运算后,纹理特征方差统计量灰度图像如图4所示。
步骤106、将形态学滤波后的所述纹理特征方差统计量灰度图像与所述原始灰度图像进行波段组合生成新的两波段合成图像。
步骤107、对所述两波段合成图像进行独立成分分析(ICA);对所述独立成分分析的结果中的独立分量IC2求反,将暗像元转换为亮像元,进行归一化处理,然后将其扩展至8比特(0~255)的数值范围,并统计其直方图,取直方图中累计百分比达到预设值时对应的像元值为船只信息分离阈值K。
本步骤中,所述两波段合成图像经过独立成分分析后,在独立成分分析的结果中,船只信息位于第二个独立分量IC2中,表现为暗像元,如图5所示。然后经过求反、归一化处理、扩展至8比特等处理后,获得的结果如图6所示。进一步的,统计其直方图,取直方图中累计百分比达到预设值时对应的像元值为船只信息分离阈值K。参考表1,以图1至图5中使用的图像为例,取预设值为98.5%,则相应的能够获得船只信息分离阈值K为254。
表1求反后8比特的第二个独立分量SAR灰度图像统计特征
步骤108、将像元值大于所述船只信息分离阈值K的像元判定为船只像元,基于所述船只像元统计得到船只的数量、大小及空间位置。
本步骤中,根据船只信息分离阈值K判定图像中的船只像元。在本实施例的示例中,像元值大于阈值K(K=254)的像元即为船只像元,如图7所示;根据合成孔径雷达图像包含的雷达遥感数据,统计可得到船只的数量为76只、其面积大小及其中心位置,如表2所示。
表2 SAR雷达探测的船只信息统计表
由上述实施例可见,本发明提出的水面船只雷达遥感检测方法,通过对水面船只雷达影像纹理特征的形态学滤波构建船只雷达影像异常体,并利用独立成分分析(ICA)方法实现船只与背景水域的分离,具有易实现,高效率和高精度的技术特点。利用多个海域30多幅合成孔径雷达图像的水面船只遥感探测结果表明,平均精度为98.5%,且其计算流程简单,算法高效,容易嵌入到各类遥感图像处理软件中。
基于同一发明构思,本发明实施例还提供了一种电子设备。所述电子设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上述水面船只雷达遥感检测方法实施例所述的方法。
上述实施例的电子设备用于实现前述实施例中相应的方法,并且具有相应的方法实施例的有益效果,在此不再赘述。
所属领域的普通技术人员应当理解:以上任何实施例的讨论仅为示例性的,并非旨在暗示本公开的范围(包括权利要求)被限于这些例子;在本发明的思路下,以上实施例或者不同实施例中的技术特征之间也可以进行组合,步骤可以以任意顺序实现,并存在如上所述的本发明的不同方面的许多其它变化,为了简明它们没有在细节中提供。
另外,为简化说明和讨论,并且为了不会使本发明难以理解,在所提供的附图中可以示出或可以不示出与集成电路(IC)芯片和其它部件的公知的电源/接地连接。此外,可以以框图的形式示出装置,以便避免使本发明难以理解,并且这也考虑了以下事实,即关于这些框图装置的实施方式的细节是高度取决于将要实施本发明的平台的(即,这些细节应当完全处于本领域技术人员的理解范围内)。在阐述了具体细节(例如,电路)以描述本发明的示例性实施例的情况下,对本领域技术人员来说显而易见的是,可以在没有这些具体细节的情况下或者这些具体细节有变化的情况下实施本发明。因此,这些描述应被认为是说明性的而不是限制性的。
尽管已经结合了本发明的具体实施例对本发明进行了描述,但是根据前面的描述,这些实施例的很多替换、修改和变型对本领域普通技术人员来说将是显而易见的。例如,其它存储器架构(例如,动态RAM(DRAM))可以使用所讨论的实施例。
本发明的实施例旨在涵盖落入所附权利要求的宽泛范围之内的所有这样的替换、修改和变型。因此,凡在本发明的精神和原则之内,所做的任何省略、修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种水面船只雷达遥感检测方法,其特征在于,包括:
输入一幅单极化的合成孔径雷达2级产品原始灰度图像,对其进行直方图均衡化,以实现图像增强;
对图像增强后的所述原始灰度图像进行3×3窗口大小的Gamma滤波;
基于所述3×3窗口大小计算灰度共生矩阵,获得纹理特征方差统计量灰度图像;
统计所述纹理特征方差统计量灰度图像的最小值和最大值,进行归一化处理,并将其扩展至8比特的数值范围;
基于所述3×3窗口大小对扩展至8比特的所述纹理特征方差统计量灰度图像进行形态学滤波的闭运算;
将形态学滤波后的所述纹理特征方差统计量灰度图像与所述原始灰度图像进行波段组合生成新的两波段合成图像;
对所述两波段合成图像进行独立成分分析;对所述独立成分分析的结果中的独立分量IC2求反,将暗像元转换为亮像元,进行归一化处理,然后将其扩展至8比特的数值范围,并统计其直方图,取直方图中累计百分比达到预设值时对应的像元值为船只信息分离阈值K;
将像元值大于所述船只信息分离阈值K的像元判定为船只像元,基于所述船只像元统计得到船只的数量、大小及空间位置。
2.根据权利要求1所述的水面船只雷达遥感检测方法,其特征在于,所述灰度共生矩阵的计算公式为:
P(i,j)=#{[(r1,c1),(r2,c2)]∈M×M|I(r1,c1)=i,I(r2,c2)=j};
其中,#{x}表示集合x的元素个数;i和j分别为两个像素的灰度值;(r1,c1)与(r2,c2)表示图像空间上的一个点对的像素坐标。
4.根据权利要求1所述的水面船只雷达遥感检测方法,其特征在于,所述预设值为98.5%。
5.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1至4任意一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711178445.9A CN107967696B (zh) | 2017-11-23 | 2017-11-23 | 一种水面船只雷达遥感检测方法、电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711178445.9A CN107967696B (zh) | 2017-11-23 | 2017-11-23 | 一种水面船只雷达遥感检测方法、电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107967696A CN107967696A (zh) | 2018-04-27 |
CN107967696B true CN107967696B (zh) | 2021-06-25 |
Family
ID=62000754
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711178445.9A Active CN107967696B (zh) | 2017-11-23 | 2017-11-23 | 一种水面船只雷达遥感检测方法、电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107967696B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109583319B (zh) * | 2018-11-07 | 2020-12-15 | 北京市遥感信息研究所 | 一种基于海面红外遥感图像的舰船检测方法 |
CN110456387A (zh) * | 2019-08-14 | 2019-11-15 | 上海卫星工程研究所 | 主动遥感卫星建立屏障带对飞机搜寻的方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101520896B (zh) * | 2009-03-30 | 2012-05-30 | 中国电子科技集团公司第十研究所 | 利用光学遥感影像自动检测云层干扰舰船目标的方法 |
CN102831440B (zh) * | 2012-08-20 | 2015-03-18 | 中国测绘科学研究院 | 一种广域遥感影像决策树分类方法及装置 |
CN104036239B (zh) * | 2014-05-29 | 2017-05-10 | 西安电子科技大学 | 基于特征融合和聚类的高分辨sar图像快速舰船检测方法 |
CN105095846B (zh) * | 2014-09-28 | 2019-07-30 | 航天恒星科技有限公司 | 面向遥感图像海陆分割的区域生长种子点提取方法及系统 |
CN105279745B (zh) * | 2015-10-28 | 2018-09-07 | 浙江工业大学 | 遥感影像色彩归一化的多特征点约束直方图规整方法 |
CN106485722A (zh) * | 2016-09-21 | 2017-03-08 | 北京航天宏图信息技术股份有限公司 | 一种遥感影像中靠港舰船检测方法 |
CN106886747B (zh) * | 2016-12-28 | 2019-09-06 | 中国科学院长春光学精密机械与物理研究所 | 一种基于扩展小波变换的复杂背景下舰船检测方法 |
-
2017
- 2017-11-23 CN CN201711178445.9A patent/CN107967696B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107967696A (zh) | 2018-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | An improved superpixel-level CFAR detection method for ship targets in high-resolution SAR images | |
Li et al. | Haze and thin cloud removal via sphere model improved dark channel prior | |
CN107808161A (zh) | 一种基于光视觉的水下目标识别方法 | |
JP6977873B2 (ja) | 画像処理装置、画像処理方法、及び、画像処理プログラム | |
Vijay et al. | Gray scale image segmentation using OTSU Thresholding optimal approach | |
Ke et al. | SAR ship detection based on an improved faster R-CNN using deformable convolution | |
CN107967696B (zh) | 一种水面船只雷达遥感检测方法、电子设备 | |
Cui et al. | An image edge detection method based on haar wavelet transform | |
Li et al. | Automatic infrared ship target segmentation based on structure tensor and maximum histogram entropy | |
Soundrapandiyan et al. | Adaptive pedestrian detection in infrared images using fuzzy enhancement and top-hat transform | |
CN107464255B (zh) | 一种基于信息量与多尺度异常检测的船舶目标检测方法 | |
Deng et al. | Ship detection from optical satellite image using optical flow and saliency | |
CN113487642A (zh) | 一种显著性视觉毫米波体外目标检测方法 | |
CN115035350B (zh) | 一种基于边缘检测增强的对空地、地面背景小目标检测方法 | |
CN111428624A (zh) | 一种光学遥感图像在轨舰船目标检测方法 | |
Zehtabian et al. | Adaptive restoration of multispectral datasets used for SVM classification | |
Maussang et al. | Automated segmentation of SAS images using the mean-standard deviation plane for the detection of underwater mines | |
CN105551029A (zh) | 一种多光谱遥感图像舰船检测方法 | |
CN112233042B (zh) | 一种含非合作目标的大场景sar图像快速生成方法 | |
Dong-dong et al. | A feature combining spatial and structural information for SAR image classification | |
CN104268874B (zh) | 一种基于正态分布函数的非相参雷达图像背景建模方法 | |
CN110472472B (zh) | 基于sar遥感图像的机场检测方法与装置 | |
el Abbadi et al. | Single image haze removal via accurate atmosphere light | |
Gu et al. | A novel procedure for land masking in ocean-land segmentation from SAR images | |
Karvir et al. | Power spectrum weighted edge analysis for straight edge detection in images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |