CN103675814A - A circumference-SAR-based method for determining the height of a ground level of a building - Google Patents

A circumference-SAR-based method for determining the height of a ground level of a building Download PDF

Info

Publication number
CN103675814A
CN103675814A CN201310445479.5A CN201310445479A CN103675814A CN 103675814 A CN103675814 A CN 103675814A CN 201310445479 A CN201310445479 A CN 201310445479A CN 103675814 A CN103675814 A CN 103675814A
Authority
CN
China
Prior art keywords
gauge point
height
image
buildings
sar
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.)
Granted
Application number
CN201310445479.5A
Other languages
Chinese (zh)
Other versions
CN103675814B (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

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
    • 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)
  • Radar Systems Or Details Thereof (AREA)
  • Image Processing (AREA)

Abstract

The invention provides a circumference-SAR-based method for determining the height of a ground level of a building. The method comprises the following steps: a step A, a height H0 of a ground level of an object building is roughly estimated; a step B, an imaging scope with the H0 as the center is selected; a plurality of height layers are divided in the imaging scope, and circumference SAR imaging is respectively carried out on the building aiming at each height level; a step C, in a plurality of images after the height level circumference SAR imaging, secondary scattering bright lines are selected to form an enclosed image without intersection; and a step D, the height of the height layer corresponding to the enclosed image without intersection formed by the secondary scattering bright lines are determined to be the bottom height of the object building. According to the circumference-SAR-based method for determining the height of the ground level of the building of the invention, advantages of a circumstance SAR for omnibearing observation of the object are utilized, so that elevation information of the ground at the bottom of the building can be accurately obtained.

Description

