WO2019127621A1 - Ultrasonic imaging method, system and device - Google Patents

Ultrasonic imaging method, system and device Download PDF

Info

Publication number
WO2019127621A1
WO2019127621A1 PCT/CN2018/070069 CN2018070069W WO2019127621A1 WO 2019127621 A1 WO2019127621 A1 WO 2019127621A1 CN 2018070069 W CN2018070069 W CN 2018070069W WO 2019127621 A1 WO2019127621 A1 WO 2019127621A1
Authority
WO
WIPO (PCT)
Prior art keywords
coordinate system
coordinates
array elements
receiving
cartesian coordinate
Prior art date
Application number
PCT/CN2018/070069
Other languages
French (fr)
Chinese (zh)
Inventor
邱维宝
刘荣
郑海荣
Original Assignee
深圳先进技术研究院
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 深圳先进技术研究院 filed Critical 深圳先进技术研究院
Publication of WO2019127621A1 publication Critical patent/WO2019127621A1/en

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging

Definitions

  • the present invention relates to the field of ultrasound technology, and more particularly to an ultrasound imaging method, system and apparatus.
  • the working principle of the ultrasonic imaging device is that the probe receives the ultrasonic wave reflected from the target, converts it into a digital signal through A/D, and then uploads it from the front end to the back end GPU (Graphics Processing Unit) for beam. Synthesizing, the plane wave delay superposition, weighting and coherent combining under multiple angles are completed in the GPU to form a frame of the original RF RF data, which is sent to the display by the digital signal processor in the computer system.
  • the ultrasonic imaging device utilizes the ultrasonic plane wave imaging method to scan the detection area by using the plane wave scanning. The scanning of the full-width region of interest can be completed by one shot control, and the time taken for one image is much lower than that of the conventional focus-by-line scan.
  • a multi-angle plane wave coherent composite imaging method is adopted to radiate a plurality of plane waves of different angles to the region of interest, and the image quality can be improved after coherent superposition.
  • a line coordinate system is needed. By calculating the delay time of the pixel point by the coordinate value of the pixel point and the coordinate value of the array element in the linear coordinate system, the coherent superimposed image can be calculated.
  • the position of different array elements in the convex array surface involves the change of the angle between the array element and the coordinate axis, so this method using only the line coordinate system is not It is suitable for the scanning transformation under the convex angle of the convex array, so the beamforming process under the multi-angle of the convex array plane wave can not be completed by the conventional coherent composite algorithm.
  • an object of the present invention is to provide an ultrasound imaging method, system and apparatus for converting a coordinate of a receiving element in a Cartesian coordinate system into a calculation mode when a plurality of array elements are arranged in a convex array form.
  • the coordinates of the receiving array elements in the polar coordinate system that is, the coordinate values of the receiving array elements are converted from linear coordinates to angled coordinates, and the beam synthesis under the convex array multi-angle or single-angle plane wave emission can be completed.
  • an embodiment of the present invention provides an ultrasound imaging method, which is applied to an ultrasound imaging apparatus, the ultrasound imaging apparatus including a probe, the probe includes a plurality of array elements, and the plurality of array elements are used to transmit an ultrasonic pair
  • the object is scanned and receives echo data
  • the ultrasonic imaging method comprises: performing, for each pixel, the following steps: acquiring echo data received by the plurality of array elements; acquiring coordinates of the pixel points in the Cartesian coordinate system, and receiving The coordinates of the receiving array element of the pixel echo data, the coordinates of the transmitting array element transmitting the ultrasonic wave to the pixel point; converting the coordinates of the receiving array element in the Cartesian coordinate system to receiving the array element in the polar coordinate system
  • the coordinates of the pixel points are calculated according to the coordinates of the pixel points and the transmitting array elements in the Cartesian coordinate system and the coordinates of the receiving array elements in the polar coordinate system; according to the delay time and the echo data,
  • an embodiment of the present invention provides a first possible implementation of the first aspect, wherein the position of receiving an array element in a Cartesian coordinate system is converted to a position of receiving an array element in polar coordinates
  • the method includes: calculating an angle of a radius of curvature of the receiving array element and a first direction axis of the Cartesian coordinate system, wherein the Cartesian coordinate system is a center of the curved surface formed by the plurality of array elements, and the curved surface is the coordinate zero point
  • the axis where the center line is located is the first direction axis, and the axis perpendicular to the center line of the curved surface is the second direction axis; according to the angle, the position of the receiving element in the Cartesian coordinate system is converted to the polar coordinate Receive the position of the array element.
  • the embodiment of the present invention provides a second possible implementation manner of the first aspect, wherein the angle is calculated as:
  • ⁇ n is the angle of curvature of the receiving array element and the first direction axis of the Cartesian coordinate system
  • num is the number of multiple array elements
  • pitch is the separation distance of multiple array elements
  • R is the radius of curvature
  • n For the nth receiving array element, n has a value range of [1, num].
  • an embodiment of the present invention provides a third possible implementation manner of the first aspect, wherein the receiving frame is in a polar coordinate system according to coordinates of a pixel point and a transmitting array element in a Cartesian coordinate system.
  • Calculating the delay time of the pixel including: setting the point near the radius of curvature of the surface formed by the pixel and the plurality of array elements as a transmitting array element, wherein calculating the delay time of the pixel point
  • the formula is:
  • (x m , z m ) is the coordinate of the pixel point P m in the Cartesian coordinate system
  • (E nx , E nz ) is the coordinate of the receiving element E n in the polar coordinate system
  • (C x , C z ) is the coordinates of the center C of the surface in the Cartesian coordinate system.
  • the embodiment of the present invention provides a fourth possible implementation manner of the first aspect, wherein the coordinates of the center of the surface of the Cartesian coordinate system are calculated as:
  • (C x , C z ) is the coordinate of the center of the circle in the Cartesian coordinate system
  • R is the radius of curvature
  • pitch is the separation distance of multiple array elements
  • num is the number of multiple array elements
  • is Multiple array elements deflect the sweep angle.
  • the embodiment of the present invention provides a fifth possible implementation manner of the first aspect, wherein, according to the delay time and echo data, the pixel point intensity value under beamforming is obtained,
  • the method includes: searching, according to the acquisition frequency and the delay time, the echo data from the pixel points received by the plurality of array elements; performing weighted summation on the echo data to obtain the pixel intensity value under beam combining .
  • an embodiment of the present invention further provides an ultrasound imaging system, including: an acquiring module, configured to acquire echo data received by a plurality of array elements, acquire coordinates of pixel points in a Cartesian coordinate system, and receive the The coordinates of the receiving array element of the pixel echo data, the coordinates of the transmitting array element transmitting the ultrasonic wave to the pixel point; the conversion module, the conversion module being connected to the acquiring module for receiving in the Cartesian coordinate system Converting the coordinates of the array elements to the coordinates of the receiving array elements in the polar coordinate system; the computing module, the computing module being coupled to the conversion module for using the coordinates of the pixel points and the transmitting array elements in the Cartesian coordinate system, Receiving the coordinates of the array element in polar coordinates, calculating the delay time of the pixel; the imaging module, the imaging module is respectively connected to the obtaining module and the calculating module, and is configured to obtain a beam according to the delay time and the echo data.
  • the pixel intensity value under synthesis configured to acquire
  • the embodiment of the present invention provides a first possible implementation manner of the second aspect, wherein the converting module is specifically configured to: calculate a radius of curvature of the receiving array element and a first direction of the Cartesian coordinate system An angle of the axis, wherein the Cartesian coordinate system is a coordinate zero point formed by a plurality of array elements, and an axis of the curved surface is a first direction axis perpendicular to a center line of the curved surface The axis is a second direction axis; according to the angle, the position of the receiving element in the Cartesian coordinate system is converted to the position at which the element is received in polar coordinates.
  • the embodiment of the present invention provides a second possible implementation manner of the second aspect, wherein the angle is calculated as:
  • ⁇ n is the angle of curvature of the receiving array element and the first direction axis of the Cartesian coordinate system
  • num is the number of multiple array elements
  • pitch is the separation distance of multiple array elements
  • R is the radius of curvature
  • n For the nth receiving array element, n has a value range of (1, num).
  • an embodiment of the present invention further provides an ultrasound imaging apparatus, including: a probe, the probe includes a plurality of array elements, wherein the plurality of array elements are configured to emit ultrasonic waves to scan an object, and receive echo data; An ultrasound imaging system according to any of claims 7-9.
  • the embodiment of the invention brings about the following beneficial effects: by acquiring the coordinates of the pixel point in the Cartesian coordinate system, the coordinates of the receiving array element of the pixel point, the coordinates of the transmitting array element of the pixel point, and then the Cartesian coordinate system
  • the coordinates of the receiving element are converted to the coordinates of the receiving element in the polar coordinate system, and the total distance of the pixel is calculated according to the coordinate value to calculate the delay time of the pixel, which is obtained by delay time and echo data.
  • the pixel points under beamforming such that when a plurality of array elements are arranged in a convex array form, the coordinates of the receiving array elements in the Cartesian coordinate system are converted into the receiving array elements in the polar coordinate system in a computational manner.
  • the coordinates, that is, the coordinate values of the receiving array elements are converted from linear coordinates to angled coordinates, and the synthesis of beams under multi-angle or single-angle plane wave emission under convex array can be completed, and ultrafast ultrasound
  • FIG. 1 is a structural diagram of an ultrasound imaging method according to an embodiment of the present invention.
  • FIG. 2 is a schematic diagram of calculation of a pixel point of a convex array probe without deflection according to an embodiment of the present invention
  • FIG. 3 is a schematic diagram showing a relationship between coordinates of a Cartesian coordinate system and coordinates of a polar coordinate system according to an embodiment of the present invention
  • FIG. 4 is a schematic diagram of calculation of a pixel point under deflection of a convex array probe according to an embodiment of the present invention
  • FIG. 5 is a structural diagram of an ultrasound imaging system according to an embodiment of the present invention.
  • FIG. 6 is a structural diagram of an ultrasound imaging apparatus according to an embodiment of the present invention.
  • 200-ultrasound imaging system 210-acquisition module; 220-conversion module; 230-calculation module; 240-imaging module; 300-ultrasonic imaging device; 310-probe; 311-array; 320-display.
  • linear array multi-angle plane wave correlation composite imaging algorithm is disclosed for ultrasound imaging.
  • the linear array refers to the linear array structure of multiple array elements in the ultrasonic transducer, and the linear array multi-angle plane wave coherent composite imaging.
  • the method is not applicable to the method of convex array multi-angle or convex array single-angle plane wave coherent composite imaging.
  • an ultrasound imaging method, system and device provided by the embodiments of the present invention obtain pixel points in a Cartesian coordinate system.
  • the coordinates, the coordinates of the receiving element of the pixel, the coordinates of the transmitting element of the pixel, and the coordinates of the receiving element in the Cartesian coordinate system are converted to the coordinates of the receiving element in the polar coordinate system, according to
  • the coordinate value calculates the total distance of propagation of the pixel to calculate the delay time of the pixel, and obtains the pixel under beamforming through the delay time and the echo data, so that when the plurality of array elements are arranged in a convex array form,
  • the coordinates of the receiving array element in the Cartesian coordinate system are converted into the coordinates of the receiving array element in the polar coordinate system, that is, the receiving array element Standard values are converted by the linear coordinate angled coordinates, it is possible to complete the synthesis of the beam angles of the convex array or a single plane wave emission angle, to achieve ultrafast ultrasound imaging.
  • an ultrasound imaging method disclosed in the embodiment of the present invention is first introduced in detail.
  • the method is applied to an ultrasound imaging apparatus 300, and the ultrasound imaging apparatus 300 includes a probe 310.
  • the probe 310 includes a plurality of array elements 311 for transmitting ultrasonic waves to scan an object and receiving echo data.
  • the specific workflow of the ultrasound imaging apparatus 300 is: after the ultrasound imaging apparatus 300 is activated, the plurality of array elements 311 emit a plane wave to scan the measured object, and the emitted plane wave is reflected by the measured object after encountering the object.
  • the plurality of array elements 311 receive the reflected ultrasonic waves, that is, the plurality of array elements 311 receive the echo data, wherein the echo data records information about the amplitude and phase of the reflected ultrasonic wavefront reflected by the measured object.
  • the ultrasound imaging method includes:
  • S110 Acquire echo data received by multiple array elements.
  • a plurality of array elements 311 receive echo data, and a plurality of array elements 311 are uploaded through channels connected thereto to enable acquisition of the echo data for post processing.
  • S120 Acquire coordinates of a pixel point in a Cartesian coordinate system, coordinates of a receiving array element that receives pixel point echo data, and coordinates of a transmitting array element that transmits an ultrasonic wave to the pixel point.
  • a Cartesian coordinate system in which the coordinates (x, z) of the pixel point P1 are obtained, and in this Cartesian coordinate system, the coordinates of the receiving array element E n of the received pixel point P1 echo data are acquired.
  • the coordinates of the receiving array element A1 receiving the pixel point P1 echo data are acquired.
  • convex array probes referred to herein are merely exemplary, and arcuate probes such as intracavity probes and the like are all within the definition of the present invention.
  • step S120 before performing step S120, it should also include: filtering the echo data to filter out other ultrasonic waves that are not part of the echo data, so that the echo data can be purified, thereby improving the accuracy of post processing. .
  • converting the position of the receiving element in the Cartesian coordinate system to the position of receiving the element in polar coordinates comprises: calculating a radius of curvature of the receiving element and a first direction axis of the Cartesian coordinate system
  • the angle of the Cartesian coordinate system is a coordinate zero point formed by a plurality of array elements, the axis of the curved line is the first direction axis, and the axis perpendicular to the center line of the curved surface is the second direction axis;
  • the position of the receiving element in the Cartesian coordinate system is converted to the position at which the element is received in polar coordinates.
  • the Cartesian coordinate system has a curved circle C1 (0, 0) formed by a plurality of array elements 311 as a coordinate zero point, and an x-axis of a Cartesian coordinate system is perpendicular to a center line of the curved surface, and a flute
  • the z-axis of the Carl coordinate system is the axis of the centerline of the surface.
  • the angle ⁇ n between the receiving element E n and the z axis is calculated, wherein the angle is calculated as:
  • ⁇ n is the angle of curvature of the receiving array element and the first direction axis of the Cartesian coordinate system
  • num is the number of multiple array elements
  • pitch is the separation distance of multiple array elements
  • R is the radius of curvature
  • n For the nth receiving array element, n has a value range of [1, num].
  • E nx R*sin ⁇ n
  • E nz R*cos ⁇ n
  • S140 Calculate the delay time of the pixel point according to the coordinates of the pixel point and the transmitting array element in the Cartesian coordinate system and the coordinates of the receiving element in the polar coordinate system.
  • step S140 includes: setting a pixel point and a point formed by a plurality of array elements along a radius of curvature along a radius of curvature as a transmitting array element, wherein calculating a delay time of the pixel point is:
  • the delay time of the pixel is the sum of the propagation time of the transmitting array element transmitting the ultrasonic wave to the pixel point and the propagation time of the receiving ultrasonic wave receiving the reflected ultrasonic wave reflected from the pixel point.
  • the pixel P1 when the pixel P1 is calculated in the receiving array element E n to the pixel position of point P1 and the receiving array element E n on the basis of the curved surface along a curvature radius P1 and the pixel array element is formed from a plurality of The next closest point is the transmit array element.
  • the pixel P1 in a convex array probe under non-deflection, the pixel P1 receives E n in the array element, the array element receiving point P1 of the pixel A1.
  • the convex array probe deflection ⁇ when the convex array probe deflection ⁇ is scanned, the pixel point P1 in the array element E n is received, and the receiving array element of the pixel point P1 is A2.
  • the delay time of the pixel point can be calculated by confirming the coordinates of the center of the surface of the surface.
  • the coordinates of the center of the surface of the Cartesian coordinate system are calculated as:
  • (C x , C z ) is the coordinate of the center of the circle in the Cartesian coordinate system
  • R is the radius of curvature
  • pitch is the separation distance of multiple array elements
  • num is the number of multiple array elements
  • is Multiple array elements deflect the sweep angle.
  • the coordinates of the center C of the curved surface formed by the plurality of array elements 311 are:
  • the coordinates of the center C can be calculated by the deflection angle ⁇ , the number of the plurality of array elements, the separation distance of the plurality of array elements, and the radius of curvature.
  • the coordinates of the center C at any deflection angle can be calculated.
  • Step S150 includes: searching for echo data from the pixel points received by the plurality of array elements 311 according to the acquisition frequency and the delay time; and performing weighted summation on the echo data to obtain pixel points under beam synthesis.
  • the array element E 1 when the probe 310 formed by the plurality of array elements 311 collects echo data, periodic acquisition is performed using a certain frequency. Taking the array element E 1 , the acquisition frequency is 1 s, and the total acquisition time is 10 s as an example, after the acquisition is completed, the array element E 1 includes echo data of a plurality of collected pixel points 10 times. Since the time of receiving the echo data of each pixel is different, the echo time of the pixel E 1 receiving the pixel is the total time of the pixel propagation, that is, the delay time of the pixel, and the delay of the pixel is calculated. At the time, the echo data of the pixel in the element E 1 can be found.
  • a plurality of echo data of the pixel point are found in the plurality of array elements 311, and then the plurality of echo data of the pixel point are weighted and summed, that is, the converted optical signal of the pixel point is converted.
  • the signals are weighted and summed.
  • the weighting function is a window function
  • the window function may be a Hanning window, a Hamming window, a rectangular window function, or the like.
  • the imaging depth can be set to limit the range of pixels to obtain the highest quality image.
  • each pixel is processed to obtain an image under beamforming.
  • the acquisition module 210 the conversion module 220, the calculation module 230, and the imaging module 240 are included.
  • the acquiring module 210 is configured to acquire echo data received by multiple array elements, and acquire coordinates of pixel points in a Cartesian coordinate system, coordinates of receiving array elements that receive pixel point echo data, and transmit ultrasonic waves to pixel points.
  • the coordinates of the transmitting element The conversion module 220 is coupled to the acquisition module 210 for converting the coordinates of the receiving array elements in the Cartesian coordinate system to the coordinates of the receiving array elements in the polar coordinate system.
  • the calculation module 230 is connected to the conversion module 220 for calculating the delay time of the pixel point according to the coordinates of the pixel point and the transmission element in the Cartesian coordinate system and the coordinates of the receiving element in the polar coordinate.
  • the imaging module 240 is connected to the acquisition module 210 and the calculation module 230, respectively, for obtaining pixel intensity values under beamforming according to the delay time and the echo data.
  • the ultrasound imaging system 200 further includes a processing module, and the processing module is connected to the imaging module 240.
  • the processing module is configured to process the digital signal of the result of the imaging module 240 after detection, filtering, logarithmic compression, dynamic range conversion, etc., to the digital signal. Display.
  • the conversion module 220 is specifically configured to: calculate an angle of a radius of curvature of the receiving array element and a first direction axis of the Cartesian coordinate system, wherein the Cartesian coordinate system is a surface center formed by a plurality of array elements For the coordinate zero point, the axis of the curved surface is the first direction axis, and the axis perpendicular to the center line of the surface is the second direction axis; according to the angle, the position of the receiving element in the Cartesian coordinate system is converted to the pole The position of the receiving element in coordinates.
  • the angle is calculated as:
  • ⁇ n is the angle of curvature of the receiving array element and the first direction axis of the Cartesian coordinate system
  • num is the number of multiple array elements
  • pitch is the separation distance of multiple array elements
  • R is the radius of curvature
  • n For the nth receiving array element, n has a value range of (1, num).
  • an ultrasound imaging apparatus 300 includes a probe 310, an ultrasound imaging system 200, and a display screen 320.
  • the probe 310 includes a plurality of array elements 311 for transmitting ultrasonic waves to scan an object and receiving echo data.
  • the ultrasonic imaging system 200 according to any of the above embodiments.
  • Display screen 320 is coupled to ultrasound imaging system 200 for displaying imaging results of ultrasound imaging system 200. It is worth noting that a plurality of array elements 311 are connected to the ultrasound imaging system 200 through channels.
  • each block of the flowchart or block diagram can represent a module, a program segment, or a portion of code that includes one or more of the Executable instructions.
  • the functions noted in the blocks may also occur in a different order than that illustrated in the drawings. For example, two consecutive blocks may be executed substantially in parallel, and they may sometimes be executed in the reverse order, depending upon the functionality involved.
  • each block of the block diagrams and/or flowcharts, and combinations of blocks in the block diagrams and/or flowcharts can be implemented in a dedicated hardware-based system that performs the specified function or function. Or it can be implemented by a combination of dedicated hardware and computer instructions.
  • the synthesis of the beam may be performed in a processor, which may be an image processor, or an integrated circuit chip having signal processing capabilities.
  • a processor which may be an image processor, or an integrated circuit chip having signal processing capabilities.
  • each step of the above method may be completed by an integrated logic circuit of hardware in a processor or an instruction in a form of software.
  • the above processor may be a general-purpose processor, including a central processing unit (CPU), a network processor (Network Processor, NP for short), or a digital signal processor (DSP). , Application Specific Integrated Circuit (ASIC), Field-Programmable Gate Array (FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware component.
  • CPU central processing unit
  • NP Network Processor
  • DSP digital signal processor
  • ASIC Application Specific Integrated Circuit
  • FPGA Field-Programmable Gate Array
  • the general purpose processor may be a microprocessor or the processor or any conventional processor or the like.
  • the steps of the method disclosed in the embodiments of the present invention may be directly implemented by the hardware decoding processor, or may be performed by a combination of hardware and software modules in the decoding processor.
  • the software module can be located in a conventional storage medium such as random access memory, flash memory, read only memory, programmable read only memory or electrically erasable programmable memory, registers, and the like.
  • the storage medium is located in the memory, and the processor reads the information in the memory and combines the hardware to complete the steps of the above method.
  • connection shall be understood broadly, and may be, for example, a fixed connection, a detachable connection, or an integral connection, unless otherwise explicitly defined and defined; It is a mechanical connection or an electrical connection; it can be directly connected or indirectly connected through an intermediate medium, which can be the internal connection of two components.
  • intermediate medium which can be the internal connection of two components.

Abstract

Provided are an ultrasonic imaging method, a system and a device, relating to the technical field of ultrasonics. Said method comprises: performing the following steps for each pixel point: acquiring echo data received by a plurality of array elements; acquiring the coordinates of a pixel point in a Cartesian coordinate system, the coordinates of receiving array elements for receiving the echo data of the pixel point, and the coordinates of transmitting array elements for transmitting ultrasonic waves to the pixel point; converting the coordinates of the receiving array elements in the Cartesian coordinate system into the coordinates of the receiving array elements in a polar coordinate system; calculating a delay time of the pixel point according to the coordinates of the pixel point and the transmitting array elements in the Cartesian coordinate system and the coordinates of the receiving array elements in the polar coordinate system; and obtaining, according to the delay time and the echo data, an intensity value of the pixel point under beam forming. By computationally converting the coordinates of the receiving array elements in a Cartesian coordinate system into coordinates in a polar coordinate system, the present technical solution is able to complete beam forming under the transmission of convex-array plane waves at single/multiple angles, so as to realize ultrafast convex-array ultrasonic imaging.

Description

超声成像方法、系统和设备Ultrasound imaging methods, systems and devices
相关申请的交叉引用Cross-reference to related applications
本申请要求于2017年12月28日提交中国专利局的申请号为201711467316.1、名称为“超声成像方法、系统和设备”的中国专利申请的优先权,其全部内容通过引用结合在本申请中。The present application claims priority to Chinese Patent Application No. PCT Application No. No. No. No. No. No. No. No. No. No No No No No No No No No No No No No No No No No No No No No No No No No No No No No No No No No No No No No No
技术领域Technical field
本发明涉及超声技术领域,尤其是涉及一种超声成像方法、系统和设备。The present invention relates to the field of ultrasound technology, and more particularly to an ultrasound imaging method, system and apparatus.
背景技术Background technique
通常情况下,超声成像设备的工作原理为,探头接收从目标反射回的超声波,经过A/D变换成数字信号,再从前端上传到后端GPU(Graphics Processing Unit,图形处理器)中进行波束合成,在GPU中完成多角度下的平面波延时叠加、加权和相干复合,形成一帧目标原始射频RF数据,该数据经计算机系统中数字信号处理器后,送于显示器显示。超声成像设备利用超声平面波成像方法,利用平面波扫描检测区域,一次发射控制即可完成扫描全幅感兴趣区,一幅图像花费的时间远低于传统聚焦逐线扫描。Usually, the working principle of the ultrasonic imaging device is that the probe receives the ultrasonic wave reflected from the target, converts it into a digital signal through A/D, and then uploads it from the front end to the back end GPU (Graphics Processing Unit) for beam. Synthesizing, the plane wave delay superposition, weighting and coherent combining under multiple angles are completed in the GPU to form a frame of the original RF RF data, which is sent to the display by the digital signal processor in the computer system. The ultrasonic imaging device utilizes the ultrasonic plane wave imaging method to scan the detection area by using the plane wave scanning. The scanning of the full-width region of interest can be completed by one shot control, and the time taken for one image is much lower than that of the conventional focus-by-line scan.
基于平面波是非聚焦波,成像质量差,采用多角度平面波相干复合成像方法,将多个不同角度的平面波辐射感兴趣区域,相干叠加后可以提高图像质量。在计算多角度平面波的合成中,需要用到一种线坐标系,通过在线性坐标系下像素点的坐标值和阵元的坐标值计算像素点的延迟时间,从而可以计算出相干叠加的图像,然而,在凸阵平面波的波束合成过程中,凸阵曲面中的不同阵元的位置涉及到阵元与坐标轴之间的角度的变化,所以,仅用线坐标系的这种方式并不适用于凸阵偏转角度下的扫描变换,因此无法利用常规的相干复合算法完成凸阵平面波多角度下的波束合成过程。Based on the fact that the plane wave is a non-focusing wave and the imaging quality is poor, a multi-angle plane wave coherent composite imaging method is adopted to radiate a plurality of plane waves of different angles to the region of interest, and the image quality can be improved after coherent superposition. In the calculation of multi-angle plane wave synthesis, a line coordinate system is needed. By calculating the delay time of the pixel point by the coordinate value of the pixel point and the coordinate value of the array element in the linear coordinate system, the coherent superimposed image can be calculated. However, in the beamforming process of a convex array plane wave, the position of different array elements in the convex array surface involves the change of the angle between the array element and the coordinate axis, so this method using only the line coordinate system is not It is suitable for the scanning transformation under the convex angle of the convex array, so the beamforming process under the multi-angle of the convex array plane wave can not be completed by the conventional coherent composite algorithm.
发明内容Summary of the invention
有鉴于此,本发明的目的在于提供超声成像方法、系统和设备,在多个阵元以凸阵形式排列时,在计算方式上通过将在笛卡尔坐标系下的接收阵元的坐标转换为极坐标系下的接收阵元的坐标,即接收阵元的坐标值由线性坐标转换为带有角度的坐标,就能够完成凸阵多角度或者单角度平面波发射下的波束的合成,实现了超快超声成像。In view of the above, an object of the present invention is to provide an ultrasound imaging method, system and apparatus for converting a coordinate of a receiving element in a Cartesian coordinate system into a calculation mode when a plurality of array elements are arranged in a convex array form. The coordinates of the receiving array elements in the polar coordinate system, that is, the coordinate values of the receiving array elements are converted from linear coordinates to angled coordinates, and the beam synthesis under the convex array multi-angle or single-angle plane wave emission can be completed. Fast ultrasound imaging.
第一方面,本发明实施例提供了一种超声成像方法,应用于超声成像设备,所述超声成像设备包括探头,所述探头包括多个阵元,所述多个阵元用于发射超声波对物体进行扫描,并接收回波数据,所述超声成像方法包括:对每个像素点进行以下步骤:获取多个阵元接收的回波数据;获取在笛卡尔坐标系下像素点的坐标、接收所述像素点回波数据的接收阵元的坐标、发射超声波到所述像素点的发射阵元的坐标;将在笛卡尔坐标系下接收阵元的坐标转换为在极坐标系下接收阵元的坐标;根据在笛卡尔坐标系下像素点和发射阵元的坐标、在极坐标系下接收阵元的坐标,计算像素点的延时时间;根据所述延时时间和回波数据,得到波束合成下的所述像素点强度值。In a first aspect, an embodiment of the present invention provides an ultrasound imaging method, which is applied to an ultrasound imaging apparatus, the ultrasound imaging apparatus including a probe, the probe includes a plurality of array elements, and the plurality of array elements are used to transmit an ultrasonic pair The object is scanned and receives echo data, and the ultrasonic imaging method comprises: performing, for each pixel, the following steps: acquiring echo data received by the plurality of array elements; acquiring coordinates of the pixel points in the Cartesian coordinate system, and receiving The coordinates of the receiving array element of the pixel echo data, the coordinates of the transmitting array element transmitting the ultrasonic wave to the pixel point; converting the coordinates of the receiving array element in the Cartesian coordinate system to receiving the array element in the polar coordinate system The coordinates of the pixel points are calculated according to the coordinates of the pixel points and the transmitting array elements in the Cartesian coordinate system and the coordinates of the receiving array elements in the polar coordinate system; according to the delay time and the echo data, The pixel intensity value under beam synthesis.
结合第一方面,本发明实施例提供了第一方面的第一种可能的实施方式,其中,所述将在笛卡尔坐标系下接收阵元的位置转换为在极坐标下接收阵元的位置,具体包括:计算接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,其中,所述笛卡尔坐标系为以多个阵元形成的曲面圆心为坐标零点,所述曲面的中心线所在的轴为第一方向轴,与所述曲面的中心线垂直的轴为第二方向轴;根据所述角度,将在笛卡尔坐标系下接收阵元的位置转换为在极坐标下接收阵元的位置。In conjunction with the first aspect, an embodiment of the present invention provides a first possible implementation of the first aspect, wherein the position of receiving an array element in a Cartesian coordinate system is converted to a position of receiving an array element in polar coordinates Specifically, the method includes: calculating an angle of a radius of curvature of the receiving array element and a first direction axis of the Cartesian coordinate system, wherein the Cartesian coordinate system is a center of the curved surface formed by the plurality of array elements, and the curved surface is the coordinate zero point The axis where the center line is located is the first direction axis, and the axis perpendicular to the center line of the curved surface is the second direction axis; according to the angle, the position of the receiving element in the Cartesian coordinate system is converted to the polar coordinate Receive the position of the array element.
结合第一方面,本发明实施例提供了第一方面的第二种可能的实施方式,其中,所述角度的算式为:With reference to the first aspect, the embodiment of the present invention provides a second possible implementation manner of the first aspect, wherein the angle is calculated as:
Figure PCTCN2018070069-appb-000001
Figure PCTCN2018070069-appb-000001
其中,β n为接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,num为多个阵元的个数,pitch为多个阵元的间隔距离,R为曲率半径,n为第n个接收阵元,n的取值范围为[1,num]。 Where β n is the angle of curvature of the receiving array element and the first direction axis of the Cartesian coordinate system, num is the number of multiple array elements, pitch is the separation distance of multiple array elements, and R is the radius of curvature, n For the nth receiving array element, n has a value range of [1, num].
结合第一方面,本发明实施例提供了第一方面的第三种可能的实施方式,其中,所述根据在笛卡尔坐标系下像素点和发射阵元的坐标、在极坐标系下接收阵元的坐标,计算像素点的延时时间,包括:设定所述像素点与距离多个阵元形成的曲面沿曲率半径下最近的点为发射阵元,其中,计算像素点的延时时间的算式为:With reference to the first aspect, an embodiment of the present invention provides a third possible implementation manner of the first aspect, wherein the receiving frame is in a polar coordinate system according to coordinates of a pixel point and a transmitting array element in a Cartesian coordinate system. Calculating the delay time of the pixel, including: setting the point near the radius of curvature of the surface formed by the pixel and the plurality of array elements as a transmitting array element, wherein calculating the delay time of the pixel point The formula is:
τ(E n,x m,z m)*c=sqrt((x m-C x) 2+(z m-C z) 2)-R+sqrt((x m-E nx) 2+(z m+E nz) 2) τ(E n ,x m ,z m )*c=sqrt((x m -C x ) 2 +(z m -C z ) 2 )-R+sqrt((x m -E nx ) 2 +(z m +E nz ) 2 )
其中,(x m,z m)为在笛卡尔坐标系下像素点P m的坐标,(E nx,E nz)为在极坐标系下接收阵元E n的坐标,(C x,C z)为在笛卡尔坐标系曲面圆心C的坐标。 Where (x m , z m ) is the coordinate of the pixel point P m in the Cartesian coordinate system, and (E nx , E nz ) is the coordinate of the receiving element E n in the polar coordinate system, (C x , C z ) is the coordinates of the center C of the surface in the Cartesian coordinate system.
结合第一方面,本发明实施例提供了第一方面的第四种可能的实施方式,其中,所述在笛卡尔坐标系曲面圆心的坐标,计算公式为:With reference to the first aspect, the embodiment of the present invention provides a fourth possible implementation manner of the first aspect, wherein the coordinates of the center of the surface of the Cartesian coordinate system are calculated as:
C x=(-1)*(R 2+R 2-2*R*R*COSθ)*cos((num/2+0.5)*pitch/R-θ/2) C x =(-1)*(R 2 +R 2 -2*R*R*COSθ)*cos((num/2+0.5)*pitch/R-θ/2)
C z=(-1)*(R 2+R 2-2*R*R*COSθ)*sin((num/2+0.5)*pitch/R-θ/2) C z =(-1)*(R 2 +R 2 -2*R*R*COSθ)*sin((num/2+0.5)*pitch/R-θ/2)
其中,(C x,C z)为所述圆心在笛卡尔坐标系下的坐标,R为曲率半径,pitch为多个阵元的间隔距离,num为多个阵元的个数,θ为在多个阵元偏转扫射角度。 Where (C x , C z ) is the coordinate of the center of the circle in the Cartesian coordinate system, R is the radius of curvature, pitch is the separation distance of multiple array elements, num is the number of multiple array elements, and θ is Multiple array elements deflect the sweep angle.
结合第一方面,本发明实施例提供了第一方面的第五种可能的实施方式,其中,所述根据所述延时时间和回波数据,得到波束合成下的所述像素点强度值,包括:根据采集频率和延时时间,查找多个阵元接收到的来自所述像素点的回波数据;对所述回波数据进行加权求和,得到波束合成下的所述像素点强度值。With reference to the first aspect, the embodiment of the present invention provides a fifth possible implementation manner of the first aspect, wherein, according to the delay time and echo data, the pixel point intensity value under beamforming is obtained, The method includes: searching, according to the acquisition frequency and the delay time, the echo data from the pixel points received by the plurality of array elements; performing weighted summation on the echo data to obtain the pixel intensity value under beam combining .
第二方面,本发明实施例还提供一种超声成像系统,包括:获取模块,用于获取多个阵元接收的回波数据,并获取在笛卡尔坐标系下像素点的坐标、接收所述像素点回波数据的接收阵元的坐标、发射超声波到所述像素点的发射阵元的坐标;转换模块,所述转换模块与所述获取模块相连,用于将在笛卡尔坐标系下接收阵元的坐标转换为在极坐标系下接收阵元的坐标;计算模块,所述计算模块与所述转换模块相连,用于根据在笛卡尔坐标系下像素点和发射阵元的坐标、在极坐标下接收阵元的坐标,计算像素点的延时时间;成像模块,所述成像模块分别与所述获取模块和计算模块相连,用于根据所述延时时间和回波数据,得到波束合成下的所述像素点强度值。In a second aspect, an embodiment of the present invention further provides an ultrasound imaging system, including: an acquiring module, configured to acquire echo data received by a plurality of array elements, acquire coordinates of pixel points in a Cartesian coordinate system, and receive the The coordinates of the receiving array element of the pixel echo data, the coordinates of the transmitting array element transmitting the ultrasonic wave to the pixel point; the conversion module, the conversion module being connected to the acquiring module for receiving in the Cartesian coordinate system Converting the coordinates of the array elements to the coordinates of the receiving array elements in the polar coordinate system; the computing module, the computing module being coupled to the conversion module for using the coordinates of the pixel points and the transmitting array elements in the Cartesian coordinate system, Receiving the coordinates of the array element in polar coordinates, calculating the delay time of the pixel; the imaging module, the imaging module is respectively connected to the obtaining module and the calculating module, and is configured to obtain a beam according to the delay time and the echo data. The pixel intensity value under synthesis.
结合第二方面,本发明实施例提供了第二方面的第一种可能的实施方式,其中,所述转换模块,具体用于:计算接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,其中,所述笛卡尔坐标系为以多个阵元形成的曲面圆心为坐标零点,所述曲面的中心线所在的轴为第一方向轴,与所述曲面的中心线垂直的轴为第二方向轴;根据所述角度,将在笛卡尔坐标系下接收阵元的位置转换为在极坐标下接收阵元的位置。With reference to the second aspect, the embodiment of the present invention provides a first possible implementation manner of the second aspect, wherein the converting module is specifically configured to: calculate a radius of curvature of the receiving array element and a first direction of the Cartesian coordinate system An angle of the axis, wherein the Cartesian coordinate system is a coordinate zero point formed by a plurality of array elements, and an axis of the curved surface is a first direction axis perpendicular to a center line of the curved surface The axis is a second direction axis; according to the angle, the position of the receiving element in the Cartesian coordinate system is converted to the position at which the element is received in polar coordinates.
结合第二方面,本发明实施例提供了第二方面的第二种可能的实施方式,其中,所述角度的算式为:With reference to the second aspect, the embodiment of the present invention provides a second possible implementation manner of the second aspect, wherein the angle is calculated as:
Figure PCTCN2018070069-appb-000002
Figure PCTCN2018070069-appb-000002
其中,β n为接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,num为多个阵元的个数,pitch为多个阵元的间隔距离,R为曲率半径,n为第n个接收阵元,n的取值范围为(1,num)。 Where β n is the angle of curvature of the receiving array element and the first direction axis of the Cartesian coordinate system, num is the number of multiple array elements, pitch is the separation distance of multiple array elements, and R is the radius of curvature, n For the nth receiving array element, n has a value range of (1, num).
第三方面,本发明实施例还提供一种超声成像设备,包括:探头,所述探头包括多个阵元,所述多个阵元用于发射超声波对物体进行扫描,并接收回波数据;根据权利要求7-9任一项所述的超声成像系统。In a third aspect, an embodiment of the present invention further provides an ultrasound imaging apparatus, including: a probe, the probe includes a plurality of array elements, wherein the plurality of array elements are configured to emit ultrasonic waves to scan an object, and receive echo data; An ultrasound imaging system according to any of claims 7-9.
本发明实施例带来了以下有益效果:通过获取笛卡尔坐标系下的像素点的坐标、该像素点的接收阵元的坐标、该像素点的发射阵元的坐标,再将笛卡尔坐标系下的接收阵元的坐标转换为极坐标系下的接收阵元的坐标,根据坐标值计算像素点的传播总距离从而计算该像素点的延时时间,在通过延时时间和回波数据得到波束合成下的该像素点,这样在多个阵元以凸阵形式排列时,在计算方式上通过将在笛卡尔坐标系下的接收阵元的坐标转换为极坐标系下的接收阵元的坐标,即接收阵元的坐标值由线性坐标转换为带有角度的坐标,就能够完成凸阵下多角度或者单角度平面波发射下的波束的合成,实现了超快超声成像。The embodiment of the invention brings about the following beneficial effects: by acquiring the coordinates of the pixel point in the Cartesian coordinate system, the coordinates of the receiving array element of the pixel point, the coordinates of the transmitting array element of the pixel point, and then the Cartesian coordinate system The coordinates of the receiving element are converted to the coordinates of the receiving element in the polar coordinate system, and the total distance of the pixel is calculated according to the coordinate value to calculate the delay time of the pixel, which is obtained by delay time and echo data. The pixel points under beamforming, such that when a plurality of array elements are arranged in a convex array form, the coordinates of the receiving array elements in the Cartesian coordinate system are converted into the receiving array elements in the polar coordinate system in a computational manner. The coordinates, that is, the coordinate values of the receiving array elements are converted from linear coordinates to angled coordinates, and the synthesis of beams under multi-angle or single-angle plane wave emission under convex array can be completed, and ultrafast ultrasound imaging is realized.
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。Other features and advantages of the invention will be set forth in the description which follows, and The objectives and other advantages of the invention are realized and attained by the invention particularly pointed in
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。The above described objects, features and advantages of the present invention will become more apparent from the aspects of the appended claims.
附图说明DRAWINGS
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the drawings to be used in the specific embodiments or the description of the prior art will be briefly described below, and obviously, the attached in the following description The drawings are some embodiments of the present invention, and those skilled in the art can obtain other drawings based on these drawings without any creative work.
图1为本发明实施例提供的超声成像方法的结构图;1 is a structural diagram of an ultrasound imaging method according to an embodiment of the present invention;
图2为本发明实施例提供的凸阵探头无偏转下像素点的计算示意图;2 is a schematic diagram of calculation of a pixel point of a convex array probe without deflection according to an embodiment of the present invention;
图3为本发明实施例提供的笛卡尔坐标系的坐标与极坐标系的坐标的关系示意图;3 is a schematic diagram showing a relationship between coordinates of a Cartesian coordinate system and coordinates of a polar coordinate system according to an embodiment of the present invention;
图4为本发明实施例提供的凸阵探头偏转下像素点的计算示意图;4 is a schematic diagram of calculation of a pixel point under deflection of a convex array probe according to an embodiment of the present invention;
图5为本发明实施例提供的超声成像系统的结构图;FIG. 5 is a structural diagram of an ultrasound imaging system according to an embodiment of the present invention; FIG.
图6为本发明实施例提供的超声成像设备的结构图。FIG. 6 is a structural diagram of an ultrasound imaging apparatus according to an embodiment of the present invention.
图标:icon:
200-超声成像系统;210-获取模块;220-转换模块;230-计算模块;240-成像模块;300-超声成像设备;310-探头;311-阵元;320-显示屏。200-ultrasound imaging system; 210-acquisition module; 220-conversion module; 230-calculation module; 240-imaging module; 300-ultrasonic imaging device; 310-probe; 311-array; 320-display.
具体实施方式Detailed ways
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The embodiments of the present invention will be clearly and completely described in detail with reference to the accompanying drawings. An embodiment. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative efforts are within the scope of the present invention.
目前针对超声成像只是公开了线阵多角度平面波相关复合成像算法,其中,线阵指的是超声换能器中的多个阵元排列成线性的结构,而线阵多角度平面波相干复合成像的方法并不适用于凸阵多角度或者凸阵单角度平面波相干复合成像的方法,基于此,本发明实施例提供的一种超声成像方法、系统以及设备,通过获取笛卡尔坐标系下的像素点的坐标、该像素点的接收阵元的坐标、该像素点的发射阵元的坐标,再将笛卡尔坐标系下的接收阵元的坐标转换为极坐标系下的接收阵元的坐标,根据坐标值计算像素点的传播总距离从而计算该像素点的延时时间,在通过延时时间和回波数据得到波束合成下的该像素点,这样在多个阵元以凸阵形式排列时,在计算方式上通过将在笛卡尔坐标系下的接收阵元的坐标转换为极坐标系下的接收阵元的坐标,即接收阵元的坐标值由线性坐标转换为带有角度的坐标,就能够完成凸阵多角度或者单角度平面波发射下的波束的合成,实现了超快超声成像。At present, only the linear array multi-angle plane wave correlation composite imaging algorithm is disclosed for ultrasound imaging. The linear array refers to the linear array structure of multiple array elements in the ultrasonic transducer, and the linear array multi-angle plane wave coherent composite imaging. The method is not applicable to the method of convex array multi-angle or convex array single-angle plane wave coherent composite imaging. Based on this, an ultrasound imaging method, system and device provided by the embodiments of the present invention obtain pixel points in a Cartesian coordinate system. The coordinates, the coordinates of the receiving element of the pixel, the coordinates of the transmitting element of the pixel, and the coordinates of the receiving element in the Cartesian coordinate system are converted to the coordinates of the receiving element in the polar coordinate system, according to The coordinate value calculates the total distance of propagation of the pixel to calculate the delay time of the pixel, and obtains the pixel under beamforming through the delay time and the echo data, so that when the plurality of array elements are arranged in a convex array form, In the calculation method, the coordinates of the receiving array element in the Cartesian coordinate system are converted into the coordinates of the receiving array element in the polar coordinate system, that is, the receiving array element Standard values are converted by the linear coordinate angled coordinates, it is possible to complete the synthesis of the beam angles of the convex array or a single plane wave emission angle, to achieve ultrafast ultrasound imaging.
为便于对本实施例进行理解,首先对本发明实施例所公开的一种超声成像方法进行详细介绍,首先,结合图6所示,该方法应用于超声成像设备300,超声成像设备300包括探头310,探头310包括多个阵元311,多 个阵元311用于发射超声波对物体进行扫描,并接收回波数据。超声成像设备300具体的工作流程为:在启动超声成像设备300后,多个阵元311会发射平面波对所测的物体进行扫描,发射的平面波在遇到物体后,经所测的物体反射,多个阵元311会接收反射回的超声波,即多个阵元311接收回波数据,其中,回波数据记载所测的物体反射超声波波阵面的振幅和相位的信息。In order to facilitate the understanding of the present embodiment, an ultrasound imaging method disclosed in the embodiment of the present invention is first introduced in detail. First, as shown in FIG. 6, the method is applied to an ultrasound imaging apparatus 300, and the ultrasound imaging apparatus 300 includes a probe 310. The probe 310 includes a plurality of array elements 311 for transmitting ultrasonic waves to scan an object and receiving echo data. The specific workflow of the ultrasound imaging apparatus 300 is: after the ultrasound imaging apparatus 300 is activated, the plurality of array elements 311 emit a plane wave to scan the measured object, and the emitted plane wave is reflected by the measured object after encountering the object. The plurality of array elements 311 receive the reflected ultrasonic waves, that is, the plurality of array elements 311 receive the echo data, wherein the echo data records information about the amplitude and phase of the reflected ultrasonic wavefront reflected by the measured object.
参考图1所示,超声成像方法包括:Referring to Figure 1, the ultrasound imaging method includes:
对每个像素点进行以下步骤:Perform the following steps for each pixel:
S110:获取多个阵元接收的回波数据。S110: Acquire echo data received by multiple array elements.
作为一个示例,结合图6所示,多个阵元311接收回波数据,多个阵元311通过与其相连的通道进行上传,以能够获取这些回波数据进行后期处理。As an example, in conjunction with FIG. 6, a plurality of array elements 311 receive echo data, and a plurality of array elements 311 are uploaded through channels connected thereto to enable acquisition of the echo data for post processing.
S120:获取在笛卡尔坐标系下像素点的坐标、接收像素点回波数据的接收阵元的坐标、发射超声波到所述像素点的发射阵元的坐标。S120: Acquire coordinates of a pixel point in a Cartesian coordinate system, coordinates of a receiving array element that receives pixel point echo data, and coordinates of a transmitting array element that transmits an ultrasonic wave to the pixel point.
作为一个示例,以待合成的像素点P1为例,在凸阵探头下的多个阵元311进行无偏转扫描或者偏转扫描下,首先,在多个阵元311形成的曲面的平面内,设立一个笛卡尔坐标系,在这个笛卡尔坐标系下获取像素点P1的坐标(x,z),在这个笛卡尔坐标系下获取接收像素点P1回波数据的接收阵元E n的坐标,获取接收像素点P1回波数据的接收阵元A1的坐标。 As an example, taking the pixel point P1 to be synthesized as an example, under the non-deflection scanning or the deflection scanning of the plurality of array elements 311 under the convex array probe, first, in the plane of the curved surface formed by the plurality of array elements 311, a Cartesian coordinate system, in which the coordinates (x, z) of the pixel point P1 are obtained, and in this Cartesian coordinate system, the coordinates of the receiving array element E n of the received pixel point P1 echo data are acquired. The coordinates of the receiving array element A1 receiving the pixel point P1 echo data.
值得注意的是,这里所说的凸阵探头只是示例性的,例如腔内探头等等弧形探头均在本发明的限定中。It is to be noted that the convex array probes referred to herein are merely exemplary, and arcuate probes such as intracavity probes and the like are all within the definition of the present invention.
在一些实施例中,在进行步骤S120之前,还应该包括:对回波数据进行滤波,将不属于回波数据的其他超声波进行滤除,这样能够净化回波数据从而可以提高后期处理的正确率。In some embodiments, before performing step S120, it should also include: filtering the echo data to filter out other ultrasonic waves that are not part of the echo data, so that the echo data can be purified, thereby improving the accuracy of post processing. .
S130:将在笛卡尔坐标系下接收阵元的坐标转换为在极坐标系下接收阵元的坐标。S130: Convert the coordinates of the receiving element in the Cartesian coordinate system to the coordinates of the receiving element in the polar coordinate system.
在一些实施例中,将在笛卡尔坐标系下接收阵元的位置转换为在极坐标下接收阵元的位置,具体包括:计算接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,其中,笛卡尔坐标系为以多个阵元形成的曲面圆心为坐标零点,曲面的中心线所在的轴为第一方向轴,与曲面的中心线垂直的轴为第二方向轴;根据角度,将在笛卡尔坐标系下接收阵元的位置转换为在极坐标下接收阵元的位置。In some embodiments, converting the position of the receiving element in the Cartesian coordinate system to the position of receiving the element in polar coordinates comprises: calculating a radius of curvature of the receiving element and a first direction axis of the Cartesian coordinate system The angle of the Cartesian coordinate system is a coordinate zero point formed by a plurality of array elements, the axis of the curved line is the first direction axis, and the axis perpendicular to the center line of the curved surface is the second direction axis; Depending on the angle, the position of the receiving element in the Cartesian coordinate system is converted to the position at which the element is received in polar coordinates.
具体来说,结合图2所示,笛卡尔坐标系以多个阵元311形成的曲面圆形C1(0,0)为坐标零点,笛卡尔坐标系的x轴是垂直于曲面中心线,笛卡尔坐标系的z轴为曲面的中心线所在的轴。在这个坐标系下,计算接收阵元E n与z轴之间的夹角β n,其中,角度的算式为: Specifically, as shown in FIG. 2, the Cartesian coordinate system has a curved circle C1 (0, 0) formed by a plurality of array elements 311 as a coordinate zero point, and an x-axis of a Cartesian coordinate system is perpendicular to a center line of the curved surface, and a flute The z-axis of the Carl coordinate system is the axis of the centerline of the surface. In this coordinate system, the angle β n between the receiving element E n and the z axis is calculated, wherein the angle is calculated as:
Figure PCTCN2018070069-appb-000003
Figure PCTCN2018070069-appb-000003
其中,β n为接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,num为多个阵元的个数,pitch为多个阵元的间隔距离,R为曲率半径,n为第n个接收阵元,n的取值范围为[1,num]。 Where β n is the angle of curvature of the receiving array element and the first direction axis of the Cartesian coordinate system, num is the number of multiple array elements, pitch is the separation distance of multiple array elements, and R is the radius of curvature, n For the nth receiving array element, n has a value range of [1, num].
在结合图3所示,想要将在笛卡尔坐标系下的接收阵元E n的坐标转换为极坐标系下的坐标,实际上就是将具体的坐标值使用角度来表示,即, As shown in FIG. 3, in order to convert the coordinates of the receiving element E n in the Cartesian coordinate system into the coordinates in the polar coordinate system, the specific coordinate values are actually expressed by using angles, that is,
E nx=R*sinβ n,E nz=R*cosβ nE nx =R*sinβ n , E nz =R*cosβ n ,
其中,R为曲率半径,结合上面的角度的算式,可以得到:Where R is the radius of curvature, combined with the above formula, you can get:
Figure PCTCN2018070069-appb-000004
Figure PCTCN2018070069-appb-000004
Figure PCTCN2018070069-appb-000005
Figure PCTCN2018070069-appb-000005
综上,根据上述的公式,当要计算接收阵元E n的坐标时,只需要知道num、pitch、R的数据。 In summary, according to the above formula, when the coordinates of the receiving array element E n are to be calculated, only the data of num, pitch, and R need to be known.
S140:根据在笛卡尔坐标系下像素点和发射阵元的坐标、在极坐标系下接收阵元的坐标,计算像素点的延时时间。S140: Calculate the delay time of the pixel point according to the coordinates of the pixel point and the transmitting array element in the Cartesian coordinate system and the coordinates of the receiving element in the polar coordinate system.
在一些实施例中,步骤S140包括:设定像素点与距离多个阵元形成的曲面沿曲率半径下最近的点为发射阵元,其中,计算像素点的延时时间的算式为:In some embodiments, step S140 includes: setting a pixel point and a point formed by a plurality of array elements along a radius of curvature along a radius of curvature as a transmitting array element, wherein calculating a delay time of the pixel point is:
τ(E n,x m,z m)*c=sqrt((x m-C x) 2+(z m-C z) 2)-R+sqrt((x m-E nx) 2+(z m+E nz) 2)其中,(x m,z m)为在笛卡尔坐标系下像素点P m的坐标,(E nx,E nz)为在极坐标系下接收阵元E n的坐标,(C x,C z)为在笛卡尔坐标系曲面圆心C的坐标。 τ(E n ,x m ,z m )*c=sqrt((x m -C x ) 2 +(z m -C z ) 2 )-R+sqrt((x m -E nx ) 2 +(z m + E nz ) 2 ) where (x m , z m ) is the coordinate of the pixel point P m in the Cartesian coordinate system, and (E nx , E nz ) is the coordinate of the receiving element E n in the polar coordinate system , (C x , C z ) is the coordinate of the center C of the surface of the Cartesian coordinate system.
其中,像素点的延时时间是发射阵元发射超声波到像素点的传播时间和接收阵元接收从像素点反射回来的反射超声波的传播时间的总和。The delay time of the pixel is the sum of the propagation time of the transmitting array element transmitting the ultrasonic wave to the pixel point and the propagation time of the receiving ultrasonic wave receiving the reflected ultrasonic wave reflected from the pixel point.
举例来说,当计算接收阵元E n中的像素点P1时,以像素点P1和接收阵元E n的位置为基础,设定像素点P1与距离多个阵元形成的曲面沿曲率半径下最近的点为发射阵元。结合图2所示,在凸阵探头无偏转下,接收阵元E n中的像素点P1,像素点P1的接收阵元为A1。结合图4所示,在凸阵探头偏转θ进行扫描时,接收阵元E n中的像素点P1,像素点P1的接收阵元为A2。 For example, when the pixel P1 is calculated in the receiving array element E n to the pixel position of point P1 and the receiving array element E n on the basis of the curved surface along a curvature radius P1 and the pixel array element is formed from a plurality of The next closest point is the transmit array element. In conjunction with FIG. 2, in a convex array probe under non-deflection, the pixel P1 receives E n in the array element, the array element receiving point P1 of the pixel A1. As shown in FIG. 4, when the convex array probe deflection θ is scanned, the pixel point P1 in the array element E n is received, and the receiving array element of the pixel point P1 is A2.
在确定接收阵元、发射阵元之后,根据接收阵元和发射阵元的坐标计算与像素点的延时时间,由上述延时时间的算式来看,当确定接收阵元和发射阵元的位置时,这时,只要确认曲面圆心的坐标就可以计算像素点的延时时间,其中,在笛卡尔坐标系曲面圆心的坐标,计算公式为:After determining the receiving array element and the transmitting array element, calculating the delay time of the pixel point according to the coordinates of the receiving array element and the transmitting array element, according to the calculation formula of the delay time, when determining the receiving array element and the transmitting array element At the time of position, the delay time of the pixel point can be calculated by confirming the coordinates of the center of the surface of the surface. The coordinates of the center of the surface of the Cartesian coordinate system are calculated as:
C x=(-1)*(R 2+R 2-2*R*R*COSθ)*cos((num/2+0.5)*pitch/R-θ/2) C x =(-1)*(R 2 +R 2 -2*R*R*COSθ)*cos((num/2+0.5)*pitch/R-θ/2)
C z=(-1)*(R 2+R 2-2*R*R*COSθ)*sin((num/2+0.5)*pitch/R-θ/2) C z =(-1)*(R 2 +R 2 -2*R*R*COSθ)*sin((num/2+0.5)*pitch/R-θ/2)
其中,(C x,C z)为所述圆心在笛卡尔坐标系下的坐标,R为曲率半径,pitch为多个阵元的间隔距离,num为多个阵元的个数,θ为在多个阵元偏转扫射角度。 Where (C x , C z ) is the coordinate of the center of the circle in the Cartesian coordinate system, R is the radius of curvature, pitch is the separation distance of multiple array elements, num is the number of multiple array elements, and θ is Multiple array elements deflect the sweep angle.
作为一个示例,结合图2所示,当在笛卡尔坐标系曲面圆心C的坐标为笛卡尔坐标系的坐标零点时,(x 1,z 1)为在笛卡尔坐标系下像素点P 1的坐标,圆心C(C x,C z)为C(0,0),则上述计算延时时间的算式为: As an example, as shown in FIG. 2, when the coordinates of the center C of the Cartesian coordinate system are the coordinate zeros of the Cartesian coordinate system, (x 1 , z 1 ) is the pixel point P 1 in the Cartesian coordinate system. The coordinates, the center of the circle C (C x , C z ) is C (0, 0), then the above calculation of the delay time is:
τ(E n,x 1,z 1)*c=sqrt(x 1 2+z 1 2)-R+sqrt((x 1-E nx) 2+(z 1+E nz) 2)。 τ(E n , x 1 , z 1 )*c=sqrt(x 1 2 +z 1 2 )-R+sqrt((x 1 -E nx ) 2 +(z 1 +E nz ) 2 ).
作为另一个示例,结合图4所示,多个阵元偏转扫射θ角度时,多个阵元311形成的曲面的圆心C的坐标为:As another example, in conjunction with FIG. 4, when a plurality of array elements deflect the sweep θ angle, the coordinates of the center C of the curved surface formed by the plurality of array elements 311 are:
C x=(-1)*(R 2+R 2-2*R*R*COSθ)*cos((num/2+0.5)*pitch/R-θ/2) C x =(-1)*(R 2 +R 2 -2*R*R*COSθ)*cos((num/2+0.5)*pitch/R-θ/2)
C z=(-1)*(R 2+R 2-2*R*R*COSθ)*sin((num/2+0.5)*pitch/R-θ/2) C z =(-1)*(R 2 +R 2 -2*R*R*COSθ)*sin((num/2+0.5)*pitch/R-θ/2)
可以看出,圆心C的坐标可以通过偏转角度θ,多个阵元的个数、多个阵元的间隔距离和曲率半径计算而来。当知道上述参数时,就可以计算在任何偏转角度下的圆心C的坐标。It can be seen that the coordinates of the center C can be calculated by the deflection angle θ, the number of the plurality of array elements, the separation distance of the plurality of array elements, and the radius of curvature. When the above parameters are known, the coordinates of the center C at any deflection angle can be calculated.
S150:根据延时时间和回波数据,得到波束合成下的像素点。S150: Obtain a pixel point under beamforming according to the delay time and the echo data.
步骤S150包括:根据采集频率和延时时间,查找多个阵元311接收到的来自像素点的回波数据;对回波数据进行加权求和,得到波束合成下的像素点。Step S150 includes: searching for echo data from the pixel points received by the plurality of array elements 311 according to the acquisition frequency and the delay time; and performing weighted summation on the echo data to obtain pixel points under beam synthesis.
具体来说,多个阵元311形成的探头310对回波数据采集时,采用一定的频率进行周期性的采集。以阵元E 1,采集频率为1s,总采集时间为10s为例,在采集结束后,阵元E 1中包括10次采集的多个像素点的回波数据。由于接收每个像素点的回波数据的时间不同,阵元E 1接收像素点的回波时间为该像素点传播的总时间,即为该像素点的延时时间,计算该像素点的延时时间就可以找到阵元E 1中该像素点的回波数据。利用这种方式,在多个阵元311中找到该像素点的多个回波数据,然后对该像素点的多个回波数据进行加权求和,即将该像素点的光信号转换后的电信号进行加权求和,实际上,从超声波的角度来看,就是将该像素点的所有波束进行叠加得到一个叠加后的波束。其中,加权函数为窗函数,窗函数可以为汉宁窗、汉明窗、矩形窗函数等。 Specifically, when the probe 310 formed by the plurality of array elements 311 collects echo data, periodic acquisition is performed using a certain frequency. Taking the array element E 1 , the acquisition frequency is 1 s, and the total acquisition time is 10 s as an example, after the acquisition is completed, the array element E 1 includes echo data of a plurality of collected pixel points 10 times. Since the time of receiving the echo data of each pixel is different, the echo time of the pixel E 1 receiving the pixel is the total time of the pixel propagation, that is, the delay time of the pixel, and the delay of the pixel is calculated. At the time, the echo data of the pixel in the element E 1 can be found. In this way, a plurality of echo data of the pixel point are found in the plurality of array elements 311, and then the plurality of echo data of the pixel point are weighted and summed, that is, the converted optical signal of the pixel point is converted. The signals are weighted and summed. In fact, from the perspective of the ultrasound, all the beams of the pixel are superimposed to obtain a superimposed beam. Wherein, the weighting function is a window function, and the window function may be a Hanning window, a Hamming window, a rectangular window function, or the like.
值得注意的是,在对每个像素点进行上述处理时,可以设置成像深度,限制像素点的范围,以能够得到最优质的图像。根据上述过程,对每一个像素点进行处理,获得一幅波束合成下的图像。It is worth noting that when performing the above processing for each pixel, the imaging depth can be set to limit the range of pixels to obtain the highest quality image. According to the above process, each pixel is processed to obtain an image under beamforming.
参见图5所示的超声成像系统200的结构图,包括:获取模块210、转换模块220、计算模块230、成像模块240。Referring to the structural diagram of the ultrasound imaging system 200 shown in FIG. 5, the acquisition module 210, the conversion module 220, the calculation module 230, and the imaging module 240 are included.
其中,获取模块210用于获取多个阵元接收的回波数据,并获取在笛卡尔坐标系下像素点的坐标、接收像素点回波数据的接收阵元的坐标、发射超声波到像素点的发射阵元的坐标。转换模块220与获取模块210相连,用于将在笛卡尔坐标系下接收阵元的坐标转换为在极坐标系下接收阵元的坐标。计算模块230与转换模块220相连,用于根据在笛卡尔坐标系下像素点和发射阵元的坐标、在极坐标下接收阵元的坐标,计算像素点的延时时间。成像模块240分别与获取模块210和计算模块230相连,用于根据延时时间和回波数据,得到波束合成下的像素点强度值。The acquiring module 210 is configured to acquire echo data received by multiple array elements, and acquire coordinates of pixel points in a Cartesian coordinate system, coordinates of receiving array elements that receive pixel point echo data, and transmit ultrasonic waves to pixel points. The coordinates of the transmitting element. The conversion module 220 is coupled to the acquisition module 210 for converting the coordinates of the receiving array elements in the Cartesian coordinate system to the coordinates of the receiving array elements in the polar coordinate system. The calculation module 230 is connected to the conversion module 220 for calculating the delay time of the pixel point according to the coordinates of the pixel point and the transmission element in the Cartesian coordinate system and the coordinates of the receiving element in the polar coordinate. The imaging module 240 is connected to the acquisition module 210 and the calculation module 230, respectively, for obtaining pixel intensity values under beamforming according to the delay time and the echo data.
超声成像系统200还包括:处理模块,处理模块和成像模块240相连,处理模块用于将成像模块240的结果经过检波、滤波、对数压缩、动态范围变换等等的数字信号处理后,发送至显示屏。The ultrasound imaging system 200 further includes a processing module, and the processing module is connected to the imaging module 240. The processing module is configured to process the digital signal of the result of the imaging module 240 after detection, filtering, logarithmic compression, dynamic range conversion, etc., to the digital signal. Display.
在一些实施例中,转换模块220,具体用于:计算接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,其中,笛卡尔坐标系为以多个阵元形成的曲面圆心为坐标零点,曲面的中心线所在的轴为第一方向轴,与曲面的中心线垂直的轴为第二方向轴;根据角度,将在笛卡尔坐标系下接收阵元的位置转换为在极坐标下接收阵元的位置。In some embodiments, the conversion module 220 is specifically configured to: calculate an angle of a radius of curvature of the receiving array element and a first direction axis of the Cartesian coordinate system, wherein the Cartesian coordinate system is a surface center formed by a plurality of array elements For the coordinate zero point, the axis of the curved surface is the first direction axis, and the axis perpendicular to the center line of the surface is the second direction axis; according to the angle, the position of the receiving element in the Cartesian coordinate system is converted to the pole The position of the receiving element in coordinates.
在一些实施例中,角度的算式为:In some embodiments, the angle is calculated as:
Figure PCTCN2018070069-appb-000006
Figure PCTCN2018070069-appb-000006
其中,β n为接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,num为多个阵元的个数,pitch为多个阵元的间隔距离,R为曲率半径,n为第n个接收阵元,n的取值范围为(1,num)。 Where β n is the angle of curvature of the receiving array element and the first direction axis of the Cartesian coordinate system, num is the number of multiple array elements, pitch is the separation distance of multiple array elements, and R is the radius of curvature, n For the nth receiving array element, n has a value range of (1, num).
本发明实施例所提供的装置,其实现原理及产生的技术效果和前述方法实施例相同,为简要描述,装置实施例部分未提及之处,可参考前述方法实施例中相应内容。The implementation principle and the technical effects of the device provided by the embodiments of the present invention are the same as those of the foregoing method embodiments. For a brief description, where the device embodiment is not mentioned, reference may be made to the corresponding content in the foregoing method embodiments.
参见图6所示,本发明实施例的超声成像设备300,包括:探头310、超声成像系统200、显示屏320。Referring to FIG. 6, an ultrasound imaging apparatus 300 according to an embodiment of the present invention includes a probe 310, an ultrasound imaging system 200, and a display screen 320.
其中,探头310,探头310包括多个阵元311,多个阵元311用于发射超声波对物体进行扫描,并接收回波数据;根据上述实施例任一项所述的超声成像系统200。显示屏320与超声成像系统200相连,用于显示超声成像系统200的成像结果。值得注意的是,多个阵元311通过通道与超声成像系统200相连。The probe 310 includes a plurality of array elements 311 for transmitting ultrasonic waves to scan an object and receiving echo data. The ultrasonic imaging system 200 according to any of the above embodiments. Display screen 320 is coupled to ultrasound imaging system 200 for displaying imaging results of ultrasound imaging system 200. It is worth noting that a plurality of array elements 311 are connected to the ultrasound imaging system 200 through channels.
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的设备的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。A person skilled in the art can clearly understand that for the convenience and brevity of the description, the specific working process of the device described above can refer to the corresponding process in the foregoing method embodiments, and details are not described herein again.
除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对步骤、数字表达式和数值并不限制本发明的范围。The relative steps, numerical expressions and numerical values of the components and steps set forth in the examples are not intended to limit the scope of the invention.
在这里示出和描述的所有示例中,任何具体值应被解释为仅仅是示例性的,而不是作为限制,因此,示例性实施例的其他示例可以具有不同的值。In all of the examples shown and described herein, any specific values should be construed as merely exemplary, and not as a limitation, and thus, other examples of the exemplary embodiments may have different values.
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。It should be noted that similar reference numerals and letters indicate similar items in the following figures, and therefore, once an item is defined in a drawing, it is not necessary to further define and explain it in the subsequent drawings.
附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程 图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。The flowchart and block diagrams in the Figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of the invention. In this regard, each block of the flowchart or block diagram can represent a module, a program segment, or a portion of code that includes one or more of the Executable instructions. It should also be noted that in some alternative implementations, the functions noted in the blocks may also occur in a different order than that illustrated in the drawings. For example, two consecutive blocks may be executed substantially in parallel, and they may sometimes be executed in the reverse order, depending upon the functionality involved. It is also noted that each block of the block diagrams and/or flowcharts, and combinations of blocks in the block diagrams and/or flowcharts, can be implemented in a dedicated hardware-based system that performs the specified function or function. Or it can be implemented by a combination of dedicated hardware and computer instructions.
在发明的实施例超声成像方法中,波束的合成可以在处理器中完成,该处理器可能为图像处理器,还可以是一种集成电路芯片,具有信号的处理能力。在实现过程中,上述方法的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器可以是通用处理器,包括中央处理器(Central Processing Unit,简称CPU)、网络处理器(Network Processor,简称NP)等;还可以是数字信号处理器(Digital Signal Processing,简称DSP)、专用集成电路(Application Specific Integrated Circuit,简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array,简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成,或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器,处理器读取存储器中的信息,结合其硬件完成上述方法的步骤。In the ultrasound imaging method of the embodiment of the invention, the synthesis of the beam may be performed in a processor, which may be an image processor, or an integrated circuit chip having signal processing capabilities. In the implementation process, each step of the above method may be completed by an integrated logic circuit of hardware in a processor or an instruction in a form of software. The above processor may be a general-purpose processor, including a central processing unit (CPU), a network processor (Network Processor, NP for short), or a digital signal processor (DSP). , Application Specific Integrated Circuit (ASIC), Field-Programmable Gate Array (FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware component. The methods, steps, and logical block diagrams disclosed in the embodiments of the present invention may be implemented or carried out. The general purpose processor may be a microprocessor or the processor or any conventional processor or the like. The steps of the method disclosed in the embodiments of the present invention may be directly implemented by the hardware decoding processor, or may be performed by a combination of hardware and software modules in the decoding processor. The software module can be located in a conventional storage medium such as random access memory, flash memory, read only memory, programmable read only memory or electrically erasable programmable memory, registers, and the like. The storage medium is located in the memory, and the processor reads the information in the memory and combines the hardware to complete the steps of the above method.
另外,在本发明实施例的描述中,除非另有明确的规定和限定,术语“相连”、应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一 体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。In addition, in the description of the embodiments of the present invention, the term "connected" shall be understood broadly, and may be, for example, a fixed connection, a detachable connection, or an integral connection, unless otherwise explicitly defined and defined; It is a mechanical connection or an electrical connection; it can be directly connected or indirectly connected through an intermediate medium, which can be the internal connection of two components. The specific meaning of the above terms in the present invention can be understood in a specific case by those skilled in the art.
在本发明的描述中,需要说明的是,术语“上”、“下”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。In the description of the present invention, it is to be noted that the orientation or positional relationship of the terms "upper", "lower", "inside", "outside", etc. is based on the orientation or positional relationship shown in the drawings, only for the purpose of The invention is not limited by the scope of the invention, and is not intended to be a limitation of the invention. Moreover, the terms "first," "second," and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。Finally, it should be noted that the above-mentioned embodiments are merely specific embodiments of the present invention, and are used to explain the technical solutions of the present invention, and are not limited thereto, and the scope of protection of the present invention is not limited thereto, although reference is made to the foregoing. The present invention has been described in detail, and those skilled in the art should understand that any one skilled in the art can still modify the technical solutions described in the foregoing embodiments within the technical scope disclosed by the present invention. The changes may be easily conceived, or equivalents may be substituted for some of the technical features. The modifications, variations, or substitutions of the present invention are not intended to depart from the spirit and scope of the technical solutions of the embodiments of the present invention. Within the scope of protection. Therefore, the scope of the invention should be determined by the scope of the claims.

Claims (10)

  1. 一种超声成像方法,其特征在于,应用于超声成像设备,所述超声成像设备包括探头,所述探头包括多个阵元,所述多个阵元用于发射超声波对物体进行扫描,并接收回波数据,所述超声成像方法包括:An ultrasound imaging method, characterized in that it is applied to an ultrasound imaging apparatus, the ultrasound imaging apparatus comprising a probe, the probe comprising a plurality of array elements, the plurality of array elements being used for transmitting ultrasonic waves to scan an object, and receiving Echo data, the ultrasound imaging method includes:
    对每个像素点进行以下步骤:Perform the following steps for each pixel:
    获取多个阵元接收的回波数据;Acquiring echo data received by multiple array elements;
    获取在笛卡尔坐标系下像素点的坐标、接收所述像素点回波数据的接收阵元的坐标、以及发射超声波到所述像素点的发射阵元的坐标;Obtaining coordinates of a pixel point in a Cartesian coordinate system, coordinates of a receiving array element receiving the pixel point echo data, and coordinates of a transmitting array element transmitting an ultrasonic wave to the pixel point;
    将在笛卡尔坐标系下接收阵元的坐标转换为在极坐标系下接收阵元的坐标;Converting the coordinates of the receiving element in the Cartesian coordinate system to the coordinates of the receiving element in the polar coordinate system;
    根据在笛卡尔坐标系下像素点和发射阵元的坐标、以及在极坐标系下接收阵元的坐标,计算像素点的延时时间;Calculating the delay time of the pixel point according to the coordinates of the pixel point and the transmitting element in the Cartesian coordinate system and the coordinates of the receiving element in the polar coordinate system;
    根据所述延时时间和回波数据,得到波束合成下的所述像素点强度值。According to the delay time and the echo data, the pixel intensity value under beamforming is obtained.
  2. 根据权利要求1所述的超声成像方法,其特征在于,所述将在笛卡尔坐标系下接收阵元的位置转换为在极坐标下接收阵元的位置,具体包括:The ultrasonic imaging method according to claim 1, wherein the converting the position of the receiving element in the Cartesian coordinate system to the position of receiving the element in the polar coordinate comprises:
    计算接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,其中,所述笛卡尔坐标系为以多个阵元形成的曲面圆心为坐标零点,所述曲面的中心线所在的轴为第一方向轴,与所述曲面的中心线垂直的轴为第二方向轴;Calculating an angle of a radius of curvature of the receiving array element and a first direction axis of the Cartesian coordinate system, wherein the Cartesian coordinate system is a coordinate center with a plurality of array elements as a coordinate zero point, where the center line of the curved surface is located The axis is a first direction axis, and an axis perpendicular to a center line of the curved surface is a second direction axis;
    根据所述角度,将在笛卡尔坐标系下接收阵元的位置转换为在极坐标下接收阵元的位置。Based on the angle, the position of the receiving element in the Cartesian coordinate system is converted to the position at which the element is received in polar coordinates.
  3. 根据权利要求2所述的超声成像方法,其特征在于,所述角度的算式为:The ultrasonic imaging method according to claim 2, wherein the angle is calculated as:
    Figure PCTCN2018070069-appb-100001
    Figure PCTCN2018070069-appb-100001
    其中,β n为接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,num为多个阵元的个数,pitch为多个阵元的间隔距离,R为曲率半径,n为第n个接收阵元,n的取值范围为[1,num]。 Where β n is the angle of curvature of the receiving array element and the first direction axis of the Cartesian coordinate system, num is the number of multiple array elements, pitch is the separation distance of multiple array elements, and R is the radius of curvature, n For the nth receiving array element, n has a value range of [1, num].
  4. 根据权利要求1-3任一项所述的超声成像方法,其特征在于,所述根据在笛卡尔坐标系下像素点和发射阵元的坐标、以及在极坐标系下接收阵元的坐标,计算像素点的延时时间,包括:The ultrasonic imaging method according to any one of claims 1 to 3, wherein the coordinates of the pixel points and the transmitting array elements in the Cartesian coordinate system and the coordinates of the receiving element elements in the polar coordinate system are Calculate the delay time of the pixel, including:
    设定所述像素点与距离多个阵元形成的曲面沿曲率半径下最近的点为发射阵元,其中,计算像素点的延时时间的算式为:The point that the pixel point and the surface formed by the plurality of array elements are closest to the radius of curvature is a transmitting array element, wherein the calculation formula of the delay time of the pixel point is:
    τ(E n,x m,z m)*c=sqrt((x m-C x) 2+(z m-C z) 2)-R+sqrt((x m-E nx) 2+(z m+E nz) 2)其中,(x m,z m)为在笛卡尔坐标系下像素点P m的坐标,(E nx,E nz)为在极坐标系下接收阵元E n的坐标,(C x,C z)为在笛卡尔坐标系曲面圆心C的坐标。 τ(E n ,x m ,z m )*c=sqrt((x m -C x ) 2 +(z m -C z ) 2 )-R+sqrt((x m -E nx ) 2 +(z m + E nz ) 2 ) where (x m , z m ) is the coordinate of the pixel point P m in the Cartesian coordinate system, and (E nx , E nz ) is the coordinate of the receiving element E n in the polar coordinate system , (C x , C z ) is the coordinate of the center C of the surface of the Cartesian coordinate system.
  5. 根据权利要求4所述的超声成像方法,其特征在于,所述在笛卡尔坐标系曲面圆心的坐标,计算公式为:The ultrasonic imaging method according to claim 4, wherein the coordinates of the center of the surface of the Cartesian coordinate system are calculated as:
    C x=(-1)*(R 2+R 2-2*R*R*COSθ)*cos((num/2+0.5)*pitch/R-θ/2) C x =(-1)*(R 2 +R 2 -2*R*R*COSθ)*cos((num/2+0.5)*pitch/R-θ/2)
    C z=(-1)*(R 2+R 2-2*R*R*COSθ)*sin((num/2+0.5)*pitch/R-θ/2) C z =(-1)*(R 2 +R 2 -2*R*R*COSθ)*sin((num/2+0.5)*pitch/R-θ/2)
    其中,(C x,C z)为所述圆心在笛卡尔坐标系下的坐标,R为曲率半径,pitch为多个阵元的间隔距离,num为多个阵元的个数,θ为在多个阵元偏转扫射角度。 Where (C x , C z ) is the coordinate of the center of the circle in the Cartesian coordinate system, R is the radius of curvature, pitch is the separation distance of multiple array elements, num is the number of multiple array elements, and θ is Multiple array elements deflect the sweep angle.
  6. 根据权利要求1所述的超声成像方法,其特征在于,所述根据所述延时时间和回波数据,得到波束合成下的所述像素点,包括:The ultrasound imaging method according to claim 1, wherein the obtaining the pixel points under beamforming according to the delay time and echo data comprises:
    根据采集频率和延时时间,查找多个阵元接收到的来自所述像素点的回波数据;Finding echo data from the pixel points received by the plurality of array elements according to the acquisition frequency and the delay time;
    对所述回波数据进行加权求和,得到波束合成下的所述像素点强度值。The echo data is weighted and summed to obtain the pixel intensity values under beamforming.
  7. 一种超声成像系统,其特征在于,包括:An ultrasound imaging system, comprising:
    获取模块,用于获取多个阵元接收的回波数据,并获取在笛卡尔坐标系下像素点的坐标、接收所述像素点回波数据的接收阵元的坐标、以及发射超声波到所述像素点的发射阵元的坐标;An acquiring module, configured to acquire echo data received by the plurality of array elements, acquire coordinates of pixel points in a Cartesian coordinate system, coordinates of receiving array elements that receive the echo data of the pixel points, and transmit ultrasonic waves to the The coordinates of the firing element of the pixel;
    转换模块,所述转换模块与所述获取模块相连,用于将在笛卡尔坐标系下接收阵元的坐标转换为在极坐标系下接收阵元的坐标;a conversion module, the conversion module being connected to the acquisition module, configured to convert coordinates of receiving array elements in a Cartesian coordinate system into coordinates of receiving array elements in a polar coordinate system;
    计算模块,所述计算模块与所述转换模块相连,用于根据在笛卡尔坐标系下像素点和发射阵元的坐标、以及在极坐标下接收阵元的坐标,计算像素点的延时时间;a calculation module, the calculation module being connected to the conversion module, configured to calculate a delay time of the pixel point according to the coordinates of the pixel point and the transmission element in the Cartesian coordinate system and the coordinates of the receiving element in the polar coordinate ;
    成像模块,所述成像模块分别与所述获取模块和计算模块相连,用于根据所述延时时间和回波数据,得到波束合成下的所述像素点强度值。An imaging module is respectively connected to the obtaining module and the calculating module, and configured to obtain the pixel intensity value under beam combining according to the delay time and echo data.
  8. 根据权利要求7所述的超声成像系统,其特征在于,所述转换模块,具体用于:计算接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,其中,所述笛卡尔坐标系为以多个阵元形成的曲面圆心为坐标零点,所述曲面的中心线所在的轴为第一方向轴,与所述曲面的中心线垂直的轴为第二方向轴;根据所述角度,将在笛卡尔坐标系下接收阵元的位置转换为在极坐标下接收阵元的位置。The ultrasound imaging system according to claim 7, wherein the conversion module is specifically configured to: calculate an angle of a radius of curvature of the receiving array element and a first direction axis of the Cartesian coordinate system, wherein the Cartesian The coordinate system is a coordinate zero point formed by a plurality of array elements, wherein the axis of the curved surface is the first direction axis, and the axis perpendicular to the center line of the curved surface is the second direction axis; Angle, the position of the receiving element in the Cartesian coordinate system is converted to the position of the receiving element in polar coordinates.
  9. 根据权利要求8所述的超声成像方法,其特征在于,所述角度的算式为:The ultrasonic imaging method according to claim 8, wherein the angle is calculated as:
    Figure PCTCN2018070069-appb-100002
    Figure PCTCN2018070069-appb-100002
    其中,β n为接收阵元的曲率半径与笛卡尔坐标系的第一方向轴的角度,num为多个阵元的个数,pitch为多个阵元的间隔距离,R为曲率半径,n为第n个接收阵元,n的取值范围为(1,num)。 Where β n is the angle of curvature of the receiving array element and the first direction axis of the Cartesian coordinate system, num is the number of multiple array elements, pitch is the separation distance of multiple array elements, and R is the radius of curvature, n For the nth receiving array element, n has a value range of (1, num).
  10. 一种超声成像设备,其特征在于,包括:An ultrasonic imaging apparatus, comprising:
    探头,所述探头包括多个阵元,所述多个阵元用于发射超声波对物体进行扫描,并接收回波数据;a probe comprising a plurality of array elements, wherein the plurality of array elements are configured to emit ultrasonic waves to scan an object and receive echo data;
    根据权利要求7-9任一项所述的超声成像系统;An ultrasound imaging system according to any of claims 7-9;
    显示屏,所述显示屏与所述超声成像系统相连,用于显示超声成像系统的成像结果。A display screen coupled to the ultrasound imaging system for displaying imaging results of the ultrasound imaging system.
PCT/CN2018/070069 2017-12-28 2018-01-02 Ultrasonic imaging method, system and device WO2019127621A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201711467316.1A CN109975814B (en) 2017-12-28 2017-12-28 Ultrasound imaging method, system and device
CN201711467316.1 2017-12-28

Publications (1)

Publication Number Publication Date
WO2019127621A1 true WO2019127621A1 (en) 2019-07-04

Family

ID=67062919

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2018/070069 WO2019127621A1 (en) 2017-12-28 2018-01-02 Ultrasonic imaging method, system and device

Country Status (2)

Country Link
CN (1) CN109975814B (en)
WO (1) WO2019127621A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112215829A (en) * 2020-10-21 2021-01-12 深圳度影医疗科技有限公司 Positioning method of hip joint standard tangent plane and computer equipment

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110477951B (en) * 2019-08-30 2020-08-25 浙江大学 Ultra-fast composite plane wave imaging method based on broadband acoustic metamaterial
CN111110279A (en) * 2020-01-13 2020-05-08 卓瑞姆生物技术有限公司 Ultrasonic imaging equipment and imaging method thereof
CN111551899A (en) * 2020-03-31 2020-08-18 河北科技大学 Ultrasonic positioning method, terminal device and storage medium
CN114081526A (en) * 2020-11-18 2022-02-25 武汉联影医疗科技有限公司 Ultrasonic imaging method, device, system and storage medium
CN115469015B (en) * 2022-10-31 2023-01-31 泉州装备制造研究所 Curved surface component ultrasonic self-adaptive detection method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1016875A2 (en) * 1998-12-29 2000-07-05 Diasonics Israel Ltd. Ultrasound system performing simultaneous parallel computer instructions
CN1492313A (en) * 2003-07-29 2004-04-28 深圳迈瑞生物医疗电子股份有限公司 Coordinate transformation method for digital scanning change-over device and processor
CN101683274A (en) * 2008-09-25 2010-03-31 深圳迈瑞生物医疗电子股份有限公司 Deflection convex array digital scanning conversion and convex array special compound method and imaging system
CN106651740A (en) * 2016-11-29 2017-05-10 中国科学院声学研究所 FPGA-based ultrasonic total data focusing quick imaging method and system
CN106680825A (en) * 2016-12-05 2017-05-17 中国科学院声学研究所 Acoustic array imaging system and method thereof
CN107102324A (en) * 2017-04-28 2017-08-29 华讯方舟科技有限公司 A kind of close shot microwave imaging method and system

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4344327B1 (en) * 1979-12-28 1994-05-03 Aloka Co Ltd Electronic scanning ultrasonic diagnostic system
KR100949066B1 (en) * 2006-09-08 2010-03-30 주식회사 메디슨 Ultrasound system and method for controlling scanline
CN101371789A (en) * 2007-08-22 2009-02-25 深圳迈瑞生物医疗电子股份有限公司 Real-time computing method and device for focalization delay parameters of beam-forming
CN101770650B (en) * 2009-01-07 2013-04-24 深圳迈瑞生物医疗电子股份有限公司 Method and device for three-dimensional ultrasonic real-time imaging and imaging system
US20100324418A1 (en) * 2009-06-23 2010-12-23 Essa El-Aklouk Ultrasound transducer
CN101596114A (en) * 2009-06-30 2009-12-09 深圳市蓝韵实业有限公司 A kind of digital scan converter
EP3960075A1 (en) * 2009-11-27 2022-03-02 Hologic, Inc. Systems and methods for tracking positions between imaging modalities and transforming a displayed three-dimensional image corresponding to a position and orientation of a probe
US9452447B2 (en) * 2013-12-27 2016-09-27 General Electric Company Ultrasound transducer and ultrasound imaging system with a variable thickness dematching layer
CN104739453A (en) * 2013-12-31 2015-07-01 深圳市鹏瑞智能技术应用研究院 Ultrasonic tomography system and method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1016875A2 (en) * 1998-12-29 2000-07-05 Diasonics Israel Ltd. Ultrasound system performing simultaneous parallel computer instructions
CN1492313A (en) * 2003-07-29 2004-04-28 深圳迈瑞生物医疗电子股份有限公司 Coordinate transformation method for digital scanning change-over device and processor
CN101683274A (en) * 2008-09-25 2010-03-31 深圳迈瑞生物医疗电子股份有限公司 Deflection convex array digital scanning conversion and convex array special compound method and imaging system
CN106651740A (en) * 2016-11-29 2017-05-10 中国科学院声学研究所 FPGA-based ultrasonic total data focusing quick imaging method and system
CN106680825A (en) * 2016-12-05 2017-05-17 中国科学院声学研究所 Acoustic array imaging system and method thereof
CN107102324A (en) * 2017-04-28 2017-08-29 华讯方舟科技有限公司 A kind of close shot microwave imaging method and system

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112215829A (en) * 2020-10-21 2021-01-12 深圳度影医疗科技有限公司 Positioning method of hip joint standard tangent plane and computer equipment
CN112215829B (en) * 2020-10-21 2021-12-14 深圳度影医疗科技有限公司 Positioning method of hip joint standard tangent plane and computer equipment

Also Published As

Publication number Publication date
CN109975814A (en) 2019-07-05
CN109975814B (en) 2020-09-22

Similar Documents

Publication Publication Date Title
WO2019127621A1 (en) Ultrasonic imaging method, system and device
EP2325672B1 (en) Spatial compound imaging in an ultrasound system
US8187192B2 (en) Method and apparatus for scan conversion and interpolation of ultrasonic linear array steering imaging
JP5473381B2 (en) Ultrasonic device
JP5357684B2 (en) Ultrasonic imaging apparatus and ultrasonic imaging method
US20060079780A1 (en) Ultrasonic imaging apparatus
US20210330293A1 (en) Ultrasound diagnostic device, ultrasound diagnostic method, and ultrasound diagnostic program
US10679349B2 (en) Method and system for estimating motion between images, particularly in ultrasound spatial compounding
US20070110291A1 (en) Image processing system and method for editing contours of a target object using multiple sectional images
JP2005152648A (en) Method and system for motion-adaptive type spatial synthesis
CN102697524B (en) Full-focus ultrasonic imaging method and application of method in blood flow imaging
US8343054B1 (en) Methods and apparatus for ultrasound imaging
JP3763924B2 (en) Ultrasonic diagnostic equipment
EP2466330B1 (en) Ultrasound system and method for processing beam-forming based on sampling data
US20110054323A1 (en) Ultrasound system and method for providing an ultrasound spatial compound image considering steering angle
US9140790B2 (en) Ultrasound system and method of forming ultrasound image
US7128712B2 (en) Adaptive ultrasound imaging system
JP4761999B2 (en) Ultrasonic diagnostic apparatus, image processing method thereof, and image processing program thereof
CN112998745A (en) Transmitting beam forming method and system for ultrasonic imaging and diagnostic equipment
JP2009000361A (en) Ultrasonic diagnostic device and ultrasonic image generation method
CN110613476B (en) Ultrasonic signal processing method and device
WO2020133332A1 (en) Method for determining ultrasonic emission angle, and ultrasonic device
US20200196990A1 (en) Ultrasonic imaging device and imaging method thereof
EP4338678A1 (en) Flexible ultrasound transducer
JP2013102959A (en) Ultrasonic diagnostic apparatus and method

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18894987

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 10/11/2020)

122 Ep: pct application non-entry in european phase

Ref document number: 18894987

Country of ref document: EP

Kind code of ref document: A1