CN103163521B - InSAR homonymy point automatic optical selection method based on multi-source information correlation - Google Patents

InSAR homonymy point automatic optical selection method based on multi-source information correlation Download PDF

Info

Publication number
CN103163521B
CN103163521B CN201310048502.7A CN201310048502A CN103163521B CN 103163521 B CN103163521 B CN 103163521B CN 201310048502 A CN201310048502 A CN 201310048502A CN 103163521 B CN103163521 B CN 103163521B
Authority
CN
China
Prior art keywords
same place
threshold value
homonymy
insar
points
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
CN201310048502.7A
Other languages
Chinese (zh)
Other versions
CN103163521A (en
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.)
Jigang Defense Technology Co.,Ltd.
Original Assignee
Institute of Electronics of CAS
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 Institute of Electronics of CAS filed Critical Institute of Electronics of CAS
Priority to CN201310048502.7A priority Critical patent/CN103163521B/en
Publication of CN103163521A publication Critical patent/CN103163521A/en
Application granted granted Critical
Publication of CN103163521B publication Critical patent/CN103163521B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)

Abstract

The invention provides an InSAR homonymy point automatic optical selection method based on multi-source information correlation and solves the problem that block adjustment is used for homonymy point selection produced by InSAR. A first step includes calculating correlation coefficients of two scene correlation coefficient pictures corresponding to homonymy points, and removing the group of homonymy points if the correlation coefficient which one of the homonymy points corresponds to is lower than a threshold value; a second step includes calculating the proportion of a bright spot of each homonymy point on a window in the two scene correlation coefficient pictures, and removing the homonymy points, wherein the proportion of the bright spots of the homonymy points is larger than a preset value; a third step includes calculating a gray-level co-occurrence matrix of each homonymy point in the two corresponding scene correlation coefficients, selecting the statistic of correlation, and removing the homonymy points, wherein the correlation of the homonymy points is larger than a threshold value T1; a fourth step includes calculating standard deviation under the window of the two scene DEM and removing the homonymy points, wherein the standard deviation of the homonymy points is larger than a threshold value T2; and a fifth step includes finishing optical selection of the standard deviation if the number of the left homonymy points ranges within a target quantity interval. When the number of the left homonymy points is smaller than a minimum value of the target quantity interval, the threshold value T1 and the threshold value T2 are increased. When the number of the left homonymy points is larger than a maxmum value of the target quantity interval, the threshold value T1 and the threshold value T2 are reduced.

Description