Based on circumference SAR, determine the method for buildings ground level height
Technical field
The present invention relates to a kind of buildings ground level method for determining height and profile formation method based on circumference SAR.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 in remote sensing field.In SAR image, due to the oblique distance imaging characteristics of radar sensor, make that typical urban structure is subject to folded covering, the impact of secondary reflection and shade.Secondary reflection effect is an important feature in buildings SAR image, this class secondary reflection shows to exist a buildings, and the position of secondary reflection effect is over against a side wall body that should buildings, and secondary reflection effect can be used as a key character that carries out buildings detection and reconstruction by SAR (1NSAR) signal.Yet secondary reflection effect depends on height (build higher secondary reflection stronger, vice versa) and the position angle of buildings
Figure BDA0000387064530000011
(building wall and sensor orientation between angle).If building wall is parallel with sensor heading, or angle is less than the closely parallel of 10 degree, back scattering extremely strong (rescattering effect the body of wall of direction of building sensor one side be almost parallel to SAR orientation in have very strong reflected signal, in addition, in the less scope in position angle, rescattering intensity decays rapidly, and larger place decay is delayed at position angle), for traditional employing straight line SAR buildings formation method, if the angle between sensor heading and building wall is larger, rescattering effect so is not clearly.Because the buildings position angle on SAR image is all different, therefore, adopt straight line SAR carry out buildings detection technique based on secondary reflection effect be have circumscribed.
In circumference SAR situation, all buildingss can be by the observation of 360 degree (except the buildings being blocked), therefore each metope in buildings outside can produce stronger rescattering effect, and this is very favourable to obtaining of contour of building information.Use circumference SAR to carry out imaging to buildings, if imaging plane is located to its construction bottom, the bright line that rescattering effect forms will form a closed polygon wire frame, thereby can obtain accurately position and the area information of buildings.Yet, in actual circumference SAR imaging process, in not by auxiliary elevation information (as obtained DEM, DSM by INSAR or Lidar) situation, due to and do not know the accurate height of buildings ground level, likely imaging plane is located on other height, thereby the bright line that dihedral angle effect produces can be offset its actual position.For the metope in buildings outside, it is very strong with the directivity of the rescattering that ground forms, so when different from construction bottom height as image height degree, meeting edge, the position direction vertical with metope of the bright line that rescattering produces be moved.In addition, for example, by INSAR or Lidar, obtain DEM, DSM, the altitude figures precision that LIDAR obtains is higher, but discrete spot elevation is unfavorable for accurately locating the limit of buildings, and it is higher to obtain cost; INSAR can obtain continuous ground elevation information, but for buildings, some technical barriers of SAR image self are as shade and fold the difficult solution of the problem of covering, thereby affected INSAR measurement.
In realizing process of the present invention, applicant find prior art determine buildings ground level height method need to by supplementary means for example the mode of interference SAR (INSAR) or airborne Lidar (Light Detection And Ranging) obtain the elevation information of accurate construction area, measurement in space method altimetry precision based on circumference SAR data is greater than 2m simultaneously, 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 of determining buildings ground level height based on circumference SAR.
(2) technical scheme
According to an aspect of the present invention, provide a kind of method of determining buildings ground level height based on circumference SAR.The method comprises: steps A, rough estimate goes out target buildings ground level height H 0; Step B, selects with H 0centered by areas imaging, in this areas imaging, divide a plurality of height layers, for each height layer, buildings is carried out respectively to circumference SAR imaging; Step C in image, selects rescattering bright line to form closed region and a disjoint sub-picture after several height layer circumference SAR imaging; And step D, the height that rescattering bright line is formed to closed region and height layer corresponding to disjoint image is defined as the bottom level of target buildings.
(3) beneficial effect
From technique scheme, can find out, the present invention is based on circumference SAR and determine that the method for buildings ground level height utilizes the advantage of circumference SAR to the omnibearing observation of target, without supplementary means, only from circumference SAR image by judging that whether the rescattering bright line that each exterior wall of buildings and ground form forms closed polygon, just can draw the elevation information on construction bottom ground.
Accompanying drawing explanation
Fig. 1 shows the position of rescattering bright line in circumference SAR image and the relation between imaging plane height;
Fig. 2 a-2c shows respectively that imaging plane is consistent with buildings ground level, imaging plane reflects the position of the bright line of buildings two-dimensional silhouette during higher than buildings ground level in image lower than buildings ground level, imaging plane;
Fig. 3 is for determining the process flow diagram of the method for buildings ground level height based on circumference SAR according to the embodiment of the present invention;
Fig. 4 determines the schematic diagram of the first gauge point connected region in method step C shown in Fig. 3.
Embodiment
For making the object, technical solutions and advantages of the present invention clearer, 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 description, similar or identical part is all used identical figure number.The implementation that does not illustrate in accompanying drawing or describe is form known to a person of ordinary skill in the art in affiliated technical field.In addition, although the demonstration of the parameter that comprises particular value can be provided herein, should be appreciated that, parameter is without definitely equaling corresponding value, but can in acceptable error margin or design constraint, be similar to corresponding value.The direction term of mentioning in embodiment, such as " on ", D score, 'fornt', 'back', " left side ", " right side " etc., be only the direction with reference to accompanying drawing.Therefore, the direction term of use is to be not used for limiting the scope of the invention for explanation.
For the limitation of the imaging technique based on straight line SAR in prior art, the present invention proposes to buildings, to carry out imaging with circumference SAR imaging technique.In circumference SAR situation, all buildingss can be by the observation of 360 degree (except the buildings being blocked), so each metope of buildings can produce stronger rescattering effect, and this detects and rebuild buildings is very favourable.Use circumference SAR to carry out imaging to buildings, if imaging plane is located to its construction bottom, the bright line that rescattering effect forms will form a closed figure (for example rectangle), thereby can obtain exactly position and the area information of buildings.Yet, in actual circumference SAR imaging process, and do not know and likely imaging plane is located at the accurate height of buildings ground level on other height, thereby the bright line that dihedral angle effect produces can be offset its actual position.Here it should be noted that if imaging plane height is not right so, so likely image can have secondary lobe in all directions if the scattering properties of scatterer is isotropic (as wheel is dialled ball).For building wall, the directivity of its rescattering effect is very strong, so as image height degree when different from buildings ground level height, the position of the bright line of dihedral angle effect generation only can be moved along the direction vertical with metope.
Fig. 1 shows the position of rescattering bright line in circumference SAR image and the relation between imaging plane height.In Fig. 1, O represents the position (zero doppler position) of radar detedtor, 0 is the incident angle of radar detedtor, H is that radar detedtor is apart from the vertical height on ground, the metope that A point is buildings and the angle on ground, the metope of buildings and the angle on ground of B point for detecting during lower than physical plane at the imaging plane of buildings bottom surface, C point is the metope of buildings and the angle on ground that imaging plane in buildings bottom surface detects during higher than physical plane, Δ h 1imaging plane difference in height between imaging plane and physical plane during lower than physical plane, Δ h 2be imaging plane difference in height between imaging plane and physical plane during higher than 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 )
As imaging plane Δ h > 0 (Δ h during higher than physical plane height 2), Δ X > 0 (that is, Δ X 2), bright line is just along the translation inwardly of metope normal direction, and translational movement is directly proportional to Δ h; As imaging plane Δ h < 0 (Δ h during lower than physical plane height 1), Δ X < 0 (Δ X 1), bright line is just along the translation outwardly of metope normal direction.As shown in Figure 2 a, if imaging plane is located to buildings ground level, because radar detedtor is 360 degree observations, the bright line being 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 during lower than buildings ground level, bright line will be towards translation outside metope normal normal direction, and as shown in Figure 2 b, bright line spreads towards surrounding, no longer closed.When imaging plane is during 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, in the situation that the digital elevation model (DEM) of building is unknown, can according to building wall secondary reflection bright line, whether closure judge that imaging is highly whether consistent with bottom surface, can determine particular location and the floor area of buildings 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 for determining the process flow diagram of the method for buildings ground level height based on circumference SAR according to the embodiment of the present invention.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 buildings ground level height H 0;
In this step, rough estimate goes out target buildings ground level height H 0can adopt method well known in the art to carry out the bottom level of target buildings according to a preliminary estimate, for example: the multi-angle observation of utilization circle mark SAR adopts the bottom level of the method rough estimate target buildings of measurement in spaces for how much, that is: 360 ° of Texas tower track annulus are divided into multi-section circular arc, to every section of arc fragmenting, be a plurality of sub-apertures, utilize the similarity between sub-subaperture image, by the maximum method of related coefficient, on other sub-subaperture images, find the pixel mating with center image, obtain the position skew of the two, and then the digital elevation DEM of extraction ground scene, it is target buildings ground level height H 0.The concrete steps of the method refer to the Palm below with reference to document 1:Stephan, H é lene M.Oriot, and Hubert M.Cantalloube.Radargrammetric DEM Extraction Over Urban AreaUsing Circular SAR Imagery.[D] .IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, vol.50, no.11, pp:935-952, NOV.2012.For the application being given prominence to the key points and briefly, known method being repeated no more.
Step B, selects with described H 0centered by areas imaging, in this areas imaging, divide a plurality of height layers, for each height layer, buildings is carried out to circumference SAR imaging;
In this step, the selection of areas imaging can be based on different factors, for example, and the target buildings ground level height H that rough estimate goes out 0accuracy.For example, 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 adopting is as areas imaging.Wherein, between described measurement in space precision and the poor and sub-subaperture image of orientation angles that is used for extracting between the sub-subaperture image of DEM, registration accuracy is relevant.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, can reach 3m-5m.
In the areas imaging of selecting, height can be divided according to certain intervals, thereby be obtained a plurality of height layers.Interval delta h between adjacent two height layers can meet wherein, H represents the height between radar detedtor and ground, and X represents that radar detedtor is in the subpoint on ground and the horizontal range between described target buildings, and C represents the light velocity, and B represents the transmitted bandwidth of radar detedtor.
Step C in image, selects rescattering bright line to form closed region and a disjoint sub-picture after several height layer circumference SAR imaging;
Here, judge in image rescattering bright line whether form closed region and non-intersect be known method.In addition, provide a typical method comparatively 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, by energy amplitude in image, higher than the pixel definition of predetermined threshold, is the first gauge point, by energy amplitude in image, lower than the pixel definition of described predetermined threshold, is the second gauge point;
In this sub-step, determining of predetermined threshold can adopt several different methods, for example can be based on empirical value or representative value.For example can adopt horizontal CFAR CFAR to detect and determine predetermined threshold.CFAR method is to detect a kind of common method of strong scattering target on SAR image, in SAR image man-made target context of detection, is widely used.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 and with respect to background, and it is by the object that relatively reaches detection object pixel of single pixel grey scale and a certain thresholding.The in the situation that of given false alarm rate, predetermined threshold is determined by the statistical property of clutter.
Wherein, the general process of CFAR method is: according to statistics etection theory, under given false alarm rate condition, first according to the target statistical property of background clutter around of living in, ask for adaptively detection threshold, 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 target around statistical property of background clutter is determined by the pixel in reference windows around object pixel conventionally.
Son is C1b step by step, selects at least one area-of-interest in image, described the first gauge point that wherein this area-of-interest comprises some;
If each interesting image regions number different (or the number of the second gauge point is different), such weights summation remains significant.This be due to, we are that the rule based on certain is given weights 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, marker extraction method for example, 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, also provide the system of selection of another kind of area-of-interest, wherein, in image, selected area-of-interest to comprise:
Specifically, in particular implementation, can use pixel window to scan image.The predefine size of pixel window can be P * Q, and P is selected from the round values that scope is 2-4, and Q is selected from the round values that scope is 2-4.For example the size of pixel window can be for example 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, can be called pixel window port area), if there is the first gauge point in this region, 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, on the basis of this adjacent pixel window port area, pixel window moves to the next pixel window port area that the pixel window port area adjacent with this is adjacent, if there is no the first gauge point in adjacent pixel window port area, do not carry out above-mentioned mobile processing, pixel window (for example moves to the one other pixel window area adjacent with this pixel window port area, pixel window port area can have two neighbor window areas of horizontal direction, two neighbor window areas of vertical direction, and four neighbor window areas of diagonal) and repeat aforesaid operations.
According to said process, traveling through after the pixel of image, can obtain at least one first gauge point connected region, 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 (span of N can be 10-20).Then with the polygon of area minimum, surround respectively these connected regions, surrounded region is defined as to area-of-interest.
Son is C1c step by step, rule-basedly gives weights to the second gauge point in described area-of-interest, and wherein this rule comprises if a second gauge point surrounding is surrounded by described the first gauge point, gives compared with high weight to this second gauge point; If second gauge point is partly surrounded by described the first gauge point, give lower weights to this second gauge point; If second gauge point only has a side to have described the first gauge point, give zero weights to this second gauge point;
In the present embodiment, according to following formula, give the second gauge point weighting value:
Figure BDA0000387064530000081
In formula 3, T (x i, y j) expression S m(x in region i, y j) weights of the second gauge point of position, conventionally select a 1=2, a 2=1.
Son is C1d step by step, for the weights of the second gauge point at least one area-of-interest in the image of each height layer, sues for peace, and 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 tomographic image the first gauge point connected region.
Sub-step C2, selects weights and maximum image to be judged as rescattering bright line and forms closed region and disjoint image.
It should be noted that, the method that this step can also detect with line segment is come after several height layer circumference SAR imaging in image, select rescattering bright line 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, and E.Pechersky, " Detection of linear features in SAR images:Application to road network extraction, " IEEE TRANSACTIONS ON GEO SCIENCE AND REMOTE SENSING, , vol.36, no.2, pp.434-453, Mar.1998.).
Step D, the height that the rescattering bright line judged is formed to closed region and height layer corresponding to disjoint image is defined as the bottom level of described target buildings, is shown below:
h ^ n = arg { max h n ( F ( h n ) , n = 1,2 , . . . , N ) } - - - ( 5 )
Wherein,
Figure BDA0000387064530000084
the height that represents final definite buildings ground level, h nthe height that represents n height layer, N represents the height number of plies.
So far, by reference to the accompanying drawings the present embodiment be have been described in detail.According to above, describe, those skilled in the art should determine that to the present invention is based on circumference SAR the method for buildings ground level height has had clearly understanding.
In addition, the above-mentioned definition to each element and method is not limited in various concrete structures, shape or the mode of mentioning in embodiment, and those of ordinary skill in the art can know simply and replace it.
In sum, the present invention utilizes the advantage of circumference SAR to the omnibearing observation of target, without supplementary means, only from circumference SAR image by judging whether the rescattering bright line that each exterior wall of buildings and ground form 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; institute is 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 modification of making, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.

Claims (10)

1. based on circumference SAR, determine a method for buildings ground level height, it is characterized in that, comprising:
Steps A, rough estimate goes out target buildings ground level height H 0;
Step B, selects with described H 0centered by areas imaging, in this areas imaging, divide a plurality of height layers, for each height layer, described buildings is carried out respectively to circumference SAR imaging;
Step C in image, selects rescattering bright line to form closed region and a disjoint sub-picture after several height layer circumference SAR imaging; And
Step D, the height that rescattering bright line is formed to closed region and height layer corresponding to disjoint image is defined as the bottom level of described target buildings.
2. method according to claim 1 wherein, in described step C, selects rescattering bright line to form closed region and a disjoint sub-picture in the following way in multiple image, 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, by energy amplitude in image, higher than the pixel definition of predetermined threshold, is the first gauge point, by energy amplitude in image, lower than the pixel definition of described predetermined threshold, is the second gauge point;
Son is C1b step by step, selects at least one area-of-interest in image, described the first gauge point that wherein this area-of-interest comprises some;
Son is C1c step by step, rule-basedly gives weights to the second gauge point in described area-of-interest, and wherein this rule comprises:
If a second gauge point surrounding is surrounded by described the first gauge point, give compared with high weight to this second gauge point;
If second gauge point is partly surrounded by described the first gauge point, give lower weights to this second gauge point;
If second gauge point only has a side to have described the first gauge point, give zero weights to this second gauge point;
Son is C1d step by step, for the weights of the second gauge point at least one area-of-interest in the image of each height layer, sues for peace;
Sub-step C2, selects weights and maximum image to be judged as rescattering bright line and forms closed region and disjoint image.
3. method according to claim 2, wherein, described son in C1a, adopts horizontal CFAR CFAR to detect and determines predetermined threshold step by step.
4. method according to claim 2, wherein, described son is step by step in C1b, and the step of the area-of-interest of the first gauge point of selecting to comprise some in image comprises:
Use the pixel window of predefine size 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 closing on of each the first gauge point in region exists at least one first gauge point in scanning window;
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;
With the polygon of area minimum, surround one by one described K the first gauge point connected region;
Surrounded region is defined as to described area-of-interest.
5. method according to claim 4, wherein, the predefine of described pixel window size is P * Q, and it is 2～4 round values that P is selected from scope, and it is 2～4 round values that Q is selected from scope.
6. method according to claim 4, is characterized in that:
Described son in C1d, is sued for peace for the weights of the second gauge point at least one area-of-interest in the image of each height layer according to following formula step by step:
F ( h n ) = &Sigma; x i , y j &Element; S m , m = 1 , . . . , M T ( x i , y j )
Wherein: S mrepresent m the first gauge point connected region in described area-of-interest, M represents the number of n tomographic image the first gauge point connected region, T (x i, y j) expression 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:
Figure FDA0000387064520000022
Wherein, described a 1=2, a 2=1.
7. method according to claim 2, wherein, described son is step by step in C1b, adopts the area-of-interest of the first gauge point that marker extraction method selects to comprise some.
8. according to the method described in any one in claim 1 to 7, wherein, in described steps A, the multi-angle observation of utilization circle mark SAR adopts the bottom level H of the method rough estimate target buildings of measurement in spaces for how much 0.
9. according to the method described in any one in 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
Figure FDA0000387064520000031
wherein, H represents the height between radar detedtor and ground, and X represents that radar detedtor is in the subpoint on ground and the horizontal range between described target buildings, 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 true CN103675814A (en) 2014-03-26
CN103675814B 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)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104007439A (en) * 2014-05-30 2014-08-27 中国科学院电子学研究所 Interferential circular SAR elevation estimation processing method
CN109583284A (en) * 2017-09-29 2019-04-05 中国科学院遥感与数字地球研究所 Urban skyscraper object height extracting method and device based on High Resolution SAR Images
CN110782465A (en) * 2019-12-30 2020-02-11 中智行科技有限公司 Ground segmentation method and device based on laser radar and storage medium
CN112415515A (en) * 2021-01-22 2021-02-26 南京邮电大学 Method for separating targets with different heights by airborne circular track SAR
WO2022050026A1 (en) * 2020-09-03 2022-03-10 日本電気株式会社 Measurement data processing device, measurement data processing method, and storage medium

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图像建筑物三维模型重构", 《红外与毫米波学报》 *
王金峰: "SAR层析三维成像技术研究", 《中国硕士学位论文全文数据库》 *
赵凌君: "高分辨率SAR图像建筑物提取方法与研究", 《中国博士学位论文全文数据库》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104007439A (en) * 2014-05-30 2014-08-27 中国科学院电子学研究所 Interferential circular SAR elevation estimation processing method
CN109583284A (en) * 2017-09-29 2019-04-05 中国科学院遥感与数字地球研究所 Urban skyscraper object height extracting method and device based on High Resolution SAR Images
CN109583284B (en) * 2017-09-29 2023-09-12 中国科学院空天信息创新研究院 Urban high-rise building height extraction method and device based on high-resolution SAR image
CN110782465A (en) * 2019-12-30 2020-02-11 中智行科技有限公司 Ground segmentation method and device based on laser radar and storage medium
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
JP7444268B2 (en) 2020-09-03 2024-03-06 日本電気株式会社 Measured data processing device, measured data processing method and program
CN112415515A (en) * 2021-01-22 2021-02-26 南京邮电大学 Method for separating targets with different heights by airborne circular track SAR
CN112415515B (en) * 2021-01-22 2021-04-13 南京邮电大学 Method for separating targets with different heights by airborne circular track SAR

