阅读说明：本技术 相位解包裹方法、装置、系统、计算机设备及存储介质 (Phase- un- wrapping method, apparatus, system, computer equipment and storage medium ) 是由 黄佐华 林盼明 罗都 彭力 于 2019-07-15 设计创作，主要内容包括：本发明公开了一种相位解包裹方法、装置、系统、计算机设备及存储介质,所述方法包括：获取包裹相位干涉图；确定阈值；通过阈值确定包裹相位干涉图的凹陷边界,形成翻转临界线；对包裹相位干涉图的凹陷部分进行自动翻转替换,得到第一解包裹相位矩阵；对第一解包裹相位矩阵进行滤波平滑处理,得到第二解包裹相位矩阵；计算并保存第一解包裹相位矩阵与第二解包裹相位矩阵的相对误差；当阈值未达到预设值时,返回重新确定阈值及执行后续操作；当阈值达到预设值时,输出最小相对误差对应的阈值和第二解包裹相位矩阵；根据最小相对误差对应的第二解包裹相位矩阵,输出重构的真实相位分布。本发明提高了包裹相位干涉图的真实相位还原的效率及准确性。(The invention discloses a kind of Phase- un- wrapping method, apparatus, system, computer equipment and storage mediums, which comprises obtains wrapped phase interference pattern；Threshold value；The recess boundary of wrapped phase interference pattern is determined by threshold value, forms overturning critical line；Automatic turning replacement is carried out to the depressed section of wrapped phase interference pattern, obtains the first Phase Unwrapping Algorithm matrix；Smoothing processing is filtered to the first Phase Unwrapping Algorithm matrix, obtains the second Phase Unwrapping Algorithm matrix；Calculate and save the relative error of the first Phase Unwrapping Algorithm matrix and the second Phase Unwrapping Algorithm matrix；When threshold value is not up to preset value, return redefines threshold value and executes subsequent operation；When threshold value reaches preset value, the corresponding threshold value of minimum relative error and the second Phase Unwrapping Algorithm matrix are exported；According to the corresponding second Phase Unwrapping Algorithm matrix of minimum relative error, the true phase distribution of reconstruct is exported.The present invention improves the efficiency and accuracy of the true phase reduction of wrapped phase interference pattern.)

1. a kind of phase unwrapping package method, which is characterized in that the described method includes:

Obtain wrapped phase interference pattern；

Threshold value；

The recess boundary of wrapped phase matrix is determined by threshold value, forms overturning critical line；

Automatic turning replacement is carried out to the depressed section of wrapped phase matrix, obtains the first Phase Unwrapping Algorithm matrix；

Smoothing processing is filtered to the first Phase Unwrapping Algorithm matrix, obtains the second Phase Unwrapping Algorithm matrix；

Calculate and save the relative error of the first Phase Unwrapping Algorithm matrix and the second Phase Unwrapping Algorithm matrix；

When threshold value is not up to preset value, return redefines threshold value, and executes subsequent operation；

When threshold value reaches preset value, the corresponding threshold value of minimum relative error and the second Phase Unwrapping Algorithm matrix are exported；

According to the corresponding second Phase Unwrapping Algorithm matrix of minimum relative error, the true phase distribution of reconstruct is exported.

2. phase unwrapping package method according to claim 1, which is characterized in that number of the threshold value between 0.9 and 1 Value, wherein 0.9 is first threshold value, 1 is the last one threshold value；It is described redefine threshold value and refer to present threshold value is added 0.01；The preset value is 1.

3. -2 described in any item phase unwrapping package methods according to claim 1, which is characterized in that described determined by threshold value is wrapped The recess boundary of phasing matrix is wrapped up in, overturning critical line is formed, specifically includes:

Using threshold value W as the foundation on the recess boundary for obtaining wrapped phase matrix, extracts phase value in wrapped phase matrix and be greater than The pixel of W*phi (max), and 1 is set by the pixel region, wherein phi (max) is that the phase of input phase matrix is maximum Value, that is, the pixel region for being set as 1 is depressed section, and rest of pixels is set as 0.

4. -2 described in any item phase unwrapping package methods according to claim 1, which is characterized in that described to wrapped phase matrix Depressed section carry out automatic turning replacement, obtain the first Phase Unwrapping Algorithm matrix, specifically include:

