CN106680825A - Acoustic array imaging system and method thereof - Google Patents

Acoustic array imaging system and method thereof Download PDF

Info

Publication number
CN106680825A
CN106680825A CN201611105251.1A CN201611105251A CN106680825A CN 106680825 A CN106680825 A CN 106680825A CN 201611105251 A CN201611105251 A CN 201611105251A CN 106680825 A CN106680825 A CN 106680825A
Authority
CN
China
Prior art keywords
array
array element
image
signal
imaging
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201611105251.1A
Other languages
Chinese (zh)
Other versions
CN106680825B (en
Inventor
吴文焘
张云翼
曹政
韩晓丽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Acoustics CAS
Original Assignee
Institute of Acoustics CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Acoustics CAS filed Critical Institute of Acoustics CAS
Priority to CN201611105251.1A priority Critical patent/CN106680825B/en
Publication of CN106680825A publication Critical patent/CN106680825A/en
Application granted granted Critical
Publication of CN106680825B publication Critical patent/CN106680825B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range 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
    • G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Abstract

The method relates to an acoustic array imaging system and a method thereof. The system comprises a data acquisition module 101, a first imaging module 102, a second imaging module 103 and an image composite module 104. The method comprises the following steps of successively emitting array elements in an acoustic array and collecting echo signals so as to acquire an analysis signal; according to a sound propagation relation between the emitted array elements and receiving array elements in an array, acquiring a time-delay value, and according to the analysis signal and the time-delay value, calculating images of the receiving array elements; using a cascaded parallel composite processing unit to recombine the image of all the receiving array elements; and superposing all the composite images and carrying out envelope processing so as to acquire a full-matrix data acquisition imaging result. The system and the method have advantages that when each array element data is collected, processing is performed so that an occupied storage space is small; a parallel mode is used to increase an operation speed; during an actual application, calculating time consuming is less, and high frame rate imaging can be realized; and the method and the system are suitable for real-time data processing and high resolution imaging. In addition, the system and the method are suitable for acoustic array imaging under different frequency and application scenes.

Description

A kind of acoustic array imaging system and method
Technical field
The present invention relates to acoustic array technical field of imaging, more particularly to a kind of distributed parallel complete matrix acoustic array into As System and method for.
Background technology
At present, acoustic array imaging is for industrial nondestructive testing, medical science and under water the industry such as detection all has important meaning Justice, was widely applied in the last few years.
Supersonic array detection technique is based on supersonic array transducer, using ultrasonic longitudinal wave or shear wave, in industry zero Part carries out a kind of method of Non-Destructive Testing.Supersonic array detection technique can also be used to detect such as composite complex structure Material.The method has great flexibility, can detect irregularly shaped part, relative to traditional single array element flaw detection side Formula, has the advantages that to detect that area is big, signal to noise ratio is high, testing result is directly perceived.
The conventional supersonic array detection method of field of industry detection includes that direct contact type phased array imaging and Jing voussoirs are rolled over Formula phased array imaging is penetrated, different imaging methods are applied to different application scenarios.
The content of the invention
It is an object of the present invention in solve prior art, in total focus imaging method, each reception process is entirely Aperture collection data, this process needs to take substantial amounts of memory space, and because end processing sequences calculate time-consuming big, causes Frame per second is low, does not apply to real time imagery.
For achieving the above object, on the one hand, the invention provides a kind of acoustic array imaging system, the system includes：Number According to acquisition module 101, the first image-forming module 102, the second image-forming module 103 and image composite module 104.
Wherein, data acquisition module 101, for launching ultrasonic wave successively to each array element in acoustic array, and adopt The echo-signal of each array element in collection acoustic array, analytic trnasformation process is carried out respectively to echo-signal, obtains corresponding Analytic signal.Wherein, data acquisition module 101 specifically for：The echo letter of each array element is obtained using Hilbert transform Number corresponding analytic signal, the Hilbert transform formula is：
A (t)=x (t)+iH { x (t) }
Wherein H { } represents Hilbert transform, and x (t) represents echo-signal, and Hilbert transform integrated form is：
Reilly narrow-long slot parts can also be adopted to obtain the corresponding analytic signal of echo-signal of each array element.
First image-forming module 102, for according to the physical geometry relation in acoustic array between each array element or passing through The acoustic propagation relation of each array element is calculated during medium refraction, each is obtained according to acoustic propagation relation and is received array element passage every The delay value of one pixel, obtains each reception array element corresponding according to the analytic signal and delay value calculating for obtaining Image.
Wherein, the first image-forming module 102 specifically adopts following form：According to the analytic signal and delay value that obtain, adopt Distributed Calculation obtains each and receives the corresponding image of array element.Distributed Calculation carries out pyrene process using polycaryon processor, Pyrene process includes：Realize that one or more single pass images are generated using a processor core, then using multinuclear or many The all receiving channels of the parallel calculating of processor；According to the analytic signal and delay value that obtain, obtain every using Distributed Calculation One receives the corresponding image of array element, and Distributed Calculation is processed using circuit module.
Second image-forming module 103, is combined to each corresponding image of reception array element, obtains each transmitting array element Array received combination picture.
Wherein, the second image-forming module 103 adopts tree-like tandem type Combined Processing unit to the corresponding figure of each array element As being combined, array received combination picture is obtained.
Image composite module 104, using array received composite diagram of the parallel composition processing unit to each transmitting array element As being combined, total focus image is obtained, then total focus image is taken into image envelope, obtain complete matrix data acquisition imaging knot Really.
On the other hand, the invention provides a kind of acoustic array imaging method, the method is comprised the following steps：To acoustic array Each array element in row launches successively ultrasonic wave, and gathers the echo-signal of each array element in acoustic array, to echo Signal carries out respectively analytic trnasformation process, obtains corresponding analytic signal.
Wherein, the step of carrying out analytic trnasformation process, the corresponding analytic signal of acquisition respectively to echo-signal includes：Adopt The corresponding analytic signal Hilbert transform formula of echo-signal for obtaining each array element with Hilbert transform is：
A (t)=x (t)+iH { x (t) }
Wherein H { } represents Hilbert transform, and x (t) represents echo-signal, and Hilbert transform integrated form is：
Reilly narrow-long slot parts can be also adopted to obtain the corresponding analytic signal of echo-signal of each array element.
Each is calculated during according to the physical geometry relation in acoustic array between each array element or by medium refraction The acoustic propagation relation of array element, obtains each and receives delay value of the array element passage in each pixel according to acoustic propagation relation, Analytic signal and delay value according to obtaining is calculated and obtains each corresponding image of reception array element.
Wherein, according to the analytic signal and delay value for obtaining, each is obtained using Distributed Calculation and receives array element correspondence Image.Distributed Calculation carries out pyrene process using polycaryon processor, and pyrene process includes：Realized using a processor core One or more single pass images are generated, then using all receiving channels of calculating of multinuclear or parallel multiprocessor；According to The analytic signal for obtaining and delay value, obtain each and receive the corresponding image of array element, Distributed Calculation using Distributed Calculation Processed using circuit module.
Each corresponding image of reception array element is combined, the array received composite diagram of each transmitting array element is obtained Picture.
Wherein, complex method is entered using tree-like tandem type Combined Processing unit to each corresponding image of reception array element Row is compound.
All array received combination pictures for obtaining are overlapped, total focus image is obtained；Total focus image is taken again Image envelope, obtains complete matrix data acquisition imaging result.
The invention provides a kind of acoustic array imaging system and method, the present invention has the beneficial effect that:Gathering each battle array Processed in real time during metadata, so the memory space for taking is few, improve calculating in arithmetic speed, practical application using parallel Time-consuming little, achievable high frame per second imaging, it is adaptable to real time data processing and high-resolution imaging.Additionally, the present invention is applied to not Acoustic array imaging under same frequency and application scenarios.
Description of the drawings
Fig. 1 (a) is the schematic diagram of existing supersonic array compressional wave scanning imagery detection；
Fig. 1 (b) is the schematic diagram of existing supersonic array shear wave scanning imagery detection；
Fig. 2 is existing total focus imaging method transceiver mode schematic diagram；
Fig. 3 is existing total focus imaging transmitting, the space geometry relation schematic diagram of reception；
Fig. 4 is existing total focus imaging method and SAFT wave beam simulation comparison schematic diagrames；
Fig. 5 is a kind of distributed acoustics array imaging system function structure chart provided in an embodiment of the present invention；
Fig. 6 is a kind of distributed acoustics array imaging system structure chart provided in an embodiment of the present invention；
Fig. 7 is a kind of distributed acoustics array imaging system schematic diagram provided in an embodiment of the present invention；
Fig. 8 realizes block diagram for a kind of first image-forming module provided in an embodiment of the present invention；
Fig. 9 is a kind of cascade formula list array element transmitting of second image-forming module provided in an embodiment of the present invention 103-array received Image combination module mechanism map；
Figure 10 is a kind of distributed acoustics array image-forming method flow diagram provided in an embodiment of the present invention；
Figure 11 (a) is that a kind of 1st array element provided in an embodiment of the present invention launches the image that the-the 1 array element is received；
Figure 11 (b) is that a kind of 1st array element provided in an embodiment of the present invention launches the image that the-the 16 array element is received；
Figure 11 (c) is that a kind of 1st array element provided in an embodiment of the present invention launches the image that the-the 32 array element is received；
Figure 12 (a) is a kind of image of 1st array element transmitting-array received provided in an embodiment of the present invention；
Figure 12 (b) is a kind of image of 16th array element transmitting-array received provided in an embodiment of the present invention；
Figure 12 (c) is a kind of image of 32nd array element transmitting-array received provided in an embodiment of the present invention；
Figure 13 is a kind of complete matrix data acquisition image provided in an embodiment of the present invention.
Specific embodiment
Below by drawings and Examples, technical scheme is described in further detail.
Fig. 1 sets forth direct contact type sector scan and Jing voussoirs produce shear wave S type scanning schematic diagrames.
Fig. 1 (a) is a kind of schematic diagram of supersonic array compressional wave scanning imagery detection.As shown in figure 1, electronic system is by control The time delay of the different array elements of system, so as to change the angle of incident longitudinal wave, reaches the purpose to sector region imaging, and this mode is referred to as B-mode is imaged.
Fig. 1 (b) is a kind of schematic diagram of supersonic array shear wave scanning imagery detection provided in an embodiment of the present invention.Such as Fig. 1 Shown in (b), the time delay of the different array elements of electronic system control, and incident longitudinal wave is mapped on part along voussoir inside is oblique, from And the purpose being imaged to sector region using refracted shear is reached, this mode is referred to as S mode imaging.S mode imaging needs to calculate Ultrasound is refracted into the sound path of test block by voussoir, is calculated using Fermat's theorem.
Fig. 2 is a kind of total focus imaging method transceiver mode schematic diagram.As shown in Fig. 2 in emission mode, sensor array Single array element in row is launched successively, in a receive mode, is sent using the single array element of whole array receiveds in sensor array Echo data, complete matrix data acquisition contains between any two array element transmitting and the syntagmatic for receiving, for composite diagram Picture.
Fig. 3 is a kind of total focus imaging transmitting, the space geometry relation schematic diagram of reception.As shown in figure 3, straight with the first Connect and be introduced as a example by contact imaging pattern, calculate the space geometry relation for focusing on time delay, to space point P, if n-th Array element time delay is penetrated for τn, it is τ that m-th receives array element time delaym, the signal of collection is sm,nT (), total array element number is N, then total focus Imaging method is expressed as：
Wherein
Using frequency-domain analysis method, when array element launches single-frequency signalsWhen, total focus imaging is represented by：
By τmAnd τnBring formula (4) abbreviation into, save the synthesis wave beam after constant and be expressed as：
Obtain simulation architecture such as Fig. 4.Wherein θ represents beam angle.PSF(θ) beam modes are represented.
Fig. 4 is a kind of total focus imaging method and SAFT wave beam simulation comparison schematic diagrames.As shown in Figure 4, it can be seen that complete Focal imaging method is superior to SAFT to the inhibition of graing lobe and secondary lobe.Therefore total focus imaging method is frequently as image quality " golden standard ".But every time reception process is full aperture gathered data, it is therefore desirable to take substantial amounts of memory space；Simultaneously It is time-consuming big due to calculating, cause frame per second low, it is not suitable for real-time processing.Therefore, need badly a kind of new processing method of proposition and System architecture, reduces memory space, improves arithmetic speed, solves the above problems.
Wherein, synthetic aperture imaging method (SAFT) and conventional Ultrasound imaging and the synthetic aperture imaging method of other modes Compare, it is possible to increase frame per second, signal to noise ratio and contrast.SAFT launches spherical wave first with single array element, covers whole imaging Region；Then full array signal is received, is superimposed using time delay, be calculated low resolution ultrasonoscopy.In recent years, with regard to synthesis Research in terms of the imaging of transmitting aperture concentrates on and how to reduce complexity, improves the side such as frame per second, investigation depth and lateral resolution Face.
Fig. 5 is a kind of distributed acoustics array imaging system function structure chart provided in an embodiment of the present invention.Such as Fig. 5 institutes Show, a kind of acoustic array imaging system includes：Data acquisition module 101, the first image-forming module 102, the and of the second image-forming module 103 Image composite module 104.
Data acquisition module 101, for launching ultrasonic wave successively to each array element in acoustic array, and gathers acoustics The echo-signal of each array element in array, analytic trnasformation process is carried out respectively to echo-signal, obtains corresponding parsing letter Number.
Wherein, data acquisition module 101 specifically for：The echo-signal of each array element is obtained using Hilbert transform Corresponding analytic signal, Hilbert transform formula is：
A (t)=x (t)+iH { x (t) } (6)
Wherein H { } represents Hilbert transform, and x (t) represents echo-signal, and Hilbert transform integrated form is：
Reilly narrow-long slot parts can also be adopted to obtain the corresponding analytic signal of echo-signal of each array element.
First image-forming module 102, for according to the physical geometry relation in acoustic array between each array element or passing through The acoustic propagation relation of each array element is calculated during medium refraction, each is obtained according to acoustic propagation relation and is received array element passage every The delay value of one pixel, calculates according to the analytic signal and delay value that obtain and obtains each corresponding figure of reception array element Picture.
Wherein, the first image-forming module 102 is adopted and divided specifically using following form according to the analytic signal and delay value for obtaining Cloth is calculated and obtains each corresponding image of reception array element, and Distributed Calculation carries out pyrene process using polycaryon processor, point Core process includes：Realize that one or more single pass images are generated using a processor core, then using multinuclear or many places The parallel all receiving channels of calculating of reason device；According to the analytic signal and delay value that obtain, obtain each using Distributed Calculation Corresponding image in individual reception array element, Distributed Calculation is processed using circuit module.
Second image-forming module 103, is combined to each corresponding image of reception array element, obtains each transmitting array element Array received combination picture.
Wherein, the second image-forming module 103 is corresponding to each reception array element using tree-like tandem type Combined Processing unit Image be combined, obtain array received combination picture.
Image composite module 104, using array received composite diagram of the parallel composition processing unit to each transmitting array element As being combined, total focus image is obtained, then total focus image is taken into image envelope, obtain complete matrix data acquisition imaging knot Really.
Fig. 6 is a kind of distributed acoustics array imaging system structure chart provided in an embodiment of the present invention.As shown in fig. 6, number Include supersonic array transducer according to acquisition module, the individual array elements of N (N is positive integer) of supersonic array transducer receive supersonic array and change First array element of energy device launches the echo-signal of ultrasonic wave, and processing system is narrow using Hilbert or Reilly to echo-signal Method with algorithm obtains the analytic signal of echo-signal.
Wherein, analytic signal is obtained using Hilbert transform：
A (t)=x (t)+iH { x (t) } (6)
Wherein H { } represents Hilbert transform, and x (t) represents echo-signal, and Hilbert transform integrated form is：
The real-time implementation of Hilbert transform is more difficult, and the present invention is provided using improvement Reilly Narrow bands.Setting It is bandpass filter according to the passband design of ultrasonic transducer during meter ptototype filter, then obtains band-pass form using frequency displacement Reilly wave filters.Can draw from theory analysis, due to band logical effect, reducing band and making an uproar outward, signal to noise ratio more directly designs solution Analysis filtered method has some improvement.
Calculate during according to the physical geometry relation between first array element and each array element or by medium refraction every The acoustic propagation relation of one array element, obtains each and receives time delay of the array element passage in each pixel according to acoustic propagation relation Value, obtains each and receives the corresponding image of array element according to the analytic signal and delay value for obtaining by Distributed Calculation.Utilize All of single array element transmitting-mono- array element receives image, by compound, can both obtain single array element transmitting-array received composite diagram Picture.Acquisition module and the step performed by the first image-forming module 102 is performed a plurality of times, until obtain single array element transmitting of N number of array element- Array received combination picture.
Specifically, the step of being performed a plurality of times includes：From first array element transmitting ultrasonic wave, multiple array elements receive successively first Echo-signal, according to first echo signal multiple corresponding images are produced；Second array element launches ultrasonic wave, and multiple array elements are successively Second echo signal is received, multiple corresponding images are produced according to second echo signal；By that analogy, launch to n-th array element Ultrasonic wave, multiple array elements receive successively N echo-signals, and according to N echo-signals multiple corresponding images are produced.
By first echo signal produce multiple corresponding images are combined, obtain first single array element transmitting-array Receive combination picture, by that analogy, by N echo-signals produce multiple corresponding images are combined, obtain n-th Dan Zhen First transmitting-array received combination picture.
The N number of single array element transmitting-array received combination picture for obtaining is carried out using parallel composition processing unit to it multiple Close, obtain total focus image, then total focus image is taken into image envelope, obtain complete matrix data acquisition imaging result.It is distributed Parallel complete matrix acoustic array imaging system works successively according to different transmitting array elements.For transmitting every time, imaging system is adopted Order is performed, and such as array element 1 and array element 2 are launched, and are twice emitting imaging processes, it is possible to use a set of imaging software and hardware system Complete, without the system using redundancy.
Additionally, single-shot, single collection of letters number generate parsing complex signal through Hilbert change.Complex signal is admitted to a data Buffer in array of indexes.Then table module is generated by time delay and launch-middle the transmitting of mono- array element reception, reception battle array according to single array element First position, is calculated the index value of each pixel, then reads corresponding signal, stores to two-dimensional image data structure In.
Preferably, if data directory is realized using software, completed by a read-write array；If data rope Draw table to realize using hardware, then index is completed in real time by a double port memory.Time delay generates table, can be calculated by processor Go out, in being then store in real-time processing module.Additionally, the second image-forming module is additionally operable to using parallel composition processing unit to described Each composite received array element image for receiving array element is combined, and obtains total combination picture.
Image composite module, for all combination pictures for obtaining to be overlapped, obtains total focus image；To gather entirely again Burnt image takes image envelope, obtains complete matrix data acquisition imaging result.
Additionally, perform launching each array element successively, all of single array element transmitting-array received image can be obtained.The mould Block is mainly made up of each single array element transmitting-array received frame buffer and composite module, is sent out by being combined each single array element - array received image is penetrated, complete matrix data acquisition image is obtained.By calculating image envelope as final imaging results.It is right Ultrasonic scanning line signal can obtain the envelope of scanning-line signal by quadrature demodulation, and demodulating process is using analytic signal What mode was obtained, analytic signal is particularly useful in signal analysis, is frequently utilized for obtaining the instantaneous amplitude and frequency of signal.
Fig. 7 is a kind of distributed acoustics array imaging system schematic diagram provided in an embodiment of the present invention.As shown in fig. 7, from First array element launches ultrasonic wave, and multiple array elements receive successively first echo signal, and it is multiple right to be produced according to first echo signal The image answered；Second array element launches ultrasonic wave, and multiple array elements receive successively second echo signal, is produced according to second echo signal The multiple corresponding images of life；By that analogy, ultrasonic wave is launched to n-th array element, multiple array elements receive successively N echo-signals, Multiple corresponding images are produced according to N echo-signals.
By first echo signal produce multiple corresponding images are combined, obtain first transmitting array element array connect Receive combination picture, by that analogy, by N echo-signals produce multiple corresponding images are combined, obtain n-th transmitting battle array The array received combination picture of unit.
The array received combination picture that the n-th for obtaining is launched into array element is multiple to carrying out using parallel composition processing unit Close, obtain total focus image, then total focus image is taken into image envelope, obtain complete matrix data acquisition imaging result.
Fig. 8 realizes block diagram for a kind of first image-forming module provided in an embodiment of the present invention.As shown in figure 8, sending out for single array element Penetrate-mono- array element reception image-forming module.Single-shot-mono- collection of letters number generates parsing complex signal through Hilbert change.Complex signal is sent Enter in a data array of indexes and buffer.Then by time delay generate table module launch according to single array element-mono- array element reception in Penetrate, receive element position, be calculated the index value of each pixel, then read corresponding signal, store to X-Y scheme As in data structure.
Preferably, if data directory is realized using software, completed by a read-write array；If data Concordance list is realized using hardware, is then completed in real time to index by a double port memory.Time delay generates table, can be calculated by processor Draw, in being then store in real-time processing module.
Fig. 9 is the tandem type list array element transmitting-array received in a kind of second image-forming module provided in an embodiment of the present invention Image combination module mechanism map.As shown in figure 9, single array element transmitting-mono- array element is received image and is carried out by tree-like cascade system It is compound, so as to produce single array element transmitting-array received image.By taking 8 array elements as an example, 8 images produced by previous step, by two Two combinations, produce 41 rank images；Then it is combined, produces 22 rank images；3 rank images are finally produced, as single array element transmitting- Array received image.
Figure 10 is a kind of acoustic array imaging method flow chart provided in an embodiment of the present invention.As shown in Figure 10, Yi Zhongsheng Learn array image-forming method, the method comprising the steps of S201-S204:
Step S201：Launch each array element in acoustic array ultrasonic wave successively, and gather every in acoustic array The echo-signal of one array element, analytic trnasformation process is carried out respectively to echo-signal, obtains corresponding analytic signal.
Specifically, the step of carrying out analytic trnasformation process, the corresponding analytic signal of acquisition respectively to echo-signal includes： The corresponding analytic signal of echo-signal that each array element is obtained using Hilbert transform,
Hilbert transform formula is：
A (t)=x (t)+iH { x (t) } (6)
Wherein H { } represents Hilbert transform, and x (t) represents echo-signal, and Hilbert transform integrated form is：
Reilly narrow-long slot parts can also be adopted to obtain the corresponding analytic signal of echo-signal of each array element.
Step S202：Count during according to the physical geometry relation in acoustic array between each array element or by medium refraction The acoustic propagation relation of each array element is calculated, each is obtained according to acoustic propagation relation and is received array element passage in each pixel Delay value, calculates according to the analytic signal and delay value that obtain and obtains each corresponding image for receiving array element.Wherein, according to The analytic signal for obtaining and delay value, obtain each and receive the corresponding image of array element, Distributed Calculation using Distributed Calculation Pyrene process is carried out using polycaryon processor, pyrene process includes：One or more single channels are realized using a processor core Image generate, then using all receiving channels of calculating of multinuclear or parallel multiprocessor；According to the analytic signal that obtains and Delay value, obtains each and receives the corresponding image of array element using Distributed Calculation, and Distributed Calculation is carried out using circuit module Process.
Step S203：Each corresponding image of reception array element is combined using parallel composition processing unit, is obtained The array received combination picture of each transmitting array element.Wherein, also include in this step：By tree-like tandem type Combined Processing Unit is combined to the composite received array element image of each reception array element.Additionally, every to obtaining using parallel composition circuit The image of one array element is combined, and is included the step of obtain combination picture：Using tree-like tandem type compound circuit to each The image of individual array element is combined, and the special circuit module of employing is realized, including special IC or FPGA etc., lead to The same circuits module of different parameters is crossed, Distributed Parallel Computing is realized.
Step S204：All array received combination pictures for obtaining are overlapped, total focus image is obtained；To gather entirely again Burnt image takes image envelope, obtains complete matrix data acquisition imaging result.
Figure 11 (a), (b), (c) be a kind of 1st array element provided in an embodiment of the present invention launch image that the-the 1 array element is received, 1st array element launches image, the image of the 1st array element transmitting the-the 32 array element reception that the-the 16 array element is received.As Figure 11 (a), (b), C shown in (), the transmitting calculated using simulation receives sound field as the input of algorithm.Imaging object is 12 scattering points, is divided respectively Cloth is 20mm, 30mm, 40mm and 50mm apart from transducer depth；Horizontal coordinate is respectively -5mm, 0mm and 5mm.What simulation was calculated Acoustic transducer parameter is：The phase array transducer in 32 array element 0.6mm gaps.
We calculate first first array element transmitting, and single array element transmitting-mono- array element of the 1st, 16 and 32 array element receptions connects Receive image.In order to directly perceived, image is shown.It can be seen that because single-shot list is received, information is not enough to be distinguished into picture Object.
The processing procedure that all single-shot lists due to more than are received is identical, and difference is only that array element coordinate position is different, can be with Realized using identical software or hardware module, so as to reach the purpose of distributed implementation.
Figure 12 (a), (b), the image that (c) is a kind of 1st array element transmitting-array received provided in an embodiment of the present invention, The image of 16 array element transmitting-array receiveds, the image of the 32nd array element transmitting-array received.As shown in Figure 12 (a), (b), (c), The transmitting calculated using simulation receives sound field as the input of algorithm.Imaging object is 12 scattering points, and respectively distribution distance is changed Energy device depth is 20mm, 30mm, 40mm and 50mm；Horizontal coordinate is respectively -5mm, 0mm and 5mm.The acoustic transducer that simulation is calculated Device parameter is：The phase array transducer in 32 array element 0.6mm gaps.It can be seen that single-shot-array received image is Jing can be imaged to target.In the imaging results of difference transmitting array element, only the bending direction of scattering point is different, array element Equivalent to direction of observation, different array element transmittings cause different direction of observations for transmitting, so result is similar, but not exclusively phase Together.
Figure 13 is a kind of complete matrix data acquisition image provided in an embodiment of the present invention.As shown in figure 13, in above single-shot Array is received on the basis of image, all single-shots-array received imaging is combined, you can to obtain complete matrix data acquisition Image.Figure 12 gives complete matrix data acquisition image, there it can be seen that complete matrix collection image carries out imaging ratio to target Image originally becomes apparent from.
The invention provides a kind of acoustic array imaging system and method, the present invention has the beneficial effect that:Gathering each battle array Processed in real time during metadata, so the memory space for taking is few, improve calculating in arithmetic speed, practical application using parallel Time-consuming little, achievable high frame per second imaging, it is adaptable to real time data processing and high-resolution imaging.Additionally, the present invention is applied to not Acoustic array imaging under same frequency and application scenarios.
Above-described specific embodiment, has been carried out further to the purpose of the present invention, technical scheme and beneficial effect Describe in detail, should be understood that the specific embodiment that the foregoing is only the present invention, be not intended to limit the present invention Protection domain, all any modification, equivalent substitution and improvements within the spirit and principles in the present invention, done etc. all should include Within protection scope of the present invention.

Claims (10)

1. a kind of acoustic array imaging system, it is characterised in that include：Data acquisition module 101, the first image-forming module 102, Two image-forming modules 103 and image composite module 104；Wherein,
Data acquisition module 101, for launching each array element in the acoustic array ultrasonic wave successively, and gathers described The echo-signal of each array element in acoustic array, analytic trnasformation process is carried out respectively to the echo-signal, obtains corresponding Analytic signal；
First image-forming module 102, for according to the physical geometry relation in the acoustic array between each array element or passing through The acoustic propagation relation of each array element is calculated during medium refraction, each is obtained according to acoustic propagation relation and is received array element passage every The delay value of one pixel, calculates according to the analytic signal and the delay value that obtain and obtains each reception array element pair The image answered；
Second image-forming module 103, is combined to described each corresponding image of reception array element, obtains each transmitting array element Array received combination picture；
Image composite module 104, using array received composite diagram of the parallel composition processing unit to each transmitting array element As being combined, total focus image is obtained, then the total focus image is taken into image envelope, obtain complete matrix data acquisition imaging As a result.
2. system according to claim 1, it is characterised in that the data acquisition module 101 specifically for：Using Xi Er The corresponding analytic signal of echo-signal that Bert is converted or obtains each array element using Reilly narrow-long slot parts, the Martin Hilb Special transformation for mula is：
A (t)=x (t)+iH { x (t) }
Wherein H { } represents Hilbert transform, and x (t) represents echo-signal, and Hilbert transform integrated form is：
H { x ( t ) } = 1 π ∫ - ∞ + ∞ x ( τ ) t - τ d τ .
3. system according to claim 1, it is characterised in that first image-forming module 102 specifically adopts following form： According to the analytic signal and the delay value that obtain, each is obtained using Distributed Calculation and receives the corresponding figure of array element Picture, the Distributed Calculation carries out pyrene process using polycaryon processor, and the pyrene is processed to be included：Using a processor core Realize that one or more single pass images are generated, then using all receiving channels of calculating of multinuclear or parallel multiprocessor.
4. system according to claim 3, it is characterised in that first image-forming module 102 specifically for：According to obtaining The analytic signal and the delay value, each is obtained using the Distributed Calculation and receives the corresponding image of array element, institute State Distributed Calculation to be processed using circuit module.
5. system according to claim 1, it is characterised in that second image-forming module 103 specifically adopts tree-like level Connection formula parallel composition circuit is combined to the corresponding image of described each array element, obtains array received combination picture.
6. a kind of acoustic array imaging method, it is characterised in that comprise the following steps：
Launch each array element in the acoustic array ultrasonic wave successively, and gather each battle array in the acoustic array The echo-signal of unit, analytic trnasformation process is carried out respectively to the echo-signal, obtains corresponding analytic signal；
Each is calculated during according to the physical geometry relation in the acoustic array between each array element or by medium refraction The acoustic propagation relation of array element, obtains each and receives delay value of the array element passage in each pixel according to acoustic propagation relation, The analytic signal and the delay value according to obtaining is calculated and obtains each corresponding image of reception array element；
Described each corresponding image of reception array element is combined, array received composite diagram of each transmitting array element is obtained Picture；
The array received combination picture of each transmitting array element is combined using parallel composition processing unit, obtains complete Focusedimage, then the total focus image is taken into image envelope, obtain complete matrix data acquisition imaging result.
7. method according to claim 6, it is characterised in that described that the echo-signal is carried out respectively at analytic trnasformation The step of reason, acquisition corresponding analytic signal, includes：Obtain every using Hilbert transform or using Reilly narrow-long slot parts The corresponding analytic signal of echo-signal of one array element, the Hilbert transform formula is：
A (t)=x (t)+iH { x (t) }
Wherein H { } represents Hilbert transform, and x (t) represents echo-signal, and Hilbert transform integrated form is：
H { x ( t ) } = 1 π ∫ - ∞ + ∞ x ( τ ) t - τ d τ .
8. method according to claim 6, it is characterised in that the analytic signal and the time delay that the basis is obtained Value calculates to obtain the step of each receives array element corresponding image to be included：
The Distributed Calculation carries out pyrene process using polycaryon processor, and the pyrene is processed to be included：Using a processor Examine existing one or more single pass images to generate, then led to using all reception of calculating of multinuclear or parallel multiprocessor Road.
9. system according to claim 8, it is characterised in that according to the analytic signal for obtaining and the delay value meter Calculating to obtain the step of each receives array element corresponding image includes：
According to the analytic signal and the delay value that obtain, each is obtained using the Distributed Calculation and receives array element pair The image answered, the Distributed Calculation is processed using circuit module.
10. method according to claim 6, it is characterised in that using parallel composition processing unit to it is described each send out The array received combination picture for penetrating array element is combined, obtain total focus image the step of include：
Complex method is carried out multiple using tree-like tandem type Combined Processing unit to described each corresponding image of reception array element Close, obtain array received combination picture.
CN201611105251.1A 2016-12-05 2016-12-05 A kind of acoustic array imaging system and method Active CN106680825B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611105251.1A CN106680825B (en) 2016-12-05 2016-12-05 A kind of acoustic array imaging system and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611105251.1A CN106680825B (en) 2016-12-05 2016-12-05 A kind of acoustic array imaging system and method

Publications (2)

Publication Number Publication Date
CN106680825A true CN106680825A (en) 2017-05-17
CN106680825B CN106680825B (en) 2019-06-04

Family

ID=58866449

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611105251.1A Active CN106680825B (en) 2016-12-05 2016-12-05 A kind of acoustic array imaging system and method

Country Status (1)

Country Link
CN (1) CN106680825B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019127621A1 (en) * 2017-12-28 2019-07-04 深圳先进技术研究院 Ultrasonic imaging method, system and device
CN110095779A (en) * 2018-01-29 2019-08-06 西南交通大学 A kind of method and apparatus of matrix imaging
CN111307934A (en) * 2019-12-10 2020-06-19 中国科学院声学研究所 Method for calculating response of surface acoustic wave sensor with non-uniform load distribution
CN111308395A (en) * 2020-03-24 2020-06-19 深圳供电局有限公司 Acoustic imaging device, method, device and system for detecting abnormal sound of transformer
CN111726534A (en) * 2020-07-01 2020-09-29 安徽慧可科技有限公司 Signal visualization processing system and processing method thereof
CN111766306A (en) * 2020-05-21 2020-10-13 广州多浦乐电子科技股份有限公司 Full-focus imaging method and device
CN113109443A (en) * 2021-03-08 2021-07-13 中国科学院声学研究所 Focusing acoustic array imaging method and system
CN113180730A (en) * 2021-03-31 2021-07-30 上海深至信息科技有限公司 Distributed beam imaging method
CN114460588A (en) * 2022-04-12 2022-05-10 杭州兆华电子股份有限公司 High-precision imaging method based on active acoustic imager
CN109709562B (en) * 2018-11-08 2023-08-29 中国船舶重工集团公司七五0试验场 Target resolution data processing method of detection sonar

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004031802A1 (en) * 2002-10-04 2004-04-15 Koninklijke Philips Electronics N.V. Method and apparatus for 1d array ultrasound probe
WO2009057062A1 (en) * 2007-10-29 2009-05-07 Koninklijke Philips Electronics, N.V. Systems and methods for ultrasound assembly including multiple imaging transducer arrays
CN104739448A (en) * 2015-04-03 2015-07-01 深圳先进技术研究院 Ultrasonic imaging method and device
WO2015150954A1 (en) * 2014-03-31 2015-10-08 Koninklijke Philips N.V. System and method for acoustic imaging with coherent compounding using intercostal spaces

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004031802A1 (en) * 2002-10-04 2004-04-15 Koninklijke Philips Electronics N.V. Method and apparatus for 1d array ultrasound probe
WO2009057062A1 (en) * 2007-10-29 2009-05-07 Koninklijke Philips Electronics, N.V. Systems and methods for ultrasound assembly including multiple imaging transducer arrays
WO2015150954A1 (en) * 2014-03-31 2015-10-08 Koninklijke Philips N.V. System and method for acoustic imaging with coherent compounding using intercostal spaces
CN104739448A (en) * 2015-04-03 2015-07-01 深圳先进技术研究院 Ultrasonic imaging method and device

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘冬青等: "基于超声阵列的合成孔径聚焦成像方法研究", 《声学技术》 *
刘哲等: "阵列成像系统的设计", 《微计算机应用》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019127621A1 (en) * 2017-12-28 2019-07-04 深圳先进技术研究院 Ultrasonic imaging method, system and device
CN110095779A (en) * 2018-01-29 2019-08-06 西南交通大学 A kind of method and apparatus of matrix imaging
CN109709562B (en) * 2018-11-08 2023-08-29 中国船舶重工集团公司七五0试验场 Target resolution data processing method of detection sonar
CN111307934B (en) * 2019-12-10 2020-11-06 中国科学院声学研究所 Method for calculating response of surface acoustic wave sensor with non-uniform load distribution
CN111307934A (en) * 2019-12-10 2020-06-19 中国科学院声学研究所 Method for calculating response of surface acoustic wave sensor with non-uniform load distribution
CN111308395A (en) * 2020-03-24 2020-06-19 深圳供电局有限公司 Acoustic imaging device, method, device and system for detecting abnormal sound of transformer
CN111766306A (en) * 2020-05-21 2020-10-13 广州多浦乐电子科技股份有限公司 Full-focus imaging method and device
CN111766306B (en) * 2020-05-21 2023-06-27 广州多浦乐电子科技股份有限公司 Full-focus imaging method and device
CN111726534A (en) * 2020-07-01 2020-09-29 安徽慧可科技有限公司 Signal visualization processing system and processing method thereof
CN113109443A (en) * 2021-03-08 2021-07-13 中国科学院声学研究所 Focusing acoustic array imaging method and system
CN113180730A (en) * 2021-03-31 2021-07-30 上海深至信息科技有限公司 Distributed beam imaging method
CN113180730B (en) * 2021-03-31 2022-10-18 上海深至信息科技有限公司 Distributed beam imaging method
CN114460588A (en) * 2022-04-12 2022-05-10 杭州兆华电子股份有限公司 High-precision imaging method based on active acoustic imager

Also Published As

Publication number Publication date
CN106680825B (en) 2019-06-04

Similar Documents

Publication Publication Date Title
CN106680825B (en) A kind of acoustic array imaging system and method
CN101203183B (en) Ultrasound imaging system with pixel oriented processing
Bulavinov et al. Sampling phased array a new technique for signal processing and ultrasonic imaging
CN104688271B (en) Ultrasonic imaging method and ultrasonic imaging device by synthetic focusing
CN103181779B (en) Ultrasonic energy conversion device and ultrasonic imaging system and method
CN100456016C (en) Multi-channel electronic parallel scanning photoacoustic real-time tomo graphic-imaging method and apparatus thereof
KR100742467B1 (en) Method and apparatus for visualization of motion in ultrasound flow imaging using packet data acquisition
CN102727256B (en) Dual focusing beam forming method and device based on virtual array elements
CN102697524B (en) Full-focus ultrasonic imaging method and application of method in blood flow imaging
US4011750A (en) Method and apparatus for ultrasonic examination of objects
CN113109443B (en) Focused acoustic array imaging method and system
CN104546003B (en) VFR (variable frame rate) color ultra-high speed ultrasonic imaging method based on plane wave transmission
CN106526602A (en) Ultrasonic ghost imaging method and device based on principle of calorescence ghost imaging
CN107303186A (en) Frequency multiplexed in elastogram
EP2705383A2 (en) Enhanced ultrasound image formation using qualified regions of overlapping transmit beams
Johnson et al. Digital computer simulation study of a real-time collection, post-processing synthetic focusing ultrasound cardiac camera
JP2000149015A (en) Method for edge enhancement of image and imaging device
US8235906B2 (en) System and method for accelerated focused ultrasound imaging
CN117147694A (en) Inverse problem-based ultrasonic full-focusing imaging sparse regularization reconstruction method and equipment
US6048313A (en) Method and apparatus for fractal-based enhancement of ultrasound imaging
CN105011965A (en) Focusing implementation method by adoption of aperture synthesis and ultrasonic holography for scanning
CN108573474A (en) A kind of Photoacoustic image optimization method using inverse convolution algorithm
CN105548363A (en) Multi-route identification based ultrasonic detection imaging method
CN108318883B (en) Distance direction calculation method based on natural weighting of full-spectrum transducer
Piedade et al. Compressive sensing strategy on sparse array to accelerate ultrasonic TFM imaging

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