A kind of InSAR same place automatic selection of optimal method based on multi-source information association
Technical field
The invention belongs to InSAR terrain mapping technology, relate to a kind of InSAR same place automatic selection of optimal method based on multi-source information association.
Background technology
When the SAR data scaling of Duo Tiao air strips, large region or splicing, tend to run into the inconsistent problem of three-dimensional coordinate of overlapping region.Can be by reference mark information to image calibration one by one, but need a large amount of reference mark, this surveys district for difficulty, is difficult to realize.By the area adjustment thought in photogrammetric, utilize a small amount of reference mark and same place, can set up the restriction relation between each scape image, thereby solve under the condition of sparse reference mark, combine Scaling Problem.By image processing method, can obtain a large amount of same places in overlapping region, many air strips, because the distortion of SAR oblique distance image geometry is serious, make some same place there will be position deviation, same place extractive technique is many simultaneously removes to evaluate same place quality good or not from image process angle, for InSAR, atural object characteristic and the physical meaning of same place representative is more important, therefore how from a large amount of same places, to filter out the useful same place of InSAR is become to a new problem.
By the geometric properties conversion between same place, can identify the point that error is larger, yet even if the same place position remaining is all correct, can not be directly used in and interferes calibration.In InSAR, the interferometric phase of same place present position will be for calibration, and interferometric phase has strict physical meaning, and it is the phase differential that relevant electromagnetic wave produces at same place place, as very weak to reflection of electromagnetic wave ability in river, lake due to water body, almost can not get correct phase differential.In addition because being actually of receiving of radar is wound around into [π, π] between phase place, must recover true value by phase unwrapping method, when same place is on trees, trees volume scattering effect plays a leading role, the phase place that solution winds out can not accurately represent the height of tree, when same place is on house, due to the distinctive folded effect of covering of SAR, make phase value herein cannot recover true value, when same place is in shadow region, because radar cannot receive any echo, can not get phase place at all.Finally, for interfering calibration, use as far as possible level land point under equal conditions, same place extracts can introduce site error, and due in the larger region of topographic relief, site error has amplification to vertical error, so level land point can effectively suppress this impact.
Above-mentioned these problems, can automatically identify by artificial or image processing method, and for the problem of being combined with real physical characteristics, simple image processing method have been difficult to, and need to find a kind of same place method for optimizing for InSAR for positional fault.
Summary of the invention
The object of this invention is to provide a kind of InSAR same place automatic selection of optimal method based on multi-source information association, to solve area adjustment thought, be applied to the same place On The Choice that InSAR produces.
For achieving the above object, technical solution of the present invention is as follows:
Be somebody's turn to do the InSAR same place automatic selection of optimal method based on multi-source information association, comprise the following steps:
Step 1: calculate the coefficient of coherence of the two scape coefficient of coherence figure that same place is corresponding, if having a corresponding coefficient of coherence of point in same place lower than threshold value, remove this group same place;
Step 2: by remaining same place after step 1 is rejected, calculate the bright spot proportion under its window in two scape map of magnitudes, removal ratio is greater than the same place of preset value;
Step 3: by remaining same place after step 2 is rejected, calculate its gray level co-occurrence matrixes in corresponding two scape magnitude image, select this statistic of correlativity, remove the same place that correlativity is greater than threshold value T1;
Step 4: calculation procedure 3 is rejected same place standard deviation under window in two scape DEM afterwards, removes the same place that standard deviation is greater than threshold value T2;
Step 5: if residue same place quantity falls into destination number interval, same place is preferably complete; When residue same place quantity is less than the minimum value in destination number interval, increase threshold value T1, T2, return to step 3 and continue down to carry out; When residue same place quantity is greater than the maximal value in destination number interval, reduce threshold value T1, T2, return to step 3 and continue down to carry out.
Beneficial effect of the present invention: the present invention utilizes image processing method to distinguish the residing type of ground objects of same place, in conjunction with the Physics eigenvector in InSAR, filter out to meet simultaneously and interfere the needed same place of calibration, effectively suppress to interfere the vertical error of calibration to be introduced.
Embodiment
An InSAR same place automatic selection of optimal method for multi-source information association, embodiment is as follows:
Step 1: the coefficient of coherence that calculates the two scape coefficient of coherence figure that same place is corresponding, if have a corresponding coefficient of coherence of point in same place lower than threshold value, be set to 0.4, remove this group same place, what this step was mainly removed is water body and the same place above shade.
Step 2: by remaining same place after step 1 is rejected, calculate the bright spot proportion under its certain window 101*101 in two scape map of magnitudes, 1/2 the point that definition amplitude is greater than maximum amplitude is bright spot.According to bright spot proportion, if ratio is greater than 0.5, remove this same place.This step is mainly removed near the same place strong scattering points such as house.
Step 3: by remaining same place after step 2 is rejected, calculate its gray level co-occurrence matrixes in corresponding two scape magnitude image, select this statistic of correlativity, remove the same place that correlativity is greater than threshold value T1, initial T1 gets 0.1.
While calculating gray level co-occurrence matrixes, window is set to 101*101, and angle is set to 0 ° (horizontal direction), and pixel distance is set to 1, and grey level quantization grade is 16.14 kinds of characteristic statistics that can obtain commonly using from gray level co-occurrence matrixes, select correlativity as statistic here.
Related definition is: COR = 1 σ x σ y Σ i = 0 L - 1 Σ j = 0 L - 1 ( i - u x ) ( j - u y ) p ij , Wherein u x = Σ i = 0 L - 1 i Σ j = 0 L - 1 p ij , u y = Σ j = 0 L - 1 j Σ i = 0 L - 1 p ij , σ x 2 = Σ i = 0 L - 1 ( i - u x ) 2 Σ j = 0 L - 1 p ij , σ y 2 = Σ j = 0 L - 1 ( j - u y ) 2 Σ i = 0 L - 1 p ij , Wherein L is the gray shade scale after quantizing.
The ASSOCIATE STATISTICS amount of calculating according to the gray level co-occurrence matrixes of each point, compares itself and threshold value T1, is initially set to empirical value.This step is mainly removed the point on trees.
Step 4: calculation procedure 3 is rejected same place afterwards standard deviation under certain window in two scape DEM, removes the same place that standard deviation is greater than threshold value T2.Wherein, the calculation window of standard deviation is identical with the window that calculates gray level co-occurrence matrixes, and threshold value T2 is initially set to empirical value 1.5m.This step is mainly removed the point that the first-class ratio of slope in mountain region is larger.
Step 5: if residue same place quantity falls into destination number interval (getting 40-60), same place is preferably complete; When residue same place quantity is less than the minimum value in destination number interval, increase threshold value T1, T2, return to step 3 and continue down to carry out; When residue same place quantity is greater than the maximal value in destination number interval, reduce threshold value T1, T2, return to step 3 and continue down to carry out.
Wherein the increase of T1, T2 and minimizing are independently, and such as the each change of T2 is set to 0.05m, and T1 is set to 0.01 at every turn, the interval of T1 [0.05,0.15], the interval of T2 is [0,2],, if maximum (little) value goes beyond the scope, with interval maximum (little) value, substitute.

