CN106651740B - A method and system for fast imaging of ultrasound full data focusing based on FPGA - Google Patents

A method and system for fast imaging of ultrasound full data focusing based on FPGA Download PDF

Info

Publication number
CN106651740B
CN106651740B CN201611077700.6A CN201611077700A CN106651740B CN 106651740 B CN106651740 B CN 106651740B CN 201611077700 A CN201611077700 A CN 201611077700A CN 106651740 B CN106651740 B CN 106651740B
Authority
CN
China
Prior art keywords
port
pixel point
pixel
array element
ultrasonic signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201611077700.6A
Other languages
Chinese (zh)
Other versions
CN106651740A (en
Inventor
王冲
毛捷
冷涛
安志武
陈秋颖
宋波
周文佳
廉国选
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Acoustics CAS
Original Assignee
Institute of Acoustics CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Acoustics CAS filed Critical Institute of Acoustics CAS
Priority to CN201611077700.6A priority Critical patent/CN106651740B/en
Publication of CN106651740A publication Critical patent/CN106651740A/en
Application granted granted Critical
Publication of CN106651740B publication Critical patent/CN106651740B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00General purpose image data processing
    • G06T1/20Processor architectures; Processor configuration, e.g. pipelining
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F15/00Digital computers in general; Data processing equipment in general
    • G06F15/76Architectures of general purpose stored program computers
    • G06F15/78Architectures of general purpose stored program computers comprising a single central processing unit
    • G06F15/7807System on chip, i.e. computer system on a single chip; System in package, i.e. computer system on one or more chips in a single package
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F15/00Digital computers in general; Data processing equipment in general
    • G06F15/76Architectures of general purpose stored program computers
    • G06F15/78Architectures of general purpose stored program computers comprising a single central processing unit
    • G06F15/7839Architectures of general purpose stored program computers comprising a single central processing unit with memory
    • G06F15/7842Architectures of general purpose stored program computers comprising a single central processing unit with memory on one IC chip (single chip microcontrollers)

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明涉及一种基于FPGA的超声全数据聚焦快速成像方法及系统,该基于FPGA的超声全数据聚焦快速成像方法包括:采集由至少一个第一阵元的每一个阵元发射的第一超声波信号至像素点,及由像素点反射至少一个第二超声波信号至至少一个第二阵元,并将采集后的数据存储;根据第一超声波信号和至少一个第二超声波信号,确定像素点对应的至少一个传输时延;根据像素点的至少一个传输时延检索存储的像素点对应的数据;根据像素点对应的数据合成像素点对应的像素值;对像素点对应的像素值进行变换,确定像素值的解析信号,并在获取解析信号的包络后成像。基于本发明的成像方法,提高了成像并行化程度,加快了成像效率,同时显著改善了阵元数和像素数增加对效率的影响。

The invention relates to an FPGA-based ultrasonic full-data focused fast imaging method and system. The FPGA-based ultrasonic full-data focused fast imaging method includes: collecting a first ultrasonic signal emitted by each array element of at least one first array element to the pixel point, and the pixel point reflects at least one second ultrasonic signal to at least one second array element, and stores the collected data; according to the first ultrasonic signal and at least one second ultrasonic signal, determine at least one a transmission delay; retrieve the data corresponding to the stored pixel point according to at least one transmission delay of the pixel point; synthesize the pixel value corresponding to the pixel point according to the data corresponding to the pixel point; transform the pixel value corresponding to the pixel point to determine the pixel value , and imaged after acquiring the envelope of the resolved signal. Based on the imaging method of the present invention, the degree of parallelization of imaging is improved, the imaging efficiency is accelerated, and the influence of the increase of the number of array elements and the number of pixels on the efficiency is significantly improved.

Description

一种基于FPGA的超声全数据聚焦快速成像方法及系统A method and system for fast imaging of ultrasound full data focusing based on FPGA

技术领域technical field

本发明涉及超声检测领域,特别涉及一种基于FPGA的超声全数据聚焦快速成像方法及系统。The invention relates to the field of ultrasonic detection, in particular to an FPGA-based method and system for fast imaging of ultrasonic full data focusing.

背景技术Background technique

超声全数据聚焦(TFM)成像是一种基于全矩阵捕捉(FMC)的高级相控阵检测方法,相比于常规相控阵成像方法,具有高分辨率、检测范围大等优势,具有着广阔的应用前景。Ultrasound full data focus (TFM) imaging is an advanced phased array detection method based on full matrix capture (FMC). application prospects.

TFM是基于FMC的一种成像方法,如图1所示,TFM对检测区域网格化,根据阵元与像素的空间距离计算传输延时,检索相应的FMC数据,完成像素点的叠加合成。TFM is an imaging method based on FMC. As shown in Figure 1, TFM grids the detection area, calculates the transmission delay according to the spatial distance between the array element and the pixel, retrieves the corresponding FMC data, and completes the superposition and synthesis of pixel points.

任一像素点P(x,z)的值I由公式1给出。The value I of any pixel point P(x, z) is given by Equation 1.

hij为激发阵元为i时,j阵元的接收数据,并经过了Hilbert变换。xi和xj分别是阵元i和j的横轴坐标,c是声速,N是阵元个数。 h ij is the received data of the j array element when the excitation array element is i, and has undergone the Hilbert transform. x i and x j are the abscissa coordinates of array elements i and j, respectively, c is the speed of sound, and N is the number of array elements.

长期以来,全数据聚焦只能在实验室中非实时实现,大量的数据处理和计算成为其走向应用的瓶颈。并行运算平台以其高效的并发执行能力,成为提高全聚焦成像效率的重要方法。For a long time, full data focus can only be realized in non-real time in the laboratory, and a large amount of data processing and calculation has become the bottleneck of its application. With its efficient concurrent execution capability, the parallel computing platform has become an important method to improve the efficiency of all-focus imaging.

并行运算平台主要有GPU,多核CPU和FPGA等。基于GPU+CPU的TFM算法加速是近年来主要的实现方式,但没有考虑硬件采集和数据传输带来的影响,同时限于软件处理器的并行能力,TFM计算效率仍然不足,在阵元或像素增多时,成像效率会急剧下降。Parallel computing platforms mainly include GPU, multi-core CPU and FPGA. The acceleration of the TFM algorithm based on GPU+CPU is the main implementation method in recent years, but it does not consider the impact of hardware acquisition and data transmission, and is limited by the parallel capability of the software processor, the TFM calculation efficiency is still insufficient, and the number of array elements or pixels increases , the imaging efficiency will drop sharply.

FPGA为TFM的高效计算提供了良好并行平台。2013年有研究提出基于多片FPGA的TFM成像,在16阵元、60*60像素量下达到73Hz成像速率,但像素提高到128*128时,效率大幅下降至20Hz以下。原因就在于当前的FPGA设计方法并行化程度不足。如聚焦延时采用软件计算和下传存储的方式,致使FPGA需要花费多个周期才能读取延时值、计算出一个像素,限制了计算效率;同时像素点只能串行合成,在像素分辨率较高时,成像速率明显下降,难以满足以机械扫查为代表的TFM快速检测要求。FPGA provides a good parallel platform for efficient computation of TFM. In 2013, a study proposed that TFM imaging based on multiple FPGAs can achieve a 73Hz imaging rate with 16 elements and 60*60 pixels, but when the pixels are increased to 128*128, the efficiency drops significantly below 20Hz. The reason is that the current FPGA design method is insufficiently parallelized. For example, the focusing delay is calculated by software and downloaded and stored, so that it takes multiple cycles for the FPGA to read the delay value and calculate a pixel, which limits the calculation efficiency; When the imaging rate is high, the imaging rate decreases significantly, and it is difficult to meet the requirements of TFM rapid detection represented by mechanical scanning.

现有技术TFM成像,由于计算方法的并行化不足而导致的成像速率低下、受阵元和像素数量影响明显等问题。In the prior art TFM imaging, due to insufficient parallelization of the calculation method, the imaging rate is low, and it is significantly affected by the number of array elements and pixels.

发明内容SUMMARY OF THE INVENTION

本发明的目的是针对现有技术的缺陷,提出了一种基于FPGA的超声全数据聚焦快速成像方法及系统,相比现有技术,有效提高了TFM成像效率,并显著改善阵元数和像素数增加对成像效率的影响,为以机械扫查为代表的高速TFM检测提供快速成像方案。The purpose of the present invention is to address the defects of the prior art, and propose an FPGA-based ultrasound full data focusing fast imaging method and system. Compared with the prior art, the TFM imaging efficiency is effectively improved, and the number of array elements and pixels are significantly improved. The influence of the increase of the number on the imaging efficiency provides a fast imaging solution for the high-speed TFM detection represented by mechanical scanning.

为实现上述目的,一方面,本发明提供了一种基于FPGA的超声全数据聚焦快速成像方法,该基于FPGA的超声全数据聚焦快速成像具体包括:采集由至少一个第一阵元的每一个阵元发射的第一超声波信号至像素点,及由像素点反射至少一个第二超声波信号至至少一个第二阵元,并将采集后的数据存储;根据第一超声波信号和至少一个第二超声波信号,确定像素点对应的至少一个传输时延;根据像素点的至少一个传输时延检索存储的像素点对应的数据;根据像素点对应的数据合成像素点对应的像素值;对像素点对应的像素值进行变换,确定像素值的解析信号,并在获取解析信号的包络后成像。In order to achieve the above object, on the one hand, the present invention provides an FPGA-based ultrasound full-data focused fast imaging method. The FPGA-based ultrasound full-data focused fast imaging specifically includes: collecting each array consisting of at least one first array element. The first ultrasonic signal transmitted by the element is sent to the pixel point, and at least one second ultrasonic signal is reflected by the pixel point to at least one second array element, and the collected data is stored; according to the first ultrasonic signal and at least one second ultrasonic signal , determine at least one transmission delay corresponding to the pixel point; retrieve the data corresponding to the stored pixel point according to at least one transmission delay of the pixel point; synthesize the pixel value corresponding to the pixel point according to the data corresponding to the pixel point; The values are transformed, the analytic signal of the pixel value is determined, and the envelope of the analytic signal is acquired and imaged.

优选地,并在采集后存储的步骤具体包括:将像素点对应的第一超声波信号和至少一个第二超声波信号存储在双端口内存。Preferably, the step of storing after collection specifically includes: storing the first ultrasonic signal and at least one second ultrasonic signal corresponding to the pixel in a dual-port memory.

优选地,方法还包括:当采集和检索并发执行时,由配置有第一端口和第二端口的双端口内存的第一端口作为数据输入的端口,配置有第一端口和第二端口的双端口内存的第二端口作为数据的检索输出;当采集完毕后,检索执行时,由配置有第一端口和第二端口的双端口内存的第一端口和配置有第一端口和第二端口的双端口内存的第二端口分别作为检索输出。Preferably, the method further includes: when the acquisition and retrieval are performed concurrently, the first port of the dual-port memory configured with the first port and the second port is used as the port for data input, and the dual-port memory configured with the first port and the second port is used as the port for data input. The second port of the port memory is used as the retrieval output of the data; after the collection is completed, when the retrieval is executed, the first port of the dual-port memory configured with the first port and the second port and the first port configured with the first port and the second port are configured. The second port of the dual-port memory is respectively used as the retrieval output.

优选地,根据第一超声波信号和至少一个第二超声波信号,确定像素点对应的至少一个传输时延的步骤具体包括:通过实时延时算法计算像素点对应的至少一个第一阵元的每一个阵元与至少一个第二阵元之间的声程;根据声程和声速,确定至少一个传输时延。Preferably, according to the first ultrasonic signal and the at least one second ultrasonic signal, the step of determining at least one transmission delay corresponding to a pixel specifically includes: calculating each of the at least one first array element corresponding to a pixel by a real-time delay algorithm The sound path between the array element and at least one second array element; at least one transmission delay is determined according to the sound path and the sound speed.

优选地,通过实时延时算法计算像素点对应的至少一个第一阵元的每一个阵元与至少一个第二阵元之间的声程的步骤为:Preferably, the steps of calculating the sound path between each of the at least one first array element and the at least one second array element corresponding to the pixel point by the real-time delay algorithm are:

其中,至少一个第一阵元的第i个阵元的横坐标xi,至少一个第二阵元的第j个阵元的横坐标xj,像素点坐标为(x,z)。The abscissa x i of the i-th array element of at least one first array element, the abscissa x j of the j-th array element of at least one second array element, and the pixel coordinates are (x, z).

优选地,方法还包括:通过现场可编程逻辑门阵列芯片对第一超声波信号和至少一个第二超声波信号并发计算,得到像素点对应的至少一个传输时延。Preferably, the method further includes: concurrently calculating the first ultrasonic signal and the at least one second ultrasonic signal through a field programmable logic gate array chip to obtain at least one transmission delay corresponding to a pixel point.

优选地,对像素进行变换,得到像素的解析信号的步骤具体包括:根据希尔伯特变换对像素值变换,将像素值由时域信号变为频域信号;根据时域信号的像素值和频域信号的像素值合成解析信号。Preferably, the step of transforming the pixel to obtain the analysis signal of the pixel specifically includes: transforming the pixel value according to the Hilbert transform, and changing the pixel value from a time-domain signal to a frequency-domain signal; The pixel values of the frequency domain signal synthesize the analytical signal.

优选地,并在获取解析信号的包络后成像的步骤具体包括:将解析信号的包络上传至上位机进行成像。Preferably, the step of imaging after acquiring the envelope of the analytical signal specifically includes: uploading the envelope of the analytical signal to a host computer for imaging.

优选地,希尔伯特变换还可以为:FIR滤波器。Preferably, the Hilbert transform can also be: an FIR filter.

另一方面,本发明实施例提供一种基于FPGA的超声全数据聚焦快速成像系统,该基于FPGA的超声全数据聚焦快速成像系统包括:采集模块、存储模块、计算模块、检索模块和变换模块;采集模块用于,采集由至少一个第一阵元的每一个阵元发射的第一超声波信号至像素点,及由像素点反射至少一个第二超声波信号至至少一个第二阵元;存储模块用于,将采集后的数据存储;计算模块用于,根据第一超声波信号和至少一个第二超声波信号,确定像素点对应的至少一个传输时延;检索模块用于,根据像素点的至少一个传输时延检索存储的像素点对应的数据;合成模块用于,根据像素点对应的数据合成像素点对应的像素值;变换模块用于,对像素点对应的像素值进行变换,确定像素值的解析信号。On the other hand, an embodiment of the present invention provides an FPGA-based ultrasound full-data focused fast imaging system, where the FPGA-based ultrasound full-data focused fast imaging system includes: an acquisition module, a storage module, a calculation module, a retrieval module, and a transformation module; The acquisition module is used to collect the first ultrasonic signal emitted by each array element of the at least one first array element to the pixel point, and to reflect at least one second ultrasonic signal from the pixel point to the at least one second array element; the storage module is used for is used to store the collected data; the calculation module is used to, according to the first ultrasonic signal and at least one second ultrasonic signal, determine at least one transmission delay corresponding to the pixel point; the retrieval module is used to, according to the at least one transmission time delay of the pixel point The data corresponding to the stored pixel points is retrieved by time delay; the synthesis module is used to synthesize the pixel value corresponding to the pixel point according to the data corresponding to the pixel point; the transformation module is used to transform the pixel value corresponding to the pixel point to determine the analysis of the pixel value Signal.

本发明采用FPGA片内计算,无需CPU参与,有效提高了TFM成像的并行化程度,加快了成像效率,同时显著改善了阵元数和像素数增加对效率的影响。The invention adopts FPGA on-chip calculation without CPU participation, effectively improves the parallelization degree of TFM imaging, accelerates the imaging efficiency, and at the same time significantly improves the influence of the increase of the number of array elements and the number of pixels on the efficiency.

附图说明Description of drawings

图1为现有技术中的超声全数据聚焦成像图;1 is an ultrasound full data focused imaging diagram in the prior art;

图2为本发明的实施例提供的一种基于FPGA的超声全数据聚焦快速成像的方法结构流程图;2 is a structural flowchart of a method for FPGA-based ultrasound full data focused fast imaging provided by an embodiment of the present invention;

图3为极坐标转换图;Figure 3 is a polar coordinate conversion diagram;

图4为FIR滤波器的频率、幅度、相位图;Fig. 4 is the frequency, amplitude, phase diagram of FIR filter;

图5为本发明实施例提供的一种基于FPGA的超声全数据聚焦成快速像系统的结构示意图;5 is a schematic structural diagram of an FPGA-based ultrasound full data focusing system for fast imaging provided by an embodiment of the present invention;

图6为基于图1中基于FPGA的超声全数据聚焦快速成像方法的一个具体的成像图。FIG. 6 is a specific imaging diagram based on the FPGA-based ultrasound full data focused fast imaging method in FIG. 1 .

具体实施方式Detailed ways

下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。The technical solutions of the present invention will be further described in detail below through the accompanying drawings and embodiments.

本发明设计了FPGA片内全并行计算方案,包括基于CORDIC(CoordinateRotational Digital Computer)的实时延时计算、并行像素合成以及像素合成后的Hilbert变换。分别实现聚焦延时的实时并行计算,多像素点的并发合成以及像素解析信号的获取。The invention designs a fully parallel computing scheme in the FPGA chip, including real-time delay calculation based on CORDIC (Coordinate Rotational Digital Computer), parallel pixel synthesis and Hilbert transform after pixel synthesis. Real-time parallel calculation of focus delay, concurrent synthesis of multi-pixel points, and acquisition of pixel resolution signals are realized respectively.

图2为本发明的实施例提供的一种基于FPGA的超声全数据聚焦快速成像方法的流程图。如图2所示,超声全数据聚焦快速成像的步骤具体包括:FIG. 2 is a flowchart of an FPGA-based ultrasound full data focused fast imaging method provided by an embodiment of the present invention. As shown in Figure 2, the steps of the ultrasound full data focusing and fast imaging specifically include:

步骤S100:采集由至少一个第一阵元的每一个阵元发射的第一超声波信号至像素点,及由像素点反射至少一个第二超声波信号至至少一个第二阵元,并在采集后存储;Step S100: Collect the first ultrasonic signal emitted by each of the at least one first array element to the pixel point, and reflect at least one second ultrasonic signal from the pixel point to the at least one second array element, and store after collection. ;

以具有N个阵元激发-接收组合为例,每次FMC采集仅对超声相控阵探头中的一个第一阵元激发发射第一超声波信号,N个第二阵元作为接收单元,对于N通道相控阵,一次激发就可以得到N个第二超声波信号,第二超声波都是经由N个第二阵元接收的回波信号。第一阵元是被激发的发射阵元用于发射第一超声波信号,第二阵元则是接收第二超声波信号的接收阵元。Taking the excitation-reception combination of N array elements as an example, each FMC acquisition only excites and transmits the first ultrasonic signal to one first array element in the ultrasonic phased array probe, and N second array elements are used as receiving units. In the channel phased array, N second ultrasonic signals can be obtained by one excitation, and the second ultrasonic waves are all echo signals received by the N second array elements. The first array element is the excited transmitting array element for transmitting the first ultrasonic signal, and the second array element is the receiving array element for receiving the second ultrasonic signal.

如16阵元系统为例,其过程为:阵元1发射,阵元1-16均采集发射波和采集超声反射回波;采集完后,阵元2发射,阵元1-16均采集发射波和超声反射回波;以此类推,直到阵元16发射,阵元1-16均采集超声反射回波。Take the 16-element system as an example, the process is: array element 1 transmits, and array elements 1-16 all collect the transmitted wave and ultrasonic reflected echo; after the acquisition, array element 2 transmits, and array elements 1-16 all collect and transmit Wave and ultrasonic reflection echoes; and so on, until the array element 16 transmits, the ultrasonic reflection echoes are collected by the array elements 1-16.

由于成像的清晰度不同像素点的个数也不同。在FMC采集完成后,将像素点对应的第一超声波信号和N个第二超声波信号存储双端口内存中。Due to the different imaging resolution, the number of pixels is also different. After the FMC acquisition is completed, the first ultrasonic signal and N second ultrasonic signals corresponding to the pixel points are stored in the dual-port memory.

双端口内存(True Dual-port RAM)是FPGA的片内资源,可满足多种数据读写要求。相比于片外缓存方式,利用片内资源存储可提高通道的并行化读写程度。在片内存储资源中,双端口内存相比于单端口RAM和FIFO,可实现双倍的数据读取带宽。Dual-port memory (True Dual-port RAM) is the on-chip resource of FPGA, which can meet various data read and write requirements. Compared with the off-chip cache method, the use of on-chip resource storage can improve the parallel read and write degree of the channel. Among the on-chip storage resources, dual-port memory can achieve double the data read bandwidth compared to single-port RAM and FIFO.

步骤S110:根据第一超声波信号和至少一个第二超声波信号,确定像素点对应的至少一个传输时延;Step S110: Determine at least one transmission delay corresponding to the pixel point according to the first ultrasonic signal and at least one second ultrasonic signal;

以具有N个阵元激发-接收组合为例,N个接收单元分别接收N个第二超声波信号,每一个接收阵元只接收一个第二超声波信号。Taking an excitation-reception combination with N array elements as an example, the N receiving units receive N second ultrasonic signals respectively, and each receiving array element receives only one second ultrasonic signal.

时延根据超声相控阵探头中的一个阵元激发发射第一超声波信号至像素点,及由像素点反射第二超声波信号至N个接收阵元中的每一个阵元的声程与声速的比值得到的。The time delay is based on the difference between the sound path and the sound speed of each of the N receiving array elements according to the excitation of one array element in the ultrasonic phased array probe to transmit the first ultrasonic signal to the pixel point, and the second ultrasonic signal reflected from the pixel point to each of the N receiving array elements. ratio obtained.

具体地,根据第一超声波信号和至少一个第二超声波信号,确定像素点对应的至少一个传输时延的步骤具体包括:Specifically, according to the first ultrasonic signal and at least one second ultrasonic signal, the step of determining at least one transmission delay corresponding to a pixel point specifically includes:

根据实时延时算法(CORDIC算法)计算超声相控阵探头中的一个阵元激发发射第一超声波信号至像素点,及由像素点反射第二超声波信号至N个接收阵元中的每一个阵元的声程;由声程与声速的比值,确定传输时延。According to the real-time delay algorithm (CORDIC algorithm), one array element in the ultrasonic phased array probe is excited to transmit the first ultrasonic signal to the pixel point, and the second ultrasonic signal is reflected from the pixel point to each of the N receiving array elements. The sound path of the element; the transmission delay is determined by the ratio of the sound path to the sound speed.

CORDIC算法是通过一系列固定的、与运算基数有关的角度,该角度不断偏摆迭代以逼近所需的旋转角度,适合坐标转换和三角函数等的计算。The CORDIC algorithm uses a series of fixed angles related to the operation base, and the angle is continuously swayed and iterated to approximate the required rotation angle, which is suitable for the calculation of coordinate transformation and trigonometric functions.

CORDIC算法应用于坐标转换可以将任意输入向量(X,Y)沿单位圆转动,直到该向量的Y坐标为0,得到旋转角度θ和向量的模X’,从而实现极坐标到直角坐标的转换。The CORDIC algorithm applied to coordinate transformation can rotate any input vector (X, Y) along the unit circle until the Y coordinate of the vector is 0, obtain the rotation angle θ and the modulo X' of the vector, so as to realize the conversion from polar coordinates to rectangular coordinates .

TFM的传输延时由声程/声速得出。像素与阵元间的声程可以在CORDIC的坐标转换中求得,不需要直接开方运算,极大地方便了声程计算。The propagation delay of TFM is derived from the sound path/sound velocity. The sound path between the pixel and the array element can be obtained in the CORDIC coordinate conversion, without direct square root operation, which greatly facilitates the sound path calculation.

如,阵元i与像素点P(x,z),CORDIC的输入坐标X为阵元i与像素点P(x,z)的横坐标之差,纵坐标Y为阵元i与像素点P(x,z)的Z轴坐标差,坐标旋转后的幅度结果就是像素点P(x,z)与该阵元i间的单向声程。该过程的数学表达为:For example, the array element i and the pixel point P(x, z), the input coordinate X of the CORDIC is the difference between the abscissa element i and the pixel point P(x, z), and the ordinate Y is the array element i and the pixel point P The Z-axis coordinate difference of (x, z), the amplitude result after the coordinate rotation is the one-way sound path between the pixel point P(x, z) and the array element i. The mathematical expression of this process is:

由阵元i激发-阵元j接收组合与像素P(x,z)之间的双向声程为:The bidirectional sound path between the combination of excitation by array element i and reception of array element j and pixel P(x, z) is:

一个传输时延为:A transmission delay is:

其中,T为阵元i和阵元j之间的传输时延,c为声速。Among them, T is the transmission delay between array element i and array element j, and c is the speed of sound.

即阵元i激发-N个阵元接收组合与声速的比值,即TFM中每一个像素点所需的N个传输时延。That is, the ratio between the excitation of array element i and the reception combination of N array elements and the speed of sound, that is, the N transmission delays required for each pixel in the TFM.

可选地,本发明实施例基于FPGA采用N个CORDIC计算模块,可实现N个传输时延的并发计算,满足像素点对应的像素值的快速计算要求。相比于常规的软件计算方式,本发明实施例完全由FPGA内部计算,不需要任何存储和读取,可真正实现实时并行计算,有效提高了TFM像素计算效率。Optionally, the embodiment of the present invention adopts N CORDIC calculation modules based on FPGA, which can realize the concurrent calculation of N transmission delays, and meet the requirements for fast calculation of pixel values corresponding to pixel points. Compared with the conventional software calculation method, the embodiment of the present invention is completely calculated by the FPGA without any storage and reading, which can truly realize real-time parallel calculation, and effectively improve the TFM pixel calculation efficiency.

步骤S120:根据像素点的至少一个传输时延检索存储的像素点对应的数据;Step S120: retrieving data corresponding to the stored pixel point according to at least one transmission delay of the pixel point;

步骤S130:根据像素点对应的数据合成像素点对应的像素值。Step S130 : synthesizing pixel values corresponding to the pixel points according to the data corresponding to the pixel points.

由N个阵元激发-接收组合为例,每一个像素点都会对应N个传输时延,将N个传输时延检索出存储在双端口内存中的数据,叠加检索出的数据该像素点的像素值。Taking the excitation-reception combination of N array elements as an example, each pixel corresponds to N transmission delays. The N transmission delays are retrieved to retrieve the data stored in the dual-port memory, and the retrieved data is superimposed on the pixel of the pixel. Pixel values.

由于TFM成像会涉及到很多的像素值,一般会先将像素值存储。在需要进行对像素值进行成像时,直接从存储的位置获取像素值。Since TFM imaging involves a lot of pixel values, the pixel values are generally stored first. When it is necessary to image the pixel value, the pixel value is directly obtained from the stored location.

步骤S140:对像素点对应的像素值进行变换,确定像素值的解析信号。Step S140: Transform the pixel value corresponding to the pixel point to determine the analysis signal of the pixel value.

根据Hilbert变换对像素值进行变换,将时域信号的像素值转变为频域信号的像素值;将时域信号的像素值与频域信号的像素值组合成解析信号。The pixel value is transformed according to the Hilbert transform, and the pixel value of the time domain signal is converted into the pixel value of the frequency domain signal; the pixel value of the time domain signal and the pixel value of the frequency domain signal are combined into an analytical signal.

对像素值阵列进行Hilbert变换,获取像素值的解析信号,保证成像效果的同时,符合FPGA的流水线设计,并减少了资源占用。Hilbert transform is performed on the pixel value array to obtain the analytic signal of the pixel value, which ensures the imaging effect, conforms to the FPGA pipeline design, and reduces resource occupation.

在通信领域中,Hilbert变换是对信号取包络的有效方式。在基于FMC的后处理成像中,通常先以Hilbert变换得到超声回波的解析信号,再获取解析信号的包络信号。相比于原始信号直接成像,解析信号的包络使成像更加平滑。In the field of communication, the Hilbert transform is an effective way to obtain the envelope of the signal. In post-processing imaging based on FMC, the analytical signal of ultrasonic echo is usually obtained by Hilbert transform first, and then the envelope signal of the analytical signal is obtained. Compared to direct imaging of the original signal, the envelope of the resolved signal results in smoother imaging.

在常规TFM成像中,Hilbert变换主要由软件实现,其逐点运算给软件实时性提高带来不便。Hilbert变换是线性变换过程,符合叠加原理,因此TFM成像中,信号先叠加得到像素点对应的像素值后再Hilbert变换得到解析信号与先Hilbert变换再叠加在理论上具有同样的效果。然而后者更符合FPGA的资源节约和流水线计算。In conventional TFM imaging, Hilbert transform is mainly implemented by software, and its point-by-point operation brings inconvenience to the real-time improvement of software. Hilbert transform is a linear transformation process, which conforms to the superposition principle. Therefore, in TFM imaging, the signals are first superimposed to obtain the pixel values corresponding to the pixels, and then the Hilbert transform to obtain the analytical signal has the same effect in theory as Hilbert transform and then superposition. However, the latter is more in line with the resource saving and pipeline computing of FPGAs.

步骤S130和步骤S140的执行顺序可以是先执行步骤S140再执行步骤S130。本发明实施例只是列举执行步骤S130后执行步骤S140,先执行步骤S140在此处不在论述,成像的效果是相同的。The execution sequence of step S130 and step S140 may be that step S140 is executed first and then step S130 is executed. This embodiment of the present invention only enumerates that step S130 is performed and then step S140 is performed, and step S140 is performed first, which is not discussed here, and the imaging effects are the same.

Hilbert变换是物理不可实现的,因此在FPGA中采用FIR滤波器来实现Hilbert变换相同的功能。The Hilbert transform is physically unrealizable, so an FIR filter is used in the FPGA to achieve the same function as the Hilbert transform.

本发明实施例了采用128阶FIR滤波器,其频率幅度相位响应如图3所示,其中幅度谱具有全通特性,而相位具有严格的线性,符合Hilbert变换要求。The embodiment of the present invention adopts a 128-order FIR filter, and its frequency, amplitude and phase response are shown in Figure 3, wherein the amplitude spectrum has an all-pass characteristic, and the phase has strict linearity, which meets the requirements of Hilbert transform.

步骤S150:获取解析信号的包络上传至上位机成像。Step S150 : acquiring the envelope of the analytical signal and uploading it to the host computer for imaging.

获取解析信号的包络;将包络信号上传至上位机显示图像,从而完成TFM成像。Obtain the envelope of the analytical signal; upload the envelope signal to the host computer to display the image, thus completing the TFM imaging.

由于双端口内存配置有两个可配置的第一端口和第二端口,超声全数据聚焦成像方法还包括:Since the dual-port memory configuration has two configurable first ports and second ports, the ultrasound full data focused imaging method further includes:

步骤S160:当FMC采集与检索并发执行时,由配置有第一端口和第二端口的双端口内存的第一端口作为数据输入的端口,配置有第一端口和第二端口的双端口内存的第二端口作为数据的检索输出;Step S160: When the FMC acquisition and retrieval are performed concurrently, the first port of the dual-port memory configured with the first port and the second port is used as the port for data input, and the dual-port memory configured with the first port and the second port is used. The second port is used as the retrieval output of data;

当采集完毕后,检索执行时,由配置有第一端口和第二端口的双端口内存的第一端口作为检索输出,由配置有第一端口和第二端口的双端口内存的第二端口作为检索输出。After the collection is completed, when the retrieval is executed, the first port of the dual-port memory configured with the first port and the second port is used as the retrieval output, and the second port of the dual-port memory configured with the first port and the second port is used as the retrieval output. Retrieve output.

具体地,双端口内存配置有可配置的第一端口和第二端口,第一端口为端口A,第二端口为端口B。Specifically, the dual-port memory is configured with a configurable first port and a second port, where the first port is port A and the second port is port B.

当FMC减速和检索并发执行时,端口A作为数据输入的端口,端口B作为数据检索输出;当FMC采集完毕,只执行检索时,端口A和端口B作为检索输出,端口A和端口B可以是相互独立的,可以同时支持2个像素点2N个传输时延的计算,进而由端口A和端口B分别检索出像素点对应的数据,合成2个像素点对应的像素值。When FMC deceleration and retrieval are executed concurrently, port A is used as data input port, and port B is used as data retrieval output; when FMC collection is completed and only retrieval is performed, port A and port B are used as retrieval output, and port A and port B can be Independent of each other, it can support the calculation of 2N transmission delays for 2 pixels at the same time, and then retrieve the data corresponding to the pixels from port A and port B respectively, and synthesize the pixel values corresponding to the two pixels.

针对一个像素点时,基于FPGA,利用N个CORDIC计算模块,可实现N个传输时延的计算;针对2个像素点,将N个CORDIC算法模块组成成一个延时模块,每一个像素点根据延时模块得到对应的N个传输延时,进而能得到2个像素点2N个传输延时的计算,进而合成2个像素点分别对应的像素值。For one pixel, based on FPGA, N CORDIC calculation modules can be used to calculate N transmission delays; for 2 pixels, N CORDIC algorithm modules are formed into a delay module, and each pixel is based on The delay module obtains the corresponding N transmission delays, and then can obtain the calculation of the 2N transmission delays of the two pixel points, and then synthesize the pixel values corresponding to the two pixel points respectively.

相比于片外数据缓存和片内普通缓存方式,本发明实施例可实现2个像素的并发计算,将像素计算效率提高了一倍,加快了成像速率。随着FPGA中资源的提升,在阵元数N增加时可提高双端口内存的例化,保证每个通道的FMC采集,降低N增加对TFM效率的影响。Compared with the off-chip data cache and the on-chip common cache, the embodiment of the present invention can realize the concurrent calculation of 2 pixels, double the pixel calculation efficiency, and speed up the imaging rate. With the improvement of resources in the FPGA, when the number of array elements N increases, the instantiation of dual-port memory can be improved to ensure the FMC acquisition of each channel and reduce the impact of increasing N on TFM efficiency.

需要说明的是,N个阵元的一维线性阵列探头进行FMC采集,顺序激发每个阵元,每一次激发之后所有阵元同时且各自独立接收数据。顺序激发每个阵元,可以采用N个延时时序,延时时序与阵元对应,实现阵元的顺序激发;也可以通过软件实现N个定时,到达第一时间时,激发第一个阵元至像素点,以此类推。It should be noted that the one-dimensional linear array probe of N array elements performs FMC acquisition, and sequentially excites each array element. After each excitation, all array elements simultaneously and independently receive data. To activate each array element sequentially, N delay sequences can be used, and the delay sequence corresponds to the array elements to realize the sequential excitation of the array elements; it is also possible to realize N timings through software, and when the first time is reached, the first array is excited. Elements to pixels, and so on.

本发明实施例提出了基于FPGA的全并行TFM计算架构,包括基于CORDIC(Coordinate Rotational Digital Computer)的实时延时计算、并行像素合成以及像素合成后的Hilbert变换。本发明可实现聚焦延时的实时并发计算,无需任何存储,同时可实现多个像素点的并发合成。相比于常规方法,本发明更好地发挥了FPGA在并行计算和信号处理中的优势,有效提高了TFM成像速率,可显著改善阵元数和像素数增加对效率的影响。本发明可为以机械扫查为代表的高速TFM检测提供快速成像方案。The embodiment of the present invention proposes an FPGA-based fully parallel TFM computing architecture, including CORDIC (Coordinate Rotational Digital Computer)-based real-time delay computation, parallel pixel synthesis, and Hilbert transform after pixel synthesis. The invention can realize the real-time concurrent calculation of the focusing delay without any storage, and can realize the concurrent synthesis of multiple pixel points at the same time. Compared with the conventional method, the present invention better utilizes the advantages of FPGA in parallel computing and signal processing, effectively improves the TFM imaging rate, and can significantly improve the influence of the increase in the number of array elements and pixels on the efficiency. The present invention can provide a fast imaging solution for high-speed TFM detection represented by mechanical scanning.

图3为极坐标转换图。如图3所示,CORDIC算法应用于极坐标转换将任意输入向量(Input Vector)沿单位圆转动,直到输入向量的为零,得到旋转角度(Outputphase)和输入向量的模(Output Mag),从而实现极坐标到直角坐标的转换。Figure 3 is a polar coordinate conversion diagram. As shown in Figure 3, the CORDIC algorithm is applied to polar coordinate transformation to convert any input vector (Input Vector) rotate along the unit circle until the If it is zero, the rotation angle (Outputphase) and the modulus of the input vector (Output Mag) are obtained, so as to realize the conversion from polar coordinates to rectangular coordinates.

像素点与至少一个第一阵元的每一个阵元与至少一个第二阵元之间的声程可以在CORDIC的极坐标到直角坐标的转换中求得。The sound path between the pixel point and each of the at least one first array element and the at least one second array element can be obtained in the conversion from polar coordinates to rectangular coordinates of CORDIC.

如,阵元i与像素点P(x,z),CORDIC的输入坐标X为阵元i与像素点P(x,z)的横坐标x之差,纵坐标Y为阵元i与像素点P(x,z)的纵坐标z之差,坐标旋转后的幅度结果就是像素点P(x,z)与该阵元i间的单向声程。该过程的数学表达为:For example, the array element i and the pixel point P(x, z), the input coordinate X of the CORDIC is the difference between the abscissa x of the array element i and the pixel point P(x, z), and the ordinate Y is the array element i and the pixel point. The difference between the vertical coordinates z of P(x, z) and the amplitude result after the coordinate rotation is the one-way sound path between the pixel point P(x, z) and the array element i. The mathematical expression of this process is:

其中,横坐标X为至少一个第一阵元的一个阵元的横坐标xi与像素点横坐标x之差,纵坐标Y为至少一个第一阵元的一个阵元的纵坐标zi与像素点的轴坐标z之差。Wherein, the abscissa X is the difference between the abscissa x i of one array element of at least one first array element and the abscissa x of the pixel point, and the ordinate Y is the ordinate zi and z i of one array element of at least one first array element The difference between the axis coordinates z of the pixel points.

与CORDIC算法具有相同功能都可以等效替代。The same function as the CORDIC algorithm can be replaced by equivalent.

图4为FIR滤波器的频率、幅度、相位图。如图4所示,频率(Frequency),单位为MHz;幅度(Magnitude),单位是dB(分贝);相位(Phase),单位为radians(弧度)。Figure 4 shows the frequency, amplitude and phase diagrams of the FIR filter. As shown in Figure 4, the frequency (Frequency), the unit is MHz; the amplitude (Magnitude), the unit is dB (decibel); the phase (Phase), the unit is radians (radians).

FIR滤波器为128阶的滤波器,用于将时域信号的像素值转变为频域信号的像素值。The FIR filter is a 128-order filter, which is used to convert the pixel value of the time domain signal into the pixel value of the frequency domain signal.

如图4所示,FIR滤波器的幅度谱具有全通特性,FIR滤波器的相位具有严格的线性,符合Hilbert变换要求。As shown in Figure 4, the amplitude spectrum of the FIR filter has an all-pass characteristic, and the phase of the FIR filter has strict linearity, which meets the requirements of Hilbert transform.

与Hilber变换具有相同功能都可以等效替代。It can be equivalently replaced with the same function as Hilber transform.

图5为本发明实施例提供的一种基于FPGA的超声全数据聚焦快速成像系统的结构示意图。如图5所示,超声全数据聚焦快速成像系统包括:FIG. 5 is a schematic structural diagram of an FPGA-based ultrasound full data focusing fast imaging system provided by an embodiment of the present invention. As shown in Figure 5, the ultrasound full data focusing rapid imaging system includes:

硬件部分是基于FPGA,基于FPGA的超声全数据聚焦快速成像系统包括FMC采集模块10、存储模块20、计算模块30、检索模块40、合成模块50和变换模块60。The hardware part is based on FPGA. The FPGA-based ultrasound full data focusing fast imaging system includes an FMC acquisition module 10 , a storage module 20 , a calculation module 30 , a retrieval module 40 , a synthesis module 50 and a transformation module 60 .

采集模块10用于采集由至少一个第一阵元的每一个阵元发射的第一超声波信号至像素点,及由像素点反射至少一个第二超声波信号至至少一个第二阵元;存储模块20用于将采集后的数据存储于双端口内存;计算模块30用于,根据第一超声波信号和至少一个第二超声波信号,确定像素点对应的至少一个传输时延;检索模块40用于,根据像素点的至少一个传输时延检索存储的像素点对应的数据;合成模块50用于,根据像素点对应的数据合成像素点对应的像素值;变换模块60用于,对像素点对应的像素值进行变换,确定像素值的解析信号。The acquisition module 10 is used to collect the first ultrasonic signal emitted by each of the at least one first array element to the pixel point, and to reflect at least one second ultrasonic signal from the pixel point to the at least one second array element; the storage module 20 Used to store the collected data in the dual-port memory; the calculation module 30 is used to, according to the first ultrasonic signal and at least one second ultrasonic signal, determine at least one transmission delay corresponding to the pixel point; the retrieval module 40 is used to, according to At least one transmission delay of the pixel point retrieves the data corresponding to the stored pixel point; the synthesis module 50 is used for synthesizing the pixel value corresponding to the pixel point according to the data corresponding to the pixel point; the transformation module 60 is used for the pixel value corresponding to the pixel point. A transformation is performed to determine the analytic signal of the pixel value.

具体地,存储模块20用于将采集后的数据存储与双端口内存中。Specifically, the storage module 20 is used for storing the collected data in the dual-port memory.

双端口内存配置有可配置的第一端口和第二端口,第一端口为端口A,第二端口为端口B。The dual-port memory is configured with a configurable first port and a second port, where the first port is port A and the second port is port B.

当FMC减速和检索并发执行时,端口A作为数据输入的端口,端口B作为数据检索输出;当FMC采集完毕,只执行检索时,端口A和端口B作为检索输出,端口A和端口B可以是相互独立的,可以同时支持2个像素点2N个传输时延的计算,进而由端口A和端口B分别检索出像素点对应的数据,合成2个像素点对应的像素值。When FMC deceleration and retrieval are executed concurrently, port A is used as data input port, and port B is used as data retrieval output; when FMC collection is completed and only retrieval is performed, port A and port B are used as retrieval output, and port A and port B can be Independent of each other, it can support the calculation of 2N transmission delays for 2 pixels at the same time, and then retrieve the data corresponding to the pixels from port A and port B respectively, and synthesize the pixel values corresponding to the two pixels.

针对一个像素点时,基于FPGA,利用N个CORDIC计算模块,可实现N个传输时延的计算;针对2个像素点,将N个CORDIC计算模块组成成一个延时模块,每一个像素点根据延时模块得到对应的N个传输延时,进而能得到2个像素点2N个传输延时的计算,进而合成2个像素点分别对应的像素值。For one pixel, based on FPGA, N CORDIC calculation modules can be used to calculate N transmission delays; for 2 pixels, N CORDIC calculation modules are formed into a delay module, and each pixel is based on The delay module obtains the corresponding N transmission delays, and then can obtain the calculation of the 2N transmission delays of the two pixel points, and then synthesize the pixel values corresponding to the two pixel points respectively.

计算模块30通过实时延时算法计算像素点对应的至少一个第一阵元的每一个阵元与至少一个第二阵元之间的声程;根据声程和声速,确定至少一个传输时延。The calculation module 30 calculates the sound path between each of the at least one first array element corresponding to the pixel and at least one second array element through a real-time delay algorithm; and determines at least one transmission delay according to the sound path and sound speed.

本发明实施例基于FPGA采用N个CORDIC计算模块,可实现N个传输时延的并发计算,满足像素点对应的像素值的快速计算要求。The embodiment of the present invention adopts N CORDIC calculation modules based on the FPGA, which can realize the concurrent calculation of N transmission delays, and meet the requirements for fast calculation of pixel values corresponding to pixel points.

合成模块50根据希尔伯特变换对像素值变换,将像素值由时域信号变为频域信号;根据时域信号的像素值和频域信号的像素值合成解析信号。The synthesis module 50 transforms the pixel value according to the Hilbert transform, and changes the pixel value from a time-domain signal to a frequency-domain signal; and synthesizes the analytical signal according to the pixel value of the time-domain signal and the pixel value of the frequency-domain signal.

由于TFM成像会涉及到很多的像素值,一般会先将像素值存储。基于FPGA的超声全数据快速成像系统还包括存储模块70,存储模块70用于将像素值存储,可以是普通的缓存或者内存,在此不限制,在附图5中没有标识。Since TFM imaging involves a lot of pixel values, the pixel values are generally stored first. The FPGA-based ultrasound full-data rapid imaging system further includes a storage module 70, which is used for storing pixel values, which may be a common cache or memory, which is not limited here, and is not marked in FIG. 5 .

软件部分,基于FPGA的超声全数据聚焦快速成像系统包括系统控制模块80和图像成像模块90,系统控制模块80用于顺序激发阵元,图像成像模块90用于基于获取解析信号的包络上传至上位机成像。In the software part, the FPGA-based ultrasound full-data focusing fast imaging system includes a system control module 80 and an image imaging module 90. The system control module 80 is used to sequentially excite the array elements, and the image imaging module 90 is used to obtain the envelope of the analytical signal and upload it to the system. PC imaging.

图6为基于图1中基于FPGA的超声全数据聚焦快速成像方法的一个具体的成像图。如图6所示,本发明实施例采用16/64(脉冲重复频率/阵元数),通过TFM实时成像,实现铁轨母材裂纹和缺陷的快速检测。FIG. 6 is a specific imaging diagram based on the FPGA-based ultrasound full data focused fast imaging method in FIG. 1 . As shown in FIG. 6 , the embodiment of the present invention adopts 16/64 (pulse repetition frequency/number of array elements), and realizes the rapid detection of cracks and defects in the rail parent material through TFM real-time imaging.

钢块中1mm横穿孔是常见的人工缺陷,本发明实施例为了验证TFM计算效率和成像效果,对20mm深的1mm横穿孔进行了TFM检测成像,成像速率可达312.5Hz,满足TFM的快速、实时成像检测。The 1mm transverse perforation in the steel block is a common artificial defect. In the embodiment of the present invention, in order to verify the TFM calculation efficiency and imaging effect, TFM detection and imaging are performed on the 1mm transverse perforation with a depth of 20mm, and the imaging rate can reach 312.5Hz, which satisfies the fast, fast, and reliable performance of TFM. Real-time imaging inspection.

本发明使超声全数据成像的效率得到了有效提高,成像速率为PRF/N(脉冲重复频率/阵元数),达到TFM检测的物理极限;在像素分辨率增加时,效率降低的情况被显著改善。The present invention effectively improves the efficiency of ultrasound full data imaging, and the imaging rate is PRF/N (pulse repetition frequency/number of array elements), which reaches the physical limit of TFM detection; when the pixel resolution increases, the efficiency reduction is significantly reduced. improve.

在100*100和200*200像素分辨率下,本发明实施例采用16阵元TFM成像速率可达312.5Hz,效率不随像素点的个数增加而降低;32阵元时,TFM成像速率均可达156.25Hz,在阵元数增加时,符合PRF/N的物理限制。可为以机械扫查为代表的高速TFM检测提供快速成像方案。Under the resolution of 100*100 and 200*200 pixels, the TFM imaging rate of 16 array elements in the embodiment of the present invention can reach 312.5 Hz, and the efficiency does not decrease with the increase of the number of pixels; when the number of pixels is 32, the TFM imaging rate can be Up to 156.25Hz, when the number of array elements increases, it meets the physical limit of PRF/N. It can provide a fast imaging solution for high-speed TFM inspection represented by mechanical scanning.

本发明采用FPGA片内计算,无需CPU参与,有效提高了TFM成像的并行化程度,加快了成像效率,同时显著改善了阵元数和像素数增加对效率的影响。The invention adopts FPGA on-chip calculation without CPU participation, effectively improves the parallelization degree of TFM imaging, accelerates the imaging efficiency, and at the same time significantly improves the influence of the increase of the number of array elements and the number of pixels on the efficiency.

以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The specific embodiments described above further describe the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above descriptions are only specific embodiments of the present invention, and are not intended to limit the scope of the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims (7)

1.一种基于FPGA的超声全数据聚焦快速成像方法,其特征在于,包括:1. a FPGA-based ultrasound full data focusing fast imaging method, is characterized in that, comprising: 采集由至少一个第一阵元的每一个阵元发射的第一超声波信号至像素点,及由所述像素点反射至少一个第二超声波信号至至少一个第二阵元,并将采集后的数据存储;其中,将所述像素点对应的所述第一超声波信号和所述至少一个第二超声波信号存储在双端口内存;当采集和检索并发执行时,由配置有第一端口和第二端口的双端口内存的所述第一端口作为数据输入的端口,配置有第一端口和第二端口的双端口内存的所述第二端口作为数据的检索输出;当采集完毕后,检索执行时,由配置有第一端口和第二端口的双端口内存的所述第一端口和配置有第一端口和第二端口的双端口内存的所述第二端口分别作为检索输出;Collect the first ultrasonic signal emitted by each of the at least one first array element to a pixel point, and reflect at least one second ultrasonic signal from the pixel point to the at least one second array element, and collect the collected data storage; wherein, the first ultrasonic signal and the at least one second ultrasonic signal corresponding to the pixel are stored in a dual-port memory; when the acquisition and retrieval are performed concurrently, the first port and the second port are configured with the first port and the second port. The first port of the dual-port memory is used as the port for data input, and the second port of the dual-port memory configured with the first port and the second port is used as the retrieval output of data; after the collection is completed, when the retrieval is executed, The first port of the dual-port memory configured with the first port and the second port and the second port of the dual-port memory configured with the first port and the second port are respectively used as retrieval outputs; 通过现场可编程逻辑门阵列(FPGA)芯片对所述第一超声波信号和所述至少一个第二超声波信号并发计算,得到所述像素点对应的至少一个传输时延;The first ultrasonic signal and the at least one second ultrasonic signal are concurrently calculated by a field programmable gate array (FPGA) chip to obtain at least one transmission delay corresponding to the pixel point; 根据所述像素点的至少一个传输时延检索存储的所述像素点对应的数据;Retrieve the stored data corresponding to the pixel point according to at least one transmission delay of the pixel point; 根据所述像素点对应的数据合成所述像素点对应的像素值;Synthesize the pixel value corresponding to the pixel point according to the data corresponding to the pixel point; 对所述像素点对应的像素值进行变换,确定所述像素值的解析信号,并在获取所述解析信号的包络后成像。Transforming the pixel value corresponding to the pixel point, determining the analytical signal of the pixel value, and obtaining an image after obtaining the envelope of the analytical signal. 2.根据权利要求1所述的成像方法,其特征在于,所述通过现场可编程逻辑门阵列(FPGA)芯片对所述第一超声波信号和所述至少一个第二超声波信号并发计算,得到所述像素点对应的至少一个传输时延的步骤具体包括:2 . The imaging method according to claim 1 , wherein the first ultrasonic signal and the at least one second ultrasonic signal are concurrently calculated by a Field Programmable Logic Gate Array (FPGA) chip, and the obtained result is obtained. 3 . The step of at least one transmission delay corresponding to the pixel point specifically includes: 通过实时延时算法计算所述像素点对应的所述至少一个第一阵元的每一个阵元与所述至少一个第二阵元之间的声程;Calculate the sound path between each array element of the at least one first array element corresponding to the pixel point and the at least one second array element by using a real-time delay algorithm; 根据所述声程和声速,确定所述至少一个传输时延。The at least one transmission delay is determined based on the sound path and the sound speed. 3.根据权利要求2所述的成像方法,其特征在于,所述通过实时延时算法计算所述像素点对应的所述至少一个第一阵元的每一个阵元与所述至少一个第二阵元之间的声程的步骤为:3 . The imaging method according to claim 2 , wherein the calculation of each array element of the at least one first array element corresponding to the pixel point by a real-time delay algorithm and the at least one second array element are 3 . The steps of the sound path between the elements are: 其中,xi为至少一个第一阵元的第i个阵元的横坐标,xj为至少一个第二阵元的第j个阵元的横坐标,(x,z)为像素点坐标。Wherein, x i is the abscissa of the i-th array element of at least one first array element, x j is the abscissa of the j-th array element of at least one second array element, and (x, z) is the pixel coordinate. 4.根据权利要求1所述的成像方法,其特征在于,所述对所述像素点对应的像素值进行变换,确定所述像素值的解析信号的步骤具体包括:4. The imaging method according to claim 1, wherein the step of transforming the pixel value corresponding to the pixel point and determining the analytical signal of the pixel value specifically comprises: 根据希尔伯特变换对所述像素值变换,将所述像素值由时域信号变为频域信号;Transform the pixel value according to the Hilbert transform, and change the pixel value from a time-domain signal to a frequency-domain signal; 根据所述时域信号的像素值和所述频域信号的像素值合成解析信号。An analytical signal is synthesized from pixel values of the time domain signal and pixel values of the frequency domain signal. 5.根据权利要求4所述的成像方法,其特征在于,所述并在获取所述解析信号的包络后成像的步骤具体包括:5. The imaging method according to claim 4, wherein the step of imaging after acquiring the envelope of the analytical signal specifically comprises: 将所述解析信号的包络上传至上位机进行成像。Upload the envelope of the analysis signal to the upper computer for imaging. 6.根据权利要求4所述的成像方法,其特征在于,所述希尔伯特变换还可以为:FIR滤波器。6 . The imaging method according to claim 4 , wherein the Hilbert transform can also be an FIR filter. 7 . 7.一种基于FPGA的超声全数据聚焦快速成像系统,其特征在于,包括:采集模块(10)、存储模块(20)、计算模块(30)、检索模块(40)、合成模块(50)和变换模块(60);7. An FPGA-based ultrasound full-data focused fast imaging system, characterized in that, comprising: an acquisition module (10), a storage module (20), a calculation module (30), a retrieval module (40), and a synthesis module (50) and a transformation module (60); 所述采集模块(10)用于,采集由至少一个第一阵元的每一个阵元发射的第一超声波信号至像素点,及由所述像素点反射至少一个第二超声波信号至至少一个第二阵元;The acquisition module (10) is used for acquiring a first ultrasonic signal emitted by each of the at least one first array element to a pixel point, and reflecting at least one second ultrasonic signal from the pixel point to at least one first ultrasonic signal. binary element; 所述存储模块(20)用于,将采集后的数据存储;其中,将所述像素点对应的所述第一超声波信号和所述至少一个第二超声波信号存储在双端口内存;当采集和检索并发执行时,由配置有第一端口和第二端口的双端口内存的所述第一端口作为数据输入的端口,配置有第一端口和第二端口的双端口内存的所述第二端口作为数据的检索输出;当采集完毕后,检索执行时,由配置有第一端口和第二端口的双端口内存的所述第一端口和配置有第一端口和第二端口的双端口内存的所述第二端口分别作为检索输出;The storage module (20) is used for storing the collected data; wherein, the first ultrasonic signal and the at least one second ultrasonic signal corresponding to the pixel are stored in a dual-port memory; when collecting and When the retrieval is executed concurrently, the first port of the dual-port memory configured with the first port and the second port is used as the port for data input, and the second port of the dual-port memory configured with the first port and the second port is used. As the retrieval output of the data; when the retrieval is performed after the collection is completed, the first port and the dual-port memory configured with the first port and the second port are configured with the first port and the second port. The second ports are respectively used as retrieval outputs; 所述计算模块(30)用于,通过现场可编程逻辑门阵列(FPGA)芯片对所述第一超声波信号和所述至少一个第二超声波信号并发计算,得到所述像素点对应的至少一个传输时延;The calculation module (30) is used for concurrently calculating the first ultrasonic signal and the at least one second ultrasonic signal through a field programmable logic gate array (FPGA) chip to obtain at least one transmission corresponding to the pixel point delay; 所述检索模块(40)用于,根据所述像素点的至少一个传输时延检索存储的所述像素点对应的数据;The retrieval module (40) is configured to retrieve the stored data corresponding to the pixel point according to at least one transmission delay of the pixel point; 所述合成模块(50)用于,根据所述像素点对应的数据合成所述像素点对应的像素值;The synthesizing module (50) is configured to synthesize the pixel value corresponding to the pixel point according to the data corresponding to the pixel point; 所述变换模块(60)用于,对所述像素点对应的像素值进行变换,确定所述像素值的解析信号,并在获取所述解析信号的包络后成像。The transformation module (60) is used for transforming the pixel value corresponding to the pixel point, determining an analytical signal of the pixel value, and imaging after acquiring the envelope of the analytical signal.
CN201611077700.6A 2016-11-29 2016-11-29 A method and system for fast imaging of ultrasound full data focusing based on FPGA Expired - Fee Related CN106651740B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611077700.6A CN106651740B (en) 2016-11-29 2016-11-29 A method and system for fast imaging of ultrasound full data focusing based on FPGA

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611077700.6A CN106651740B (en) 2016-11-29 2016-11-29 A method and system for fast imaging of ultrasound full data focusing based on FPGA

Publications (2)

Publication Number Publication Date
CN106651740A CN106651740A (en) 2017-05-10
CN106651740B true CN106651740B (en) 2019-09-27

Family

ID=58814108

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611077700.6A Expired - Fee Related CN106651740B (en) 2016-11-29 2016-11-29 A method and system for fast imaging of ultrasound full data focusing based on FPGA

Country Status (1)

Country Link
CN (1) CN106651740B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109975814B (en) * 2017-12-28 2020-09-22 深圳先进技术研究院 Ultrasound imaging method, system and device
CN111812205B (en) * 2020-06-18 2023-05-26 广州多浦乐电子科技股份有限公司 Full-focusing and phased array double-scanning imaging method
CN113156431B (en) * 2021-04-30 2022-04-26 珠海微度芯创科技有限责任公司 FPGA-based back projection rapid imaging architecture design method
CN117405780A (en) * 2023-12-15 2024-01-16 广州多浦乐电子科技股份有限公司 Dual-pixel ultrasonic full focusing system based on FPGA and implementation method
CN117710184B (en) * 2024-01-04 2024-08-23 浙江大学 Ultrasonic full-focusing wave number domain imaging method, system and device based on FPGA and GPU

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101783138A (en) * 2009-01-17 2010-07-21 华为技术有限公司 Image-scaling method and device thereof
CN103908304A (en) * 2014-03-14 2014-07-09 中瑞科技(常州)有限公司 Ultrasonic elastography system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050148872A1 (en) * 2003-12-23 2005-07-07 Richard William D. Apparatus and method for synthetic focus ultrasonic imaging

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101783138A (en) * 2009-01-17 2010-07-21 华为技术有限公司 Image-scaling method and device thereof
CN103908304A (en) * 2014-03-14 2014-07-09 中瑞科技(常州)有限公司 Ultrasonic elastography system

Also Published As

Publication number Publication date
CN106651740A (en) 2017-05-10

Similar Documents

Publication Publication Date Title
CN106651740B (en) A method and system for fast imaging of ultrasound full data focusing based on FPGA
CN106680825B (en) An acoustic array imaging system and method
CN109975814B (en) Ultrasound imaging method, system and device
CN108490077A (en) A kind of real-time 3D total focus phased array supersonic imaging method, device and storage medium
CN102928844B (en) Underwater sub-wavelength resolution ratio three-dimensional imaging method
CN102621223A (en) Ultrasonic scanning and detection method based on sonic time reversal method
CN102508251B (en) A Fast Realization Method of Sector Transformation in Multibeam Image Sonar
CN110531338A (en) Multimode SAR self-focusing immediate processing method and system based on FPGA
CN106596736A (en) Real-time ultrasonic phased array total focusing imaging method
JP2000028589A (en) Three-dimensional ultrasonic imaging device
CN109471108A (en) Three-dimensional imaging method, device and electronic device for precessing vertebral body target
US20200150013A1 (en) Method and device for quantifying viscoelasticity of a medium
JP2006508729A (en) High frame rate 3D ultrasound imager
JP5456367B2 (en) Phased array aperture synthesis processing method
Wang et al. Efficient acceleration for total focusing method based on advanced parallel computing in FPGA
CN114185046A (en) Video SAR imaging system and method based on adjustable transmit waveform parameters
CN105068076B (en) The device of the weather radar image interpolation method based on Fourier spectrum analysis
CN114415181B (en) Method and device for generating original echo of synthetic aperture radar
Njiki et al. Total focusing method for non destructive evaluation: Toward real-time imaging systems
Njiki et al. A multi-FPGA implementation of real-time reconstruction using Total Focusing Method
CN117092218A (en) Ultrasonic phased array surface leakage wave full focusing imaging method based on root mean square speed
Dziewierz et al. A design methodology for 2D sparse NDE arrays using an efficient implementation of refracted-ray TFM
Romero-Laorden et al. Strategies for hardware reduction on the design of portable ultrasound imaging systems
Gyöngy et al. Experimental validation of a convolution-based ultrasound image formation model using a planar arrangement of micrometer-scale scatterers
CN111965257A (en) Space weighting optimized rapid ultrasonic plane wave imaging detection method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190927

Termination date: 20211129