CN103675814B - The method of buildings ground level height is determined based on circumference SAR - Google Patents

The method of buildings ground level height is determined based on circumference SAR Download PDF

Info

Publication number
CN103675814B
CN103675814B CN201310445479.5A CN201310445479A CN103675814B CN 103675814 B CN103675814 B CN 103675814B CN 201310445479 A CN201310445479 A CN 201310445479A CN 103675814 B CN103675814 B CN 103675814B
Authority
CN
China
Prior art keywords
gauge point
image
height
buildings
area
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
CN201310445479.5A
Other languages
Chinese (zh)
Other versions
CN103675814A (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.)
Institute of Electronics of CAS
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 CN201310445479.5A priority Critical patent/CN103675814B/en
Publication of CN103675814A publication Critical patent/CN103675814A/en
Application granted granted Critical
Publication of CN103675814B publication Critical patent/CN103675814B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01S13/904—SAR modes
    • G01S13/9088—Circular SAR [CSAR, C-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/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/9021—SAR image post-processing techniques
    • G01S13/9023—SAR image post-processing techniques combined with interferometric 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

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

Abstract

The invention provides a kind of method determining buildings ground level height based on circumference SAR.The method comprises: steps A, and rough estimate goes out target structures thing ground level height H 0; Step B, selects with H 0centered by areas imaging, in this areas imaging, divide multiple height layer, for each height layer, circumference SAR imaging carried out respectively to buildings; Step C, after several height layer circumference SAR imaging in image, selects rescattering bright line to form closed region and a disjoint sub-picture; And step D, rescattering bright line is formed closed region and the height of height layer corresponding to disjoint image is defined as the bottom level of target structures thing.The present invention is based on circumference SAR and determine that the method for buildings ground level height utilizes circumference SAR to the advantage of the omnibearing observation of target, accurately can draw the elevation information on construction bottom ground.

Description

The method of buildings ground level height is determined based on circumference SAR
Technical field
The present invention relates to a kind of buildings ground level method for determining height based on circumference SAR and profile formation method.Wherein, buildings ground level refers to the ground of buildings surrounding.And buildings ground level height refers to the hypo height on buildings surrounding ground.
Background technology
Synthetic-aperture radar (SAR) is widely used at remote sensing fields.In SAR image, due to the oblique distance imaging characteristics of radar sensor, typical urban structure is made to be subject to folded covering, the impact of secondary reflection and shade.Secondary reflection effect is an important feature in buildings SAR image, this kind of secondary reflection shows existence buildings, and the position of secondary reflection effect is just to should the side body of wall of buildings, and secondary reflection effect can as the key character being carried out buildings detection and reconstruction by SAR (1NSAR) signal.But secondary reflection effect depends on height (build higher secondary reflection stronger, vice versa) and the position angle of buildings (building wall and sensor orientation between angle).If building wall is parallel with sensor heading, or angle is less than 10 degree closely parallel, then back scattering extremely strong (rescattering effect the body of wall of direction of building sensor side be almost parallel to SAR orientation in have very strong reflected signal, in addition, in the scope that position angle is less, rescattering intensity decays rapidly, and the place decay larger at position angle is slow), for traditional employing straight line SAR buildings formation method, if the angle between sensor heading and building wall is larger, so rescattering effect is not clearly.Because the buildings position angle in SAR image is all different, therefore, adopt straight line SAR carry out based on the buildings detection technique of secondary reflection effect be have circumscribed.
In circumference SAR situation, all buildingss can by the observation (except the buildings be blocked) of 360 degree, therefore each metope outside buildings can produce stronger rescattering effect, and this is very favourable to the acquisition of contour of building information.Use circumference SAR to carry out imaging to buildings, if imaging plane is located at its construction bottom, the bright line that rescattering effect is formed by polygon wire frame closed for formation one, thus can obtain position and the area information of buildings accurately.But, in actual circumference SAR imaging process, not by auxiliary elevation information (as obtained DEM, DSM by INSAR or Lidar) situation, due to and do not know the exact height of buildings ground level, likely imaging plane is located on other height, thus the bright line that dihedral angle effect produces can offset its actual position.For the metope outside buildings, the directivity of the rescattering that itself and ground are formed is very strong, so when image height is different from construction bottom height, the position of the bright line that rescattering produces can be moved along the direction vertical with metope.In addition, such as, obtain DEM, DSM by INSAR or Lidar, the altitude figures precision that LIDAR obtains is higher, but discrete spot elevation is unfavorable for the limit of accurately locating buildings, and it is higher to obtain cost; INSAR can obtain continuous print ground elevation information, but for buildings, some technical barriers of SAR image self are as shade and foldedly cover problem and not easily solve, thus have impact on INSAR and measure.
Realizing in process of the present invention, applicant finds that the method for prior art determination buildings ground level height needs to be obtained by the mode of supplementary means such as interference SAR (INSAR) or airborne Lidar (LightDetectionAndRanging) elevation information of accurate construction area, stereoscopic measurement method altimetry precision simultaneously based on circumference SAR data is greater than 2m, and precision is also limited.
Summary of the invention
(1) technical matters that will solve
In view of above-mentioned technical matters, the invention provides a kind of method determining buildings ground level height based on circumference SAR.
(2) technical scheme
According to an aspect of the present invention, a kind of method determining buildings ground level height based on circumference SAR is provided.The method comprises: steps A, and rough estimate goes out target structures thing ground level height H 0; Step B, selects with H 0centered by areas imaging, in this areas imaging, divide multiple height layer, for each height layer, circumference SAR imaging carried out respectively to buildings; Step C, after several height layer circumference SAR imaging in image, selects rescattering bright line to form closed region and a disjoint sub-picture; And step D, rescattering bright line is formed closed region and the height of height layer corresponding to disjoint image is defined as the bottom level of target structures thing.
(3) beneficial effect
As can be seen from technique scheme, the present invention is based on circumference SAR and determine that the method for buildings ground level height utilizes circumference SAR to the advantage of the omnibearing observation of target, without the need to supplementary means, whether the rescattering bright line only formed by each exterior wall and ground of judging buildings from circumference SAR image forms closed polygon, just can draw the elevation information on construction bottom ground.
Accompanying drawing explanation
Fig. 1 shows the relation between the position of rescattering bright line in circumference SAR image and imaging plane height;
Fig. 2 a-2c respectively illustrates that imaging plane is consistent with buildings ground level, imaging plane lower than buildings ground level, imaging plane higher than the position reflecting the bright line of buildings two-dimensional silhouette in image during buildings ground level;
Fig. 3 is the process flow diagram of the method determining buildings ground level height according to the embodiment of the present invention based on circumference SAR;
Fig. 4 is for determining the schematic diagram of the first gauge point connected region in the C of method step shown in Fig. 3.
Embodiment
For making the object, technical solutions and advantages of the present invention clearly understand, below in conjunction with specific embodiment, and with reference to accompanying drawing, the present invention is described in more detail.It should be noted that, in accompanying drawing or instructions describe, similar or identical part all uses identical figure number.The implementation not illustrating in accompanying drawing or describe is form known to a person of ordinary skill in the art in art.In addition, although herein can providing package containing the demonstration of the parameter of particular value, should be appreciated that, parameter without the need to definitely equaling corresponding value, but can be similar to corresponding value in acceptable error margin or design constraint.The direction term mentioned in embodiment, such as " on ", D score, "front", "rear", "left", "right" etc., be only the direction with reference to accompanying drawing.Therefore, the direction term of use is used to illustrate and is not used for limiting the scope of the invention.
For in prior art based on the limitation of the imaging technique of straight line SAR, the present invention propose circumference SAR imaging technique can be used to carry out imaging to buildings.In circumference SAR situation, all buildingss can by the observation (except the buildings be blocked) of 360 degree, and therefore each metope of buildings can produce stronger rescattering effect, and this detects buildings and rebuild is very favourable.Use circumference SAR to carry out imaging to buildings, if imaging plane is located at its construction bottom, the bright line that rescattering effect is formed by figure (such as rectangle) closed for formation one, thus can obtain position and the area information of buildings exactly.But, in actual circumference SAR imaging process, and do not know the exact height of buildings ground level, likely imaging plane is located on other height, thus the bright line that dihedral angle effect produces can offset its actual position.Here it should be noted that if so imaging plane height is not right, so likely image can have secondary lobe in all directions if the scattering properties of scatterer is isotropic (as wheel dials ball).For building wall, the directivity of its rescattering effect is very strong, so when image height is different from buildings ground level height, the position of the bright line that dihedral angle effect produces only can be moved along the direction vertical with metope.
Fig. 1 shows the relation between the position of rescattering bright line in circumference SAR image and imaging plane height.In FIG, O represents the position (zero doppler position) of radar detedtor, 0 is the incident angle of radar detedtor, H is the vertical height on radar detedtor distance ground, A point is the metope of buildings and the angle on ground, B point is the metope of imaging plane lower than the buildings detected during physical plane in buildings bottom surface and the angle on ground, and C point is the metope of imaging plane higher than the buildings detected during physical plane in buildings bottom surface and the angle on ground, Δ h 1be imaging plane lower than the difference in height between imaging plane during physical plane and physical plane, Δ h 2be imaging plane higher than the difference in height between imaging plane during physical plane and physical plane, X is the horizontal range between O point and building wall, Δ X 1the horizontal range between A point and B point, Δ X 2it is the horizontal range between A point and C point.
With reference to figure 1, due to | OA|=|OB|=|OC|, can obtain
(H-Δh) 2+(X+ΔX) 2＝H 2+X 2(1)
ΔX = H 2 + X 2 - ( H - Δh ) 2 - X = X 2 + 2 HΔh - Δh 2 - X ≈ H X Δh - - - ( 2 )
When imaging plane is higher than Δ h > 0 (Δ h during the height of physical plane 2), Δ X > 0 (that is, Δ X 2), namely bright line is just along the translation inwardly of metope normal direction, and translational movement is directly proportional to Δ h; When imaging plane is lower than Δ h < 0 (Δ h during the height of physical plane 1), Δ X < 0 (Δ X 1), namely bright line is just along the translation outwardly of metope normal direction.As shown in Figure 2 a, if imaging plane to be located at buildings ground level, because radar detedtor is 360 degree of observations, the bright line produced by secondary reflection is closed, and the rectangular wire frame that can form according to bright line obtains two-dimensional position and the profile information (supposing that buildings is rectangle along the tangent plane of height here) of buildings accurately.When imaging plane is lower than buildings ground level, bright line will towards along the outer translation of metope normal normal direction, and as shown in Figure 2 b, bright line spreads towards surrounding, no longer closed.When imaging plane is higher than buildings ground level, bright line is just along the translation inwardly of metope normal direction, and each bright line can intersect, as shown in Figure 2 c.
According to the secondary reflection characteristic of circumference SAR building wall, in actual imaging process, when the digital elevation model (DEM) built is unknown, whether can close according to building wall secondary reflection bright line and judge that whether image height is consistent with bottom surface, particular location and the floor area of buildings can be determined simultaneously.The range value of secondary reflection bright line is the most remarkable in buildings scattering signatures, and based on this characteristic, embodiments of the present invention provide a kind of circumference SAR buildings formation method.
Fig. 3 is the process flow diagram of the method determining buildings ground level height according to the embodiment of the present invention based on circumference SAR.As shown in Figure 3, according to an embodiment of the invention, provide a kind of buildings ground level method for determining height based on circumferential synthetic aperture radar (SAR), the method comprises:
Steps A, rough estimate goes out target structures thing ground level height H 0;
In this step, rough estimate goes out target structures thing ground level height H 0method well known in the art can be adopted to carry out the bottom level of target structures thing according to a preliminary estimate, such as: utilize the multi-angle observation geometry of circle mark SAR to adopt the bottom level of the method rough estimate target structures thing of measurement in space, that is: 360 ° of Texas tower track annulus are divided into multi-section circular arc, be multiple sub-aperture to every section of arc fragmenting, utilize the similarity between sub-aperture image, by the maximum method of related coefficient, other sub-aperture images are found the pixel mated with center image, obtain the position skew of the two, and then extract the digital elevation DEM of ground scene, i.e. target structures thing ground level height H 0.The concrete steps of the method refer to below with reference to document 1:StephanPalm, H é leneM.Oriot, andHubertM.Cantalloube.RadargrammetricDEMExtractionOverU rbanAreaUsingCircularSARImagery. [D] .IEEETRANSACTIONSONGEOSCIENCEANDREMOTESENSING, vol.50, no.11, pp:935-952, NOV.2012.Give prominence to the key points for making the application and briefly, known method repeated no more.
Step B, selects with described H 0centered by areas imaging, in this areas imaging, divide multiple height layer, for each height layer, circumference SAR imaging carried out to buildings;
In this step, the selection of areas imaging can based on Different factor, such as, and the target structures thing ground level height H that rough estimate goes out 0accuracy.Such as, can select with described H 0centered by highly for rough estimate H 0the scope of 4 to 8 times of the precision of the measurement in space adopted is as areas imaging.Wherein, described measurement in space precision is relevant with registration accuracy between the orientation angles difference being used for extracting between the sub-aperture image of DEM and sub-aperture image.Its measurement in space precision of different-waveband is different, and at X-band, circumference SAR measurement in space precision can reach 1m-2m; At pattern-band, 3m-5m can be reached.
Height can be divided according to certain intervals in the areas imaging selected, thus obtain multiple height layer.Interval delta h between adjacent two height layers can meet wherein, H represents the height between radar detedtor and ground, and X represents the horizontal range between the subpoint of radar detedtor on ground and described target structures thing, and C represents the light velocity, and B represents the transmitted bandwidth of radar detedtor.
Step C, after several height layer circumference SAR imaging in image, selects rescattering bright line to form closed region and a disjoint sub-picture;
Here, to judge in image rescattering bright line whether form closed region and non-intersect be known method.In addition, provide a comparatively typical method herein, the method comprises:
Sub-step C1, obtain the image after the circumference SAR imaging of each height layer weights and, wherein, obtain image weights and step comprise:
Son is C1a step by step, is the first gauge point by energy magnitude in image higher than the pixel definition of predetermined threshold, is the second gauge point by energy magnitude in image lower than the pixel definition of described predetermined threshold;
In this sub-step, the determination of predetermined threshold can adopt multiple method, such as can based on empirical value or representative value.Horizontal CFAR CFAR such as can be adopted to detect and to determine predetermined threshold.CFAR method detects a kind of common method of strong scattering target in SAR image, is widely used in SAR image man-made target context of detection.CFAR detection method is a kind of object detection method of Pixel-level, and its prerequisite is that object pixel gray scale has stronger contrast higher than background relative to background, and it reaches by the comparison of single pixel grey scale and a certain thresholding the object detecting object pixel.When given false alarm rate, predetermined threshold is determined by the statistical property of clutter.
Wherein, the general process of CFAR method is: theoretical according to statistic mixed-state, under given false alarm rate condition, first residing for target, the statistical property of ambient background clutter asks for detection threshold adaptively, then pixel to be detected and adaptive threshold are compared, judge whether it is impact point.By reference to the slip of window, realize the self-adapting detecting to all pixels.The statistical property of target ambient background clutter is determined by the pixel in object pixel surrounding reference window usually.
Son is C1b step by step, selects at least one area-of-interest in the picture, and wherein this area-of-interest comprises described first gauge point of some;
If each interesting image regions number difference (or the number of the second gauge point is different), such weights summation remains significant.This is because, we give weights based on certain rule to the second gauge point in area-of-interest, how many second mark tally objects can not determine the size of weights.
The system of selection of area-of-interest can adopt existing method in prior art, such as marker extraction method, its particular content refers to reference paper 2 (High Resolution SAR Images buildings Study on Extraction Method [D], Hunan, graduate school of the National University of Defense technology, 2009).It will be appreciated by those skilled in the art that this Threshold is exemplary, is not restrictive.
In an embodiment of the invention, additionally provide the system of selection of another kind of area-of-interest, wherein, select area-of-interest to comprise in the picture:
Specifically, in particular implementation, pixel window can be used to scan image.The predefine size of pixel window can be the round values that P × Q, P are selected from that scope is 2-4, and Q is selected from the round values that scope is 2-4.The size of such as pixel window can be such as 2 × 2,2 × 3,2 × 4,3 × 2,3 × 3,3 × 4,4 × 2,4 × 3,4 × 4.(this area size equals pixel window size to suppose to scan certain region by pixel window, pixel window port area can be called), if there is the first gauge point in this region, then this pixel window is respectively along horizontal direction, vertical direction and diagonal move to the pixel window port area adjacent with this region (this pixel window port area and adjacent window area can be limit is adjacent or angle is adjacent and not overlapping), if there is the first gauge point in adjacent pixel window port area, then on the basis of this adjacent pixel window port area, pixel window moves to the adjacent next pixel window port area of the pixel window port area adjacent with this, if there is no the first gauge point in adjacent pixel window port area, then do not carry out above-mentioned mobile process, pixel window moves to the one other pixel window area adjacent with this pixel window port area (such as, pixel window port area can have two neighbor window areas of horizontal direction, two neighbor window areas of vertical direction, and four of diagonal neighbor window areas) and repeat aforesaid operations.
After the pixel having traveled through image according to said process, at least one first gauge point connected region can be obtained, as shown in white portion in Fig. 4.From these connected regions, select number of pixels to be greater than the connected region of N number of (span of N can be 10-20).Then surround these connected regions respectively with the polygon that area is minimum, surrounded region is defined as area-of-interest.
Son is C1c step by step, and rule-based the second gauge point imparting weights given in described area-of-interest, if wherein this rule comprises the first gauge point encirclement described in second gauge point four periderm, then give comparatively high weight to this second gauge point; If second gauge point is surrounded by described first gauge point part, then give lower weights to this second gauge point; If second gauge point only has side to have described first gauge point, then give zero weights to this second gauge point;
In the present embodiment, the second gauge point weighting value is given according to following formula:
In formula 3, T (x i, y j) represent S m(x in region i, y j) weights of the second gauge point of position, usually select a 1=2, a 2=1.
Son is C1d step by step, and the weights for the second gauge point at least one area-of-interest in the image of each height layer are sued for peace, weights and F (h n) be shown below:
F ( h n ) = &Sigma; x i , y j &Element; S m , m = 1 , . . . , M T ( x i , y j ) - - - ( 4 )
Wherein: S mrepresent m the first gauge point connected region, M represents the number of n-th layer image first gauge point connected region.
Sub-step C2, rescattering bright line forms closed region and disjoint image to select weights and maximum image to be judged as.
It should be noted that, this step can also to be come after several height layer circumference SAR imaging in image by the method for Line segment detection, rescattering bright line is selected to form closed region and a disjoint sub-picture, its particular content refers to list of references 3 (F.Tupin, H.Maitre, J.-F.Mangin, J.-M.Nicolas, andE.Pechersky, " DetectionoflinearfeaturesinSARimages:Applicationtoroadne tworkextraction, " IEEETRANSACTIONSONGEOSCIENCEANDREMOTESENSING, , vol.36, no.2, pp.434-453, Mar.1998.).
Step D, forms closed region by the rescattering bright line judged and the height of height layer corresponding to disjoint image is defined as the bottom level of described target structures thing, is namely shown below:
h ^ n = arg { max h n ( F ( h n ) , n = 1,2 , . . . , N ) } - - - ( 5 )
Wherein, represent the height of the buildings ground level finally determined, h nrepresent the height of the n-th height layer, N represents the height number of plies.
So far, by reference to the accompanying drawings the present embodiment has been described in detail.Describe according to above, those skilled in the art should determine that the method for buildings ground level height has had and are clearly familiar with the present invention is based on circumference SAR.
In addition, the above-mentioned definition to each element and method is not limited in various concrete structures, shape or the mode mentioned in embodiment, and those of ordinary skill in the art can replace it with knowing simply.
In sum, the present invention utilizes circumference SAR to the advantage of the omnibearing observation of target, without the need to supplementary means, whether the rescattering bright line only formed by each exterior wall and ground of judging buildings from circumference SAR image forms closed polygon, just can draw the elevation information on construction bottom ground, compared with prior art there is plurality of advantages.
Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only specific embodiments of the invention; be not limited to the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (10)

1. determine a method for buildings ground level height based on circumference SAR, it is characterized in that, comprising:
Steps A, rough estimate goes out target structures thing ground level height H 0;
Step B, selects with described H 0centered by areas imaging, in this areas imaging, divide multiple height layer, for each height layer, circumference SAR imaging carried out respectively to described buildings;
Step C, after several height layer circumference SAR imaging in image, selects rescattering bright line to form closed region and disjoint piece image; And
Step D, forms closed region by rescattering bright line and the height of height layer corresponding to disjoint image is defined as the bottom level of described target structures thing.
2. method according to claim 1, wherein, in described step C, in multiple image, select rescattering bright line to form closed region and disjoint piece image in the following way, comprising:
Sub-step C1, obtain the image after the circumference SAR imaging of each height layer weights and, wherein, obtain an image weights and step comprise:
Son is C1a step by step, is the first gauge point by energy magnitude in image higher than the pixel definition of predetermined threshold, is the second gauge point by energy magnitude in image lower than the pixel definition of described predetermined threshold;
Son is C1b step by step, selects at least one area-of-interest in the picture, and wherein this area-of-interest comprises described first gauge point of some;
Son is C1c step by step, rule-based the second gauge point imparting weights given in described area-of-interest, and wherein this rule comprises:
If the first gauge point surrounds described in second gauge point four periderm, then give comparatively high weight to this second gauge point;
If second gauge point is surrounded by described first gauge point part, then give lower weights to this second gauge point;
If second gauge point only has side to have described first gauge point, then give zero weights to this second gauge point;
Son is C1d step by step, and the weights for the second gauge point at least one area-of-interest in the image of each height layer are sued for peace;
Sub-step C2, rescattering bright line forms closed region and disjoint image to select weights and maximum image to be judged as.
3. method according to claim 2, wherein, described son, step by step in C1a, adopts horizontal CFAR CFAR to detect and determines predetermined threshold.
4. method according to claim 2, wherein, described son, step by step in C1b, selects the step comprising the area-of-interest of the first gauge point of some to comprise in the picture:
The pixel window of predefine size is used to scan to obtain J the first gauge point connected region to the pixel in described image, wherein, this first gauge point connected region refers to that in region, closing in scanning window of each first gauge point exists at least one first gauge point;
From described J the first gauge point connected region, select number of pixels to be greater than K the first gauge point connected region of N, wherein N is selected from the round values that scope is 10-20, K < J;
Described K the first gauge point connected region is surrounded one by one with the polygon that area is minimum;
Surrounded region is defined as described area-of-interest.
5. method according to claim 4, wherein, the predefine size of described pixel window is the round values that P × Q, P are selected from that scope is 2 ~ 4, and Q is selected from the round values that scope is 2 ~ 4.
6. method according to claim 4, is characterized in that:
Described son, step by step in C1d, is sued for peace according to the weights of following formula for the second gauge point at least one area-of-interest in the image of each height layer:
Wherein: S mrepresent m the first gauge point connected region in described area-of-interest, M represents the number of n-th layer image first gauge point connected region, T (x i, y j) represent S m(x in region i, y j) weights of the second gauge point of position;
Described son is step by step in C1c, according to following rule to the second gauge point assignment:
Wherein, described a 1=2, a 2=1.
7. method according to claim 2, wherein, described son, step by step in C1b, adopts marker extraction method to select the area-of-interest of the first gauge point comprising some.
8. method according to any one of claim 1 to 7, wherein, in described steps A, utilizes the multi-angle observation geometry of circle mark SAR to adopt the bottom level H of the method rough estimate target structures thing of measurement in space 0.
9. method according to any one of claim 1 to 7, wherein, in described step B, areas imaging is with described H 0centered by height, the scope of 4 to 8 times of measurement in space precision, wherein, at X-band, described measurement in space precision is between 1m-2m; At pattern-band, described measurement in space precision is between 3m-5m.
10. method according to claim 9, wherein, in described step B, the interval delta h between adjacent two described height layers meets wherein, H represents the height between radar detedtor and ground, and X represents the horizontal range between the subpoint of radar detedtor on ground and described target structures thing, and C represents the light velocity, and B represents the transmitted bandwidth of radar detedtor.
CN201310445479.5A 2013-09-25 2013-09-25 The method of buildings ground level height is determined based on circumference SAR Active CN103675814B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310445479.5A CN103675814B (en) 2013-09-25 2013-09-25 The method of buildings ground level height is determined based on circumference SAR

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310445479.5A CN103675814B (en) 2013-09-25 2013-09-25 The method of buildings ground level height is determined based on circumference SAR

Publications (2)

Publication Number Publication Date
CN103675814A CN103675814A (en) 2014-03-26
CN103675814B true CN103675814B (en) 2016-03-23

Family

ID=50313951

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310445479.5A Active CN103675814B (en) 2013-09-25 2013-09-25 The method of buildings ground level height is determined based on circumference SAR

Country Status (1)

Country Link
CN (1) CN103675814B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104007439B (en) * 2014-05-30 2017-02-08 中国科学院电子学研究所 Interferential circular SAR elevation estimation processing method
CN109583284B (en) * 2017-09-29 2023-09-12 中国科学院空天信息创新研究院 Urban high-rise building height extraction method and device based on high-resolution SAR image
CN110782465B (en) * 2019-12-30 2020-03-27 中智行科技有限公司 Ground segmentation method and device based on laser radar and storage medium
WO2022050026A1 (en) * 2020-09-03 2022-03-10 日本電気株式会社 Measurement data processing device, measurement data processing method, and storage medium
CN112415515B (en) * 2021-01-22 2021-04-13 南京邮电大学 Method for separating targets with different heights by airborne circular track SAR

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101261323A (en) * 2008-04-29 2008-09-10 中国科学院上海光学精密机械研究所 Synthetic aperture laser image-forming radar non-linear chirp phase error overcoming method and device
WO2008125929A2 (en) * 2007-02-14 2008-10-23 Universita' Di Pisa Method for processing multi-pass radar data for sensing and analysing multiple components of non-stationary scatterers
CN101344594A (en) * 2008-05-14 2009-01-14 中国科学院上海光学精密机械研究所 Scanning synthetic pore diameter laser imaging radar
CN102129067A (en) * 2011-03-15 2011-07-20 北京理工大学 Ground moving target signal modeling and imaging method for forward scattering radar (FSR)

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008125929A2 (en) * 2007-02-14 2008-10-23 Universita' Di Pisa Method for processing multi-pass radar data for sensing and analysing multiple components of non-stationary scatterers
CN101261323A (en) * 2008-04-29 2008-09-10 中国科学院上海光学精密机械研究所 Synthetic aperture laser image-forming radar non-linear chirp phase error overcoming method and device
CN101344594A (en) * 2008-05-14 2009-01-14 中国科学院上海光学精密机械研究所 Scanning synthetic pore diameter laser imaging radar
CN102129067A (en) * 2011-03-15 2011-07-20 北京理工大学 Ground moving target signal modeling and imaging method for forward scattering radar (FSR)

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
SAR层析三维成像技术研究;王金峰;《中国硕士学位论文全文数据库》;20101215;正文第2，4，14，37，65，66，85，90，116页 *
单幅高分辨率SAR图像建筑物三维模型重构;傅兴玉;《红外与毫米波学报》;20121231;正文569-576页 *
高分辨率SAR图像建筑物提取方法与研究;赵凌君;《中国博士学位论文全文数据库》;20100415;正文第16，86页 *

Also Published As

Publication number Publication date
CN103675814A (en) 2014-03-26

Similar Documents

Publication Publication Date Title
CN103675814B (en) The method of buildings ground level height is determined based on circumference SAR
CN104361590B (en) High-resolution remote sensing image registration method with control points distributed in adaptive manner
Bürgmann et al. Matching of TerraSAR-X derived ground control points to optical image patches using deep learning
CN103065307B (en) The area merges detection method of SAR/SPOT image under non-precision registration
CN104376330A (en) Polarization SAR image ship target detection method based on superpixel scattering mechanism
CN103578110A (en) Multi-band high-resolution remote sensing image segmentation method based on gray scale co-occurrence matrix
CN104200521A (en) Three-dimensional reconstruction method of building objects based on model prior high-resolution SAR images
CN113239953B (en) SAR image rotation ship detection method based on directed Gaussian function
CN103729853A (en) Three-dimensional GIS assisted high-resolution remote sensing image building collapse-damage detecting method
Shimoni et al. Detection of vehicles in shadow areas using combined hyperspectral and lidar data
CN107607951A (en) A kind of SAR image rescattering characteristic simulation method
CN103268496A (en) Target identification method of SAR (synthetic aperture radar) images
CN116597365A (en) Underground disease object identification method based on neural network
Zhang et al. The characteristics of the multipath scattering and the application for geometry extraction in high-resolution SAR images
CN104166838A (en) Multi-layer CFAR-based layered high-resolution SAR image ship detection method
Liu et al. Building height detection from high-resolution TerraSAR-X imagery and GIS data
CN110276240A (en) A kind of SAR image building wall window information extracting method
Cao et al. Detecting the number of buildings in a single high-resolution SAR image
Tang et al. A knowledge-based 3-D building reconstruction from single very high resolution SAR images
Wegner et al. Building extraction in urban scenes from high-resolution InSAR data and optical imagery
Palm et al. DEM extraction over urban area using circular SAR imagery
Dubois et al. Parameterisation of layover detector for urban analysis in InSAR data
Wegner et al. Fusion of optical and InSAR features for building recognition in urban areas
Liu et al. Building footprint extraction by fusing dual-aspect SAR images
Kodors et al. Simple method of LiDAR point density definition for automatic building recognition

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent for invention or patent application
CB03 Change of inventor or designer information

Inventor after: Tan Weixian

Inventor after: Liu Yan

Inventor after: Lin Bin

Inventor after: Wang Yanping

Inventor after: Hong Wen

Inventor after: Huang Pingping

Inventor before: Tan Weixian

Inventor before: Liu Yan

Inventor before: Lin Bin

Inventor before: Wang Yanping

Inventor before: Hong Wen

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: TAN WEIXIAN LIU YAN LIN YUN WANG YANPING HONG WEN TO: TAN WEIXIAN LIU YAN LIN YUN WANG YANPING HONG WEN HUANG PINGPING

C14 Grant of patent or utility model
GR01 Patent grant
CB03 Change of inventor or designer information

Inventor after: Tan Weixian

Inventor after: Liu Yan

Inventor after: Huang Pingping

Inventor after: Lin Bin

Inventor after: Wang Yanping

Inventor after: Hong Wen

Inventor before: Tan Weixian

Inventor before: Liu Yan

Inventor before: Lin Bin

Inventor before: Wang Yanping

Inventor before: Hong Wen

Inventor before: Huang Pingping

COR Change of bibliographic data