CN104463895B - 一种基于sar的地表监测图像处理方法 - Google Patents

一种基于sar的地表监测图像处理方法 Download PDF

Info

Publication number
CN104463895B
CN104463895B CN201410826625.3A CN201410826625A CN104463895B CN 104463895 B CN104463895 B CN 104463895B CN 201410826625 A CN201410826625 A CN 201410826625A CN 104463895 B CN104463895 B CN 104463895B
Authority
CN
China
Prior art keywords
mrow
msub
sar
earth
sar image
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
CN201410826625.3A
Other languages
English (en)
Other versions
CN104463895A (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.)
Suzhou Dark Blue Space Remote Sensing Technology Co ltd
Original Assignee
Aoshima Hirokoyasu Information Technology Co Ltd
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 Aoshima Hirokoyasu Information Technology Co Ltd filed Critical Aoshima Hirokoyasu Information Technology Co Ltd
Priority to CN201410826625.3A priority Critical patent/CN104463895B/zh
Publication of CN104463895A publication Critical patent/CN104463895A/zh
Application granted granted Critical
Publication of CN104463895B publication Critical patent/CN104463895B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing
    • G06T2207/10044Radar image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth observation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Image Processing (AREA)

Abstract

一种基于SAR的地表监测图像处理方法,包括分别获取同一时刻和同一地表监测区域的SAR图像,进行边缘轮廓提取,进行滤波预处理,像素处理后得到平均像素,利用最小二乘法选出预处理SAR图像,将预处理SAR图像中的所有像素对应求和后取平均值,作为处理后SAR图像中的对应像素,重构后得到处理后的SAR图像,再获取不同时刻和同一地表监测区域SAR图像,得到所有时刻的经过处理后的SAR图像,本方法处理方式简单,计算量小,匹配对准容易,并且经过处理后的地表监测图像更加精确,对于土地情况的地表监测更加准确。

Description