Claims (3)

1. the InSAR same place automatic selection of optimal method based on multi-source information association, the method is first calculated the coefficient of coherence of the two scape coefficient of coherence figure that same place is corresponding, if having a corresponding coefficient of coherence of point in same place lower than threshold value, removes this group same place; It is characterized in that, the method is further comprising the steps of:
Step 1: by remaining same place after rejecting, calculate the bright spot proportion under its window in two scape map of magnitudes, removal ratio is greater than the same place of preset value;
Step 2: by remaining same place after step 1 is rejected, calculate its gray level co-occurrence matrixes in corresponding two scape magnitude image, select this statistic of correlativity, remove the same place that correlativity is greater than threshold value T1;
Step 3: calculation procedure 2 is rejected same place standard deviation under window in two scape DEM afterwards, removes the same place that standard deviation is greater than threshold value T2;
Step 4: if residue same place quantity falls into destination number interval, same place is preferably complete; When residue same place quantity is less than the minimum value in destination number interval, increase threshold value T1, T2, return to step 2 and continue down to carry out; When residue same place quantity is greater than the maximal value in destination number interval, reduce threshold value T1, T2, return to step 2 and continue down to carry out.
2. the InSAR same place automatic selection of optimal method based on multi-source information association as claimed in claim 1, is characterized in that, wherein threshold value T1 initially gets 0.1.
3. the InSAR same place automatic selection of optimal method based on multi-source information association as claimed in claim 1, is characterized in that, wherein threshold value T2 initial setting up 1.5.
CN201310048502.7A 2013-02-06 2013-02-06 InSAR homonymy point automatic optical selection method based on multi-source information correlation Active CN103163521B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310048502.7A CN103163521B (en) 2013-02-06 2013-02-06 InSAR homonymy point automatic optical selection method based on multi-source information correlation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310048502.7A CN103163521B (en) 2013-02-06 2013-02-06 InSAR homonymy point automatic optical selection method based on multi-source information correlation

Publications (2)

Publication Number Publication Date
CN103163521A CN103163521A (en) 2013-06-19
CN103163521B true CN103163521B (en) 2014-02-12

Family

ID=48586765

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310048502.7A Active CN103163521B (en) 2013-02-06 2013-02-06 InSAR homonymy point automatic optical selection method based on multi-source information correlation

Country Status (1)

Country Link
CN (1) CN103163521B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108762688A (en) * 2018-06-06 2018-11-06 北京世纪东方通讯设备有限公司 A kind of media stream storage method and device for video monitoring system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102213762A (en) * 2011-04-12 2011-10-12 中交第二公路勘察设计研究院有限公司 Method for automatically matching multisource space-borne SAR (Synthetic Aperture Radar) images based on RFM (Rational Function Model)
CN102680972A (en) * 2012-06-04 2012-09-19 中国神华能源股份有限公司 Method and device for monitoring surface deformation and data processing equipment
JP2012208034A (en) * 2011-03-30 2012-10-25 Remote Sensing Technology Center Of Japan Radar testing device

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012208034A (en) * 2011-03-30 2012-10-25 Remote Sensing Technology Center Of Japan Radar testing device
CN102213762A (en) * 2011-04-12 2011-10-12 中交第二公路勘察设计研究院有限公司 Method for automatically matching multisource space-borne SAR (Synthetic Aperture Radar) images based on RFM (Rational Function Model)
CN102680972A (en) * 2012-06-04 2012-09-19 中国神华能源股份有限公司 Method and device for monitoring surface deformation and data processing equipment

Non-Patent Citations (15)