The depressed section of wrapped phase matrix is overturn, the recess matrix obtained after overturning is replaced in wrapped phase matrix The position of recess obtains new phasing matrix；

When the matrix that is recessed is non-zero matrix, which is overturn, the recess matrix obtained after overturning replacement is worked as The position of recess in preceding phasing matrix, obtains new phasing matrix；

When the matrix that is recessed is null matrix, using current phasing matrix as the first Phase Unwrapping Algorithm matrix, and exported.

5. -2 described in any item phase unwrapping package methods according to claim 1, which is characterized in that the acquisition wrapped phase is dry After relating to figure, further includes:

Wrapped phase interference pattern is pre-processed.

6. a kind of Phase- un- wrapping device, which is characterized in that described device includes:

Module is obtained, for obtaining wrapped phase interference pattern；

First determining module is used for threshold value；

Second determining module forms overturning critical line for determining the recess boundary of wrapped phase matrix by threshold value；

Replacement module is overturn, automatic turning replacement is carried out for the depressed section to wrapped phase matrix, obtains the first unpacking Phasing matrix；

Processing module obtains the second Phase Unwrapping Algorithm matrix for being filtered smoothing processing to the first Phase Unwrapping Algorithm matrix；

Computing module, for calculating and saving the relative error of the first Phase Unwrapping Algorithm matrix and the second Phase Unwrapping Algorithm matrix；

Return module, for when threshold value is not up to preset value, return to redefine threshold value, and executes subsequent operation；

First output module is unpacked for when threshold value reaches preset value, exporting the corresponding threshold value of minimum relative error and second Wrap up in phasing matrix；

Second output module, for exporting the true of reconstruct according to the corresponding second Phase Unwrapping Algorithm matrix of minimum relative error Phase distribution.

7. Phase- un- wrapping device according to claim 6, which is characterized in that the overturning replacement module specifically includes:

First overturning replacement unit, overturns, the recess that will be obtained after overturning for the depressed section to wrapped phase matrix Matrix replaces the position of recess in wrapped phase matrix, obtains new phasing matrix；

Second overturning replacement unit, for being overturn to the recess matrix, after overturning when the matrix that is recessed is non-zero matrix Obtained recess matrix replaces the position of recess in current phasing matrix, obtains new phasing matrix；

Output unit is used for when the matrix that is recessed is null matrix, using current phasing matrix as the first Phase Unwrapping Algorithm matrix, and It is exported.

8. a kind of Phase- un- wrapping system, which is characterized in that the system comprises the semiconductor laser being sequentially connected, adjustable decline Subtract device, beam expanding lens, pinhole filter, the first collimation lens, the second collimation lens, phase-plate, third collimation lens, video camera and meter Calculation machine, drop target object between first collimation lens and the second collimation lens, the target object are located at the second collimation At the front focal plane of lens, the phase-plate is located at the back focal plane of the second collimation lens, and the position for video camera is in third collimation lens Back focal plane at；

It is saturating that the laser beam that the semiconductor laser issues passes sequentially through attenuator, beam expanding lens, pinhole filter, the first collimation Plane light wave is formed after mirror, to irradiate target object；The frequency spectrum of Object light wave, phase-plate are obtained on the back focal plane of the second collimation lens The zero-frequency region for making phase-modulation region alignment Object light wave frequency spectrum, by the light wave in zero-frequency region with pass through the light outside zero-frequency region Wave experienced different phase-modulations, wrapped phase interference pattern be generated on the back focal plane of third collimation lens, by video camera pair Wrapped phase interference pattern is acquired, and inputs computer；

The computer requires the described in any item phase unwrapping package methods of 1-5 for perform claim.

9. a kind of computer equipment, including processor and for the memory of storage processor executable program, feature exists In, when the processor executes the program of memory storage, the realization described in any item Phase- un- wrapping sides claim 1-5 Method.

10. a kind of storage medium, is stored with program, which is characterized in that when described program is executed by processor, realize claim The described in any item phase unwrapping package methods of 1-5.

Technical field

The present invention relates to a kind of Phase- un- wrapping method, apparatus, system, computer equipment and storage mediums, belong to image Treatment research field.

Background technique

Phase- un- wrapping calculating is widely used in interference synthetic aperture radar, the three-dimensional measurement based on structure light, optics In the measurement methods such as interferometry and imaging.During obtaining phase from package interference pattern, arc tangent norm is generally used Extract phase, cause continuous phase to be truncated, codomain be limited at (- π, π], formed wrapped phase figure.Therefore it needs to package Phase carries out unpacking, phase reconstruction will be truncated, in the hope of true continuous phase.But under actual conditions, wrapped phase due to Lack sampling contains the influence factors such as noise, local shades, often will appear residue points region, influences the precision of unpacking, to unpacking It wraps up in and increases difficulty, also influence the accuracy of true phase reconstruct indirectly.

Existing phase unwrapping algorithm and method are numerous, be summed up there are mainly two types of: path tracking algorithm and path Unrelated algorithm.Wherein, path tracking algorithm mainly has Branch cut, Quality Map guiding method and exposure mask secant method etc.；The unrelated calculation in path Method is mainly to seek to meet the solution of minimum norm, least square method such as based on discrete cosine transform, based on fast Fourier change The least square method etc. changed.Serious lack sampling and be two problems for unpacking compared with very noisy, Qian Xiaofan etc. proposes cross To shearing least square method (LSBLS), Guo Yuan etc. proposes four-way lateral shear least square method (FSLBS).These algorithms can To restore lack sampling phase diagram, but it still can not effectively solve the wrapped phase figure of serious lack sampling.

In recent years, adaptive threshold image segmentation, denoising, defogging, in terms of be widely used.Foreign scholar Eser Sert et al. is based on structured-light system and proposes that adaptive threshold phase transition three-level calibration algorithm carries out three-dimension object survey Amount.In terms of image procossing, there is researcher to propose the edge projection profilometry for removing background based on adaptive threshold, also have Adaptive threshold is used for cell image processing and image denoising sound etc..Currently, adaptive threshold is gone back in terms of image unpacking It is not widely applied.

Summary of the invention

In view of this, the present invention provides a kind of Phase- un- wrapping method, apparatus, system, computer equipment and storages to be situated between The precision of matter, Phase- un- wrapping is high, improves the efficiency and accuracy of the true phase reduction of wrapped phase interference pattern, solves Existing phase unwrapping algorithm there are the problem of.

The first purpose of this invention is to provide a kind of phase unwrapping package method.

Second object of the present invention is to provide a kind of Phase- un- wrapping device.

Third object of the present invention is to provide a kind of Phase- un- wrapping system.

Fourth object of the present invention is to provide a kind of computer equipment.

Of the invention the 5th is designed to provide a kind of storage medium.

The first purpose of this invention can be reached by adopting the following technical scheme that:

A kind of phase unwrapping package method, which comprises

Obtain wrapped phase interference pattern；

Threshold value；

The recess boundary of wrapped phase matrix is determined by threshold value, forms overturning critical line；

Automatic turning replacement is carried out to the depressed section of wrapped phase matrix, obtains the first Phase Unwrapping Algorithm matrix；

Smoothing processing is filtered to the first Phase Unwrapping Algorithm matrix, obtains the second Phase Unwrapping Algorithm matrix；

Calculate and save the relative error of the first Phase Unwrapping Algorithm matrix and the second Phase Unwrapping Algorithm matrix；

When threshold value is not up to preset value, return redefines threshold value, and executes subsequent operation；

When threshold value reaches preset value, the corresponding threshold value of minimum relative error and the second Phase Unwrapping Algorithm matrix are exported；

According to the corresponding second Phase Unwrapping Algorithm matrix of minimum relative error, the true phase distribution of reconstruct is exported.

Further, numerical value of the threshold value between 0.9 and 1, wherein 0.9 is first threshold value, 1 is the last one threshold Value；The threshold value that redefines refers to present threshold value plus 0.01；The preset value is 1.

Further, the recess boundary that wrapped phase matrix is determined by threshold value forms overturning critical line, specific to wrap It includes:

Using threshold value W as the foundation on the recess boundary for obtaining wrapped phase matrix, phase value in wrapped phase matrix is extracted Be set as 1 greater than the pixel of W*phi (max), and by the pixel region, wherein phi (max) be input phase matrix phase most Big value, that is, the pixel region for being set as 1 is depressed section, and rest of pixels is set as 0.

Further, the depressed section to wrapped phase matrix carries out automatic turning replacement, obtains the first unpacking Phasing matrix specifically includes:

The depressed section of wrapped phase matrix is overturn, the recess matrix obtained after overturning is replaced into wrapped phase square The position of recess, obtains new phasing matrix in battle array；

When the matrix that is recessed is non-zero matrix, which is overturn, the recess matrix obtained after overturning is replaced The position for changing recess in current phasing matrix obtains new phasing matrix；

When the matrix that is recessed is null matrix, using current phasing matrix as the first Phase Unwrapping Algorithm matrix, and exported.

Further, after the acquisition wrapped phase interference pattern, further includes:

Wrapped phase interference pattern is pre-processed, it whenever necessary can be by anticosine wrapped phase interference pattern in real process It is converted into arc tangent wrapped phase interference pattern.

Second object of the present invention can be reached by adopting the following technical scheme that:

A kind of Phase- un- wrapping device, described device include:

Module is obtained, for obtaining wrapped phase interference pattern；

First determining module is used for threshold value；

Second determining module forms overturning critical line for determining the recess boundary of wrapped phase matrix by threshold value；

Replacement module is overturn, automatic turning replacement is carried out for the depressed section to wrapped phase matrix, obtains the first solution Wrapped phase matrix；

Processing module obtains the second Phase Unwrapping Algorithm for being filtered smoothing processing to the first Phase Unwrapping Algorithm matrix Matrix；

Computing module is missed for calculating and saving the first Phase Unwrapping Algorithm matrix and the opposite of the second Phase Unwrapping Algorithm matrix Difference；

Return module, for when threshold value is not up to preset value, return to redefine threshold value, and executes subsequent operation；

First output module, for when threshold value reaches preset value, exporting the corresponding threshold value of minimum relative error and second Phase Unwrapping Algorithm matrix；

Second output module, for exporting reconstruct according to the corresponding second Phase Unwrapping Algorithm matrix of minimum relative error True phase distribution.

Further, the overturning replacement module, specifically includes:

First overturning replacement unit, overturns for the depressed section to wrapped phase matrix, by what is obtained after overturning The matrix that is recessed replaces the position of recess in wrapped phase matrix, obtains new phasing matrix；

Second overturning replacement unit, for overturning, turning over to the recess matrix when the matrix that is recessed is non-zero matrix The recess matrix obtained after turning replaces the position of recess in current phasing matrix, obtains new phasing matrix；

Output unit is used for when the matrix that is recessed is null matrix, using current phasing matrix as the first Phase Unwrapping Algorithm square Battle array, and exported.

Further, after the acquisition module, further includes:

Preprocessing module can incite somebody to action anti-remaining whenever necessary for pre-processing to wrapped phase interference pattern, in real process String wrapped phase interference pattern is converted into arc tangent wrapped phase interference pattern.

Third object of the present invention can be reached by adopting the following technical scheme that:

A kind of Phase- un- wrapping system the system comprises the semiconductor laser being sequentially connected, adjustable attenuator, expands Mirror, pinhole filter, the first collimation lens, the second collimation lens, phase-plate, third collimation lens, video camera and computer, it is described Drop target object between first collimation lens and the second collimation lens, the target object are located at the preceding coke of the second collimation lens At face, the phase-plate is located at the back focal plane of the second collimation lens, and the position for video camera is at the back focal plane of third collimation lens；

The laser beam that the semiconductor laser issues passes sequentially through attenuator, beam expanding lens, pinhole filter, the first standard Plane light wave is formed after straight lens, to irradiate target object；The frequency spectrum of Object light wave is obtained on the back focal plane of the second collimation lens, Phase-plate makes the zero-frequency region of phase-modulation region alignment Object light wave frequency spectrum, by the light wave in zero-frequency region and by outside zero-frequency region Light wave experienced different phase-modulations, on the back focal plane of third collimation lens generate wrapped phase interference pattern, by imaging Machine is acquired wrapped phase interference pattern, and inputs computer；

The computer, for executing above-mentioned phase unwrapping package method.

Fourth object of the present invention can be reached by adopting the following technical scheme that:

A kind of computer equipment, including processor and for the memory of storage processor executable program, the place When managing the program of device execution memory storage, above-mentioned phase unwrapping package method is realized.

5th purpose of the invention can be reached by adopting the following technical scheme that:

A kind of storage medium is stored with program, when described program is executed by processor, realizes above-mentioned Phase- un- wrapping side Method.

The present invention have compared with the existing technology it is following the utility model has the advantages that

1, for the present invention after obtaining wrapped phase interference pattern, threshold value determines the recessed of wrapped phase matrix by threshold value Boundary is fallen into, overturning critical line is formed, automatic turning replacement is carried out to the depressed section of wrapped phase matrix, automatic turning is replaced Phasing matrix afterwards is filtered smoothing processing, calculate and save filtering before and after the processing two phasing matrixs it is opposite accidentally Difference redefines threshold value, executes the subsequent operation of threshold value, when threshold value reaches preset value, just terminates circulation, leads at this time The relative error that more each threshold value saves is crossed, minimum relative error is obtained, exports the corresponding threshold value of minimum relative error With filtering treated phasing matrix, restore the method for phase using phase bit flipping in Phase- un- wrapping, it is maximum excellent Point is that overturning critical line is determined using threshold value, complete true phase can be solved with the smallest relative error, to improve The precision of Phase- un- wrapping, and by way of Automatic-searching threshold value, it can effectively avoid manual setting threshold value, provide correct Threshold value selection, improves the efficiency and accuracy of phase reduction.

2, present invention employs depressed section of the automatic turning algorithm to wrapped phase matrix to handle, to wrapped phase The depressed section of matrix is overturn, and the recess matrix obtained after overturning is replaced to the position of recess in wrapped phase matrix, Obtain new phasing matrix, at this time can to recess matrix judge, if recess matrix be non-zero matrix, to the recess matrix into The recess matrix obtained after overturning, is replaced the position of recess in current phasing matrix, obtains new phasing matrix by row overturning, When the matrix that is recessed is null matrix, circulation is just terminated, in such a way that automatic turning is replaced, so that in Phase Unwrapping Algorithm matrix It can not find recess.

3, the present invention can accurately restore true phase compared with very noisy, have certain noiseproof feature, solve The problem of noise will affect the precision of unpacking.

4, the error for the reconstructed phase and true phase that the present invention obtains is smaller, can be to the biggish wrapped phase of size of mesh opening Interference pattern (such as 1501 × 1501) realizes Phase- un- wrapping, improves the precision of Phase- un- wrapping, is high-precision, intensive Unpacking provides new reference method.

5, the present invention is suitable for the unpacking for the various wrapped phase interference patterns that wrapped cycle is π, 2 π, 3 π etc., using model It encloses and field is extensive.

Detailed description of the invention

In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with The structure shown according to these attached drawings obtains other attached drawings.

Fig. 1 is the index path of the Phase- un- wrapping system of the embodiment of the present invention 1.

Fig. 2 is the background plot of light intensity of the embodiment of the present invention 1.

Fig. 3 is the wrapped phase interference pattern of the embodiment of the present invention 1.

Fig. 4 is the thinking flow chart of the phase unwrapping package method of the embodiment of the present invention 1.

Fig. 5 is the entire flow figure of the phase unwrapping package method of the embodiment of the present invention 1.

Fig. 6 is that the automatic turning of the embodiment of the present invention 1 replaces the flow chart of algorithm.

Fig. 7 is the true phase point of the embodiment of the present invention 1 obtained using adaptive threshold overturning algorithm progress unpacking Butut.

Fig. 8 is to carry out the true phase distribution map that unpacking obtains using existing lateral shear least square solution.

Fig. 9 is the structural block diagram of the Phase- un- wrapping device of the embodiment of the present invention 3.

Figure 10 is the structural block diagram of the overturning replacement module of the embodiment of the present invention 3.

Figure 11 is the structural block diagram of the computer equipment of the embodiment of the present invention 4.

Specific embodiment

In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is A part of the embodiments of the present invention, instead of all the embodiments, based on the embodiments of the present invention, ordinary skill people Member's every other embodiment obtained without making creative work, shall fall within the protection scope of the present invention.