Also Published As

Publication number Publication date
CN103675814B (en) 2016-03-23

Similar Documents

Publication Publication Date Title
CN103675814B (en) The method of buildings ground level height is determined based on circumference SAR
Fieber et al. Analysis of full-waveform LiDAR data for classification of an orange orchard scene
CN104200521B (en) High Resolution SAR Images building target three-dimensional rebuilding method based on model priori
CN104361590B (en) High-resolution remote sensing image registration method with control points distributed in adaptive manner
CN103578110B (en) Multiband high-resolution remote sensing image dividing method based on gray level co-occurrence matrixes
CN103729853B (en) High score remote sensing image building under three-dimension GIS auxiliary damages detection method
CN103065307B (en) The area merges detection method of SAR/SPOT image under non-precision registration
Bürgmann et al. Matching of TerraSAR-X derived ground control points to optical image patches using deep learning
CN104463164A (en) Tree canopy structure information extraction method based on rib method and crown height ratio
CN106443593B (en) Sweep the adaptive oil spilling information extracting method of enhancing slowly based on coherent radar
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
CN103400137A (en) Method for extracting geometrical building parameters of synthetic aperture radar (SAR) image
CN116597365A (en) Underground disease object identification method based on neural network
Sun et al. Building height estimation in single SAR image using OSM building footprints
CN105023254A (en) Height reconstruction method of synthetic aperture radar image
Liu et al. Building height detection from high-resolution TerraSAR-X imagery and GIS data
Cao et al. Detecting the number of buildings in a single high-resolution SAR image
CN105069771A (en) Active-contour-based SAR image after-calamity water body information extraction method
Tang et al. A knowledge-based 3-D building reconstruction from single very high resolution SAR images
CN106897985A (en) A kind of multi-angle SAR image fusion method based on visibility classification
Ardekani et al. Recovering the structure of a layered soil, including layer thickness and dielectric permittivity, using the interfaces and objects backscatter detected in GPR B-scans
CN112084943A (en) High-voltage line detection method combining depth semantic segmentation with Hough transform
Wegner et al. Fusion of optical and InSAR features for building recognition in urban areas
Michelin et al. Building edge detection using small-footprint airborne full-waveform lidar data

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 of 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