* Cited by examiner, † Cited by third party
Title
付琨.基于大尺度双边SIFT的SAR图像同名点自动提取方法.《电子与信息学报》.2012,第34卷(第2期),
向茂生.面向InSAR稀疏控制点测图的同名点提取方法.《电子与信息学报》.2011,第33卷(第12期),
基于大尺度双边SIFT的SAR图像同名点自动提取方法;王山虎; 尤红建; 付琨;《电子与信息学报》;20120229;第34卷(第2期);287-293 *
基于控制点的SAR图像亚像素自动配准;易成龙; 王正明; 郭良益;《计算机仿真》;20100831;第27卷(第8期);243-246 *
姜丽敏
尤红建
易成龙
机载干涉雷达数据联合定位及拼接技术研究;马婧; 尤红建; 胡东辉;《武汉大学学报(信息科学版)》;20120430;第37卷(第4期);410-414 *
王山虎
王正明
胡东辉.机载干涉雷达数据联合定位及拼接技术研究.《武汉大学学报(信息科学版)》.2012,第37卷(第4期),
郭良益.基于控制点的SAR图像亚像素自动配准.《计算机仿真》.2010,第27卷(第8期),
陈曙暄
面向InSAR稀疏控制点测图的同名点提取方法;姜丽敏; 陈曙暄; 向茂生;《电子与信息学报》;20111231;第33卷(第12期);2837-2845 *
马婧

Also Published As

Publication number Publication date
CN103163521A (en) 2013-06-19

Similar Documents

Publication Publication Date Title
Wanders et al. Correction of real-time satellite precipitation with multi-sensor satellite observations of land surface variables
KR101483617B1 (en) Quantitative precipitation estimation system and method thereof
Gallant Adaptive smoothing for noisy DEMs
CN106093890A (en) The method setting up residual clutter suppression radar residual spur
CN104361590A (en) High-resolution remote sensing image registration method with control points distributed in adaptive manner
Meijer et al. Measuring erosion in long-term tillage plots using ground-based lidar
CN105204079B (en) A kind of method using TanDEM-X dual station InSAR extraction Earthquake-landslide volumes
CN113960595A (en) Surface deformation monitoring method and system
CN110780297A (en) Method and system for monitoring soil erosion by adopting satellite-borne synthetic aperture radar
WO2013121340A1 (en) Digital elevation model
CN109186561B (en) Method for estimating volume change of under-ice lake
CN103440420A (en) Agricultural drought disaster-inducing factor risk assessment method
Olika et al. Assessment of soil erosion using RUSLE and GIS techniques: A case of Fincha’a Watershed, Western Ethiopia
Piman et al. Development of a window correlation matching method for improved radar rainfall estimation
CN103163521B (en) InSAR homonymy point automatic optical selection method based on multi-source information correlation
CN110703252A (en) Digital elevation model correction method for interferometric synthetic aperture radar shadow area
CN114202698A (en) Water body accurate extraction method combining active and passive remote sensing data
Gaur et al. Identification of rice crop phenology using Scatsat-1 Ku-band scatterometer in Punjab and Haryana
KR101571438B1 (en) Method of adjusting radar precipitation estimation, recording medium and device for performing the method
CN116029162B (en) Flood disaster inundation range monitoring method and system by using satellite-borne GNSS-R data
CN107067384B (en) A kind of loess plateau terraced fields extracting method based on ray method
Zhang et al. Comparative strengths of SCaMPR satellite QPEs with and without TRMM ingest versus gridded gauge-only analyses
Tian et al. Beijiang water body information extraction based on ENVISAT-ASAR
Trisakti et al. DEM generation from stereo ALOS PRISM and its quality improvement
CN109143188A (en) - 1 data ionosphere corrections method of TOPS sentry

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20201224

Address after: 250101 No.9, Kuangyuan Road, Gongye North Road, Wangsheren street, Licheng District, Jinan City, Shandong Province

Patentee after: Jigang Defense Technology Co.,Ltd.

Address before: 100190 No. 19 West North Fourth Ring Road, Haidian District, Beijing

Patentee before: Research Institute of aerospace information innovation, Chinese Academy of Sciences

Effective date of registration: 20201224

Address after: 100190 No. 19 West North Fourth Ring Road, Haidian District, Beijing

Patentee after: Research Institute of aerospace information innovation, Chinese Academy of Sciences

Address before: 100080, No. 19 West Fourth Ring Road, Beijing, Haidian District

Patentee before: Institute of Electronics, Chinese Academy of Sciences

TR01 Transfer of patent right