CN115474965A - Ultrasonic CT array probe array element directivity self-adaptive evaluation method and device - Google Patents
Ultrasonic CT array probe array element directivity self-adaptive evaluation method and device Download PDFInfo
- Publication number
- CN115474965A CN115474965A CN202211332333.5A CN202211332333A CN115474965A CN 115474965 A CN115474965 A CN 115474965A CN 202211332333 A CN202211332333 A CN 202211332333A CN 115474965 A CN115474965 A CN 115474965A
- Authority
- CN
- China
- Prior art keywords
- array
- directivity
- virtual
- array element
- receiving
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 54
- 239000000523 sample Substances 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 claims abstract description 29
- 238000004088 simulation Methods 0.000 claims abstract description 10
- 238000003491 array Methods 0.000 claims description 10
- 238000005457 optimization Methods 0.000 claims description 8
- 238000002604 ultrasonography Methods 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000002939 conjugate gradient method Methods 0.000 claims description 3
- 238000013215 result calculation Methods 0.000 claims description 3
- 238000002945 steepest descent method Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 3
- 238000011478 gradient descent method Methods 0.000 abstract description 2
- 239000011159 matrix material Substances 0.000 abstract description 2
- 238000004613 tight binding model Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 11
- 238000003384 imaging method Methods 0.000 description 9
- 238000013461 design Methods 0.000 description 8
- 230000003044 adaptive effect Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000013170 computed tomography imaging Methods 0.000 description 4
- 230000005284 excitation Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 238000007654 immersion Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4416—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to combined acquisition of different diagnostic modalities, e.g. combination of ultrasound and X-ray acquisitions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4411—Constructional features of apparatus for radiation diagnosis the apparatus being modular
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5247—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from an ionising-radiation diagnostic technique and a non-ionising radiation diagnostic technique, e.g. X-ray and ultrasound
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4411—Device being modular
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4444—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5261—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from different diagnostic modalities, e.g. ultrasound and X-ray
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/58—Testing, adjusting or calibrating the diagnostic device
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The application provides an ultrasonic CT array probe array element directivity self-adaptive evaluation method and device, which are used for acquiring transmission data of one transmission and multiple receptions in water by utilizing a multi-array element probe in a full-matrix acquisition mode. The transmitting array elements are replaced by the virtual array with variable weight, the rest receiving adopts receiving point receiving, and array receiving data at the moment is calculated by a finite difference method or a Green function analysis method. And comparing the difference between the calculated received data and the actual observed data, and updating the weight of the virtual array by using a gradient descent method. After the updating is finished, the virtual array weight is fixed, and sound pressures surrounding the transmitting array element at the same distance and different angles are calculated in a simulation mode, so that the directivity of the array element can be obtained.
Description
Technical Field
The invention relates to the technical field of medical ultrasound, in particular to an ultrasonic CT array probe array element directivity self-adaptive evaluation method and device.
Background
A probe in the field of medical ultrasonic imaging transmits ultrasonic waves to a target body and then receives reflected echoes of human tissues to perform acoustic impedance structural parameter imaging, or receives transmitted waves penetrating through the human tissues to perform acoustic velocity and acoustic attenuation parameter imaging. The human tissue image reconstructed by the ultrasonic can be used as an imaging basis for auxiliary diagnosis. At present, a medical ultrasonic CT probe is generally an array (ring, arc, cylindrical surface and the like) consisting of a plurality of array elements so as to be beneficial to imaging human tissues by utilizing ultrasonic waves from multiple angles.
The size of a single array element of the medical ultrasonic probe, especially a high-frequency probe, is often difficult to ignore compared with the central wavelength of an excitation ultrasonic signal, and the excitation sound pressure loaded on the array element has non-uniformity, so that the ultrasonic signal emitted by the array element has directivity. The ideal state of the ultrasonic signal emitted by the array element is that the array element sends ultrasonic waves in a mode of a point sound source, and the array element in reality has directivity instead of sending the ultrasonic signal in the ideal state.
For the ultrasonic CT imaging technology, the knowledge of probe array element directivity is the key for improving the structural imaging resolution and increasing the parameter imaging accuracy. According to the traditional array element directivity evaluation method, sound pressure data on a circle of position with a probe to be measured as a circle center needs to be measured one by a hydrophone, the workload is huge for a multi-probe array, and the change of the array element directivity after the equipment is used for multiple times is difficult to evaluate.
Disclosure of Invention
The purpose of this application is to solve the defect that prior art exists.
The application provides a rapid, accurate and automatic ultrasonic array element directivity self-adaptive evaluation method based on array observation data, and the method can be used in the fields of ultrasonic CT imaging and the like.
In a first aspect, the present application provides an ultrasonic CT array probe array element directivity adaptive evaluation method, where the array probe includes N array elements, where N is an integer greater than 1, and the evaluation method includes: determining array elements to be tested in the N array elements; transmitting ultrasonic waves by using the array elements to be detected, and receiving the ultrasonic waves by using the other M array elements to form a first observation data set, wherein M is an integer which is larger than zero and smaller than or equal to N; simulating the array element to be measured by using a virtual array consisting of a plurality of variable-weight virtual points to be measured, simulating the ultrasonic waves emitted by the array element to be measured, simulating the M array elements by using M receiving point arrays, and receiving the ultrasonic waves emitted by simulation to form a second observation data set; obtaining a target function of variable weights of the multiple virtual points to be measured based on the first observation data set and the second observation data set; obtaining directional weights corresponding to the variable weights of the multiple virtual points to be detected based on the target function; and obtaining a directivity evaluation result of the virtual array based on the directivity weight, namely the directivity evaluation result of the array element to be tested.
In a feasible embodiment, the directivity evaluation results of all array elements are measured by using a virtual array method and compared with a non-directional virtual array, so that the directivity evaluation of the transmitting probe is realized successfully through the observation data of the array by using the virtual array sound source method. The accurate sound source modeling simulation method has important significance for ultrasonic CT image reconstruction with millimeter-scale and even submillimeter-scale resolution.
In a possible design, after obtaining a directivity estimation result of the virtual array based on the directivity weight, that is, the directivity estimation result of the array element to be measured, the method further includes: and sequentially taking the N array elements of the array probe as the array elements to be detected, and further obtaining the directivity evaluation result of the N array elements, namely the directivity evaluation result of the array probe.
In another possible design, the simulating the M array elements using M receive point arrays and receiving ultrasonic waves to form a second observation data set includes: simulating the M array elements by using M receiving point arrays; and calculating the receiving data of the M receiving points based on the propagation function to form a second observation data set.
In another possible design, the propagation function includes: finite difference or green's function.
In another possible design, the objective function includes: a two-norm comparison function, or variance function, or difference-by-difference function, or least-squares function.
In another possible design, the obtaining, based on the objective function, directional weights corresponding to the variable weights of the multiple virtual points to be measured includes: and optimizing the objective function based on the variable-weight objective function of the multiple virtual points to be measured, so as to obtain the corresponding variable weights of the multiple virtual points to be measured, namely the directional weights.
In another possible design, the optimizing the objective function includes: and optimizing the objective function by a gradient descent optimization algorithm.
In another possible design, the gradient descent optimization algorithm includes: steepest descent method, or conjugate gradient method, or L-BFGS method.
In another possible design, the obtaining a directivity estimation result of the virtual array based on the directivity weight includes: setting the variable weight of the plurality of virtual points to be measured as the directional weight; setting new receiving points at the same distance and different angle positions by taking the virtual array as a circle center; when the virtual array is simulated to transmit ultrasonic waves, the maximum value of the received signals of the new receiving points is taken; and searching to obtain the maximum values of a plurality of receiving signals with different distances, and obtaining the directivity evaluation result of the virtual array element to be tested.
In a second aspect, the present application provides an ultrasonic CT array probe array element directivity adaptive evaluation device, the array probe includes N array elements, where N is an integer greater than 1, and includes: the array element determining module is used for determining array elements to be detected in the N array elements; a first observation data set forming module, configured to transmit ultrasonic waves by using the array element to be detected, and receive the ultrasonic waves by using the remaining M array elements to form a first observation data set, where M is an integer greater than zero and less than or equal to N; the virtual array module is used for simulating the array element to be tested by using a virtual array consisting of a plurality of variable-weight virtual points to be tested and simulating the ultrasonic waves emitted by the array element to be tested; the second observation data set forming module is used for simulating the M array elements by using M receiving point arrays and receiving the ultrasonic waves transmitted by simulation to form a second observation data set; a directional weight calculation module, configured to obtain a target function of variable weights of the multiple virtual points to be measured based on the first observation data set and the second observation data set, and obtain directional weights corresponding to the variable weights of the multiple virtual points to be measured based on the target function; and the directivity evaluation result calculation module is used for obtaining the directivity evaluation result of the virtual array based on the directivity weight, namely the directivity evaluation result of the array element to be tested.
The application provides an ultrasonic CT array probe array element directivity self-adaptive evaluation method and device, which is a method for inverting ultrasonic array element directivity by utilizing array self-observation data through an optimization method based on gradient descent and can solve the problems of ultrasonic CT probe directivity evaluation and correction. The observation data of the ultrasonic array can be used for reducing extra directivity measurement work, and an automatic and high-precision array element directivity evaluation scheme is formed.
Drawings
Fig. 1 is a schematic flow chart of an ultrasonic CT array probe array element directivity adaptive evaluation method according to an embodiment of the present invention;
fig. 2 is a schematic composition diagram of an ultrasonic CT array probe array element directivity adaptive evaluation device according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a 512-element ultrasonic ring array system according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of waveform signals of a first observation data set of array elements # 384 in accordance with an embodiment of the present invention;
FIG. 5 is a schematic view of the sound wave propagation process of the loading point sound source according to the embodiment of the present invention;
fig. 6 is a schematic view illustrating a sound wave propagation process of a virtual array sound source according to an embodiment of the present invention;
fig. 7 is a schematic diagram of directivity estimation data of a load point sound source according to an embodiment of the present invention;
fig. 8 is a schematic diagram of directivity estimation data of a virtual array sound source according to an embodiment of the present invention;
FIG. 9 is a schematic directional diagram of a loading point sound source according to an embodiment of the present invention;
fig. 10 is a schematic directional diagram of a virtual array sound source according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the accompanying drawings and embodiments, and it is apparent that the described embodiments are some, but not all embodiments of the present invention. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without any creative effort, shall fall within the protection scope of the present invention.
In the following detailed description of the present invention, certain specific details are set forth in order to provide a better understanding of the present invention. It will be apparent to one skilled in the art that the present invention may be practiced without these specific details.
In the practical application of medical ultrasound CT, if there is no effective analysis on the directivity of the sound source, the resolution of the acquired structural image is reduced, and the error of the parametric image is large. However, an automatic and self-adaptive array element directivity evaluation method is lacked in the traditional practice. The invention provides an array element directivity analysis method based on array self-observation data, which aims at an ultrasonic CT multi-array element array probe system, does not need to use a hydrophone to carry out repeated additional measurement, and automatically determines the array element directivity with high precision.
For the ultrasonic CT imaging technology, the knowledge of probe array element directivity is the key for improving the structural imaging resolution and increasing the parameter imaging accuracy. Therefore, the method for rapidly, accurately and automatically evaluating the directivity of the ultrasonic array element based on the array observation data is provided, and can be used in the fields of ultrasonic CT imaging and the like.
As shown in fig. 1, the present application provides an adaptive evaluation method for directivity of an array element of an ultrasonic CT array probe, which mainly comprises the following steps:
step S110, the array probe comprises N array elements, wherein N is an integer greater than 1, and array elements to be detected in the N array elements are determined.
And step S120, transmitting ultrasonic waves by using the array elements to be detected, and receiving the ultrasonic waves by using the rest M array elements to form a first observation data set, wherein M is an integer which is larger than zero and smaller than or equal to N.
Step S130, simulating the array element to be measured by using a virtual array composed of a plurality of variable-weight virtual points to be measured, simulating the ultrasonic waves emitted by the array element to be measured, simulating the M array elements by using M receiving point arrays, and receiving the ultrasonic waves emitted by simulation to form a second observation data set.
Step S140, obtaining a variable weight objective function of the multiple virtual points to be measured based on the first observation data set and the second observation data set.
And step S150, obtaining directional weights corresponding to the variable weights of the virtual points to be measured based on the objective function.
And step S160, obtaining a directivity evaluation result of the virtual array based on the directivity weight, namely the directivity evaluation result of the array element to be tested.
In steps S110 and S120, in a possible embodiment, under the water immersion condition, each array element of the multi-array-element ultrasonic transducer is utilized to transmit ultrasonic waves one by one, and each time of transmission, all the array elements are received, or received at intervals of several array elements, where when all the array elements are received, the directivity is higher in accuracy, the precision is higher, but the simulation process is slightly complex and the calculation time is longer. If the j array element is transmitted, the wave trains received by all the array elements are arranged together in sequence to form a first observation data set of the array pair j array element
In step S130, in one possible embodiment, the j array element is virtualizedThe array is replaced by a receiving array element which adopts receiving point analog receiving, wherein a virtual array can be circular, square or linear, and the effect is best when the j number array element is taken as the center. Virtual array initial weight w 0 Can be set to the same value, and calculates the receiving data set of the array when the array element with the number j transmits through finite difference or Green function
In a possible embodiment, the size of a single array element of the medical ultrasonic probe, especially the high-frequency probe, is often difficult to ignore compared with the central wavelength of the excitation ultrasonic signal, and the excitation sound pressure loaded on the array element has non-uniformity, so that the emitted ultrasonic signal has directivity. Therefore, the virtual array is used for replacing the j array element, so that the real state of the j array element when the ultrasonic signal is sent can be better simulated, and the directivity of the ultrasonic signal can be more truly simulated. The ideal state of a single array element is that the array element transmits ultrasonic waves in a point sound source mode, and the array element in reality has directivity instead of transmitting ultrasonic signals in the ideal state.
In step S140, in a possible embodiment, an objective function of variable weight of the virtual array element to be measured is obtained based on the first observation data set and the second observation data set. Designing an objective function
Where E is a function for measuring the difference between the received data of the two methods, and a two-norm may be used as a more conventional case, but any other optional comparison function, such as a variance function, or a difference-by-difference function, or a least-squares function, is not excluded. E.g. two norms
Where w is the virtual array weight, i is the number of the receiving array element, and t is the time series. Such as in a possibleIn the embodiment of (1), if 42 virtual points are used to simulate the transmitted j array elements, w is a matrix or vector consisting of 42 unequal elements.
In step S150, based on the objective function, directional weights corresponding to the variable weights of the multiple virtual points to be measured are obtained. In one possible embodiment, the optimization problem is solved by solving it
The optimal virtual array weight can be obtained to calculate the array element directivity. The optimization problem can be solved by a gradient descent optimization algorithm (such as a steepest descent method, a conjugate gradient method and an L-BFGS method).
In step S160, in a possible embodiment, the virtual array weight w obtained in the previous step is substituted into the virtual array, the received signals at the same distance and different angle positions with the transmitting probe as the center of the circle are calculated, the maximum value is taken, the maximum values of the received signals at multiple different distances are found, and a polar coordinate diagram related to the angle is drawn, so as to obtain a directivity diagram of the j array element.
Repeating the process from step S110 to step S160, the directivity of all array elements of the whole array can be obtained, and the directivity evaluation of the whole ultrasound CT array probe can be performed.
As shown in fig. 2, an embodiment of the present invention provides an ultrasonic CT array probe array element directivity adaptive evaluation apparatus, where the array probe includes N array elements, where N is an integer greater than 1, and mainly includes the following modules:
a to-be-detected array element determining module 210, configured to determine to-be-detected array elements in the N array elements;
a first observation data set forming module 220, configured to transmit ultrasonic waves by using the array element to be detected, and receive the ultrasonic waves by using the remaining M array elements to form a first observation data set, where M is an integer greater than zero and less than or equal to N;
a virtual array module 230, configured to simulate the array element to be tested by using a virtual array composed of multiple virtual points to be tested with variable weights, and simulate the ultrasonic waves emitted by the array element to be tested;
a second observation data set forming module 240, configured to simulate the M array elements by using M receiving point arrays, and receive the ultrasound waves emitted by the simulation to form a second observation data set;
a directional weight calculation module 250, configured to obtain an objective function of the variable weights of the multiple virtual points to be measured based on the first observation data set and the second observation data set, and obtain directional weights corresponding to the variable weights of the multiple virtual points to be measured based on the objective function;
and a directivity evaluation result calculation module 260, configured to obtain a directivity evaluation result of the virtual array based on the directivity weight, where the directivity evaluation result is a directivity evaluation result of the array element to be detected.
Example 1
As shown in fig. 3, the embodiment of the present invention is directed to a set of 512-element ultrasound ring array system, where the ring is exemplary and may be in other shapes, such as an ellipse, a polygon, etc. The invention provides a specific embodiment, which explains the specific application method of the invention, and the directivity of the ultrasonic annular array system can be well simulated by the virtual array obtained by using the point sound source to replace the transmitting array element and using the virtual array to replace the transmitting array element for comparison. Note that this embodiment is only for illustrating the application method of the present invention, and is not intended to limit a specific application scenario. The specific steps of this example are as follows:
the method comprises the following steps: under the condition of water immersion, each array element of the annular multi-array element ultrasonic transducer is utilized to transmit ultrasonic waves one by one, and all the array elements (including the transmitting array elements) receive each time of transmission; taking 384 array element as an example, the wave trains received by all the array elements are arranged together in sequence, and the wave signals are
As shown in fig. 4.
Step two: the transmitting array element is replaced by a virtual array of 42 virtual points (which can be understood as 42 array elements), the weight of each virtual point of the initial virtual array is set to be 1/42, and 51 is used2 receiving points simulate 512 array elements to receive signals, and a received data set when 384 array elements are transmitted is calculated
And similarly, a point sound source is used for replacing a transmitting array element to transmit signals to obtain a receiving data set.
Step three: designing an objective function
Solving by gradient descent method
Step four: the sound field propagation processes of the nondirectional point sound source shown in fig. 5 and the optimized virtual array sound source shown in fig. 6 are calculated, and it can be seen that after the method of the embodiment of the present invention is optimized, the loading of the directional sound source is successfully realized in the simulation.
Step five: a directivity estimation array of 80 array elements was set in an angle of 108 ° inside the circular array of 10 cm from the center of the virtual array, and estimation data of a point sound source as shown in fig. 7 and estimation data of a virtual array sound source as shown in fig. 8 were calculated.
Step six: the connecting line of the circle center of the probe and the center of the circular array is 0 DEG, the angle is increased in the counterclockwise direction, the maximum value of the evaluation data is taken, the forward directivity diagram of the point sound source shown in figure 9 is drawn, and the forward directivity diagram of the virtual array sound source shown in figure 10 is drawn.
Compared with a method for transmitting signals by using a point sound source instead of a transmitting array element, the virtual array sound source method successfully realizes directivity evaluation of the transmitting probe through the observation data of the array. The accurate sound source modeling has important significance for ultrasonic CT image reconstruction seeking millimeter-level or even submillimeter-level resolution.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment may be implemented by software plus a necessary general hardware platform, and may also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, so that it should be understood that the above-mentioned embodiments are only one of the specific embodiments of the present invention, and should not be used to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. An ultrasonic CT array probe array element directivity self-adaptive evaluation method is disclosed, wherein the array probe comprises N array elements, wherein N is an integer greater than 1, and the evaluation method comprises the following steps:
determining array elements to be tested in the N array elements;
transmitting ultrasonic waves by using the array elements to be detected, and receiving the ultrasonic waves by using the other M array elements to form a first observation data set, wherein M is an integer which is greater than zero and less than or equal to N;
simulating the array element to be measured by using a virtual array consisting of a plurality of variable-weight virtual points to be measured, simulating the ultrasonic waves emitted by the array element to be measured, simulating the M array elements by using M receiving point arrays, and receiving the ultrasonic waves emitted by simulation to form a second observation data set;
obtaining a target function of variable weights of the multiple virtual points to be measured based on the first observation data set and the second observation data set;
obtaining directional weights corresponding to the variable weights of the multiple virtual points to be detected based on the target function;
and obtaining a directivity evaluation result of the virtual array based on the directivity weight, namely the directivity evaluation result of the array element to be tested.
2. The evaluation method according to claim 1, wherein after obtaining the directivity evaluation result of the virtual array based on the directivity weight, that is, the directivity evaluation result of the array element to be tested, the method further comprises:
and sequentially taking the N array elements of the array probe as array elements to be detected, and further obtaining the directivity evaluation result of the N array elements, namely the directivity evaluation result of the array probe.
3. The method of claim 1, wherein said simulating the M array elements using the M receive spot arrays and receiving ultrasound waves to form a second observation data set comprises:
simulating the M array elements by using M receiving dot arrays;
and calculating the receiving data of the M receiving points based on the propagation function to form a second observation data set.
4. The method of claim 3, wherein the propagation function comprises:
finite difference or green's function.
5. The evaluation method of claim 1, wherein the objective function comprises:
a two-norm comparison function, or variance function, or a difference-by-difference function, or a least-squares function.
6. The evaluation method according to claim 1, wherein the obtaining of the directional weights corresponding to the variable weights of the plurality of virtual points to be measured based on the objective function comprises:
and optimizing the objective function based on the variable-weight objective function of the multiple virtual points to be measured, so as to obtain the corresponding variable weights of the multiple virtual points to be measured, namely the directional weights.
7. The evaluation method according to claim 6, wherein the optimizing the objective function comprises:
and optimizing the objective function through a gradient descent optimization algorithm.
8. The evaluation method of claim 7, wherein the gradient descent optimization algorithm comprises:
steepest descent method, or conjugate gradient method, or L-BFGS method.
9. The method according to claim 1, wherein the obtaining a directivity estimation result of the virtual array based on the directivity weight comprises:
setting the variable weight of the plurality of virtual points to be measured as the directional weight;
setting new receiving points at the same distance and different angle positions by taking the virtual array as a circle center;
when the virtual array is simulated to transmit ultrasonic waves, the maximum value of the received signals of the new receiving point is taken;
and searching to obtain the maximum values of a plurality of receiving signals with different distances, and obtaining the directivity evaluation result of the virtual array element to be tested.
10. The utility model provides an supersound CT array probe array element directive property self-adaptation evaluation device, the array probe includes N array elements, and wherein, N is for being greater than 1 integer, its characterized in that includes:
the array element determining module is used for determining array elements to be detected in the N array elements;
the first observation data set forming module is used for transmitting ultrasonic waves by using the array elements to be detected, and receiving the ultrasonic waves by using the rest M array elements to form a first observation data set, wherein M is an integer which is larger than zero and is less than or equal to N;
the virtual array module is used for simulating the array element to be tested by using a virtual array consisting of a plurality of variable-weight virtual points to be tested and simulating the ultrasonic waves emitted by the array element to be tested;
the second observation data set forming module is used for simulating the M array elements by using M receiving point arrays and receiving the ultrasonic waves transmitted by simulation to form a second observation data set;
a directional weight calculation module, configured to obtain a target function of variable weights of the multiple virtual points to be measured based on the first observation data set and the second observation data set, and obtain directional weights corresponding to the variable weights of the multiple virtual points to be measured based on the target function;
and the directivity evaluation result calculation module is used for obtaining the directivity evaluation result of the virtual array based on the directivity weight, namely the directivity evaluation result of the array element to be tested.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211332333.5A CN115474965B (en) | 2022-10-28 | 2022-10-28 | Ultrasonic CT array probe array element directivity self-adaptive evaluation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211332333.5A CN115474965B (en) | 2022-10-28 | 2022-10-28 | Ultrasonic CT array probe array element directivity self-adaptive evaluation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115474965A true CN115474965A (en) | 2022-12-16 |
| CN115474965B CN115474965B (en) | 2023-03-28 |
Family
ID=84396242
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211332333.5A Active CN115474965B (en) | 2022-10-28 | 2022-10-28 | Ultrasonic CT array probe array element directivity self-adaptive evaluation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115474965B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002048853A (en) * | 2000-08-02 | 2002-02-15 | Matsushita Electric Ind Co Ltd | Radio wave arrival direction estimator and directivity variable transmitter-receiver |
| CN101238992A (en) * | 2008-02-01 | 2008-08-13 | 哈尔滨工业大学 | Ultrasonic imaging system self-adaption beam former based on correlation analysis |
| CN107789008A (en) * | 2017-11-29 | 2018-03-13 | 声泰特（成都）科技有限公司 | A kind of self-adapting ultrasonic beam synthetic method and system based on channel data |
| CN110554358A (en) * | 2019-09-25 | 2019-12-10 | 哈尔滨工程大学 | noise source positioning and identifying method based on virtual ball array expansion technology |
| CN114886469A (en) * | 2022-05-11 | 2022-08-12 | 中国科学院声学研究所 | Array element positioning method and device of ultrasonic CT array probe |
-
2022
- 2022-10-28 CN CN202211332333.5A patent/CN115474965B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002048853A (en) * | 2000-08-02 | 2002-02-15 | Matsushita Electric Ind Co Ltd | Radio wave arrival direction estimator and directivity variable transmitter-receiver |
| CN101238992A (en) * | 2008-02-01 | 2008-08-13 | 哈尔滨工业大学 | Ultrasonic imaging system self-adaption beam former based on correlation analysis |
| CN107789008A (en) * | 2017-11-29 | 2018-03-13 | 声泰特（成都）科技有限公司 | A kind of self-adapting ultrasonic beam synthetic method and system based on channel data |
| CN110554358A (en) * | 2019-09-25 | 2019-12-10 | 哈尔滨工程大学 | noise source positioning and identifying method based on virtual ball array expansion technology |
| CN114886469A (en) * | 2022-05-11 | 2022-08-12 | 中国科学院声学研究所 | Array element positioning method and device of ultrasonic CT array probe |
Non-Patent Citations (2)
| Title |
|---|
| 王波,万明习,王素品,陈忠民: "医用超声阵列换能器波束容差分析与变迹处理" * |
| 郑跃滨;武湛君;雷振坤;高东岳;周凯;张佳奇;杨正岩;邹建超;: "基于超声导波的航空航天结构损伤诊断成像技术研究进展" * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115474965B (en) | 2023-03-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20220167949A1 (en) | Calibration of multiple aperture ultrasound probes | |
| CN109900256B (en) | Self-adaptive ocean mobile acoustic tomography system and method | |
| Karaman et al. | A phase aberration correction method for ultrasound imaging | |
| CN102590804B (en) | Overland testing system of Doppler sonar and testing method thereof | |
| CN112754527B (en) | Data processing method for low-frequency ultrasonic thoracic imaging | |
| CN106500820A (en) | A kind of sonic velocity measurement method and device for estimating two-dimensional direction-of-arrival | |
| CN113252779A (en) | Method for improving defect inversion imaging quality | |
| Liu et al. | Estimation and correction of ultrasonic wavefront distortion using pulse-echo data received in a two-dimensional aperture | |
| Noda et al. | Ultrasound imaging with a flexible probe based on element array geometry estimation using deep neural network | |
| Arcondoulis et al. | Experimental validation and performance analysis of deep learning acoustic source imaging methods | |
| CN115474965B (en) | Ultrasonic CT array probe array element directivity self-adaptive evaluation method and device | |
| Omidvar et al. | An intrinsic shape estimation algorithm for flexible ultrasound probes intended for clinical applications | |
| CN111308421B (en) | Method for acquiring acoustic radiation of target free field in shallow sea | |
| CN215340310U (en) | Underwater vertical array element position measuring system | |
| CN110274613B (en) | Laboratory detection device and method suitable for acoustic wave and tide measuring instrument | |
| Li | Position and time-delay calibration of transducer elements in a sparse array for underwater ultrasound imaging | |
| Filipik et al. | Time-of-flight based calibration of an ultrasonic computed tomography system | |
| RU2798390C1 (en) | Method for passive determination of coordinates of a noise-generating object | |
| CN114777696B (en) | Ultrasonic thickness measuring method | |
| Burov et al. | Increasing Resolution Capability of Two-Dimensional Tomograph over Third Coordinate: Separating Reconstruction of c (r)-, ρ (r)-m α (r, ω)-Scatterer Characteristics | |
| Maisano et al. | An empirical regressive model to improve the metrological performance of mobile spatial coordinate measuring systems | |
| Li et al. | The cross algorithm for phase-aberration correction in medical ultrasound images formed with two-dimensional arrays | |
| Lou et al. | Ultrasound Planar Array Imaging Metric Analysis | |
| Shen et al. | Delay Calibration for Ultrasound Computed Tomography System Using Neural Network | |
| CN117991242A (en) | Spherical shell array beam forming method for correcting scattering sound field |
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 |