CN110068823A - A kind of block parallel quality guidance fast phase solution twines algorithm - Google Patents
A kind of block parallel quality guidance fast phase solution twines algorithm Download PDFInfo
- Publication number
- CN110068823A CN110068823A CN201910263071.3A CN201910263071A CN110068823A CN 110068823 A CN110068823 A CN 110068823A CN 201910263071 A CN201910263071 A CN 201910263071A CN 110068823 A CN110068823 A CN 110068823A
- Authority
- CN
- China
- Prior art keywords
- phase
- block
- point
- pixel
- writhing number
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
-
- 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
Abstract
The invention discloses a kind of block parallel quality guidance fast phase solutions to twine algorithm, the specific steps are that: obtain original interference phase diagram data, original interference phase diagram is decomposed into several phase blocks according to preset phase block piecemeal parameter, the local writhing number of the pixel in each phase block is solved using quality bootstrap algorithm；The local writhing number gradient between adjacent phase block is calculated, number of hops and the corresponding mass value of each phase block between phase block is calculated according to the local writhing number between all adjacent phase blocks, the block writhing number of each phase block is solved using quality bootstrap algorithm；The corresponding solution of the pixel is obtained according to the corresponding local writhing number of each pixel, block writhing number and winding phase and twines phase, to realize that fast and accurately obtaining the corresponding solution of original interference phase diagram all pixels point twines phase.
Description
Technical field
The invention belongs to interference signal process fields, and in particular to a kind of block parallel quality guidance fast phase solution twines calculation
Method.
Background technique
Synthetic aperture radar interferometry technology by synthetic aperture radar image-forming technology and interferometry technology successfully into
Having gone combination can be accurate using the geometrical relationship between sensor height, radar wavelength, wave beam visual direction and antenna baseline distance
Three positions for measuring every bit on image and change information.
Interference synthetic aperture sonar is because secondly peacekeeping three-dimensional high definition characteristic, in bathymetric surveying, underwater archaeology, seabed
Object is lost in oil exploration and search under water, and underwater Small object identification aspect (such as submarine mine) plays a significant role.In three-dimensional figure
As in reconstruction process, phase unwrapping is step of crucial importance in signal processing, and solution twines performance and is directly related to interference
The service performance of synthetic aperture sonar system.
As application demand is constantly promoted, resolution requirement is higher and higher, in the unit time not for the data volume of processing
It is disconnected to increase, so that conventional phase solution twines algorithm and is difficult to meet real-time demand.It is existing one that quality, which guides phase unwrapping algorithm,
The kind stronger phase unwrapping method of robustness, solution twine the reliability that precision depends entirely on phase masses, however due to quality
Guidance phase unwrapping algorithm is serial computing mode, and for phase interference figure biggish for data volume, quality guides phase solution
It twines algorithm seriously and limits the promotion of phase unwrapping efficiency, simultaneously as there are solutions to twine error propagation in low quality region, thus,
Quality guides the accuracy of phase unwrapping algorithm not high.
Summary of the invention
Aiming at the above defects or improvement requirements of the prior art, the present invention provides a kind of guidance of block parallel quality is quick
Phase unwrapping algorithm, by after original interference phase diagram piecemeal, will calculate the local writhing number of pixel in single phase block,
The block writhing number of single phase block is somebody's turn to do according to the corresponding local writhing number of each pixel, block writhing number and winding phase
The corresponding solution of pixel twines phase, to realize that fast and accurately obtaining the corresponding solution of original interference phase diagram all pixels point twines
Phase.
To achieve the above object, according to one aspect of the present invention, a kind of quick phase of block parallel quality guidance is provided
Position solution twines algorithm, specific steps are as follows:
S1. original interference phase diagram data is obtained, divides original interference phase diagram according to preset phase block piecemeal parameter
Solution is several phase blocks, and the local writhing number of the pixel in each phase block is solved using quality bootstrap algorithm；
S2. the local writhing number gradient between adjacent phase block is calculated, is twined according to the part between all adjacent phase blocks
Winding number gradient calculates number of hops and the corresponding mass value of each phase block between adjacent phase block, is asked using quality bootstrap algorithm
Solve the block writhing number of each phase block；
S3. it is corresponding the pixel to be obtained according to the corresponding local writhing number of each pixel, block writhing number and winding phase
Solution twine phase, to realize that fast and accurately obtaining the corresponding solution of original interference phase diagram all pixels point twines phase.
As a further improvement of the present invention, the part of the pixel in each phase block is solved using quality bootstrap algorithm
Writhing number specifically:
The corresponding phase masses of each pixel in single phase block are calculated using original interference phase diagram data matrix
Value；
It is initial calculation point that the maximum pixel of phase masses value in phase block, which is arranged, and part corresponding to the pixel twines
Winding number is 0；
Selecting the maximum pixel of phase masses value in a calculating vertex neighborhood is current calculating point, and the current point that calculates corresponds to
Local writhing number be that the upper one local writhing number for calculating point and upper one calculate point and current the sum of the number of hops calculated between point.
As a further improvement of the present invention, the corresponding phase masses value of each pixel can utilize puppet in single phase block
Correlation coefficient process, phase gradient variance method or maximum phase gradient method are calculated.
As a further improvement of the present invention, the calculation method of the current number of hops for calculating point are as follows: relatively currently calculate point
Winding phase and upper one calculate point winding phase difference, difference is greater than the first preset value, then it is current calculate it is corresponding
Number of hops is -1；Difference is less than the second preset value, then the current corresponding number of hops of point that calculates is 1；Otherwise, the current point that calculates corresponds to
Number of hops be 0.
As a further improvement of the present invention, the jump number calculating method between adjacent phase block are as follows:
Number of corresponding pixels is N, the phase of one of phase block at corresponding adjacent boundary between adjacent phase block
The corresponding local writhing number Cr of adjacent side circle ith pixel pointiIt indicates, the adjacent boundary ith pixel point of another phase block is corresponding
Local writhing number CliIt indicates, then the number of hops between adjacent phase block isWherein wiFor normalization
Weight afterwards.
As a further improvement of the present invention, the weight after normalization is wi=1/N.
As a further improvement of the present invention, the calculation method of the number of hops between adjacent phase block are as follows: count adjacent phase
The difference Cl of local writhing number corresponding to the adjacent boundary ith pixel point of position blocki-Cri, the most difference of the frequency of occurrences is made
For the number of hops between adjacent phase block.
As a further improvement of the present invention, the corresponding mass value of phase block is the side of phase block phase block adjacent thereto
The variance of the difference of local writhing number corresponding to pixel at boundary determines.
As a further improvement of the present invention, the block writhing number for solving each phase block using quality bootstrap algorithm is specific
Are as follows: it is initial calculation point that the corresponding the smallest phase block of mass value, which is arranged, and block writhing number corresponding to initial calculation point is 0, choosing
Selecting the smallest phase block of phase masses value in a calculating vertex neighborhood is current calculating point, current to calculate corresponding piece of writhing number of point
The local writhing number and upper one for calculating point for upper one calculate point and current the sum of the number of hops calculated between point.
As a further improvement of the present invention, the corresponding solution of the i-th row jth column pixel of original interference phase diagram twines phase
Are as follows:
Wherein,For the corresponding local writhing number of the i-th row jth column pixel,For phase where the i-th row jth column pixel
Position corresponding piece of writhing number of block.
In general, through the invention it is contemplated above technical scheme is compared with the prior art, have below beneficial to effect
Fruit:
A kind of block parallel quality guidance fast phase solution of the invention twines algorithm, and original interference phase diagram is carried out piecemeal
Afterwards, parallel processing, and then the computational efficiency of boosting algorithm can be carried out to multiple phase blocks simultaneously, since out of phase block is using only
Vertical method is handled, and the higher method of reliability can be selected in low quality region, and the solution that can effectively avoid low quality region twines mistake
Difference is propagated, thus the accuracy of boosting algorithm.
A kind of block parallel quality guidance fast phase solution of the invention twines algorithm, by the corresponding local writhing number of phase block
Merging process be changed into block writhing number determined according to local writhing number, and then extrapolate the process that solution twines phase, wherein part
The local writhing number that facet phases need to be only utilized when writhing number merges, greatly simplifies the acquisition process for understanding and twining phase, thus fastly
Solution is obtained fastly twine phase and promote the solution of phase interference figure and twine efficiency, accelerate interference signal processing speed, so that real-time high quality
Digital elevation model reconstruction becomes possible.
Detailed description of the invention
Fig. 1 is that a kind of block parallel quality guidance fast phase solution of the embodiment of the present invention twines the schematic diagram of algorithm；
Fig. 2 is the schematic diagram of the emulation interferometric phase image of the embodiment of the present invention；
Fig. 3 corresponding writhing number schematic diagram when being the non-piecemeal of emulation interferometric phase image of the embodiment of the present invention；
The solution that Fig. 4 is solved when being the non-piecemeal of emulation interferometric phase image of the embodiment of the present invention twines the schematic diagram of phase；
Fig. 5, which is the emulation interferometric phase image block count of the embodiment of the present invention, solves the local writhing number of pixel when being 4 × 4
Schematic diagram；
Fig. 6, which is the emulation interferometric phase image block count of the embodiment of the present invention, solves the block writhing number of block of pixels when being 4 × 4
Schematic diagram；
Fig. 7, which is the emulation interferometric phase image block count of the embodiment of the present invention, to be solved the writhing number of pixel when being 4 × 4 and shows
It is intended to.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
In addition, as long as technical characteristic involved in the various embodiments of the present invention described below is each other not
Constituting conflict can be combined with each other.The present invention is described in more detail With reference to embodiment.
Fig. 1 is that a kind of block parallel quality guidance fast phase solution of the embodiment of the present invention twines the schematic diagram of algorithm.Such as Fig. 1
It is shown, the specific steps of the algorithm are as follows:
S1. original interference phase diagram data is obtained, divides original interference phase diagram according to preset phase block piecemeal parameter
Solution is several phase blocks, and the local writhing number of the pixel in each phase block is solved using quality bootstrap algorithm；
Specifically, original interference phase diagram data matrix is obtainedIts ranks number is respectively i=NaAnd j=Nr, NaAnd Nr
Respectively natural number,The corresponding winding phase of the i-th row jth column pixel of original interference phase diagram is represented, winding phase is
The phase value that original non-solution twines.
Original interference phase diagram is decomposed into several phase blocks, wherein preset phase block piecemeal parameter is single phase
Block corresponding ranks number in position is respectively BaAnd Br, calculated using original interference phase diagram data matrix each in single phase block
The corresponding phase masses value M of pixelij, can use spurious correlation Y-factor method Y, phase gradient variance method or maximum phase gradient method
It carries out the corresponding phase masses value of pixel to calculate, it is initial calculation that the maximum pixel of phase masses value in phase block, which is arranged,
Point, local writhing number corresponding to the pixel are 0, select one calculate in vertex neighborhood the maximum pixel of phase masses value for
It is current to calculate point, wherein the current corresponding local writhing number of point that calculates is that the upper one local writhing number for calculating point and upper one calculate
Point and current the sum of the number of hops calculated between point, the current number of hops for calculating point are as follows: the relatively more current winding phase for calculating point
The difference for the winding phase for calculating point with upper one, difference are greater than the first preset value, then the current corresponding number of hops of point that calculates is -1；
Difference is less than the second preset value, then the current corresponding number of hops of point that calculates is 1；Otherwise, the current corresponding number of hops of point that calculates is
0.As an example,For the winding phase of a upper pixel,For the winding phase of current pixel point, ifThen number of hops cbaIt is -1, ifThen number of hops cbaIt is 1, c in other situationsbaIt is 0, according to
The secondary all pixels point part writhing number for solving all phase blocks.
S2. the local writhing number gradient between adjacent phase block is calculated, is twined according to the part between all adjacent phase blocks
Winding number gradient calculates number of hops and the corresponding mass value of each phase block between adjacent phase block, is asked using quality bootstrap algorithm
Solve the block writhing number of each phase block；
The calculation method of number of hops between adjacent phase block are as follows: set A and B as left and right adjacent phase block, at adjacent boundary
Number of corresponding pixels is N, the corresponding local writhing number Cr of the adjacent boundary ith pixel point of data block AiIt indicates, data
The corresponding local writhing number Cl of the adjacent boundary ith pixel point of block BiIt indicates, then the number of hops between data block A and B isWherein wiFor the weight after normalization, in the case of taking uniform weight as an example, wi=1/N；
The calculation method of number of hops between adjacent phase block can be with are as follows: the adjacent boundary i-th of statistics adjacent phase block
The difference Cl of local writhing number corresponding to pixeli-Cri, using the most difference of the frequency of occurrences as between adjacent phase block
Number of hops；
The corresponding mass value of each phase block, specifically, phase are calculated according to the number of hops between all adjacent phase blocks
The corresponding mass value of block is that the average value of the boundary writhing number gradient variance of phase block phase block adjacent thereto determines, phase block
Corresponding mass value is smaller, shows that the quality of the phase block is higher, can be used as the object of phase masses guidance next time；The phase
The average value of the boundary writhing number gradient variance of block phase block adjacent thereto, the as boundary of phase block phase block adjacent thereto
Pixel corresponding to local writhing number difference.
The block writhing number of each phase block is solved using quality bootstrap algorithm specifically: it is the smallest that corresponding mass value is set
Phase block is initial calculation point, and block writhing number corresponding to the phase block is 0, selects upper one to calculate phase masses value in vertex neighborhood
The smallest phase block is current calculating point, wherein current to calculate the part winding that corresponding piece of writhing number of point is upper calculating point
It is several to calculate point and current the sum of the number of hops calculated between point with upper one, successively solve the block writhing number of all phase blocks.
S3. it is corresponding the pixel to be obtained according to the corresponding local writhing number of each pixel, block writhing number and winding phase
Solution twine phase, to realize that fast and accurately obtaining the corresponding solution of original interference phase diagram all pixels point twines phase.
Specifically, the corresponding solution of the i-th row jth column pixel of original interference phase diagram twines phase are as follows:
Wherein,For the corresponding local writhing number of the i-th row jth column pixel,For phase where the i-th row jth column pixel
Position corresponding piece of writhing number of block.
Fig. 2 is the schematic diagram of the emulation interferometric phase image of the embodiment of the present invention.As shown in Fig. 2, utilizing peaks in Matlab
The emulation interferometric phase image that function generates, size are 2000 × 2000 pixels, including 3 low quality phase ranges are located at
Left side, middle part and lower right, the corresponding related coefficient in low quality region is 0.9 in the emulation interferometric phase image.
The schematic diagram of Fig. 3 corresponding writhing number when being the emulation non-piecemeal of interferometric phase image of the embodiment of the present invention.Such as Fig. 3 institute
Show, the corresponding writhing number of pixel solved to whole interference pattern using quantization quality bootstrap algorithm, pixel it is corresponding minimum and
Maximum writhing number is respectively -6 and 6.Contour line position in writhing number figure is high-visible, small by low quality regional effect.
The solution that Fig. 4 is solved when being the non-piecemeal of emulation interferometric phase image of the embodiment of the present invention twines the schematic diagram of phase.Such as Fig. 4
It is shown, according to writhing number and solution twine phase relationship can directly calculated solution twine phase, phase is integrally smooth, does not tire out
Product error.
Fig. 5, which is the emulation interferometric phase image block count of the embodiment of the present invention, solves the local writhing number of pixel when being 4 × 4
Schematic diagram.As shown in figure 5, solving pixel using quantization quality guidance to each piecemeal when block count is set as 4 × 4
Local writhing number schematic diagram, since the part winding numerical value in the corresponding initial phase site of each phase block is different, directly spelling
There is significant jump in block margin in local writhing number after connecing, the writhing number gradient referred to as between adjacent block；It solves
The range of the local writhing number of the pixel arrived is between -6 and 8.
Fig. 6, which is the emulation interferometric phase image block count of the embodiment of the present invention, solves the block writhing number of block of pixels when being 4 × 4
Schematic diagram.As shown in fig. 6, the block writhing number acquired using quantization quality bootstrap technique, it is -4 that smallest blocks, which wind numerical value, largest block
Winding numerical value is 6.
Fig. 7, which is the emulation interferometric phase image block count of the embodiment of the present invention, to be solved the writhing number of pixel when being 4 × 4 and shows
It is intended to.As shown in fig. 7, block writhing number is added the global writhing number obtained with corresponding local writhing number, and then find out pixel
The solution of point twines phase, and because gradient caused by piecemeal is eliminated, contour becomes continuously, when the result of solution and not piecemeal
Unanimously, range is also between -6 and 6.The relationship that phase is twined according to writhing number and solution is closed in the identical situation of writhing number figure
And it is also identical that rear resulting solution, which twines phase,.
Table 1 is that the different piecemeal condition solutions of the embodiment of the present invention twine the comparison schematic table of efficiency.As shown in table 1, using imitative
True interference pattern has carried out under different partitioned modes the solution that i.e. block count is respectively 1 × 1,2 × 2,4 × 4 and 8 × 8 and has twined calculating, piecemeal
Number is existing quality bootstrap algorithm for 1 × 1, is merged without carrying out the solution of block writhing number and writhing number, when block count is distinguished
It is respectively 1212,409,346 and 317ms that total solution when being 1 × 1,2 × 2,4 × 4 and 8 × 8, which twines the time, with the increasing of block count
Add, total solution twines that the time is accordingly reduced and its downward trend is more and more slow.As can be seen from the table, total solution twine the reduction of time be due to
The reduction of local writhing number acquisition time.Due to using partition strategy, multiple computing resource calculating office of CPU can be utilized simultaneously
Portion's writhing number is to promote computational efficiency.
The different piecemeal condition solutions of 1. embodiment of the present invention of table twine the comparison schematic table of efficiency
As it will be easily appreciated by one skilled in the art that the foregoing is merely illustrative of the preferred embodiments of the present invention, not to
The limitation present invention, any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should all include
Within protection scope of the present invention.
Claims (10)
1. a kind of block parallel quality guidance fast phase solution twines algorithm, which is characterized in that specific steps are as follows:
S1. original interference phase diagram data is obtained, is decomposed into original interference phase diagram according to preset phase block piecemeal parameter
Several phase blocks solve the local writhing number of the pixel in each phase block using quality bootstrap algorithm；
S2. the local writhing number gradient between adjacent phase block is calculated, according to the local writhing number between all adjacent phase blocks
Gradient calculates number of hops and the corresponding mass value of each phase block between phase block, is solved using quality bootstrap algorithm each
The block writhing number of phase block；
S3. the corresponding solution of the pixel is obtained according to the corresponding local writhing number of each pixel, block writhing number and winding phase
Phase is twined, to realize that fast and accurately obtaining the corresponding solution of original interference phase diagram all pixels point twines phase.
2. a kind of block parallel quality guidance fast phase solution according to claim 1 twines algorithm, which is characterized in that described
The local writhing number of the pixel in each phase block is solved using quality bootstrap algorithm specifically:
The corresponding phase masses value of each pixel in single phase block is calculated using original interference phase diagram data matrix；
It is initial calculation point, local writhing number corresponding to the pixel that the maximum pixel of phase masses value in phase block, which is arranged,
It is 0；
Selecting the maximum pixel of phase masses value in a calculating vertex neighborhood is current calculating point, current to calculate the corresponding office of point
Portion's writhing number is that the upper one local writhing number for calculating point and upper one calculate point and current the sum of the number of hops for calculating point.
3. a kind of block parallel quality guidance fast phase solution according to claim 2 twines algorithm, which is characterized in that described
The corresponding phase masses value of each pixel can utilize spurious correlation Y-factor method Y, phase gradient variance method or maximum in single phase block
Phase gradient method is calculated.
4. a kind of block parallel quality guidance fast phase solution according to claim 1 or 2 twines algorithm, which is characterized in that
The calculation method of the number of hops of current calculating point are as follows: the winding phase that the relatively more current winding phase for calculating point calculates point with upper one
Difference, difference is greater than the first preset value, then the current corresponding number of hops of point that calculates is -1；Difference is less than the second preset value, then
The current corresponding number of hops of point that calculates is 1；Otherwise, the current corresponding number of hops of point that calculates is 0.
5. a kind of guidance of block parallel quality described in any one of -4 fast phase solution twines algorithm according to claim 1, special
Sign is, the calculation method of the number of hops between the adjacent phase block are as follows:
Number of corresponding pixels is N, the adjacent edge of one of phase block at corresponding adjacent boundary between adjacent phase block
The corresponding local writhing number Cr of boundary's ith pixel pointiIt indicates, the corresponding office of adjacent boundary ith pixel point of another phase block
Portion writhing number CliIt indicates, then the number of hops between adjacent phase block isWherein wiAfter normalization
Weight.
6. a kind of block parallel quality guidance fast phase solution according to claim 5 twines algorithm, which is characterized in that normalizing
Weight after change is wi=1/N.
7. a kind of guidance of block parallel quality described in any one of -4 fast phase solution twines algorithm according to claim 1, special
Sign is, the calculation method of the number of hops between the adjacent phase block are as follows: the adjacent boundary ith pixel of statistics adjacent phase block
The difference Cl of the corresponding local writhing number of pointi-Cri, using the most difference of the frequency of occurrences as the jump between adjacent phase block
Number.
8. a kind of guidance of the block parallel quality according to claim 5 or 7 fast phase solution twines algorithm, which is characterized in that
The corresponding mass value of the phase block is that part corresponding to the pixel of the boundary of phase block phase block adjacent thereto twines
The variance of the difference of winding number determines.
9. a kind of block parallel quality guidance fast phase solution according to claim 8 twines algorithm, which is characterized in that utilize
Quality bootstrap algorithm solves the block writhing number of each phase block specifically: it is initial that the corresponding the smallest phase block of mass value, which is arranged,
Point is calculated, block writhing number corresponding to initial calculation point is 0, selects upper one to calculate the smallest phase of phase masses value in vertex neighborhood
Block be it is current calculate point, it is current calculate corresponding piece of writhing number of point be the upper one local writhing number for calculating point and upper one calculate point with
Current the sum of the number of hops calculated between point.
10. a kind of block parallel quality guidance fast phase solution according to claim 1 to 9 twines algorithm, special
Sign is that the corresponding solution of the i-th row jth column pixel of original interference phase diagram twines phase are as follows:
Wherein,For the corresponding local writhing number of the i-th row jth column pixel,For phase block where the i-th row jth column pixel
Corresponding piece of writhing number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910263071.3A CN110068823B (en) | 2019-04-02 | 2019-04-02 | Block parallel quality guide rapid phase unwrapping algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910263071.3A CN110068823B (en) | 2019-04-02 | 2019-04-02 | Block parallel quality guide rapid phase unwrapping algorithm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110068823A true CN110068823A (en) | 2019-07-30 |
| CN110068823B CN110068823B (en) | 2021-03-09 |
Family
ID=67367031
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910263071.3A Active CN110068823B (en) | 2019-04-02 | 2019-04-02 | Block parallel quality guide rapid phase unwrapping algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110068823B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117724101A (en) * | 2024-02-18 | 2024-03-19 | 中国特种设备检测研究院 | Kalman filtering phase unwrapping method, system and equipment based on phase block |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002025474A2 (en) * | 2000-09-19 | 2002-03-28 | California Institute Of Technology | Efficient method of identifying non-solution or non-optimal regions of the domain of a function |
| WO2008118739A2 (en) * | 2007-03-22 | 2008-10-02 | Harris Corporation | Method and apparatus for compression of sar images |
| US20120278005A1 (en) * | 2005-01-21 | 2012-11-01 | Chikayoshi Sumi | Clinical apparatuses |
| CN103279945A (en) * | 2013-04-26 | 2013-09-04 | 北京理工大学 | Quality map guide method and branch cutting method-based interferogram unwrapping method |
| CN103325092A (en) * | 2013-03-13 | 2013-09-25 | 中国科学院电子学研究所 | Method and device for generating two-dimensional phase disentanglement quality picture |
| CN103913128A (en) * | 2014-03-12 | 2014-07-09 | 哈尔滨工程大学 | Rapid phase unwrapping method based on quality diagram region segmentation |
| CN103942095A (en) * | 2014-03-18 | 2014-07-23 | 中国科学院软件研究所 | Two-dimensional phase position unwrapping method based on heterogeneous accelerating platform |
| CN104316922A (en) * | 2014-10-11 | 2015-01-28 | 南京邮电大学 | Regional division based multi-strategy InSAR (Interferometric synthetic aperture) radar phase unwrapping method |
| CN105844626A (en) * | 2016-03-18 | 2016-08-10 | 电子科技大学 | Magnetic resonance phase unwrapping method based on wrapping identification and local curved surface fitting |
| CN107092022A (en) * | 2017-04-21 | 2017-08-25 | 哈尔滨工业大学 | Region filter quality guiding phase unwrapping method based on InSAL |
| CN107730491A (en) * | 2017-10-13 | 2018-02-23 | 北京工业大学 | A kind of phase unwrapping package method based on Quality Map |
| CN107783079A (en) * | 2017-09-28 | 2018-03-09 | 淮海工学院 | One kind uses L0The quick unwrapping method of two-dimensional phase of norm cost function |
| CN108254729A (en) * | 2018-04-09 | 2018-07-06 | 北京环境特性研究所 | A kind of double fit phase unwrapping methods and double fit phase unwrapping devices |
| CN108523890A (en) * | 2018-03-06 | 2018-09-14 | 厦门大学 | Local field method of estimation based on magnetic resonance dipole field spatial distribution |
| CN108986120A (en) * | 2018-07-13 | 2018-12-11 | 长春师范大学 | A kind of quick and reliable phase unwrapping algorithm |
| CN109490895A (en) * | 2018-10-25 | 2019-03-19 | 中国人民解放军海军工程大学 | A kind of interference synthetic aperture signal processing system based on blade server |
-
2019
- 2019-04-02 CN CN201910263071.3A patent/CN110068823B/en active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002025474A2 (en) * | 2000-09-19 | 2002-03-28 | California Institute Of Technology | Efficient method of identifying non-solution or non-optimal regions of the domain of a function |
| US20120278005A1 (en) * | 2005-01-21 | 2012-11-01 | Chikayoshi Sumi | Clinical apparatuses |
| WO2008118739A2 (en) * | 2007-03-22 | 2008-10-02 | Harris Corporation | Method and apparatus for compression of sar images |
| CN103325092A (en) * | 2013-03-13 | 2013-09-25 | 中国科学院电子学研究所 | Method and device for generating two-dimensional phase disentanglement quality picture |
| CN103279945A (en) * | 2013-04-26 | 2013-09-04 | 北京理工大学 | Quality map guide method and branch cutting method-based interferogram unwrapping method |
| CN103913128A (en) * | 2014-03-12 | 2014-07-09 | 哈尔滨工程大学 | Rapid phase unwrapping method based on quality diagram region segmentation |
| CN103942095A (en) * | 2014-03-18 | 2014-07-23 | 中国科学院软件研究所 | Two-dimensional phase position unwrapping method based on heterogeneous accelerating platform |
| CN104316922A (en) * | 2014-10-11 | 2015-01-28 | 南京邮电大学 | Regional division based multi-strategy InSAR (Interferometric synthetic aperture) radar phase unwrapping method |
| CN105844626A (en) * | 2016-03-18 | 2016-08-10 | 电子科技大学 | Magnetic resonance phase unwrapping method based on wrapping identification and local curved surface fitting |
| CN107092022A (en) * | 2017-04-21 | 2017-08-25 | 哈尔滨工业大学 | Region filter quality guiding phase unwrapping method based on InSAL |
| CN107783079A (en) * | 2017-09-28 | 2018-03-09 | 淮海工学院 | One kind uses L0The quick unwrapping method of two-dimensional phase of norm cost function |
| CN107730491A (en) * | 2017-10-13 | 2018-02-23 | 北京工业大学 | A kind of phase unwrapping package method based on Quality Map |
| CN108523890A (en) * | 2018-03-06 | 2018-09-14 | 厦门大学 | Local field method of estimation based on magnetic resonance dipole field spatial distribution |
| CN108254729A (en) * | 2018-04-09 | 2018-07-06 | 北京环境特性研究所 | A kind of double fit phase unwrapping methods and double fit phase unwrapping devices |
| CN108986120A (en) * | 2018-07-13 | 2018-12-11 | 长春师范大学 | A kind of quick and reliable phase unwrapping algorithm |
| CN109490895A (en) * | 2018-10-25 | 2019-03-19 | 中国人民解放军海军工程大学 | A kind of interference synthetic aperture signal processing system based on blade server |
Non-Patent Citations (7)
| Title |
|---|
| HEPING ZHONG 等: "A Quality-Guided and Local Minimum Discontinuity Based Phase Unwrapping Algorithm for InSAR/InSAS Interferograms", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
| T.J. FLYNN 等: "phase unwrapping,local wrap count,wrap number", 《IGARSS "96. 1996 INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM》 * |
| ZHANG YAN 等: "A Hybrid Phase Unwrapping Algorithm Based on Quality-Guided and Surface-Fitting", 《2018 IEEE INTERNATIONAL WORKSHOP ON ELECTROMAGNETICS:APPLICATIONS AND STUDENT INNOVATION COMPETITION (IWEM)》 * |
| 张斌 等: "多基线InSAR图割相位解缠算法研究", 《多基线INSAR图割相位解缠算法研究》 * |
| 郝红星 等: "质量图引导的干涉合成孔径雷达图像相位解缠算法", 《装备学院学报》 * |
| 钟何平 等: "一种分区域的复合相位解缠算法", 《武汉大学学报 信息科学版》 * |
| 钟何平 等: "基于最小不连续的分块相位解缠算法", 《系统工程与电子技术》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117724101A (en) * | 2024-02-18 | 2024-03-19 | 中国特种设备检测研究院 | Kalman filtering phase unwrapping method, system and equipment based on phase block |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110068823B (en) | 2021-03-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10223802B2 (en) | Method for determining depth maps from stereo images with improved depth resolution in a range | |
| CN106093849B (en) | A kind of Underwater Navigation method based on ranging and neural network algorithm | |
| CN102884397B (en) | Structured-light measuring method and system | |
| CN110675418A (en) | Target track optimization method based on DS evidence theory | |
| CN110427807A (en) | A kind of temporal events motion detection method | |
| CN110348356A (en) | A kind of successive frame RD images steganalysis method based on depth light stream network | |
| CN109308713B (en) | Improved nuclear correlation filtering underwater target tracking method based on forward-looking sonar | |
| CN105487061A (en) | Multi-characteristic information fusion method for target data correlation | |
| CN110909591A (en) | Self-adaptive non-maximum value inhibition processing method for pedestrian image detection by using coding vector | |
| CN111414931B (en) | Multi-branch multi-scale small target detection method based on image depth | |
| CN104316922A (en) | Regional division based multi-strategy InSAR (Interferometric synthetic aperture) radar phase unwrapping method | |
| CN109839619A (en) | Based on radar signal rough segmentation choosing method, system and the storage medium for adaptively dividing bucket | |
| CN113311433B (en) | InSAR interferometric phase two-step unwrapping method combining quality map and minimum cost flow | |
| CN111126116A (en) | Unmanned ship river channel garbage identification method and system | |
| CN110688905A (en) | Three-dimensional object detection and tracking method based on key frame | |
| US7602966B2 (en) | Image processing method, image processing apparatus, program and recording medium | |
| CN106067172A (en) | A kind of underwater topography image based on suitability analysis slightly mates and mates, with essence, the method combined | |
| CN108629809B (en) | Accurate and efficient stereo matching method | |
| CN110068823A (en) | A kind of block parallel quality guidance fast phase solution twines algorithm | |
| CN105469392A (en) | High spectral image significance detection method based on regional spectrum gradient characteristic comparison | |
| KR100930286B1 (en) | Cost function calculation method and hidden region processing method for variance estimation | |
| CN110598711A (en) | Target segmentation method combined with classification task | |
| Li et al. | Pillar-based 3D object detection from point cloud with multiattention mechanism | |
| CN104881670A (en) | Rapid target extraction method used for SAR azimuth estimation | |
| CN110060282B (en) | Point cloud registration method and device |
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 |