一种基于SAR的地表监测图像处理方法
技术领域
本发明涉及遥感图像处理领域,更具体地说涉及到基于SAR的土地资源监测图像处理方法。
背景技术
近年来城市发展迅速,深化土地管理制度改革,实行统一规划、统一征地、统一开发、统一出让、统一管理的“五统一”土地使用制度,做到集约用地、节约用地,对于促进城乡经济协调发展,加快城市化进程具有极为重要作用,而对于土地情况的地表监测也变的尤为重要。
星载SAR可以全天候、全天时、大面积获取高分辨率的地标图像数据,并且不受烟尘、云雾限制,因此有着广泛的应用,可应用于军事、农业、导航、地理监视等诸多领域。SAR数据具有固定的重访周期,相对于光学遥感而言,更加适合地表监视和动态的监测。SAR图像的分析和处理要远比一股光学图像的分析和处理困难得多,SAR图像不能像光学图像那样能直观地被理解。作为一种成像雷达,SAR系统工作在微波波段,靠相干成像,因此,其成像机理要比传统的光学遥感成像机制复杂得多,并且SAR图像具有特殊的辐射和几何畸变,其信息形成的机理和信息提取的方法也和传统的光学遥感图像有很大的不同。SAR图像是地面目标平均散射系数的反应,由于种种内部或者外部的原因,使得SAR图像中具有强烈的斑点噪声,图像的失真度较高,因此需要对图像进行处理。然而,目前对于地表监测图像的处理局限于同一时间的单幅SAR图像,或者不同时间的不同图像,处理方式复杂,计算量大,对准困难,精度较低。
发明内容
本发明的目的在于克服现有技术的不足,提供了一种基于SAR的地表监测图像处理方法,其处理方式简单,计算量小,匹配对准容易,并且经过处理后的地表监测图像更加精确,对于土地情况的地表监测更加准确。
基于SAR的地表监测图像处理方法,依次包括以下步骤:
(1)分别获取同一时刻和同一地表监测区域的分辨率为N×M的5幅SAR图像,采用边缘检测算法分别对5幅SAR图像进行边缘轮廓提取,将边缘轮廓提取部分,分别进行滤波预处理;
(2)将分别进行滤波预处理过的边缘轮廓提取部分中的每一个像素与其邻近的8个像素构成3×3的像素块,按照如下公式分别得到处理后的每一个像素:
其中,Pij为对应处理前第i行、第j列的像素值,PIJ为对应经过处理后的第i行、第j列的像素值,其中i=1,2,…,N,j=1,2,…,M;
(3)将与5幅SAR图像对应的经过处理后的第i行、第j列的像素值PIJ进行平均值计算,得到平均像素
(4)利用最小二乘法,分别将与5幅SAR图像对应的经过处理后的第i行、第j列的像素值PIJ与平均像素相减后作平方,求得每一幅SAR图像对应的误差平方和,选出误差平方和最小的2幅作为预处理SAR图像;
(5)将所述预处理SAR图像中的所有像素对应求和后取平均值,作为处理后SAR图像中的对应像素,重构后得到处理后的SAR图像;
(6)分别获取不同时刻和同一地表监测区域的分辨率为N×M的5幅SAR图像,重复步骤(1)-(5)得到所有时刻的经过处理后的SAR图像。
附图说明
图1SAR的地表监测图像处理方法流程图
具体实施方式
基于SAR的地表监测图像处理方法,依次包括以下步骤:
(1)分别获取同一时刻和同一地表监测区域的分辨率为N×M的5幅SAR图像,采用边缘检测算法分别对5幅SAR图像进行边缘轮廓提取,将边缘轮廓提取部分,分别进行平滑滤波预处理,经过预处理后的图像更加清晰,对比度更高;
(2)将分别进行滤波预处理过的边缘轮廓提取部分中的每一个像素与其邻近的8个像素构成3×3的像素块,按照如下公式分别得到处理后的每一个像素:
其中,Pij为对应处理前第i行、第j列的像素值,PIJ为对应经过处理后的第i行、第j列的像素值,其中i=1,2,…,N,j=1,2,…,M;
(3)将与5幅SAR图像对应的经过处理后的第i行、第j列的像素值PIJ进行平均值计算,得到平均像素
(4)利用最小二乘法,分别将与5幅SAR图像对应的经过处理后的第i行、第j列的像素值PIJ与平均像素相减后作平方,求得每一幅SAR图像对应的误差平方和,选出误差平方和最小的2幅作为预处理SAR图像;
(5)将所述预处理SAR图像中的所有像素对应求和后取平均值,作为处理后SAR图像中的对应像素,重构后得到处理后的SAR图像;
(6)分别获取不同时刻和同一地表监测区域的分辨率为N×M的5幅SAR图像,重复步骤(1)-(5)得到所有时刻的经过处理后的SAR图像。
优选地,还包括步骤(7)对所有时刻的经过处理后的SAR图像进行配准。
优选地,所述步骤(7)中对所有时刻的经过处理后的SAR图像进行配准是通过对每一幅处理后的SAR图像进行均匀分块,在每一个均匀分块中提取一个特征窗口,基于第一幅处理后的SAR图像的特征窗口的中心位置,进行配准。
优选地,步骤(7)中的配准采用归一化互相关匹配模板法进行匹配。
优选地,还包括步骤(8),建立云数据库,将步骤(6)中所有经过处理后的SAR图像以及步骤(7)中经过配准的所有SAR图像上传至建立云数据库进行存储,作为备份数据,供远程客户端查询。
优选地,远程客户端为手机、平板电脑、PC和/或其他终端。
优选地,所述边缘检测算法为Canny边缘检测算法或Touzi边缘检测算法。
本方法处理后的SAR图像更加真实的反应了监测地表情况,处理方式简单,计算量小,匹配对准容易,对于土地情况的地表监测更加准确,对于集约用地、节约用地,对于促进城乡经济协调发展,加快城市化进程具有极为重要作用。
尽管为了说明的目的,已描述了本发明的示例性实施方式,但是本领域的技术人员将理解,不脱离所附权利要求中公开的发明的范围和精神的情况下,可以在形式和细节上进行各种修改、添加和替换等的改变,而所有这些改变都应属于本发明所附权利要求的保护范围,并且本发明要求保护的产品各个部门和方法中的各个步骤,可以以任意组合的形式组合在一起。因此,对本发明中所公开的实施方式的描述并非为了限制本发明的范围,而是用于描述本发明。相应地,本发明的范围不受以上实施方式的限制,而是由权利要求或其等同物进行限定。

Claims (6)

1.一种基于SAR的地表监测图像处理方法,其特征在于,依次包括以下步骤:
(1)分别获取同一时刻和同一地表监测区域的分辨率为N×M的5幅SAR图像,采用边缘检测算法分别对5幅SAR图像进行边缘轮廓提取,将边缘轮廓提取部分,分别进行滤波预处理;
(2)将分别进行滤波预处理过的边缘轮廓提取部分中的每一个像素与其邻近的8个像素构成3×3的像素块,按照如下公式分别得到处理后的每一个像素:
<mrow> <msub> <mi>P</mi> <mrow> <mi>I</mi> <mi>J</mi> </mrow> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mfrac> <mrow> <msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <mn>2</mn> <msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mi>j</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> <mn>4</mn> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <mn>2</mn> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mo>+</mo> <msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> <mn>4</mn> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <mn>2</mn> <msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mi>j</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> <mn>4</mn> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <mn>2</mn> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mo>+</mo> <msub> <mi>P</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> <mo>(</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </msub> </mrow> <mn>4</mn> </mfrac> </mrow> <mn>4</mn> </mfrac> <mo>+</mo> <msub> <mi>P</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> </mrow>
其中,Pij为对应处理前第i行、第j列的像素值,PIJ为对应经过处理后的第i行、第j列的像素值,其中i=1,2,…,N,j=1,2,…,M;
(3)将与5幅SAR图像对应的经过处理后的第i行、第j列的像素值PIJ进行平均值计算,得到平均像素
(4)利用最小二乘法,分别将与5幅SAR图像对应的经过处理后的第i行、第j列的像素值PIJ与平均像素相减后作平方,求得每一幅SAR图像对应的误差平方和,选出误差平方和最小的2幅作为预处理SAR图像;
(5)将所述预处理SAR图像中的所有像素对应求和后取平均值,作为处理后SAR图像中的对应像素,重构后得到处理后的SAR图像;
(6)分别获取不同时刻和同一地表监测区域的分辨率为N×M的5幅SAR图像,重复步骤(1)-(5)得到所有时刻的经过处理后的SAR图像。
2.如权利要求1所述的基于SAR的地表监测图像处理方法,其特征在于:还包括步骤(7)对所有时刻的经过处理后的SAR图像进行配准。
3.如权利要求2所述的基于SAR的地表监测图像处理方法,其特征在于:所述步骤(7)中对所有时刻的经过处理后的SAR图像进行配准是通过对每一幅处理后的SAR图像进行均匀分块,在每一个均匀分块中提取一个特征窗口,基于第一幅处理后的SAR图像的特征窗口的中心位置,进行配准。
4.如权利要求3所述的基于SAR的地表监测图像处理方法,其特征在于:所述步骤(7)中的配准采用归一化互相关匹配模板法进行匹配。
5.如权利要求4所述的基于SAR的地表监测图像处理方法,其特征在于:还包括步骤(8),建立云数据库,将步骤(6)中所有经过处理后的SAR图像以及步骤(7)中经过配准的所有SAR图像上传至云数据库进行存储,供远程客户端查询。
6.如权利要求5所述的基于SAR的地表监测图像处理方法,其特征在于:所述边缘检测算法为Canny边缘检测算法或Touzi边缘检测算法。
CN201410826625.3A 2014-12-26 2014-12-26 一种基于sar的地表监测图像处理方法 Active CN104463895B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410826625.3A CN104463895B (zh) 2014-12-26 2014-12-26 一种基于sar的地表监测图像处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410826625.3A CN104463895B (zh) 2014-12-26 2014-12-26 一种基于sar的地表监测图像处理方法

Publications (2)

Publication Number Publication Date
CN104463895A CN104463895A (zh) 2015-03-25
CN104463895B true CN104463895B (zh) 2017-10-24

Family

ID=52909877

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410826625.3A Active CN104463895B (zh) 2014-12-26 2014-12-26 一种基于sar的地表监测图像处理方法

Country Status (1)

Country Link
CN (1) CN104463895B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111223172B (zh) * 2020-01-10 2023-06-30 成都中科合迅科技有限公司 一种雷达探测范围绘制和融合方法、设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101179745A (zh) * 2007-12-05 2008-05-14 宁波大学 一种多视点图像的预处理方法
JP2012185712A (ja) * 2011-03-07 2012-09-27 Mitsubishi Electric Corp 画像照合装置または画像照合方法
CN104185981A (zh) * 2013-10-23 2014-12-03 华为终端有限公司 从连拍图像中选择图像的方法和终端

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101179745A (zh) * 2007-12-05 2008-05-14 宁波大学 一种多视点图像的预处理方法
JP2012185712A (ja) * 2011-03-07 2012-09-27 Mitsubishi Electric Corp 画像照合装置または画像照合方法
CN104185981A (zh) * 2013-10-23 2014-12-03 华为终端有限公司 从连拍图像中选择图像的方法和终端

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
一种基于区域特性与加权平均相结合的遥感图像融合方法;马艳军 等;《中国科技论文在线》;20060626;1-5 *
基于影像匹配技术的地震形变监测研究;董晓燕;《中国优秀硕士学位论文全文数据库 信息科技辑》;20110115;I140-253 *
基于边缘分割的多光谱图像融合方法;兴旺 等;《微计算机信息》;20100515;第26卷(第5-2期);198,204-205 *
多幅非均匀噪声图像的联合去噪法的应用及仿真;杨红丽 等;《现代电子技术》;20061115;121-133 *

Also Published As

Publication number Publication date
CN104463895A (zh) 2015-03-25

Similar Documents

Publication Publication Date Title
Zhao et al. A robust adaptive spatial and temporal image fusion model for complex land surface changes
Wang et al. Fusion of Landsat 8 OLI and Sentinel-2 MSI data
Bosch et al. A multiple view stereo benchmark for satellite imagery
Brigot et al. Adaptation and evaluation of an optical flow method applied to coregistration of forest remote sensing images
Cimoli et al. Application of low-cost UASs and digital photogrammetry for high-resolution snow depth mapping in the Arctic
Huang et al. Feature-based image registration using the shape context
CN102609918A (zh) 基于图像特征配准的航空多光谱遥感影像几何精校正方法
CN108230375B (zh) 基于结构相似性快速鲁棒的可见光图像与sar图像配准方法
Grohmann et al. Dune migration and volume change from airborne LiDAR, terrestrial LiDAR and Structure from Motion-Multi View Stereo
CN107274441B (zh) 一种高光谱图像的波段校准方法和系统
Zhu et al. Robust registration of aerial images and LiDAR data using spatial constraints and Gabor structural features
Pehani et al. Automatic geometric processing for very high resolution optical satellite data based on vector roads and orthophotos
Zhou Urban High-Resolution Remote Sensing: Algorithms and Modeling
CN107341781A (zh) 基于改进相位一致性特征矢量底图匹配的sar影像校正方法
Mozgovoy et al. Automated recognition of vegetation and water bodies on the territory of megacities in satellite images of visible and IR bands
CN106920235A (zh) 基于矢量底图匹配的星载光学遥感影像自动校正方法
Liu et al. Phase correlation pixel‐to‐pixel image co‐registration based on optical flow and median shift propagation
Ehrlich et al. Quantifying the building stock from optical high-resolution satellite imagery for assessing disaster risk
Mukherjee et al. Downscaling MODIS spectral bands using deep learning
Parmehr et al. Automatic registration of optical imagery with 3d lidar data using local combined mutual information
Yoo et al. True orthoimage generation by mutual recovery of occlusion areas
Dai et al. Assessment of karst rocky desertification from the local to regional scale based on unmanned aerial vehicle images: A case‐study of Shilin County, Yunnan Province, China
Sun et al. Identifying terraces in the hilly and gully regions of the Loess Plateau in China
Xiao et al. Detection and tracking method of maritime moving targets based on geosynchronous orbit satellite optical images
Schneider et al. Matching of high-resolution optical data to a shaded DEM

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20170914

Address after: 266555, room 716, Wuyi Mountain Road, Huangdao District, Shandong, Qingdao 697, China

Applicant after: Aoshima Hirokoyasu Information Technology Co.,Ltd.

Address before: 250100, building 11, Qilu cultural and creative base, No. 101 Shun Feng Road, hi tech Zone, Shandong, Ji'nan, A501

Applicant before: Han Zhe

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20190930

Address after: 233010 Incubation Workshop, 414 Changle Road, Yuhui District, Bengbu City, Anhui Province

Patentee after: Bengbu Xingshi Intellectual Property Operations Co.,Ltd.

Address before: 697 room 716, Wuyi Mountain Road, Huangdao District, Shandong, Qingdao, 266555

Patentee before: Aoshima Hirokoyasu Information Technology Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220727

Address after: 215500 3rd floor, R & D building, No.95, Liuzhou Road, Changfu street, Changshu City, Suzhou City, Jiangsu Province

Patentee after: Suzhou dark blue space remote sensing technology Co.,Ltd.

Address before: 233010 1st floor, incubation workshop, Industrial Automation Research Institute, 414 Changle Road, Yuhui District, Bengbu City, Anhui Province

Patentee before: Bengbu Xingshi Intellectual Property Operations Co.,Ltd.

TR01 Transfer of patent right