CN109239710A - Acquisition methods and device, the computer readable storage medium of radar elevation information - Google Patents
Acquisition methods and device, the computer readable storage medium of radar elevation information Download PDFInfo
- Publication number
- CN109239710A CN109239710A CN201811009149.0A CN201811009149A CN109239710A CN 109239710 A CN109239710 A CN 109239710A CN 201811009149 A CN201811009149 A CN 201811009149A CN 109239710 A CN109239710 A CN 109239710A
- Authority
- CN
- China
- Prior art keywords
- information
- dimensional coordinate
- obtains
- image points
- elevation
- 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
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/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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/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 embodiment of the invention discloses a kind of acquisition methods of radar elevation information and devices, computer readable storage medium, this method comprises: obtaining the interferometric phase information of the target area according to the interferometer radar image of target area；Using radar stereo measuring technology of image, the initial elevation information of the target area is obtained；Level land is carried out to the interferometric phase information according to the initial elevation information to handle, and obtains differential phase information；Phase unwrapping and elevation inverting are successively carried out to the differential phase information, and in conjunction with the initial elevation information, obtain absolute elevation information.
Description
Technical field
The present invention relates to data processing technique more particularly to the acquisition methods and device, calculating of a kind of radar elevation information
Machine readable storage medium storing program for executing.
Background technique
Digital elevation model (Digital Elevation Model, DEM) is the 3-dimensional digital for describing earth surface shape
Model, by it is a series of include that the data set of geographic plane coordinate and elevation forms, in scientific research, economic construction and military neck
Domain all has an important application value, high-resolution, high-precision especially in the special scenes such as earthquake deformation detection and Geography monitor
The DEM of degree is just more important.
Generally use interference synthetic aperture radar (Interferometric Synthetic Aperture Radar,
InSAR) technology, the technology utilize synthetic aperture radar (Synthetic Aperture Radar, SAR) image pair, obtain compared with
High-precision interferometric phase image calculates elevation information further according to interferometric phase image, and then rebuilds high-precision DEM, wherein
It is the principal element for influencing DEM precision, however the interference fringe packet in interferometric phase image that interferometric phase solution in InSAR technology, which twines,
The level land phase and InSAR technical application included all can cause interference fringe intensive in the biggish area of relief, increase
The difficulty that interferometric phase solution twines, also, the removal of level land phase is easy by satellite orbital error or the shadow of other factors again
It rings, so that the removal effect of level land phase is poor, the difficulty of phase unwrapping is still larger, and then affects the precision of DEM.
Summary of the invention
It is a primary object of the present invention to propose a kind of radar height reconstruction method and device, computer-readable storage medium
Matter, it is intended to solve the problems, such as the level land phase poor removal effect of the height reconstruction method of existing interference synthetic aperture radar.
The technical scheme of the present invention is realized as follows:
The embodiment of the present invention provides a kind of acquisition methods of radar elevation information, which comprises
According to the interferometer radar image of target area, the interferometric phase information of the target area is obtained；
Using radar stereo measuring technology of image, the initial elevation information of the target area is obtained；
Level land is carried out to the interferometric phase information according to the initial elevation information to handle, and obtains differential phase letter
Breath；
Phase unwrapping and elevation inverting successively carried out to the differential phase information, and in conjunction with the initial elevation information,
Obtain absolute elevation information.
Optionally, described to utilize radar stereo measuring technology of image, obtain the initial elevation information of the target area, comprising:
Using the radar stereo measuring technology of image, the corresponding stereogram in the target area is obtained；
According to the corresponding conformation model of the stereogram and the target area, picture of the same name in the stereogram is obtained
The three-dimensional coordinate of point；
Coordinate conversion is carried out to the three-dimensional coordinate of the corresponding image points, obtains the initial elevation information.
Optionally, the three-dimensional coordinate to the corresponding image points carries out coordinate conversion, obtains the initial elevation information,
Include:
Using the method for geocoding, the three-dimensional coordinate of the corresponding image points is converted, obtains the corresponding image points
Conversion after three-dimensional coordinate, the three-dimensional coordinate after the conversion is the coordinate in earth coordinates；
According to the Z axis coordinate value of the three-dimensional coordinate after the conversion, the initial elevation information is obtained.
Optionally, described according to the corresponding conformation model of the stereogram and the target area, obtain the solid
As the three-dimensional coordinate of centering corresponding image points, comprising:
Using Range-Doppler equations, the conformation model is established；
According to the corresponding imaging geometry of the radar stereo measuring technology of image, picture of the same name in the stereogram is determined
The two-dimensional coordinate of point；
According to the two-dimensional coordinate of the conformation model and the corresponding image points, the corresponding three-dimensional seat of the corresponding image points is obtained
Mark.
Optionally, the two-dimensional coordinate according to the conformation model and the corresponding image points, obtains the corresponding image points
Corresponding three-dimensional coordinate, comprising:
Obtain the image grayscale information of the stereogram, be based on described image grayscale information, to the stereogram into
Row window search obtains the parallax of the corresponding image points；
According to the conformation model, the parallax of the two-dimensional coordinate of the corresponding image points and the corresponding image points, obtain described
The corresponding three-dimensional coordinate of corresponding image points.
Optionally, the two-dimensional coordinate according to the conformation model and the corresponding image points, obtains the corresponding image points
Corresponding three-dimensional coordinate, comprising:
It carries out that abscissa makes the difference and ordinate makes the difference respectively to the two-dimensional coordinate of the corresponding image points, obtains the picture of the same name
The parallax of point；
According to the conformation model, the parallax of the two-dimensional coordinate of the corresponding image points and the corresponding image points, obtain described
The corresponding three-dimensional coordinate of corresponding image points.
Optionally, described level land is carried out to the interferometric phase information according to the initial elevation information to handle, it obtains
Differential phase information, comprising:
According to the imaging geometry of the initial elevation information and interference synthetic aperture radar, initial phase letter is obtained
Breath；
Level land is carried out to the interferometric phase information using the original phase information to handle, and obtains differential phase letter
Breath.
Optionally, described that phase unwrapping and elevation inverting are successively carried out to the differential phase information, and in conjunction with described first
Beginning elevation information, obtains absolute elevation information, comprising:
Phase unwrapping is carried out to the differential phase information, obtains the differential phase information after solution twines；
Differential phase information after twining to the solution carries out elevation inverting, obtains difference elevation information；
According to the initial elevation information and the difference elevation information, absolute elevation information is obtained.
The embodiment of the invention provides a kind of acquisition device of radar elevation information, which is characterized in that described device includes:
First information acquiring unit, the second information acquisition unit and information process unit；Wherein
The first information acquiring unit obtains the target area for the interferometer radar image according to target area
Interferometric phase information；
Second information acquisition unit obtains the initial of the target area for utilizing radar stereo measuring technology of image
Elevation information；
The information process unit, for carrying out level land to the interferometric phase information according to the initial elevation information
Processing, obtains differential phase information；Phase unwrapping and elevation inverting are successively carried out to the differential phase information, and in conjunction with described
Initial elevation information, obtains absolute elevation information.
Optionally, second information acquisition unit is specifically used for utilizing the radar stereo measuring technology of image, described in acquisition
The corresponding stereogram in target area；And according to the corresponding conformation model of the stereogram and the target area, obtain institute
State the three-dimensional coordinate of corresponding image points in stereogram；And coordinate conversion is carried out to the three-dimensional coordinate of the corresponding image points, it obtains
The initial elevation information.
Optionally, second information acquisition unit, specifically for the method using geocoding, to the corresponding image points
Three-dimensional coordinate converted, the three-dimensional coordinate after obtaining the conversion of the corresponding image points, the three-dimensional coordinate after the conversion is
Coordinate in earth coordinates；And the Z axis coordinate value according to the three-dimensional coordinate after the conversion, obtain the initial elevation letter
Breath.
Optionally, second information acquisition unit is specifically used for utilizing Range-Doppler equations, establishes the conformation
Model；And according to the corresponding imaging geometry of the radar stereo measuring technology of image, determine picture of the same name in the stereogram
The two-dimensional coordinate of point；And the two-dimensional coordinate according to the conformation model and the corresponding image points, obtain the corresponding image points pair
The three-dimensional coordinate answered.
Optionally, second information acquisition unit, specifically for obtaining the image grayscale information of the stereogram, base
In described image grayscale information, window search is carried out to the stereogram, obtains the parallax of the corresponding image points；And according to institute
The parallax for stating conformation model, the two-dimensional coordinate of the corresponding image points and the corresponding image points, it is corresponding to obtain the corresponding image points
Three-dimensional coordinate.
Optionally, second information acquisition unit carries out respectively specifically for the two-dimensional coordinate to the corresponding image points
Abscissa makes the difference to be made the difference with ordinate, obtains the parallax of the corresponding image points；And according to the conformation model, the corresponding image points
Two-dimensional coordinate and the corresponding image points parallax, obtain the corresponding three-dimensional coordinate of the corresponding image points.
Optionally, the information process unit is specifically used for according to the initial elevation information and interference synthetic aperture thunder
The imaging geometry reached, obtains original phase information；And using the original phase information to the interferometric phase information into
Row goes to level land to handle, and obtains differential phase information.
Optionally, the information process unit is specifically used for carrying out phase unwrapping to the differential phase information, be solved
Differential phase information after twining；And the differential phase information after twining to the solution carries out elevation inverting, obtains difference elevation information；
And according to the initial elevation information and the difference elevation information, obtain absolute elevation information.
The embodiment of the invention provides a kind of acquisition device of radar elevation information, described device includes: processor, storage
Device and communication bus, the memory are communicated by the communication bus with the processor, and the memory stores institute
One or more executable program of processor is stated, when one or more of programs are performed, passes through the processing
Device executes the step of acquisition methods of any one of the above radar elevation information.
The embodiment of the invention provides a kind of computer readable storage medium, the computer-readable recording medium storage has
One or more program, one or more of programs can be executed by one or more processor, to realize above-mentioned
It anticipates the acquisition methods of radar elevation information a kind of.
The embodiment of the invention provides a kind of acquisition methods of radar elevation information and devices, computer-readable storage medium
Matter after obtaining initial elevation information by radar stereo measuring technology of image, utilizes the initial elevation using above-mentioned technic relization scheme
Information carries out level land to interferometric phase information and handles, and treated that phase carries out that solution twines obtains absolute elevation letter to this
Breath, since treated, earth's surface topography phase that phase includes is few, and the solution for the phase that reduces that treated twines difficulty, therefore, mentions
The precision of high absolute elevation information improves.
Detailed description of the invention
Fig. 1 is a kind of flow chart one of the acquisition methods of radar elevation information provided in an embodiment of the present invention；
Fig. 2 is the imaging geometry figure of InSAR system provided in an embodiment of the present invention and target area；
Fig. 3 is the imaging geometry schematic diagram provided in an embodiment of the present invention using StereoSAR technology；
Fig. 4 is a kind of flowchart 2 of the acquisition methods of radar elevation information provided in an embodiment of the present invention；
Fig. 5 is the schematic diagram of the initial elevation information provided in an embodiment of the present invention obtained using StereoSAR technology；
Fig. 6 is obtained to be provided in an embodiment of the present invention using the biradical data of TanDEM-X/TerraSAR-X and InSAR technology
Interferometric phase image；
Fig. 7 is the differential phase information provided in an embodiment of the present invention obtained using StereoSAR technology and InSAR technology
Schematic diagram；
Fig. 8 is the radar image coordinate provided in an embodiment of the present invention obtained using StereoSAR technology and InSAR technology
The schematic diagram of absolute elevation information under system；
Fig. 9 is under the earth coordinates provided in an embodiment of the present invention obtained using StereoSAR technology and InSAR technology
Absolute elevation information schematic diagram；
Figure 10 is the mesh that execution Space Shuttle Radar Topographic Mapping Mission SRTM task provided in an embodiment of the present invention obtains
Mark the schematic diagram of first elevation information in region；
Figure 11 is a kind of structural schematic diagram one of the acquisition device of radar elevation information provided in an embodiment of the present invention；
Figure 12 is a kind of structural schematic diagram two of the acquisition device of radar elevation information provided in an embodiment of the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description.
Embodiment one
The embodiment of the present invention provides a kind of acquisition methods of radar elevation information, as shown in Figure 1, this method comprises:
S101: according to the interferometer radar image of target area, the interferometric phase information of target area is obtained.
In embodiments of the present invention, in order to obtain the elevation information of target area, a kind of work of radar surveying method is former
Reason is to obtain two width SAR images of target area, and establish three-dimensional coordinate model, according to the corresponding three-dimensional seat of two width SAR images
Mark, obtains the elevation information of target area, for example, using the method for synthetic aperture radar measurement in space (StereoSAR) technology；
The working principle of another radar surveying method is that two width for obtaining target area have the SAR image of coherence, and are done
Processing is related to, interferometric phase information is obtained, according to interferometric phase information acquisition elevation information, for example, using the side of InSAR technology
Method；Wherein, the phase that InSAR technology obtains, it is higher than the phase accuracy that StereoSAR technology obtains, therefore, use InSAR
Technology obtains the elevation information of target area, however, when according to interferometric phase information acquisition elevation information, due to interferometric phase
It include level land phase information in information, level land phase is exactly the line caused by Image of Flat Ground in the interference treatment process of image
Property variation phase, its presence increases the difficulty of phase unwrapping, and then affects the precision of the elevation information of acquisition, thus need
Level land is carried out to interferometric phase information to handle, and improve the clean level of removal level land phase as far as possible, to reduce phase solution
The difficulty twined, to improve the precision of elevation information.
Optionally, aiming at the problem that phase unwrapping, can by execute Space Shuttle Radar Topographic Mapping Mission (SRTM,
Shuttle Radar Topography Mission) task, to generate the elevation information of global most area, and utilize
The corresponding phase of elevation information that SRTM is generated is carried out auxiliary phase solution and is twined, but the elevation information that SRTM is generated results from 2000
Year, timeliness is poor, and vertical error is higher, and error amount is approximately equal to 16m；Can also use by TanDEM-X satellite and
The spaceborne double star distributed system of TerraSAR-X satellite composition obtains global Precision Elevation information, but effective lacking
Ground control point in the case where, spaceborne double star distributed system is difficult to obtain elevation information due to lacking external auxiliary data,
Also the level land phase poor removal effect of InSAR can not be solved the problems, such as；In order to obtain, timeliness is good, the lesser elevation letter of error
Breath, the elevation information that can be obtained using StereoSAR technology, according to the corresponding phase of the elevation information of StereoSAR, to dry
It relates to phase information and carries out level land and handle, in this way, the phase information of the smooth horizontal plane of the earth can not only be removed, additionally it is possible to go
Except the phase information that earth surface topography generates, the difficulty of phase unwrapping is thus greatly reduced, improves the elevation understood after twining
The accuracy of information.
Illustratively, the interferometric phase image that target area can be obtained using InSAR technology, be specifically included: using having
Two SAR antennas of interference imaging ability emit microwave to target area, then receive the echo of target area reflection, obtain
Take same target area two width with certain subtense angle have coherence haplopia complex pattern, to two width haplopia complex patterns into
Row interference processing, obtains a width interferometric phase image；Wherein, interference processing is carried out to two width haplopia complex patterns, specifically included: to two
Width haplopia complex pattern is subtracted each other, alternatively, carrying out conjugate multiplication to two width haplopia complex patterns；With any one in target area
For target point P, there are certain imaging geometries by two SAR antennas and target point P, as shown in Fig. 2, the B in Fig. 2 is indicated
Two SAR antenna A1And A2Between parallax range, α indicates that the baseline angle between baseline and reference planes, β indicate A2To A1's
Straight line and A1Angle between the straight line of P, θ indicate radar antenna A1Downwards angle of visibility of the phase center to target point P, R, R ' point
It Biao Shi not radar antenna A1、A2To the oblique distance of target point P, Δ R indicates that the oblique distance between two oblique distances is poor, and H indicates radar platform phase
For the height of the earth reference planes, h indicates height above sea level of the target point P relative to the earth reference planes.
S102: radar stereo measuring technology of image is utilized, the initial elevation information of target area is obtained.
In embodiments of the present invention, it can use radar stereo measuring technology of image, obtain the corresponding stereogram in target area；
According to the corresponding three-dimensional coordinate model of stereogram and target area, the two-dimensional coordinate of the corresponding image points of stereogram is obtained；Root
According to the two-dimensional coordinate and conformation model of corresponding image points, the corresponding three-dimensional coordinate of corresponding image points is obtained；Coordinate is carried out to three-dimensional coordinate
Conversion, obtains initial elevation information.
Illustratively, radar stereo measuring technology of image is that areal is imaged at two different positions, is obtained a pair of
Photo, referred to as stereogram, picture referred to as corresponding image points, the picture point elevation that same object point is constituted in stereogram is to pass through measurement
What the position of corresponding image points and coordinate value obtained；In order to obtain the initial elevation information of target area, it is possible, firstly, to using two
The SAR antenna of different location is respectively imaged target area, or right respectively on different location using a SAR antenna
Target area is imaged, and to obtain the stereogram of target area, stereogram is made of left photo and right photo；Wherein,
The imaging geometry of StereoSAR technology is as shown in figure 3, Fig. 3 expression is with SAR satellite track initiation position on ground
Projected position is the origin of X-coordinate and Y coordinate, with the direction of motion (direction of dotted arrow meaning in Fig. 3) of SAR satellite for X
Axis positive direction, using the vertical line perpendicular to the earth as Z axis, and with X-coordinate origin be directed toward SAR satellite track initiation position direction be
Z axis positive direction, the rectangular coordinate system in space of foundation, the Ps in Fig. 3 indicate SAR antenna instantaneous position, Ps0Indicate that SAR antenna rises
Beginning position, V indicate the flying speed of SAR antenna, and t indicates the flight time of SAR antenna, MslantIndicate slant range resolution, R0Table
Show most short oblique distance when SAR imaging, Rs indicates the oblique distance of a target point on radar illumination ground, R0It is indicated with the relationship of Rs
Are as follows: Rs^2=(R0+Mslant* x) ^2, x indicate the coordinate value of the X-coordinate of target point.
Further, after the stereogram for obtaining target area, with the X-axis and Y-axis in Fig. 3, stereogram is built
Vertical rectangular coordinate system, that is, can be used its abscissa (X-axis) value and ordinate (Y-axis) value to indicate position of each picture point on photo
It sets, obtains the two-dimensional coordinate of each picture point in each pair of corresponding image points, and the Diff E of each pair of corresponding image points is found out with this and is indulged
Coordinate difference, the also known as direction x parallax and the direction y parallax, the two are referred to as the parallax of corresponding image points；On the other hand, it can also obtain
The image grayscale information of stereogram is taken, the image grayscale information is based on, obtains the parallax of each pair of corresponding image points, specifically include:
Some control points are chosen in left photo and right photo, establish centered on control point on left photo and right photo respectively
With window and search window, for being registrated to left photo and right photo；Shiding matching window acquires the correlated measure function of window, into
And offset is obtained, in order to improve search precision, to interpolation is carried out at correlated measure function maxima, obtain more accurate inclined
Shifting amount；Using search window, the search offset of other pixels is obtained according to the corresponding offset in control point, to obtain picture of the same name
The parallax information of point.
Further, after obtaining the parallax of two-dimensional coordinate and corresponding image points of corresponding image points, according to StereoSAR technology
The parallax information for surveying the geometrical relationship of high model, the two-dimensional coordinate of corresponding image points and corresponding image points obtains the three-dimensional of corresponding image points
Coordinate specifically includes: according to the imaging geometry of StereoSAR technology, can use distance-Doppler (Range-
Doppler, RD) equation establishes the conformation model of left picture point and the conformation model of right picture point respectively, each conformation model by
Distance condition equation and Doppler's conditional equation two parts composition, the conformation model of left picture point can indicate are as follows:
Wherein,It is the coordinate value of left picture point imaging moment aerial position,It is left
Three coordinate components of piece picture point imaging moment aerial position,Doppler centroid when being imaged for left picture point, λ table
Show radar wavelength, xLIndicate the coordinate value of X-coordinate of the left picture point in left photo；
The conformation model of right picture point can indicate are as follows:
Wherein,It is the coordinate value of right corresponding image points imaging moment aerial position；It is
Three coordinate components of right corresponding image points imaging moment aerial position；In Doppler when being imaged for right corresponding image points
Frequency of heart, xRIndicate the coordinate value of X-coordinate of the corresponding right picture point of left picture point in right photo；It can be according to corresponding image points
Parallax and SAR flight parameter, obtain With
Numerical value, then by the x of each pair of corresponding image pointsLAnd xRFormula (1) to (4) are substituted into, the three-dimensional coordinate of each pair of corresponding image points is acquired.
Further, after the three-dimensional coordinate for obtaining corresponding image points, it can carry out ground to the three-dimensional coordinate of each corresponding image points
Reason coding, i.e., by the three-dimensional coordinate in the rectangular coordinate system in space of the earth's core, be converted to the address coordinate in earth coordinates, then
According to the Z coordinate value of the three-dimensional coordinate after conversion, the elevation of each corresponding image points is obtained, and then obtains the initial height of target area
Journey information.
S103: level land is carried out to interferometric phase information according to initial elevation information and is handled, differential phase information is obtained.
In embodiments of the present invention, in target area any one target point P absolute phaseIt can indicate are as follows:
Wherein, absolute phaseIt can be interferometric phase,Phase is twined for Baseline solution,For earth surface phase
Position, it is generally the case thatIt is generated, but that there are acquisition times is more early for other elevation informations, and is used by other elevation informations
Time unmatched problem, and the precision of other elevation informations is lower, the embodiment of the present invention is according to initial elevation information and radar
Imaging geometry, obtains original phase information, and initial phase is equal to earth surface phaseCompared to other methods
Obtained level land phase, initial phase further comprise the phase information of ground topography generation；Recycle original phase information to dry
It relates to phase information and carries out level land and handle, differential phase information is obtained, in this way, to the ground topography phase in interferometric phase information
What is removed is more clean, and the solution for thus greatly reducing differential phase twines difficulty.
Illustratively, any one target point in target area can be obtained according to imaging geometry as shown in Figure 2
The initial elevation h of P is indicated are as follows:
Oblique distance difference Δ R is indicated are as follows:
Wherein, Δ φ indicates initial phase, and formula (7) are substituted into formula (6), it is available about initial elevation h and
The relational expression of earth surface phase delta phi can find out initial phase Δ φ, i.e., in the case where known initial elevation h
Earth surface phaseInitial phase is subtracted with the interferometric phase of target point P again, obtains the differential phase of target point P.
S104: phase unwrapping and elevation inverting are successively carried out to differential phase information, and combine initial elevation information, is obtained
Absolute elevation information.
In embodiments of the present invention, since the value range of the interferometric phase obtained according to two width SAR images is (- π, π),
It is lost the phase information of 2K π, needs to add the correct 2K π of itself to interferometric phase by phase unwrapping, is obtained each dry
Relate to the corresponding true phase of phase, wherein K indicates the corresponding lane issue of each true phase, and value is 1 to N, and N is positive
Integer；Similarly, for the differential phase obtained according to interferometric phase, it is also desirable to phase unwrapping is carried out, firstly, to differential phase
Information carries out phase unwrapping, then the differential phase information after twining to solution carries out elevation inverting, obtains difference elevation information, last root
According to initial elevation information and difference elevation information, absolute elevation information is obtained；It should be noted that the absolute elevation information is thunder
The corresponding transverse and longitudinal coordinate value of elevation information under up to image coordinate system, i.e. each absolute elevation information is in radar image (example
Such as, stereogram) in transverse and longitudinal coordinate value, need to absolute elevation information carry out geocoding, obtain geocoding after it is exhausted
To elevation information, the absolute elevation information after geocoding is the elevation information under earth coordinates, i.e., exhausted after geocoding
Transverse and longitudinal coordinate to elevation information is the coordinate indicated with the earth longitude and latitude.
Illustratively, after can be using Chinese remainder theorem CRT, maximal possibility estimation MLE and total variational maximum is based on
The method for testing estimation TV-MAP carries out phase unwrapping to each differential phase information, obtains the differential phase information after solution twines, then
Differential phase after twining to solution carries out elevation inverting, difference elevation is obtained, finally with each difference elevation and corresponding initial height
Cheng Xiangjia obtains absolute elevation.
It should be noted that step S101 to step S104 can be realized by processor.
It is understood that being obtained using synthetic aperture radar (Synthetic Aperture Radar, SAR) image pair
The interferometric phase image for obtaining degree of precision calculates elevation information further according to interferometric phase image, and then rebuilds high-precision DEM,
In, it is the principal element for influencing DEM precision, however the interference item in interferometric phase image that the interferometric phase solution in InSAR technology, which twines,
The level land phase and InSAR technical application that line includes all can cause interference fringe intensive in the biggish area of relief,
Increase the difficulty that interferometric phase solution twines, however, the removal to level land phase is easy by satellite orbital error or other factors again
Influence increase phase unwrapping difficulty so that the removal effect of level land phase is poor, and then affect the precision of DEM.
Embodiment two
In order to more embody the purpose of the present invention, on the basis of the above embodiments, progress is further illustrated
It is bright.
Second embodiment of the present invention provides a kind of acquisition methods of radar elevation information, as shown in figure 4, this method comprises:
S401: according to the interferometer radar image of target area, the interferometric phase information of target area is obtained.
S402: radar stereo measuring technology of image is utilized, the initial elevation information of target area is obtained.
In embodiments of the present invention, using Envisat satellite, it is based on StereoSAR technology, absorbs the solid of target area
As right, the three-dimensional coordinate of corresponding image points in stereogram is calculated, geographical volume is carried out to the three-dimensional coordinate of each corresponding image points
Code, obtains the initial elevation information of target area, the corresponding initial number elevation model (Digital of initial elevation information
Elevation Model, DEM), as shown in figure 5, left side numerical value indicates the pixel number on the direction y, lower section numerical value in Fig. 5
Indicate the pixel number on the direction x, right scale indicates elevation, and unit is m.
S403: according to initial elevation information and radar imagery geometrical relationship, original phase information is obtained.
S404: level land is carried out to interferometric phase information using original phase information and is handled, differential phase information is obtained.
In embodiments of the present invention, it can be produced in conjunction with the biradical data of TanDEM-X/TerraSAR-X using InSAR technology
The interferometric phase information of raw target area, the corresponding interferometric phase image of interferometric phase information are as shown in Figure 6；Interferometric phase is used again
Information subtracts the original phase information obtained using StereoSAR technology, obtains differential phase information, differential phase information pair
The differential phase figure answered is as shown in fig. 7, left side numerical value indicates the pixel number on the direction y, lower section numerical tabular in Fig. 6 and Fig. 7
Show the pixel number on the direction x, right scale indicates phase, and unit is rad；As can be seen that the differential phase figure of Fig. 7, phase
Interferometric phase image compared with Fig. 6, striped is more sparse, i.e., the earth's surface phase for including in differential phase information is less, reduces phase
Position solution twines difficulty.
S405: carrying out phase unwrapping to differential phase information, and the differential phase information after twining to solution carries out elevation inverting,
Obtain difference elevation information.
S406: according to initial elevation information and difference elevation information, absolute elevation information is obtained.
In embodiments of the present invention, initial elevation information is added with differential phase information, obtains absolute elevation letter
Breath, absolute elevation information are the elevation information under radar image coordinate system, the absolute elevation information pair under radar image coordinate system
The absolute DEM answered as shown in figure 8, carrying out geocoding to absolute elevation information again, believe by the absolute elevation after obtaining geocoding
It ceases, the absolute elevation information after geocoding is the elevation information under earth coordinates, the absolute elevation information after geocoding
Corresponding absolute DEM is as shown in Figure 9；Furthermore it is also possible to the first elevation information of target area is obtained by executing SRTM task,
Corresponding first DEM of first elevation information is as shown in Figure 10, and left side numerical value indicates the pixel number on the direction y in Fig. 8-10,
Lower section numerical value indicates the pixel number on the direction x, and right scale indicates elevation, and unit is m；As can be seen that Fig. 8 and Fig. 9 is
It is finer than the texture of Figure 10, i.e., the first DEM corresponding compared to SRTM task, point of the absolute DEM before and after geocoding
Resolution is higher, i.e. the practical elevation information that more meets target area of absolute elevation information before and after geocoding, after geocoding
Absolute elevation information and the error mean of practical elevation information be 1.6586m, error mean square root is 5.06m, lower than SRTM
The corresponding vertical error of task (16m), also, the acquisition time of corresponding first elevation information of SRTM task is too early, timeliness
Difference, in comparison, the acquisition time of the corresponding initial elevation information of StereoSAR technology depend on the SAR image of user's selection
Time, closer to InSAR technology generation radar image time, so that it may avoid time interval too long, topography variation band
The influence come.
It should be noted that step S401 to step S406 can be realized by processor.
Embodiment three
Based on technical concept identical with previous embodiment, the embodiment of the present invention three provides a kind of radar elevation information
Acquisition device 1, as shown in figure 11, the device 1 include: first information acquiring unit 11, the second information acquisition unit 12 and information
Processing unit 13；Wherein
First information acquiring unit 11 obtains the interference of target area for the interferometer radar image according to target area
Phase information；
Second information acquisition unit 12 obtains the initial elevation letter of target area for utilizing radar stereo measuring technology of image
Breath；
Information process unit 13 is handled for carrying out level land to interferometric phase information according to initial elevation information, is obtained
Differential phase information；Phase unwrapping and elevation inverting are successively carried out to differential phase information, and combine initial elevation information, is obtained
Absolute elevation information.
Optionally, the second information acquisition unit 12 is specifically used for utilizing radar stereo measuring technology of image, obtains target area pair
The stereogram answered；And according to the corresponding conformation model of stereogram and target area, obtain corresponding image points in stereogram
Three-dimensional coordinate；And coordinate conversion is carried out to the three-dimensional coordinate of corresponding image points, obtain initial elevation information.
Optionally, the second information acquisition unit 12, specifically for the method using geocoding, to the three-dimensional of corresponding image points
Coordinate is converted, and the three-dimensional coordinate after obtaining the conversion of corresponding image points, the three-dimensional coordinate after conversion is in earth coordinates
Coordinate；And the Z axis coordinate value according to the three-dimensional coordinate after conversion, obtain initial elevation information.
Optionally, the second information acquisition unit 12 is specifically used for utilizing Range-Doppler equations, establishes conformation model；And
According to the corresponding imaging geometry of radar stereo measuring technology of image, the two-dimensional coordinate of corresponding image points in stereogram is determined；With
And the two-dimensional coordinate according to conformation model and corresponding image points, obtain the corresponding three-dimensional coordinate of corresponding image points.
Optionally, the second information acquisition unit 12 is based on image specifically for obtaining the image grayscale information of stereogram
Grayscale information carries out window search to stereogram, obtains the parallax of corresponding image points；And according to conformation model, corresponding image points
The parallax of two-dimensional coordinate and corresponding image points obtains the corresponding three-dimensional coordinate of corresponding image points.
Optionally, the second information acquisition unit 12 carries out abscissa specifically for the two-dimensional coordinate to corresponding image points respectively
It makes the difference and is made the difference with ordinate, obtain the parallax of corresponding image points；And according to conformation model, the two-dimensional coordinate and picture of the same name of corresponding image points
The parallax of point, obtains the corresponding three-dimensional coordinate of corresponding image points.
Optionally, information process unit 13, specifically for according to initial elevation information and interference synthetic aperture radar at
As geometrical relationship, original phase information is obtained；And level land is carried out to interferometric phase information using original phase information and is handled, it obtains
To differential phase information.
Optionally, information process unit 13 are specifically used for carrying out phase unwrapping to differential phase information, obtain after solution twines
Differential phase information；And the differential phase information after twining to solution carries out elevation inverting, obtains difference elevation information；And according to first
Beginning elevation information and difference elevation information, obtain absolute elevation information.
It should be noted that in practical applications, above-mentioned first information acquiring unit 11,12 and of the second information acquisition unit
Information process unit 13 can be realized, specially CPU by the processor 14 on the acquisition device 1 of radar elevation information
(Central Processing Unit, central processing unit), MPU (Microprocessor Unit, microprocessor), DSP
(Digital Signal Processing, digital signal processor) or field programmable gate array (FPGA, Field
Programmable Gate Array) etc. realize.
The embodiment of the invention also provides a kind of acquisition device 1 of radar elevation information, as shown in figure 12, the device 1 packet
Include: processor 14, memory 15 and communication bus 110, memory 15 are communicated by communication bus 110 with processor 14,
One or more executable program of 15 storage processor 14 of memory passes through when one or more program is performed
Processor 14 executes the acquisition methods of the radar elevation information as described in embodiment one.
The embodiment of the invention provides a kind of computer readable storage medium, the computer-readable recording medium storage has
One or more program, one or more of programs can be executed by one or more processor, and described program is located
Manage the acquisition methods that the radar elevation information as described in embodiment one, two is realized when device 14 executes.
It should be understood by those skilled in the art that, the embodiment of the present invention can provide as method, system or computer program
Product.Therefore, the shape of hardware embodiment, software implementation or embodiment combining software and hardware aspects can be used in the present invention
Formula.Moreover, the present invention, which can be used, can use storage in the computer that one or more wherein includes computer usable program code
The form for the computer program product implemented on medium (including but not limited to magnetic disk storage and optical memory etc.).
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
The foregoing is only a preferred embodiment of the present invention, is not intended to limit the scope of the present invention.
Claims (11)
1. a kind of acquisition methods of radar elevation information, which is characterized in that the described method includes:
According to the interferometer radar image of target area, the interferometric phase information of the target area is obtained；
Using radar stereo measuring technology of image, the initial elevation information of the target area is obtained；
Level land is carried out to the interferometric phase information according to the initial elevation information to handle, and obtains differential phase information；
Phase unwrapping and elevation inverting are successively carried out to the differential phase information, and in conjunction with the initial elevation information, obtained
Absolute elevation information.
2. obtaining the mesh the method according to claim 1, wherein described utilize radar stereo measuring technology of image
Mark the initial elevation information in region, comprising:
Using the radar stereo measuring technology of image, the corresponding stereogram in the target area is obtained；
According to the corresponding conformation model of the stereogram and the target area, corresponding image points in the stereogram is obtained
Three-dimensional coordinate；
Coordinate conversion is carried out to the three-dimensional coordinate of the corresponding image points, obtains the initial elevation information.
3. according to the method described in claim 2, it is characterized in that, the three-dimensional coordinate to the corresponding image points carries out coordinate
Conversion, obtains the initial elevation information, comprising:
Using the method for geocoding, the three-dimensional coordinate of the corresponding image points is converted, obtains turning for the corresponding image points
Three-dimensional coordinate after changing, the three-dimensional coordinate after the conversion are the coordinates in earth coordinates；
According to the Z axis coordinate value of the three-dimensional coordinate after the conversion, the initial elevation information is obtained.
4. according to the method described in claim 2, it is characterized in that, described according to the stereogram and the target area pair
The conformation model answered, obtains the three-dimensional coordinate of corresponding image points in the stereogram, comprising:
Using Range-Doppler equations, the conformation model is established；
According to the corresponding imaging geometry of the radar stereo measuring technology of image, corresponding image points in the stereogram is determined
Two-dimensional coordinate；
According to the two-dimensional coordinate of the conformation model and the corresponding image points, the corresponding three-dimensional coordinate of the corresponding image points is obtained.
5. according to the method described in claim 4, it is characterized in that, described according to the conformation model and the corresponding image points
Two-dimensional coordinate obtains the corresponding three-dimensional coordinate of the corresponding image points, comprising:
The image grayscale information of the stereogram is obtained, described image grayscale information is based on, window is carried out to the stereogram
Mouth search, obtains the parallax of the corresponding image points；
According to the conformation model, the parallax of the two-dimensional coordinate of the corresponding image points and the corresponding image points, obtain described of the same name
The corresponding three-dimensional coordinate of picture point.
6. according to the method described in claim 4, it is characterized in that, described according to the conformation model and the corresponding image points
Two-dimensional coordinate obtains the corresponding three-dimensional coordinate of the corresponding image points, comprising:
It carries out that abscissa makes the difference and ordinate makes the difference respectively to the two-dimensional coordinate of the corresponding image points, obtains the corresponding image points
Parallax；
According to the conformation model, the parallax of the two-dimensional coordinate of the corresponding image points and the corresponding image points, obtain described of the same name
The corresponding three-dimensional coordinate of picture point.
7. the method according to claim 1, wherein it is described according to the initial elevation information to the interference phase
Position information carries out level land and handles, and obtains differential phase information, comprising:
According to the imaging geometry of the initial elevation information and interference synthetic aperture radar, original phase information is obtained；
Level land is carried out to the interferometric phase information using the original phase information to handle, and obtains differential phase information.
8. the method according to claim 1, wherein described successively carry out phase solution to the differential phase information
It twines and obtains absolute elevation information and in conjunction with the initial elevation information with elevation inverting, comprising:
Phase unwrapping is carried out to the differential phase information, obtains the differential phase information after solution twines；
Differential phase information after twining to the solution carries out elevation inverting, obtains difference elevation information；
According to the initial elevation information and the difference elevation information, absolute elevation information is obtained.
9. a kind of acquisition device of radar elevation information, which is characterized in that described device includes: first information acquiring unit,
Two information acquisition units and information process unit；Wherein
The first information acquiring unit obtains the dry of the target area for the interferometer radar image according to target area
Relate to phase information；
Second information acquisition unit obtains the initial elevation of the target area for utilizing radar stereo measuring technology of image
Information；
The information process unit, for being carried out at level land according to the initial elevation information to the interferometric phase information
Reason, obtains differential phase information；Phase unwrapping and elevation inverting are successively carried out to the differential phase information, and in conjunction with described first
Beginning elevation information, obtains absolute elevation information.
10. a kind of acquisition device of radar elevation information, which is characterized in that described device includes: processor, memory and communication
Bus, the memory are communicated by the communication bus with the processor, and the memory stores the processor
One or more executable program is executed such as by the processor when one or more of programs are performed
The described in any item methods of claim 1-8.
11. a kind of computer readable storage medium, which is characterized in that the computer-readable recording medium storage have one or
Multiple programs, one or more of programs can be executed by one or more processor, to realize that claim 1-8 such as appoints
Method described in one.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811009149.0A CN109239710B (en) | 2018-08-31 | 2018-08-31 | Method and device for acquiring radar elevation information and computer-readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811009149.0A CN109239710B (en) | 2018-08-31 | 2018-08-31 | Method and device for acquiring radar elevation information and computer-readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109239710A true CN109239710A (en) | 2019-01-18 |
| CN109239710B CN109239710B (en) | 2020-05-08 |
Family
ID=65069055
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811009149.0A Active CN109239710B (en) | 2018-08-31 | 2018-08-31 | Method and device for acquiring radar elevation information and computer-readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109239710B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109839635A (en) * | 2019-03-13 | 2019-06-04 | 武汉大学 | A method of it is extracted by L1b grades of Wave datas of CryoSat-2 SARIn mode and surveys high pin point elevation |
| CN110072189A (en) * | 2019-04-24 | 2019-07-30 | 北京邮电大学 | A kind of determination method and device of floor locating for terminal device |
| WO2020154903A1 (en) * | 2019-01-29 | 2020-08-06 | 华为技术有限公司 | Method and device for determining elevation, and radar |
| CN115128609A (en) * | 2022-09-01 | 2022-09-30 | 中国科学院空天信息创新研究院 | Satellite-borne SAR three-dimensional product generation method and device |
| CN117665809A (en) * | 2023-12-21 | 2024-03-08 | 西南林业大学 | Method for inverting forest canopy height |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101706577A (en) * | 2009-12-01 | 2010-05-12 | 中南大学 | Method for monitoring roadbed subsidence of express way by InSAR |
| US8154435B2 (en) * | 2008-08-22 | 2012-04-10 | Microsoft Corporation | Stability monitoring using synthetic aperture radar |
| CN108445457A (en) * | 2018-02-12 | 2018-08-24 | 中国人民解放军61540部队 | Spaceborne distribution interference synthetic aperture radar threedimensional baseline determines calibration method |
-
2018
- 2018-08-31 CN CN201811009149.0A patent/CN109239710B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8154435B2 (en) * | 2008-08-22 | 2012-04-10 | Microsoft Corporation | Stability monitoring using synthetic aperture radar |
| CN101706577A (en) * | 2009-12-01 | 2010-05-12 | 中南大学 | Method for monitoring roadbed subsidence of express way by InSAR |
| CN108445457A (en) * | 2018-02-12 | 2018-08-24 | 中国人民解放军61540部队 | Spaceborne distribution interference synthetic aperture radar threedimensional baseline determines calibration method |
Non-Patent Citations (1)
| Title |
|---|
| 赵争: "地形复杂区域InSAR高精度DEM提取方法", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020154903A1 (en) * | 2019-01-29 | 2020-08-06 | 华为技术有限公司 | Method and device for determining elevation, and radar |
| CN113302519A (en) * | 2019-01-29 | 2021-08-24 | 华为技术有限公司 | Method and device for determining elevation and radar |
| CN109839635A (en) * | 2019-03-13 | 2019-06-04 | 武汉大学 | A method of it is extracted by L1b grades of Wave datas of CryoSat-2 SARIn mode and surveys high pin point elevation |
| CN109839635B (en) * | 2019-03-13 | 2022-12-27 | 武汉大学 | Method for extracting elevation of height measurement foot points through Cryosat-2 SARIn mode L1 b-level waveform data |
| CN110072189A (en) * | 2019-04-24 | 2019-07-30 | 北京邮电大学 | A kind of determination method and device of floor locating for terminal device |
| CN110072189B (en) * | 2019-04-24 | 2020-05-29 | 北京邮电大学 | Method and device for determining floor where terminal equipment is located |
| CN115128609A (en) * | 2022-09-01 | 2022-09-30 | 中国科学院空天信息创新研究院 | Satellite-borne SAR three-dimensional product generation method and device |
| CN117665809A (en) * | 2023-12-21 | 2024-03-08 | 西南林业大学 | Method for inverting forest canopy height |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109239710B (en) | 2020-05-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109239710A (en) | Acquisition methods and device, the computer readable storage medium of radar elevation information | |
| Massonnet et al. | Radar interferometry: limits and potential | |
| CN107102333B (en) | Satellite-borne InSAR long and short baseline fusion unwrapping method | |
| CN102269813B (en) | Interference processing technology of airborne non-vertical dual-antenna InSAR system | |
| KR100437890B1 (en) | Three dimensional map producing device and method | |
| CN101551455B (en) | 3D terrain imaging system of interferometric synthetic aperture radar and elevation mapping method thereof | |
| CN103885059B (en) | A kind of multi-baseline interference synthetic aperture radar three-dimensional rebuilding method | |
| JP6421395B2 (en) | 3D topographic map formation method from SAR map | |
| CN105677942B (en) | A kind of spaceborne natural scene SAR complex image data rapid simulation method of repeat track | |
| CN106249236B (en) | A kind of spaceborne InSAR long-short baselines image joint method for registering | |
| CN107991676B (en) | Troposphere error correction method of satellite-borne single-navigation-pass InSAR system | |
| JP6320577B2 (en) | Synthetic aperture radar signal processor | |
| CN109031301A (en) | Alpine terrain deformation extracting method based on PSInSAR technology | |
| CN111208512A (en) | Interferometric measurement method based on video synthetic aperture radar | |
| CN109242872A (en) | Interference baseline estimation method based on SRTM DEM | |
| CN107631700B (en) | The three-dimensional vision information method that spatial digitizer is combined with total station | |
| CN109270527B (en) | Circular SAR sub-aperture image sequence combined correlation DEM extraction method | |
| CN108562900B (en) | SAR image geometric registration method based on elevation correction | |
| CN107504919B (en) | Wrapped phase three-dimension digital imaging method and device based on phase mapping | |
| CN109116351B (en) | Spaceborne InSAR positioning and analyzing algorithm | |
| CN109886910A (en) | External digital elevation model DEM modification method and device | |
| CN109946682A (en) | GF3 data baseline estimation method based on ICESat/GLAS | |
| CN108333562A (en) | A kind of dimensionality reduction method for registering images that landform altitude is adaptive | |
| CN110133653B (en) | Satellite-borne SAR image rapid indirect positioning method based on DSM data | |
| CN108957454B (en) | Interference phase simulation method of airborne interference synthetic aperture radar |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |