CN102608205A - Multilayer-object nondestructive testing ultrasonic imaging method based on variable wave speed phase shifting - Google Patents
Multilayer-object nondestructive testing ultrasonic imaging method based on variable wave speed phase shifting Download PDFInfo
- Publication number
- CN102608205A CN102608205A CN2012100461511A CN201210046151A CN102608205A CN 102608205 A CN102608205 A CN 102608205A CN 2012100461511 A CN2012100461511 A CN 2012100461511A CN 201210046151 A CN201210046151 A CN 201210046151A CN 102608205 A CN102608205 A CN 102608205A
- Authority
- CN
- China
- Prior art keywords
- max
- value
- medium
- delta
- omega
- 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
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention discloses a multilayer-object nondestructive testing ultrasonic imaging method based on variable wave speed phase shifting, which belongs to the technical field of fast imaging of an inner structure of a multilayer heterogeneous object. The method is characterized in that under a condition that equal-interval moving detection can be carried out in the length direction of the multilayer heterogeneous object based on an ultrasonic transducer, in a homogeneous area, images in each row in the homogeneous area are calculated by a constant wave speed phase shifting ultrasonic imaging technology and a fast Fourier transform technology, in a heterogeneous area, a bounding box bounding medium boundaries is built, and for the heterogeneous areas positioned above and below the boundaries in the bounding box, images in each row in the heterogeneous area are calculated by a variable wave speed phase shifting method and an N-time single output fast Fourier transform technology. The multilayer-object nondestructive testing ultrasonic imaging method is suitable for the fast ultrasonic imaging of the multilayer heterogeneous object with heterogeneous media in the depth direction and the horizontal direction, and a respective treating method adopting variable wave speed as a characteristic can be adopted close to bonding surfaces of the media, so that the accuracy is high.
Description
Technical field
The present invention relates to Ultrasonic Nondestructive technology, synthetic aperture focusing imaging technique and Phase shift technology, realize inner structure fast imaging the non-regular layering object of complicacy.
Background technology
At present; In the Ultrasonic Nondestructive field; (comprise contact ultrasonic imaging and immersion method ultrasonic imaging for the ultrasonic imaging that contains interfacial layering object; If the couplant layer is regarded as the part of testee; Then the immersion method ultrasonic imaging also can be considered the contact ultrasonic imaging of layering object; Therefore among the present invention with both unified ultrasonic imagings that is called the layering object), currently mainly contain two kinds of schemes: the one, the method that adopts time domain ultrasonic imaging by synthetic aperture focusing technology to combine with ray trace (Ray Tracing) technology, the 2nd, adopt based on Phase shift (Phase Shift Migration) technology and the technological Phase shift ultrasonic imaging method of frequency domain ultrasonic imaging by synthetic aperture focusing.
The basic thought of ultrasonic imaging by synthetic aperture focusing technology (SAFT) is to utilize a single transducer to simulate array of apertures; Testee is carried out orderly scanning, and adopt time-delay stack (DAS) method (time delay or phase delay) that the pulse echo signal that scanning obtains is carried out focal imaging.The SAFT working model is as shown in Figure 1, and wherein, Z is to being depth direction, and X is to being the transducer direction of scanning, when transducer is positioned at u
nDuring the place, the distance of target reflection thing (x ', z ') to transducer is r
nThe SAFT technology has the branch of time domain and frequency domain:
Time domain SAFT calculates different time delays based on DAS principle and multiple spot dynamic focus technology to the focus point on the different depth of target imaging zone.Like Fig. 1, for (x ', z ') focuses at target reflection thing place, the echoed signal that time domain SAFT technology obtains transducer each scanning position place in the length of synthetic aperture L overlap-add procedure of delaying time: establish s (u
n, be that transducer is at u t)
nThe echoed signal that the place receives, t is sampling instant, u
nThe place about the time-delay of target reflection thing (x ', z ') does
Wherein, c is ultrasonic velocity of propagation in medium, then is imaged as what (x ', z ') located
Wherein, ω
nBe apodizing function, the δ function is the unit sampling sequence.
Frequency domain SAFT is called as wave number algorithm (wavenumber algorithm), mainly is based on the method that phase delay generates wave beam, and it stems from the wave equation inversion theory.Ultimate principle is that Green function (Green function) is carried out Fourier decomposition.Algorithmic procedure is for carrying out two-dimensional Fourier transform to ultrasound data earlier; Obtain 2-d spectrum; Utilize Nonlinear Mapping (Stolt mapping) and interpolation to realize coordinate transform then; At last the frequency spectrum after the conversion is done Filtering Processing, and carry out two-dimensional inverse Fourier transform, the reconstructed image under the span-time domain.
Phase shift (Phase Shift Migration) ultrasonic imaging method is that the migrating technology (Migration Technique) in the reflection seismology (Reflection Seismology) is combined with frequency domain SAFT, and a kind of ultrasonic imaging method that obtains.The Phase shift technology is proposed in 1978 by Gazdag the earliest, is used for earthquake/radar imagery always, and the Phase shift ultrasonic imaging method that proposes has recently just applied to the Ultrasonic NDT field with this technology.
The Phase shift ultrasonic imaging method is regarded as the explosive reflection model with the supersonic sounding model, supposes that the reverberation in the object under test explodes at t=0 constantly simultaneously, and the bursting strength of each reverberation is proportional to its reflectivity, and whole field intensity is measured with one group of transducer.Its main thought is to calculate the sound field of other positions of depth direction according to the sound field extrapolation that observes from horizontal level (being depth direction first row).Specific algorithm comprises two main steps: the first step is carried out two-dimensional Fourier transform to time domain data, obtains 2-d spectrum,
S(k,ω)=2D-FT(s(u
n,t))
Second step was the circulation at depth direction: the 2-d spectrum that earlier the last time circulation is obtained is done phase shift, makes two-dimensional inverse Fourier transform then and gets t=0, obtains delegation's image,
s(u
n,t=0)=2D-IFT(S(k,ω)·α(Δz,k,ω))
Wherein,
is the pairing Phase shift amount of Δ z step-length on the depth direction; C is ultrasonic velocity of propagation in medium; Its value is constant, and therefore, what this ultrasonic imaging method used is the Phase shift technology of constant velocity of wave.
Can know that from the principle of DAS time domain SAFT algorithm need be confirmed time-delay to relative distance between the reverberation and velocity of wave according to transducer.And acoustic characteristics such as for multi-layer body, sound wave can reflect at the interphase place of layer and layer, reflection, thereby travel path can change.In addition, the medium of each layer of object is different, and sound wave is also unequal usually in the velocity of propagation of each layer; Therefore; The calculating of time delay can not be adopted the simple computation method of relative distance divided by velocity of wave again, but needs to obtain earlier the travel path of wave beam, obtains the travel-time then piecemeal.Its emphasis and difficult point are to find out quickly and accurately the travel path of wave beam, and the ray trace technology realizes the best approach of this purpose just.The principle of ray trace mainly based on Snell theorem or Fermat principle, is found out the shortest time-consuming acoustic wave propagation path through iterative computation.To be this technology do not have special requirement to the interphase shape of multi-layer body to the advantage of ray trace technology, can be applicable to the ultrasonic imaging of arbitrary shaped body.But its important disadvantages is a ray tracing algorithm comprises interative computation, and time complexity is high, and simultaneously, the DAS principle is actually a kind of convolution algorithm, and its calculated amount is also bigger, and imaging time is long.Shown in Fig. 2 (b), this scheme generates this image (5cm*5cm) just needs 2 hours.This shortcoming seriously limited should technology promote the use.
For the imaging of multi-layer body, the Phase shift ultrasonic imaging method does not have demand raypath and convolution algorithm, and image taking speed can significantly improve, and only needs 4 minutes as generating the image shown in Fig. 2 (c).And the Phase shift technology is a kind of image reconstructing method based on frequency domain, than the SAFT technology under the time domain have higher distance to the orientation of resolution and Geng Gao to resolution and lower secondary lobe, its image quality is higher.But; Because the velocity of wave of hypothesis sound wave is constant in the Phase shift technology, make this method to be applicable to have foreign medium on the depth direction and horizontal direction must be medium of the same race that promptly layer is the regular situation on level straight line/plane perhaps parallel to each other and so on interphase between the layer medium; And for the object that all has foreign medium on depth direction and the horizontal direction; The object that promptly comprises non-regular interphase (curve/curved surface) since sound wave the velocity of wave in the medium is different throughout, and cause testing result accurate inadequately; Shown in Fig. 2 (c), grave error has appearred in the interfacial imaging results of the second layer.
Therefore, aspect the ultrasonic imaging that contains the interfacial complicated layering object of non-rule, also lack effective method fast and accurately at present.
Summary of the invention
The objective of the invention is to propose a kind of ultrasonic imaging method, realize containing quick, accurate, the effectively imaging of the interfacial complicated layering object of non-rule.
The invention is characterized in, can or only contain the object image-forming of single medium regular layering object, simultaneously also can be to containing non-level and non-interfacial non-regular layering object image-forming parallel to each other, and image taking speed is fast.
The invention is characterized in; Adopt the Phase shift ultrasonic imaging technique; And the Phase shift technology improved, on the depth direction of testee and horizontal direction, considered to exist the situation of foreign medium respectively, in the Phase shift formula, introduce different velocities of wave; Propose to become the Phase shift technology of velocity of wave, so that this ultrasonic imaging method is applicable to non-regular layering object.One-piece construction is as shown in Figure 3; In imaging process to the separatrix individual processing; Outside the depth range at place, separatrix; Use the Phase shift formula and the Fast Fourier Transform (FFT) technology of constant velocity of wave, and in this scope, use Phase shift technology and N the single output Fast Fourier Transform (FFT) technology that becomes velocity of wave, to reduce calculated amount.
The invention is characterized in, contain following steps successively:
Step (1): make up one forms by a computing machine, ultrasonic transducer, a cover register control and analog to digital converter one being used for based on change velocity of wave Phase shift on the vertical section that the degree of depth and horizontal both direction form, make the system of not damaged ultrasonic imaging at horizontal and vertical three layers of heterogeneous object that all have a foreign medium, wherein:
Said ultrasonic transducer is provided with: the pulse signal input terminal that links to each other with the output terminal of said register control; The input end of said register control links to each other with the corresponding positioning control signal output part of said computing machine; Said ultrasonic transducer also is provided with: the echoed signal output terminal that links to each other with the input end of said analog to digital converter, the output terminal of said analog to digital converter links to each other with the echo samples signal input part of said computing machine.Said transducer is controlled by said register control, moves with 1 step-length/ms fixed rate at three layers of heterogeneous body surface, and said register control is the gearing of the said ultrasonic transducer of control shift position, and its parameter is imported by said computing machine,
Three layers of heterogeneous object are L along the axial horizontal length of x, and it is interval to be divided into x of L/ Δ, and Δ x is a burst length, and the point 0 that also is said ultrasonic transducer from the x axle is to terminal point N
x-1 ends each step-length that moves, N
x=1+L/ Δ x, said ultrasonic transducer move the point that is reached at every turn and are called sensing point, total N
xIndividual, sequence number n
x=0,1 ..., N
x-1; Said register control produces a TTL pulse at each sensing point place; Trigger said ultrasonic transducer driving pulse of depth direction z emission perpendicular to the x axle to three layers of heterogeneous object; Transducer transfers receiving mode to and picks up counting subsequently, receive from the echoed signal of testee reflection, said analog to digital converter to said transducer at sensing point n
xThe echoed signal that the place receives is carried out N
tInferior sampling is also stored in the computing machine sampling sequence number n into
t=0,1 ..., N
t-1, SF is f
s, its value is preset for analog to digital converter, and (z, x are to be x=n at horizontal ordinate t) to note s
xThe echoed signal that Δ x place receives is at t=n
t/ f
sSampled value constantly;
Step (2): said computing machine is from n
tSensing point n is read in=0 beginning in regular turn
xThe sampled value at=0 place, and write down the sampling sequence number that sampled value fluctuates for the first time
Subscript 0 expression sensing point sequence number then, repeats above-mentioned steps and reads n successively
x=1 ..., N
xThe sampled value at-1 each sensing point place, and write down the sequence number that each point place sampled value fluctuates for the first time
Obtain N
xN altogether on the individual sensing point
xThe sequence number that individual sampled value fluctuates for the first time
From
In choose minimum sequence number
And calculate this sequence number pairing depth coordinate value on the z of said three layers of heterogeneous object direction
With from
In choose maximum sequence number
And calculate this sequence number pairing depth coordinate value on the z of said three layers of heterogeneous object direction
V wherein
1Be the velocity of wave in the ground floor medium; Make two horizontal lines that are parallel to the x axle through these two depth coordinate values respectively; Be surrounded the separatrix between ground floor medium and second layer medium, form one and be positioned at first bounding box of crossing over first, second two layer medium on said three layers of heterogeneous object x-z vertical section simultaneously;
Step (3): (z, x t) carry out one dimension Fast Fourier Transform (FFT) 1D-FFT to the sampled value s at each sensing point place one by one along time shaft t
t(s (z, x, t)) gets the initial value z of z
0=0, obtain x-ω territory 2D signal:
S(z
0,x,ω)=1D-FFT
t(s(z
0,x,t))
Wherein, x is n for the sensing point sequence number
xThe level at place is to coordinate figure, and ω is the sampling angular frequency, then, and to above-mentioned x-ω territory 2D signal S (z
0, x ω) carries out one dimension Fast Fourier Transform (FFT) 1D-FFT along transverse axis x successively
x(S (z
0, x, ω)), obtain z
0The k-ω territory 2-d spectrum S (z at place
0, k, ω):
S(z
0,k,ω)=1D-FFT
x(1D-FFT
t(s(z
0,x,t)))
Wherein, k be the x of wave-number vector to component, sequence number is n
k, said wave-number vector is represented the spatial gradient of wave phase,
Wherein, n
ωBe the sequence number of sampling angular frequency, value is at 0~(N
t-1) between, gets t=t
0=0, one dimension sampled data s (z
0, x, t
0) be z=z on three layers of heterogeneous object vertical section
0The pixel value of=0 each pixel of row;
Step (4): generate z direction homogeneity district z successively according to the following steps
0<z≤z
MinInterior each row image that forms by each row pixel respectively:
Step (4.1): utilize following Phase shift formula to calculate said vertical section at z
Min2-d spectrum S (the z at place
Min, k, ω):
Wherein,
Wherein, i is an imaginary unit, z
0=0, ω is sampling angular frequency, v
zFor testee acoustic velocity in the capable medium of z on depth direction, at z≤z
MinThe time, v
z=v
1, v
1Be the velocity of wave of sound wave in the ground floor medium, Δ z is the spacing distance between last adjacent two pixels of depth direction z, α (Δ z, k, ω, v
z) be the Phase shift amount under the constant velocity of wave,
Step (4.2): carry out following steps successively:
Step (4.2.1): get t=t
0=0, the S (z that said step (4.1) is obtained
Min, k ω) makes the Gauss integration of independent variable ω,
Step (4.2.2): the Gauss integration result to said step (4.2.1) gets conjugation, again k is done the one dimension Fast Fourier Transform (FFT), multiply by 1/N after result of calculation is got conjugation again
x, obtain in the ground floor medium that coordinate figure is z on the depth direction
MinOne-row pixels put formed image:
At z direction ground floor homogeneity district z
0<z≤z
MinIn, each row image is all identical:
s(z,x,t
0)=s(z
min,x,t
0),z
0<z≤z
min;
Step (5): get the v in the step (4)
z=v
1, be that following step (5.1), step (5.2) are carried out in the step-length circulation to the z value with Δ z, up to z>=z
MaxEnd, generate depth direction z in first bounding box
Min<z≤z
MaxBy formed each the row image of each row pixel, its step is following respectively in the interval:
Step (5.1): utilize the 2-d spectrum at the Phase shift formula calculating z+ Δ z place under the said constant velocity of wave, z
Min<z+ Δ z≤z
Max:
S(z+Δz,k,ω)=S(z,k,ω)·α(Δz,k,ω,v
z),
Step (5.2): carry out following steps successively:
Step (5.2.1): get t=t
0=0, the S that said step (5.1) is obtained (z+ Δ z, k ω) make the Gauss integration of independent variable ω,
Step (5.2.2): the Gauss integration result to said step (5.2.1) gets conjugation, again k is done the one dimension Fast Fourier Transform (FFT), multiply by 1/N after result of calculation is got conjugation again
x, the one-row pixels that obtains coordinate figure on the depth direction and be z+ Δ z is put formed image, z
Min<z+ Δ z≤z
Max:
Step (6): revise depth direction z in said first bounding box according to the following steps
MinTo z
MaxInterval z
Min<z≤z
MaxIn each row image, the error that causes to eliminate the not homogeneity between the second layer and ground floor medium:
Step (6.1): the z that obtains with step (5)
MinTo z
MaxInterval image block uses Canny operator edge extracting algorithm to extract said z as input quantity
Min, z
MaxBetween ground floor medium and the separatrix c of second layer medium
1(x, z),
Step (6.2): from z=z
MinBeginning is up to z>=z
MaxOnly, the z that obtains with said step (4.1)
Min2-d spectrum S (the z at place
Min, k ω) makes initial value, is step-length circulation execution in step (6.2.1) and step (6.2.2) with Δ z:
Step (6.2.1): use single output Fast Fourier Transform Inverse (FFTI) method N time, be calculated as follows the x-ω territory 2D signal that obtains z+ Δ z place in said first bounding box:
Wherein, the value of N equals N
x, n
x=0,1 ..., N
x-1, n
k=0,1 ..., N
x-1, α (Δ z, k, ω, v
z(x)) for becoming the Phase shift amount of velocity of wave:
Wherein, the same step of the value of k, ω (3), v
z(x) be coordinate (x, the velocity of wave of sound wave in the medium of z) locating:
v
1For being positioned at separatrix c
1(x, the velocity of propagation of sound wave in the ground floor medium on z), v
2For being positioned at separatrix c
1(x, the velocity of propagation of sound wave in the second layer medium under z),
Calculation procedure is following:
Step (a): respectively to different n
x, n
k(n
x=0,1 ..., N
x-1; n
k=0,1 ..., N
x-1) calculates
Step (b): to different n
x(n
x=0,1 ..., N
x-1), with G (n
x, n
k), n
k=0,1 ..., N
x-1 sequence is as the input quantity of cutting output fft algorithm, to n
kMake G (n
x, n
k) the discrete fourier inverse transformation, but only get n
xThe value of this time domain point is the x-ω territory 2D signal at z+ Δ z place as coordinate figure on the depth direction:
Step (6.2.2): get t=t
0=0, (z+ Δ z, x ω) make the Gauss integration of independent variable ω to the S that said step (6.2.1) is obtained, and the one-row pixels that obtains coordinate figure on the depth direction and be z+ Δ z is put formed image, z
Min<z+ Δ z<z
Max:
Simultaneously, preserve S (z+ Δ z, x, ω) and do one-dimensional Fourier transform obtain S (z+ Δ z, k, ω), as the initial value of next round iteration;
Step (7): (5) described method set by step, get v
z=v
2, be that said step (5.1), step (5.2) are carried out in the step-length circulation to the z value with Δ z, up to z>=z
DepthEnd, generate depth direction z
Max<z≤z
DepthFormed each the row image of interval each row pixel, v
2Be the velocity of propagation of sound wave in the second layer medium, z
DepthRepresent the maximum coordinates value of three layers of heterogeneous object on depth direction z, be the systemic presupposition value;
Step (8): the z that obtains with step (7)
MaxTo z
DepthInterval image block uses Canny operator edge extracting algorithm to extract said z as input quantity
Max, z
DepthBetween second layer medium and the separatrix c of the 3rd layer of medium
2(x z), asks separatrix c
2(x, z) the min coordinates value z ' on depth direction z
MinWith maximum coordinates value z '
Max, respectively through said z '
MinAnd z '
MaxTwo depth coordinate values are made two horizontal lines that are parallel to the x axle; Be surrounded the separatrix between second layer medium and the 3rd layer of medium, form one and be positioned at second bounding box crossing over second, third two layer medium on said three layers of heterogeneous object x-z vertical section simultaneously;
Step (9): according to depth direction z ' in described method correction second bounding box of step (6.2)
MinTo z '
MaxInterval z '
Min<z≤z '
MaxIn each row image, the error that causes to eliminate the not homogeneity between the 3rd layer and second layer medium;
Step (10): (5) described method set by step, get v
z=v
3, to the z value with z '
MaxBeing initial value, is that said step (5.1), step (5.2) are carried out in the step-length circulation with Δ z, up to z>=z
DepthEnd, generate in the 3rd layer of medium except said second bounding box with formed each the row image of each row pixel of exterior portions, v
3It is the velocity of propagation of sound wave in the 3rd layer of medium.
The invention has the advantages that imaging precision is high, for the testee shown in Fig. 2 (a); If transducer diameter is 0.5mm, the transducer moving step length is 0.5mm, and transducer ultrasonic waves transmitted centre frequency is 5MHz; SF 100MHz; The step delta z of depth direction is taken as 0.05mm, and then ground floor interphase error can be controlled in the 0.4mm, and second layer interphase can be controlled at 1.0mm with interior (like Fig. 4).
The present invention has following advantage compared with prior art:
1. can carry out fast ultrasonic imaging to containing the interfacial object of non-rule;
2. the resolution and the precision of imaging are higher.
Description of drawings
Fig. 1 is the working model of SAFT: X to being the transducer direction of scanning, and Z is to being depth direction, u
nFor transducer with etc. the position of step-length when moving, also be the position of sensing point, (x ', z ') is the coordinate of target reflection thing, r
nDistance for target reflection thing to transducer.
The comparison diagram of Fig. 2 right and wrong rule layering object ultrasonic imaging method under the identical experiment environment: 2 (a) are the sectional view of testee; 2 (b) are that time domain SAFT and ray trace technology combines image that 2 (a) are generated, and imaging time is more than 2 hours, and to the unusual sensitivity of the parameter of algorithm, robustness is not fine; 2 (c) are the image that the Phase shift ultrasonic imaging technique of constant velocity of wave generates, imaging time 4 minutes, and its error is bigger; The image that 2 (d) generate for this formation method, imaging time 6 minutes.
Fig. 3 is this formation method overall process synoptic diagram.
Fig. 4 is the analysis of experimental data figure to Fig. 2 (d): 4 (a) are ground floor separatrix Error Graph; Solid-line curve is the marginal standard value of testee ground floor; The ground floor separatrix that the imaginary point curve obtains for this formation method, dashed curve are the marginal graph of errors of ground floor that this formation method obtained; 4 (b) are second layer separatrix Error Graph, and solid-line curve is the marginal standard value of the testee second layer, the second layer separatrix that the imaginary point curve obtains for this formation method, and dashed curve is the marginal graph of errors of the second layer that this formation method obtained.
Fig. 5 is this ultrasonic image-forming system schematic flow sheet.
Fig. 6 is this ultrasonic imaging hardware system structure figure.
Fig. 7 is a ultrasonic transducer work synoptic diagram.
Fig. 8 is based on the multi-layer body Non-Destructive Testing ultrasonic imaging algorithm flow chart that becomes velocity of wave Phase shift technology.
Fig. 9 asks the bounding box synoptic diagram in this ultrasonic imaging method.
FFT butterfly computation synoptic diagram when Figure 10 is N=8.
Embodiment
Practical implementation process of the present invention comprises three parts (like Fig. 5): ultrasound data obtains, forms images and calculates and the image demonstration.The hardware platform system structural drawing is as shown in Figure 6; Ultrasonic image-forming system is made up of a computing machine, ultrasonic transducer, a cover register control and an analog to digital converter; The pulse signal input terminal of ultrasonic transducer links to each other with the output terminal of register control, and the input end of register control links to each other with the positioning control signal output part of computing machine.The echoed signal output terminal of ultrasonic transducer links to each other with the input end of analog to digital converter, and the output terminal of analog to digital converter links to each other with the echo samples signal input part of computing machine.
The mode that ultrasound data obtains can be selected: contact, the i.e. direct contact measured object surfaces of ultrasonic transducer; Perhaps liquid immersion type, promptly object under test places liquid, and ultrasonic transducer is measured at liquid surface.According to the difference of detection mode, transducer can use contact probe or liquid immersion type probe.System uses single emission/receiving transducer, and transducer moves with about 1 step/ms fixed rate along x direction (like Fig. 7) with uniform step delta x in testee surface or liquid through register control.A moment generation TTL (transistor-transistor logic level) pulse that controller is stable in each target location; This pulse is used for triggering ultrasonic transducer to testee and the perpendicular driving pulse of depth direction emission of x direction; Transducer transfers receiving mode to and picks up counting subsequently, receives from the echo of testee reflection.Each position of transducer transponder pulse and reception echo is a sensing point.The echoed signal that transducer receives is by the analog to digital converter collection and be stored in the storer.Δ x needs to require to confirm that its value is more little that the image of generation is accurate more based on the actual size of object under test and imaging precision, but computing time is also long more.
Imaging is calculated and is imported as computing machine with the sampled data at testee each sensing point place on a vertical section exactly; Calculate the skiagraph picture of testee then by following image-forming step (process flow diagram is referring to Fig. 8); S (z wherein; X t) is illustrated in the signals sampling value that the transducer at horizontal direction x place receives at t constantly; Z representes the coordinate figure of depth direction; Δ z representes the step-length (relevant with precision) of depth direction, and the image that is promptly generated is in the spacing of vertical (that is depth direction) last adjacent two pixels; z
DepthBe preset image maximum depth value in the vertical:
Step (1): the ground floor medium is the isomorphism medium, and its ultrasonic reflection signal is identical, and at ground floor and second layer boundary, ultrasonic reflection signal can produce fluctuation.Then, search for the moment of each sensing point place generation signals fluctuation for the first time, and from these fluctuations constantly, find out the minimum t that constantly is worth according to the sampled data of ultrasound echo signal
MinBe worth t constantly with maximum
Max, then according to the velocity of wave v in the ground floor medium
1Ask marginal bounding box between ground floor and second layer medium, obtain two boundary line: the zs of bounding box on depth direction (be z to)
Min=v
1T
MinAnd z
Max=v
1T
Max(like Fig. 9);
Step (2): get depth direction initial value z=z
0=0, parameter x and t are carried out two-dimensional Fourier transform, original ultrasound echo signal is transformed to k-ω territory signal:
S(z
0,k,ω)=∫∫s(z
0,x,t)e
-ikxe
-iωtdxdt (1)
Wherein, ω is that the x of wave-number vector (wavenumber vector) is to component (wave-number vector refers to the spatial gradient of wave phase) for sampling angular frequency, k;
Step (3): generate depth direction z
0To z
MinInterval (z
0≤z≤z
Min) image.Because the medium of testee in this section interval is identical, so this section image is also just the same:
Step (3.1): utilize the Phase shift formula to calculate z
MinThe k-ω territory signal at place:
Wherein,
I is an imaginary unit, v
zFor testee acoustic velocity in the capable medium of z on depth direction, as z≤z
MinThe time, v
z=v
1, v
1Velocity of wave for sound wave in the ground floor medium;
Step (3.2): make two-dimensional inverse Fourier transform, make t=t
0=0, obtaining the along slope coordinate value is delegation's image of z:
Step (4): generate depth direction z
MinTo z
MaxInterval (z
Min<z≤z
Max) image.Get v
z=v
1, to the z value with Δ z step-length circulation execution in step (4.1) and step (4.2) up to z>=z
Max:
Step (4.1): utilize the Phase shift formula to calculate the k-ω territory signal at z+ Δ z place:
S(z+Δz,k,ω)=S(z,k,ω)α(Δz,k,ω,v
z) (4)
Step (4.2): make z ← z+ Δ z, t=t
0=0, make two-dimensional inverse Fourier transform, obtaining the along slope coordinate value is delegation's image of z:
Step (5): revise z
MinTo z
MaxInterval (z
Min<z≤z
Max) image.Because the velocity of wave of ground floor has only been used in the calculating in the step (4), soon should still be regarded as medium of the same race in the interval, the result of calculation that obtains certainly exists error, so this step needs to revise.Specifically comprised for two steps:
Step (5.1): the z that step (4) is obtained
MinTo z
MaxInterval image uses Canny edge extracting operator extraction z
MinWith z
MaxBetween separatrix c (x, z), z
Min≤z≤z
MaxThe computation process of Canny edge extracting operator was divided into for four steps:
Step (a): image smoothing.To original image with the two-dimensional Gaussian function smothing filtering obtain image I (x, z);
Step (b): computed image I (x, z) the gradient M of each pixel and direction Q.Adopt the first approximation of 2 * 2 templates conduct, promptly to the partial differential of x direction and z direction
Then gradient magnitude M and direction Q do
Step (c): the non-maximum value of gradient image suppresses.To each pixel I (x, z), if its Grad M (x, z) less than along its gradient direction Q (x, the gradient magnitude of two consecutive point z) explain that this point is not a marginal point, then with I (x, gray scale z) is made as 0;
Step (d): the dual threashold value is handled.Configure the dual threshold method detects and adjoining edge needs low threshold value Low and high threshold High, gradient image is carried out the dual threashold value handle.Gradient magnitude is the edge greater than high threshold High's; Gradient magnitude is not the edge less than low threshold value Low's; Gradient magnitude is marginal, then judges in eight neighbors of this pixel whether have the edge pixel greater than high threshold High, then is edge pixel if exist, otherwise is not;
Step (5.2): to the z value from z=z
MinBeginning, with Δ z step-length circulation execution in step (5.2.1) and step (5.2.2) up to z>=z
Max:
Step (5.2.1): utilize the Phase shift formula that becomes velocity of wave to calculate the k-ω territory signal at z+ Δ z place:
S(z+Δz,k,ω)=S(z,k,ω)α(Δz,k,ω,v
z(x)) (6)
Wherein, α (Δ z, k, ω, v
z(x)) for becoming the Phase shift amount of velocity of wave:
The value of vz (x) is following:
v
1, v
2Equal separatrix c (x, z) velocity of propagation of sound wave in the two layer medium up and down at place respectively;
Step (5.2.2): make z ← z+ Δ z, t=t
0=0, make two-dimensional inverse Fourier transform, obtaining the along slope coordinate value is delegation's image of z:
Step (6): generate depth direction z
MaxTo z
DepthInterval (z
Max<z≤z
Depth) image.Make v
z=v
2, the method for (4) is calculated z set by step
MaxTo z
DepthEach interval row image;
Step (7): handle follow-up separatrix.Image section (the z that step (6) is generated
Max<z≤z
Depth) (x z), makes v to use the new separatrix c of Canny edge extracting operator extraction
1, v
2Be respectively separatrix c (x, the z) velocity of propagation of sound wave in the two layer medium up and down at place, and make z
MinAnd z
Max(x, the z) up-and-down boundary on depth direction form new bounding box, if z to be respectively new separatrix c
Min==z
Max==z
Depth, then finish; Otherwise execution in step (5.2) is revised z
MinTo z
MaxInterval image;
Step (8): if z
Max==z
Depth, then finish; Otherwise, repeated execution of steps (6) and step (7).
Wherein, be the two-dimensional Fourier transform in the performing step (2) (like formula (1)) that the present invention adopts the method that parameter x and t are carried out twice one dimension Fast Fourier Transform (FFT) (1D-FFT) in order, with the raising computing velocity, promptly formula (1) is embodied as:
S(z
0,k,ω)=1D-FFT
x(1D-FFT
t(s(z
0,x,t)))
Need make two-dimensional inverse Fourier transform and get t=t in step (3.2) and the step (4.2)
0=0 value then can shilling t=t as view data
0=0, and ω made integration, and then k is made integration, to reduce the data volume of the second dimension integral operation.Formula (3) (5) turns to so:
Wherein, the integration
to ω can utilize numerical integration method (like Gauss Gaussian integration) to try to achieve fast.And the Fourier inversion of k is realized through one dimension Fast Fourier Transform (FFT) (1D-FFT) still the data that just earlier numerical integration obtained are got conjugation ([Integration
ω(S (z, k, ω))] *), directly utilize the 1D-FFT subroutine to carry out computing then, again operation result is got conjugation one time, and multiply by 1/N
x, can obtain.N
xFor transducer moves step number, i.e. N
x=1+L
x/ Δ x, wherein L
xBe the length of object under test in the x direction.Therefore, the concrete implementation procedure of formula (3) (5) is expressed as:
Phase shift amount α (Δ z, k, ω, v
z) need to each speed v
z, respectively different k, ω being calculated, its result is kept in the three-dimensional array in advance.Wherein, the value of k, ω is respectively:
Wherein, Δ x is the step-length that transducer moves along the x direction, f
sBe SF (its value is systemic presupposition), N
x=1+L
x/ Δ x is that transducer is counted L along the detection of x direction
x/ Δ x is that transducer moves step number, L
xBe the length (x direction) of testee, N
tFor transducer at each sensing point place the sampling number to echoed signal.Because the Phase shift ultrasonic imaging technique is regarded as the explosive reflection model with its working model, only consider that sound wave is from reverberation (along z axle negative direction) part of propagating upwards, so at Phase shift amount α (Δ z, k, ω, v
z) calculating in only get the value of ω<0 part, be 0 and make the Phase shift amount of ω>=0 part, simultaneously if b<0 then get
Know by the periodicity of ω value, when requiring ω<0, desirable N in the calculating
t/ 2≤n
ω<N
tPart.
In the step (5.2), become velocity of wave Phase shift amount α (Δ z, k, ω, v
z(x)) be the function of x, velocity of wave v
zNon-constant so step (5.2.1) can not directly be calculated again, and need merge in the Fourier inversion of step (5.2.2), promptly
Outer integration in the formula (8) can utilize numerical integration method (like Gauss Gaussian integration) to try to achieve, and internal layer integral operation (suc as formula (9)) then needs to try to achieve (suc as formula (10)) through revising discrete fourier inverse transformation computing formula,
In the formula (10), n
xBe the sequence number of sensing point on the x direction at transducer place, span is same as n
k, and α (Δ z, k, ω, v (
zX)) be that independent variable is n
kAnd n
xFunction, can not directly use discrete fourier inverse transformation (IDFT) or Fast Fourier Transform Inverse (FFTI) (IFFT) method.Therefore, use single output Fast Fourier Transform Inverse (FFTI) method (NsIFFTPruning) among the present invention N time.This method is based on the thought of butterfly computation in the Fast Fourier Transform (FFT), but a value in the transform domain is only exported in conversion each time.The NsIFFTPruning specific algorithm is following:
For?n
x=0?to?N
x-1
For?n
k=0?to?N
x-1
end
S(z+Δz,x=n
x·Δx,ω)=IFFTPruning(G,N
x,n
x)
end
Wherein, work as n
xIn the time of fixedly, function IFFTPruning (G, N
x, n
x) calculating G (n
x, n
k), n
k=0,1 ..., N
xThe discrete fourier inverse transformation of-1 sequence, but only get n
xThe value of this time domain point, promptly
Formula (11) can use IDFT or IFFT algorithm to calculate, but different with general IDFT problem is: formula (11) only need be calculated G (n
x, n
k) at n
xThe value at this time domain point place of=h if directly use IDFT or IFFT algorithm, then certainly exists a large amount of redundant computation, for example, and in Figure 10, if only need to calculate n
x=2 point, then the butterfly computation of dotted portion is redundant among the figure.In order to reduce calculated amount, the present invention uses the fft algorithm (IFFTPruning) of cutting output.This algorithm adopts the butterfly computation method still based on fft algorithm, and its principle is following: according to butterfly group technology in the fft algorithm, calculating is divided into H step (Stage), wherein H satisfies 2
H=N
xStep h (h=1 ..., H) in, with N
xPoint list entries G (n
x, n
k), n
k=0,1 ..., N
x-1 is equally divided into 2
H-hIndividual group, increase progressively 1 numbering b since 0 respectively in each group
h(b
h=0,1 ..., 2
h-1), in every group of butterfly computation except being numbered b
h=b
H+1Mod 2
H-1Point outside, remaining the point all be the redundancy.In each calculation procedure, only calculating is respectively organized and is numbered b in the butterfly computation
hThe point.Shown in figure 10, N
x=8, n
x=2 o'clock, calculate in every group of butterfly computation in the step (1) and be numbered b
1Calculate in every group of butterfly computation in=0 the point, step (2) and be numbered b
2Calculate in every group of butterfly computation in=2 the point, step (3) and be numbered b
3=2 point.
The picture traverse that calculates according to aforementioned image-forming step is N
x, but not the wide L of object under test
x, therefore, need carry out x to interpolation arithmetic in the final step of imaging calculation stages.Detailed process is: at x linear z-1 point of inst=Δ x/ Δ that insert between every adjacent two points.
Image shows that the two-dimensional image data that the calculation stages that is about to form images obtains is presented on the display, as required display gray scale image or coloured image.
To liking the situation of regular layering object, (x z) is straight line to the separatrix c of use Canny edge extracting operator extraction, and the imaging calculation procedure is constant in the image-forming step (5.1) for to be measured.
For the situation that only contains single medium in the object to be measured (be illustrated as contact ultrasonic measurement mode, otherwise should be regarded as regular layering object), then z
Min==z
Max==z
Depth, imaging promptly stops after calculating execution of step (3).
Claims (1)
1. based on the multi-layer body Non-Destructive Testing ultrasonic imaging method that becomes the velocity of wave Phase shift, it is characterized in that, contain following steps successively:
Step (1): make up one forms by a computing machine, ultrasonic transducer, a cover register control and analog to digital converter one being used for based on change velocity of wave Phase shift on the vertical section that the degree of depth and horizontal both direction form, make the system of not damaged ultrasonic imaging at horizontal and vertical three layers of heterogeneous object that all have a foreign medium, wherein:
Said ultrasonic transducer is provided with: the pulse signal input terminal that links to each other with the output terminal of said register control; The input end of said register control links to each other with the corresponding positioning control signal output part of said computing machine; Said ultrasonic transducer also is provided with: the echoed signal output terminal that links to each other with the input end of said analog to digital converter, the output terminal of said analog to digital converter links to each other with the echo samples signal input part of said computing machine.Said transducer is controlled by said register control, moves with 1 step-length/ms fixed rate at three layers of heterogeneous body surface, and said register control is the gearing of the said ultrasonic transducer of control shift position, and its parameter is imported by said computing machine,
Three layers of heterogeneous object are L along the axial horizontal length of x, and it is interval to be divided into x of L/ Δ, and Δ x is a burst length, and the point 0 that also is said ultrasonic transducer from the x axle is to terminal point N
x-1 ends each step-length that moves, N
x=1+L/ Δ x, said ultrasonic transducer move the point that is reached at every turn and are called sensing point, total N
xIndividual, sequence number n
x=0,1 ..., N
x-1; Said register control produces a TTL pulse at each sensing point place; Trigger said ultrasonic transducer driving pulse of depth direction z emission perpendicular to the x axle to three layers of heterogeneous object; Transducer transfers receiving mode to and picks up counting subsequently, receive from the echoed signal of testee reflection, said analog to digital converter to said transducer at sensing point n
xThe echoed signal that the place receives is carried out N
tInferior sampling is also stored in the computing machine sampling sequence number n into
t=0,1 ..., N
t-1, SF is f
s, its value is preset for analog to digital converter, and (z, x are to be x=n at horizontal ordinate t) to note s
xThe echoed signal that Δ x place receives is at t=n
t/ f
sSampled value constantly;
Step (2): said computing machine is from n
tSensing point n is read in=0 beginning in regular turn
xThe sampled value at=0 place, and write down the sampling sequence number that sampled value fluctuates for the first time
Subscript 0 expression sensing point sequence number then, repeats above-mentioned steps and reads n successively
x=1 ..., N
xThe sampled value at-1 each sensing point place, and write down the sequence number that each point place sampled value fluctuates for the first time
Obtain N
xN altogether on the individual sensing point
xThe sequence number that individual sampled value fluctuates for the first time
From
In choose minimum sequence number
And calculate this sequence number pairing depth coordinate value on the z of said three layers of heterogeneous object direction
With from
In choose maximum sequence number
And calculate this sequence number pairing depth coordinate value on the z of said three layers of heterogeneous object direction
V wherein
1Be the velocity of wave in the ground floor medium; Make two horizontal lines that are parallel to the x axle through these two depth coordinate values respectively; Be surrounded the separatrix between ground floor medium and second layer medium, form one and be positioned at first bounding box of crossing over first, second two layer medium on said three layers of heterogeneous object x-z vertical section simultaneously;
Step (3): (z, x t) carry out one dimension Fast Fourier Transform (FFT) 1D-FFT to the sampled value s at each sensing point place one by one along time shaft t
t(s (z, x, t)) gets the initial value z of z
0=0, obtain x-ω territory 2D signal:
S(z
0,x,ω)=1D-FFT
t(s(z
0,x,t))
Wherein, x is n for the sensing point sequence number
xThe level at place is to coordinate figure, and ω is the sampling angular frequency, then, and to above-mentioned x-ω territory 2D signal S (z
0, x ω) carries out one dimension Fast Fourier Transform (FFT) 1D-FFT along transverse axis x successively
x(S (z
0, x, ω)), obtain z
0The k-ω territory 2-d spectrum S (z at place
0, k, ω):
S(z
0,k,ω)=1D-FFT
x(1D-FFT
t(s(z
0,x,t)))
Wherein, k be the x of wave-number vector to component, sequence number is n
k, said wave-number vector is represented the spatial gradient of wave phase,
Wherein, n
ωBe the sequence number of sampling angular frequency, value is at 0~(N
t-1) between, gets t=t
0=0, one dimension sampled data s (z
0, x, t
0) be z=z on three layers of heterogeneous object vertical section
0The pixel value of=0 each pixel of row;
Step (4): generate z direction homogeneity district z successively according to the following steps
0<z≤z
MinInterior each row image that forms by each row pixel respectively:
Step (4.1): utilize following Phase shift formula to calculate said vertical section at z
Min2-d spectrum S (the z at place
Min, k, ω):
Wherein,
Wherein, i is an imaginary unit, z
0=0, ω is sampling angular frequency, v
zFor testee acoustic velocity in the capable medium of z on depth direction, at z≤z
MinThe time, v
z=v
1, v
1Be the velocity of wave of sound wave in the ground floor medium, Δ z is the spacing distance between last adjacent two pixels of depth direction z, α (Δ z, k, ω, v
z) be the Phase shift amount under the constant velocity of wave,
Step (4.2): carry out following steps successively:
Step (4.2.1): get t=t
0=0, the S (z that said step (4.1) is obtained
Min, k ω) makes the Gauss integration of independent variable ω,
Step (4.2.2): the Gauss integration result to said step (4.2.1) gets conjugation, again k is done the one dimension Fast Fourier Transform (FFT), multiply by 1/N after result of calculation is got conjugation again
x, obtain in the ground floor medium that coordinate figure is z on the depth direction
MinOne-row pixels put formed image:
At z direction ground floor homogeneity district z
0<z≤z
MinIn, each row image is all identical:
s(z,x,t
0)=s(z
min,x,t
0),z
0<z≤z
min;
Step (5): get the v in the step (4)
z=v
1, be that following step (5.1), step (5.2) are carried out in the step-length circulation to the z value with Δ z, up to z>=z
MaxEnd, generate depth direction z in first bounding box
Min<z≤z
MaxBy formed each the row image of each row pixel, its step is following respectively in the interval:
Step (5.1): utilize the 2-d spectrum at the Phase shift formula calculating z+ Δ z place under the said constant velocity of wave, z
Min<z+ Δ z≤z
Max:
S(z+Δz,k,ω)=S(z,k,ω)·α(Δz,k,ω,v
z),
Step (5.2): carry out following steps successively:
Step (5.2.1): get t=t
0=0, the S that said step (5.1) is obtained (z+ Δ z, k ω) make the Gauss integration of independent variable ω,
Step (5.2.2): the Gauss integration result to said step (5.2.1) gets conjugation, again k is done the one dimension Fast Fourier Transform (FFT), multiply by 1/N after result of calculation is got conjugation again
x, the one-row pixels that obtains coordinate figure on the depth direction and be z+ Δ z is put formed image, z
Min<z+ Δ z≤z
Max:
Step (6): revise depth direction z in said first bounding box according to the following steps
MinTo z
MaxInterval z
Min<z≤z
MaxIn each row image, the error that causes to eliminate the not homogeneity between the second layer and ground floor medium:
Step (6.1): the z that obtains with step (5)
MinTo z
MaxInterval image block uses Canny operator edge extracting algorithm to extract said z as input quantity
Min, z
MaxBetween ground floor medium and the separatrix c of second layer medium
1(x, z),
Step (6.2): from z=z
MinBeginning is up to z>=z
MaxOnly, the z that obtains with said step (4.1)
Min2-d spectrum S (the z at place
Min, k ω) makes initial value, is step-length circulation execution in step (6.2.1) and step (6.2.2) with Δ z:
Step (6.2.1): use single output Fast Fourier Transform Inverse (FFTI) method N time, be calculated as follows the x-ω territory 2D signal that obtains z+ Δ z place in said first bounding box:
Wherein, the value of N equals N
x, n
x=0,1 ..., N
x-1, n
k=0,1 ..., N
x-1, α (Δ z, k, ω, v
z(x)) for becoming the Phase shift amount of velocity of wave:
Wherein, the same step of the value of k, ω (3), v
z(x) be coordinate (x, the velocity of wave of sound wave in the medium of z) locating:
v
1For being positioned at separatrix c
1(x, the velocity of propagation of sound wave in the ground floor medium on z), v
2For being positioned at separatrix c
1(x, the velocity of propagation of sound wave in the second layer medium under z),
Calculation procedure is following:
Step (a): respectively to different n
x, n
k(n
x=0,1 ..., N
x-1; n
k=0,1 ..., N
x-1) calculates
Step (b): to different n
x(n
x=0,1 ..., N
x-1), with G (n
x, n
k), n
k=0,1 ..., N
x-1 sequence is as the input quantity of cutting output fft algorithm, to n
kMake G (n
x, n
k) the discrete fourier inverse transformation, but only get n
xThe value of this time domain point is the x-ω territory 2D signal at z+ Δ z place as coordinate figure on the depth direction:
Step (6.2.2): get t=t
0=0, (z+ Δ z, x ω) make the Gauss integration of independent variable ω to the S that said step (6.2.1) is obtained, and the one-row pixels that obtains coordinate figure on the depth direction and be z+ Δ z is put formed image, z
Min<z+ Δ z<z
Max:
Simultaneously, preserve S (z+ Δ z, x, ω) and do one-dimensional Fourier transform obtain S (z+ Δ z, k, ω), as the initial value of next round iteration;
Step (7): (5) described method set by step, get v
z=v
2, be that said step (5.1), step (5.2) are carried out in the step-length circulation to the z value with Δ z, up to z>=z
DepthEnd, generate depth direction z
Max<z≤z
DepthFormed each the row image of interval each row pixel, v
2Be the velocity of propagation of sound wave in the second layer medium, z
DepthRepresent the maximum coordinates value of three layers of heterogeneous object on depth direction z, be the systemic presupposition value;
Step (8): the z that obtains with step (7)
MaxTo z
DepthInterval image block uses Canny operator edge extracting algorithm to extract said z as input quantity
Max, z
DepthBetween second layer medium and the separatrix c of the 3rd layer of medium
2(x z), asks separatrix c
2(x, z) the min coordinates value z ' on depth direction z
MinWith maximum coordinates value z '
Max, respectively through said z '
MinAnd z '
MaxTwo depth coordinate values are made two horizontal lines that are parallel to the x axle; Be surrounded the separatrix between second layer medium and the 3rd layer of medium, form one and be positioned at second bounding box crossing over second, third two layer medium on said three layers of heterogeneous object x-z vertical section simultaneously;
Step (9): according to depth direction z ' in described method correction second bounding box of step (6.2)
MinTo z '
MaxInterval z '
Min<z≤z '
MaxIn each row image, the error that causes to eliminate the not homogeneity between the 3rd layer and second layer medium;
Step (10): (5) described method set by step, get v
z=v
3, to the z value with z '
MaxBeing initial value, is that said step (5.1), step (5.2) are carried out in the step-length circulation with Δ z, up to z>=z
DepthEnd, generate in the 3rd layer of medium except said second bounding box with formed each the row image of each row pixel of exterior portions, v
3It is the velocity of propagation of sound wave in the 3rd layer of medium.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210046151.1A CN102608205B (en) | 2012-02-24 | 2012-02-24 | Multilayer-object nondestructive testing ultrasonic imaging method based on variable wave speed phase shifting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210046151.1A CN102608205B (en) | 2012-02-24 | 2012-02-24 | Multilayer-object nondestructive testing ultrasonic imaging method based on variable wave speed phase shifting |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102608205A true CN102608205A (en) | 2012-07-25 |
CN102608205B CN102608205B (en) | 2014-10-22 |
Family
ID=46525747
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210046151.1A Expired - Fee Related CN102608205B (en) | 2012-02-24 | 2012-02-24 | Multilayer-object nondestructive testing ultrasonic imaging method based on variable wave speed phase shifting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102608205B (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103018333A (en) * | 2012-12-07 | 2013-04-03 | 清华大学 | Synthetic aperture focused ultrasonic imaging method of layered object |
CN103033816A (en) * | 2012-12-07 | 2013-04-10 | 清华大学 | Synthetic aperture focused ultrasonic imaging implementation method based on arc scanning transition |
CN103126721A (en) * | 2013-03-05 | 2013-06-05 | 飞依诺科技(苏州)有限公司 | Method and device for compound ultrasonic imaging based on multiple focuses and multiple ultrasonic velocities |
WO2014086322A1 (en) * | 2012-12-07 | 2014-06-12 | Tsinghua University | Ultrasonic imaging method and system for detecting multi-layered object with synthetic aperture focusing technique |
CN105157741A (en) * | 2015-08-31 | 2015-12-16 | 西安科技大学 | Rapid determination method of spatial pulse response of circular transducer |
CN108020268A (en) * | 2018-01-19 | 2018-05-11 | 河海大学常州校区 | Transceiver ultrasonic probe dielectric stratifying property detection system |
CN108606811A (en) * | 2018-04-12 | 2018-10-02 | 上海交通大学医学院附属上海儿童医学中心 | A kind of ultrasound stone age detecting system and its method |
CN113177992A (en) * | 2021-05-18 | 2021-07-27 | 清华大学 | Efficient synthetic aperture ultrasonic imaging method |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6519532B2 (en) * | 2001-01-31 | 2003-02-11 | Phillips Petroleum Company | Method and apparatus for 3D depth migration |
-
2012
- 2012-02-24 CN CN201210046151.1A patent/CN102608205B/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6519532B2 (en) * | 2001-01-31 | 2003-02-11 | Phillips Petroleum Company | Method and apparatus for 3D depth migration |
Non-Patent Citations (3)
Title |
---|
LIANJIE HUANG ET AL.: "Ultrasonic breast imaging using a wave-equation migration method", 《MEDICAL IMAGING 2003: ULTRASONIC IMAGING AND SIGNAL PROCESSING》 * |
MARTIN H. SKJE1VAREID ET AL.: "Ultrasonic Imaging of Pitting using Multilayer Synthetic Aperture Focusing", 《2011 IEEE INTERNATIONAL ULTRASONICS SYMPOSIUM PROCEEDINGS》 * |
TOMAS OLOFSSON: "Phase Shift Migration for Imaging Layered Objects and Objects Immersed in Water", 《IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103018333A (en) * | 2012-12-07 | 2013-04-03 | 清华大学 | Synthetic aperture focused ultrasonic imaging method of layered object |
CN103033816A (en) * | 2012-12-07 | 2013-04-10 | 清华大学 | Synthetic aperture focused ultrasonic imaging implementation method based on arc scanning transition |
CN103033816B (en) * | 2012-12-07 | 2014-06-04 | 清华大学 | Synthetic aperture focused ultrasonic imaging implementation method based on arc scanning transition |
WO2014086322A1 (en) * | 2012-12-07 | 2014-06-12 | Tsinghua University | Ultrasonic imaging method and system for detecting multi-layered object with synthetic aperture focusing technique |
CN103018333B (en) * | 2012-12-07 | 2014-10-29 | 清华大学 | Synthetic aperture focused ultrasonic imaging method of layered object |
CN103126721A (en) * | 2013-03-05 | 2013-06-05 | 飞依诺科技(苏州)有限公司 | Method and device for compound ultrasonic imaging based on multiple focuses and multiple ultrasonic velocities |
CN105157741A (en) * | 2015-08-31 | 2015-12-16 | 西安科技大学 | Rapid determination method of spatial pulse response of circular transducer |
CN108020268A (en) * | 2018-01-19 | 2018-05-11 | 河海大学常州校区 | Transceiver ultrasonic probe dielectric stratifying property detection system |
CN108606811A (en) * | 2018-04-12 | 2018-10-02 | 上海交通大学医学院附属上海儿童医学中心 | A kind of ultrasound stone age detecting system and its method |
CN113177992A (en) * | 2021-05-18 | 2021-07-27 | 清华大学 | Efficient synthetic aperture ultrasonic imaging method |
CN113177992B (en) * | 2021-05-18 | 2022-06-10 | 清华大学 | Efficient synthetic aperture ultrasonic imaging method |
Also Published As
Publication number | Publication date |
---|---|
CN102608205B (en) | 2014-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102608205B (en) | Multilayer-object nondestructive testing ultrasonic imaging method based on variable wave speed phase shifting | |
CN103033816B (en) | Synthetic aperture focused ultrasonic imaging implementation method based on arc scanning transition | |
Huthwaite | Evaluation of inversion approaches for guided wave thickness mapping | |
Skjelvareid et al. | Synthetic aperture focusing of ultrasonic data from multilayered media using an omega-k algorithm | |
CN106501367B (en) | Implementation method is imaged in phased array supersonic based on elliptic arc scan transformation | |
CN102770079B (en) | Ultrasonic imaging apparatus and method of controlling delay | |
EP2115386B1 (en) | Ultrasonic surface monitoring | |
KR101581369B1 (en) | Imaging method and apparatus using shear waves | |
JP4776707B2 (en) | Ultrasonic imaging device | |
CN111122700B (en) | Method for improving laser ultrasonic SAFT defect positioning speed | |
CN103018333B (en) | Synthetic aperture focused ultrasonic imaging method of layered object | |
Merabet et al. | 2-D and 3-D reconstruction algorithms in the Fourier domain for plane-wave imaging in nondestructive testing | |
JP2009532089A (en) | Delay controller for ultrasonic receiving beamformer | |
CN104688224B (en) | One kind is applied to the non-homogeneous medium magnetosonic coupling imaging method for reconstructing of acoustics | |
Skjelvareid | Synthetic aperture ultrasound imaging with application to interior pipe inspection | |
CN103654732B (en) | A kind of Photoacoustic image optimization method based on linear delay compensation | |
CN113109443A (en) | Focusing acoustic array imaging method and system | |
EP2279407B1 (en) | Ultrasonic modelling | |
CN113552571B (en) | Underwater laser induced acoustic SAFT imaging method based on PSM algorithm | |
Yang et al. | Real-time ultrasonic imaging for multi-layered objects with synthetic aperture focusing technique | |
JP2019033822A (en) | Ultrasound signal processor, ultrasound diagnostic apparatus, and ultrasound signal processing method | |
Cui et al. | Real-time total focusing method imaging for ultrasonic inspection of three-dimensional multilayered media | |
Dziewierz et al. | A design methodology for 2D sparse NDE arrays using an efficient implementation of refracted-ray TFM | |
JP2017500553A (en) | How to rebuild the surface of a fragment | |
CN116171382A (en) | Method for detecting a discontinuity and system for implementing said method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141022 Termination date: 20190224 |
|
CF01 | Termination of patent right due to non-payment of annual fee |