CN108205134B - 一种极化合成孔径雷达影像的次地表信息增强方法 - Google Patents

一种极化合成孔径雷达影像的次地表信息增强方法 Download PDF

Info

Publication number
CN108205134B
CN108205134B CN201611164449.7A CN201611164449A CN108205134B CN 108205134 B CN108205134 B CN 108205134B CN 201611164449 A CN201611164449 A CN 201611164449A CN 108205134 B CN108205134 B CN 108205134B
Authority
CN
China
Prior art keywords
phase
image
point
pixel
synthetic aperture
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
Application number
CN201611164449.7A
Other languages
English (en)
Other versions
CN108205134A (zh
Inventor
朱黎江
叶发旺
李美玉
孟树
刘洪成
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Research Institute of Uranium Geology
Original Assignee
Beijing Research Institute of Uranium Geology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Research Institute of Uranium Geology filed Critical Beijing Research Institute of Uranium Geology
Priority to CN201611164449.7A priority Critical patent/CN108205134B/zh
Publication of CN108205134A publication Critical patent/CN108205134A/zh
Application granted granted Critical
Publication of CN108205134B publication Critical patent/CN108205134B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904—SAR modes
    • G01S13/9076—Polarimetric features in SAR
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明属于雷达影像技术领域，具体涉及一种极化合成孔径雷达影像的次地表信息增强方法；本发明的目的，针对现有技术不足，提供一种基于两个交叉极化图像，通过相位参考处理，抑制了地表信息，提高了次地表信息的表达能力的一种极化合成孔径雷达影像的次地表信息增强方法；步骤一，设定极化合成孔径雷达图像两个交叉极化通道的复数图像分别为：Shv和Svh，针对每一像素点的复数数值分离幅度和相位，相应形成的相位图像分别记录为Phasehv和Phasevh；步骤二，分别计算相位图像Phasehv和Phasevh的均值，分别记录为Mhv和Mvh；步骤三，逐点计算。

Description

一种极化合成孔径雷达影像的次地表信息增强方法
技术领域
本发明属于雷达影像技术领域，具体涉及一种极化合成孔径雷达影像的次地表信息增强方法。
背景技术
合成孔径雷达采用微波作为探测手段，对地表具有一定的探测深度。通常，雷达图像中的像素数值不但包含地表目标表面的雷达回波信息，也包含地表目标一定深度次地表下的雷达回波信息；并且次地表回波强度远小于地表雷达回波强度。
在获取目标区域全极化雷达图像的过程中，雷达传感器需要先后两次发射和接收雷达信号，每次接收形成一个同极化图像和一个交叉极化图像；两次信号发射接收的间隔时间非常短，因此一般被近似为同时发射和接收，并认为两个交叉极化图像包含的信息是相同的。
本方法基于两个交叉极化图像，通过相位参考处理，抑制了地表信息，提高了次地表信息的表达能力。
发明内容
本发明的目的，针对现有技术不足，提供一种基于两个交叉极化图像，通过相位参考处理，抑制了地表信息，提高了次地表信息的表达能力的一种极化合成孔径雷达影像的次地表信息增强方法。
本发明的技术方案是：
包括以下步骤：
步骤一，设定极化合成孔径雷达图像两个交叉极化通道的复数图像分别为：Shv和Svh，针对每一像素点的复数数值分离幅度和相位，相应形成的相位图像分别记录为Phasehv和Phasevh；
步骤二，分别计算相位图像Phasehv和Phasevh的均值，分别记录为Mhv和Mvh；
每一幅图像都可以认为是一个二维矩阵，有固定的长和宽；图像中的每个像素都有一个具体的位置和像素数值，所谓均值的计算即是将所有像素数值相加，然后除以像素总个数的结果，设图像的长宽分别为M和N，每个像素表达为xi,j，x是其像素数值，下标对(i，j)是像素在图像中的位置，则均值可以表达为
Figure GDA0003142532290000021
步骤三，逐点计算如果(Phasehv－Mhv)＜0，
则Phasehv＝Phasehv－Mhv+π；
逐点计算，如果(Phasevh－Mvh)＜0，则Phasevh＝Phasevh－Mvh+π；
逐点计算Phasehv和Phasevh差的绝对值，即|Phasevh－Phasehv|记录为
Figure GDA0003142532290000022
逐点计算，如果
Figure GDA0003142532290000023
则
Figure GDA0003142532290000024
计算的结果记录为R，R即为目标区域次地表信息增强的计算结果；
逐点计算，即针对图像中的每一个像素数值进行具体的数学运算。
本发明的有益效果是：
由于雷达传感器两次发射接收信号的时间间隔非常短，两次的观测几何完全一致，利用这一特性，本方法巧妙的削弱了地表回波信息，通过交叉极化通道图像的相位参考，有效突出了目标区的次地表信息，提高了遥感解译人员对目标性质的认识，有利于提高极化雷达图像分类、分析的精度和准确度。
附图说明
图1是Shv通道实部(左)和虚部(右)图像；
图2 Shv通道幅度(左)和相位(右)图像；
图3 Svh通道实部(左)和虚部(右)图像；
图4 Svh通道幅度(左)和相位(右)图像；
图5对Phasehv的计算结果，仍然记录为Phasehv
图6对Phasevh的计算结果，仍然记录为Phasevh
图7新的计算结果，记录为
Figure GDA0003142532290000031
图8最终计算结果，记录为R。
具体实施方式
下面结合附图与实施例对本发进行进一步的介绍：
包括以下步骤：
步骤一，设定极化合成孔径雷达图像两个交叉极化通道的复数图像分别为：Shv和Svh，针对每一像素点的复数数值分离幅度和相位，相应形成的相位图像分别记录为Phasehv和Phasevh；
步骤二，分别计算相位图像Phasehv和Phasevh的均值，分别记录为Mhv和Mvh；
每一幅图像都可以认为是一个二维矩阵，有固定的长和宽；图像中的每个像素都有一个具体的位置和像素数值，所谓均值的计算即是将所有像素数值相加，然后除以像素总个数的结果。设图像的长宽分别为M和N，每个像素表达为xi,j，x是其像素数值，下标对(i，j)是像素在图像中的位置，则均值可以表达为
步骤三，逐点计算如果(Phasehv－Mhv)＜0，
则Phasehv＝Phasehv－Mhv+π；
逐点计算，如果(Phasevh－Mvh)＜0，则Phasevh＝Phasevh－Mvh+π；
逐点计算Phasehv和Phasevh差的绝对值，即|Phasevh－Phasehv|记录为
Figure GDA0003142532290000041
逐点计算，如果
Figure GDA0003142532290000042
则
Figure GDA0003142532290000043
计算的结果记录为R，R即为目标区域次地表信息增强的计算结果；
逐点计算，即针对图像中的每一个像素数值进行具体的数学运算。
实施例
全极化雷达传感器开展对地遥感探测工作时，通常可以同时获得地面同一场景的4幅雷达图像，并分别记录为：Shh、Shv、Svh和Svv，其中Shh和Svv称为同极化图像，而Shv和Svh称为交叉极化图像。
每一个图像都是复图像，所谓复图像，即图像中的每一个像素的数值不是一般意义图像中的实数值，而是复数(复数既可以用实部和虚部表示，也可以用幅度和相位角度来表示，即
Figure GDA0003142532290000044
本文计算中用到的是相位)，所以每一景图像由一般意义上的两景图像构成，分别称为实部和虚部图像，也可以表示成一个幅度图像和一个相位图像。
下面应用一实际计算例子来说明算法中每一步计算过程和结果。本实例用到的数据为新疆某地区radarsat-2全极化数据。
1、雷达复图像的幅度相位分离
求图像中每一像素的复数模和相位(标准的数据变换，复数的指数表达方式)，实现幅度和相位分离。
2、分别计算相位图像Phasehv和Phasevh的均值，分别记录为Mhv和Mvh。
将两幅相位图像的像素值各自求和，除以图像的像素总个数，即可以得到各自的均值，这一步也可以通过常规的图像处理工具软件获得，比较简单和方便。
经过计算，两幅相位图像的均值分别为：
Mhv＝-0.001608和Mvh＝0.014770
3、逐点计算，如果(Phasehv－Mhv)＜0，则Phasehv＝Phasehv－Mhv+π。
对图2中的相位图像的每一个像素数值按上述公式计算，并用计算结果代替原图像的像素数值。
4、逐点计算，如果(Phasevh－Mvh)＜0，则Phasevh＝Phasevh－Mvh+π。
对图4中的相位图像的每一个像素数值按上述公式计算，并用计算结果代替原图像的像素数值。
5、逐点计算Phasehv和Phasevh差的绝对值，即|Phasevh－Phasehv|记录为
Figure GDA0003142532290000051
对第3步和第4步的计算结果，应用上式求差的绝对值。
6、逐点计算，如果
Figure GDA0003142532290000052
则
Figure GDA0003142532290000053
计算的结果记录为R，R即为目标区域次地表信息增强的计算结果。
对图7的计算结果应用上式进行计算，结果记录为R，R是最后的增强结果。

Claims (1)

1.一种极化合成孔径雷达影像的次地表信息增强方法，其特征在于：包括以下步骤：
步骤一，设定极化合成孔径雷达图像两个交叉极化通道的复数图像分别为：Shv和Svh，针对每一像素点的复数数值分离幅度和相位，相应形成的相位图像分别记录为Phasehv和Phasevh；
步骤二，分别计算相位图像Phasehv和Phasevh的均值，分别记录为Mhv和Mvh；
每一幅图像都可以认为是一个二维矩阵，有固定的长和宽；图像中的每个像素都有一个具体的位置和像素数值，所谓均值的计算即是将所有像素数值相加，然后除以像素总个数的结果，设图像的长宽分别为M和N，每个像素表达为xi,j，x是其像素数值，下标对(i,j)是像素在图像中的位置，则均值可以表达为
Figure FDA0003142532280000011
步骤三，逐点计算如果(Phasehv－Mhv)＜0，
则Phasehv＝Phasehv－Mhv+π；
逐点计算，如果(Phasevh－Mvh)＜0，则Phasevh＝Phasevh－Mvh+π；
逐点计算Phasehv和Phasevh差的绝对值，即|Phasevh－Phasehv|记录为
Figure FDA0003142532280000012
逐点计算，如果
Figure FDA0003142532280000013
则
Figure FDA0003142532280000014
计算的结果记录为R，R即为目标区域次地表信息增强的计算结果；
逐点计算，即针对图像中的每一个像素数值进行具体的数学运算。
CN201611164449.7A 2016-12-16 2016-12-16 一种极化合成孔径雷达影像的次地表信息增强方法 Active CN108205134B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611164449.7A CN108205134B (zh) 2016-12-16 2016-12-16 一种极化合成孔径雷达影像的次地表信息增强方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611164449.7A CN108205134B (zh) 2016-12-16 2016-12-16 一种极化合成孔径雷达影像的次地表信息增强方法

Publications (2)

Publication Number Publication Date
CN108205134A CN108205134A (zh) 2018-06-26
CN108205134B true CN108205134B (zh) 2021-09-17

Family

ID=62601325

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611164449.7A Active CN108205134B (zh) 2016-12-16 2016-12-16 一种极化合成孔径雷达影像的次地表信息增强方法

Country Status (1)

Country Link
CN (1) CN108205134B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101078769A (zh) * 2006-05-25 2007-11-28 中国科学院中国遥感卫星地面站 单次全极化合成孔径雷达图像反演数字高程模型的方法
CN101419284A (zh) * 2008-08-08 2009-04-29 哈尔滨工业大学 由森林覆盖下目标参数反演模型获得人造目标信息的方法
CN101482616A (zh) * 2008-08-13 2009-07-15 中国科学院电子学研究所 一种地形测量方法
CN101498789A (zh) * 2009-02-25 2009-08-05 中国测绘科学研究院 一种基于全极化合成孔径雷达的地物目标分类方法和装置
CN101685155A (zh) * 2008-09-27 2010-03-31 中国科学院电子学研究所 基于极化干涉合成孔径雷达数据优化干涉相干系数的方法
CN101706577A (zh) * 2009-12-01 2010-05-12 中南大学 InSAR监测高速公路路面沉降方法
CN103940834A (zh) * 2014-05-09 2014-07-23 中国科学院电子学研究所 采用合成孔径雷达技术测量土壤湿度的方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103529447A (zh) * 2012-07-02 2014-01-22 中国科学院对地观测与数字地球科学中心 面向城市区域的极化四分量分解方法
US9164170B2 (en) * 2013-01-14 2015-10-20 Honeywell International Inc. Systems and methods for autotilting a ground-mapping radar
JP2015094657A (ja) * 2013-11-12 2015-05-18 三菱電機株式会社 合成開口レーダ信号処理装置
FR3029045B1 (fr) * 2014-11-26 2018-03-02 Thales Architecture optimisee d'une station sol secondaire de generation d'un signal sbas secondaire dans un systeme sbas et procede de generation d'un signal sbas secondaire
CN105785363B (zh) * 2014-12-23 2018-02-13 核工业北京地质研究院 一种极化合成孔径雷达图像的span增强方法
CN105372631B (zh) * 2015-10-29 2017-12-05 中国科学院遥感与数字地球研究所 基于遗传‑偏最小二乘算法的极化雷达反演方法及其应用

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101078769A (zh) * 2006-05-25 2007-11-28 中国科学院中国遥感卫星地面站 单次全极化合成孔径雷达图像反演数字高程模型的方法
CN101419284A (zh) * 2008-08-08 2009-04-29 哈尔滨工业大学 由森林覆盖下目标参数反演模型获得人造目标信息的方法
CN101482616A (zh) * 2008-08-13 2009-07-15 中国科学院电子学研究所 一种地形测量方法
CN101685155A (zh) * 2008-09-27 2010-03-31 中国科学院电子学研究所 基于极化干涉合成孔径雷达数据优化干涉相干系数的方法
CN101498789A (zh) * 2009-02-25 2009-08-05 中国测绘科学研究院 一种基于全极化合成孔径雷达的地物目标分类方法和装置
CN101706577A (zh) * 2009-12-01 2010-05-12 中南大学 InSAR监测高速公路路面沉降方法
CN103940834A (zh) * 2014-05-09 2014-07-23 中国科学院电子学研究所 采用合成孔径雷达技术测量土壤湿度的方法

Also Published As

Publication number Publication date
CN108205134A (zh) 2018-06-26

Similar Documents

Publication Publication Date Title
JP5376814B2 (ja) 目標識別装置
CN111766577A (zh) 一种基于三阶段算法p波段的输电线路通道树木高度反演方法
Nagashree et al. Detection and classification of ground penetrating radar image using textrual features
Wang et al. Using multibeam backscatter strength to analyze the distribution of manganese nodules: A case study of seamounts in the Western Pacific Ocean
CN108205134B (zh) 一种极化合成孔径雷达影像的次地表信息增强方法
Zhu et al. Robust moving targets detection and velocity estimation using multi-channel and multi-look SAR images
CN112630741B (zh) 一种全极化合成孔径雷达图像目标补偿peoc方法
CN106680794B (zh) 基于空间建模技术的雷达数据快速去噪方法
CN110297241B (zh) 一种用于图像处理的上下文协方差矩阵构建方法
CN113723483A (zh) 一种基于鲁棒主成分分析的图像融合方法及系统
CN110865391B (zh) 用于目标增强的毫米波太赫兹多极化成像方法及成像系统
JP2010230462A (ja) ポラリメトリｓａｒ装置及びポラリメトリｓａｒデータを用いる目標物同定方法
Pastina et al. Super-resolution of polarimetric SAR images of ship targets
Liu et al. Sparsity-driven distributed array imaging
CN110045374B (zh) 一种基于Chirp Scaling的多角度超长聚束SAR超分辨率成像算法
CN109584255B (zh) 一种基于全色与高光谱图像融合的目标轮廓提取方法
JP3659577B2 (ja) クラッタ抑圧装置およびクラッタ抑圧方法
CN110335214B (zh) 结合上下文协方差矩阵的全极化sar图像相干斑滤波方法
Wang et al. Seafloor classification based on deep-sea multibeam data—Application to the southwest Indian Ridge at 50.47° E
CN105785363A (zh) 一种极化合成孔径雷达图像的span增强方法
Jo et al. Rapid generation of flood maps using dual-polarimetric synthetic aperture radar imagery
Kuhn et al. Development of Methods and Equipment for the Exploration of Manganese Nodules in the German License Area in the Central Equatorial Pacific
Xu et al. Angular response classification of multibeam sonar based on multi-angle interval division
CN117269908A (zh) 基于极化响应及散射特性的极化定标效能评估方法
Nevstad Use of different imaging systems for ROV-based mapping of complex benthic habitats

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