WO2018090654A1 - Fpga-based space compound imaging method and device - Google Patents
Fpga-based space compound imaging method and device Download PDFInfo
- Publication number
- WO2018090654A1 WO2018090654A1 PCT/CN2017/092936 CN2017092936W WO2018090654A1 WO 2018090654 A1 WO2018090654 A1 WO 2018090654A1 CN 2017092936 W CN2017092936 W CN 2017092936W WO 2018090654 A1 WO2018090654 A1 WO 2018090654A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- image data
- spatial composite
- composite image
- processing
- ultrasonic non
- Prior art date
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 57
- 150000001875 compounds Chemical class 0.000 title abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 105
- 238000000034 method Methods 0.000 claims abstract description 44
- 239000002131 composite material Substances 0.000 claims description 173
- 238000012937 correction Methods 0.000 claims description 37
- 238000002604 ultrasonography Methods 0.000 claims description 33
- 238000009499 grossing Methods 0.000 claims description 23
- 230000003313 weakening effect Effects 0.000 claims description 15
- 238000007781 pre-processing Methods 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- 239000000872 buffer Substances 0.000 description 14
- 238000010586 diagram Methods 0.000 description 6
- 238000001914 filtration Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 239000000523 sample Substances 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000003139 buffering effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000005215 recombination Methods 0.000 description 2
- 230000006798 recombination Effects 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000013329 compounding Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/60—Rotation of whole images or parts thereof
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T9/00—Image coding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Definitions
- the present invention relates to the field of medical ultrasonic imaging technology, and in particular to an FPGA-based spatial composite imaging method and apparatus.
- Spatial composite imaging in medical ultrasound imaging systems is an imaging method that scans scanned objects at different angles and then superimposes the corresponding pixels of these different angles to form an image.
- Space composite imaging technology can improve the detail performance of the image, reduce the influence of speckle noise, clutter and other ultrasonic artifacts on the image quality, enhance the resolving power between tissues, and obviously improve the low contrast tissue in the ultrasound image. The sharpness of the tiny lesions clearly shows the boundaries between the tissues.
- the existing spatial composite imaging technology scans the target tissue from different angles by the deflection of the scan line, and obtains images of these different scan angles, and then weights the pixels corresponding to the images to form an image.
- spatial compounding algorithms are implemented in the CPU through software programming.
- the spatial composite imaging it is not necessary to wait for all angles to be scanned before performing a composite image processing, but a scrolling process is adopted.
- the angle A/B/C the angle A/B/C.
- the previous cycle scan obtains the image of three angles of A0/B0/C0 to obtain a spatial composite image Y0.
- an FPGA-based spatial composite imaging method including:
- the pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving are superimposed to generate a spatial composite image.
- the step of performing interleaving processing on the ultrasound non-spatial composite image data includes:
- the interleaved image data is stored to an external storage device.
- the step of performing coordinate correction on the interlaced ultrasonic non-spatial composite image data comprises:
- Reading parameter information required for coordinate correction from the external storage device
- the coordinate corrected image data is stored to the external storage device.
- the step of performing coordinate correction on the interlaced ultrasonic non-spatial composite image data comprises:
- the interpolated image data is registered by a linear interpolation algorithm
- the step of performing deinterleaving processing on the coordinate-corrected image data includes:
- the image data after deinterleaving is stored to the external storage device.
- the pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving are superimposed, and generating the spatial composite image includes:
- the pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving are superimposed to generate a spatial composite image.
- the method further includes:
- the ultrasonic non-spatial composite image data is subjected to line smoothing processing and/or edge weakening processing.
- the step of performing line smoothing on the ultrasound non-spatial composite image data comprises:
- the segmentation method is used for line smoothing preprocessing, and the original value is maintained at the edge of the line, and the third stage FIR filter is used for smoothing the middle segment line;
- the edge of the line is defined as n ⁇ 2, n ⁇ N-2, and the middle segment line is positioned as 3 ⁇ n ⁇ n-3.
- the invention also provides an FPGA-based spatial composite imaging device, comprising:
- An interleaving module configured to perform interleaving processing on the ultrasound non-spatial composite image data
- a coordinate correction module configured to perform coordinate correction on the interlaced ultrasonic non-spatial composite image data to restore the angularly deflected image data to the image data with zero angle deflection
- a deinterleaving module configured to perform deinterleaving processing on the coordinate corrected image data
- Composite module for ultrasonic non-spatial composite image with different deflection angles after deinterleaving A spatial composite image is generated by superimposing corresponding pixel points.
- it also includes:
- a pre-processing module configured to: after acquiring the ultrasonic non-spatial composite image data of different deflection angles, the ultrasonic non-spatial composite image data before the step of performing coordinate correction on the interlaced ultrasonic non-spatial composite image data Perform line smoothing and/or edge weakening.
- the FPGA-based spatial composite imaging method and device obtain ultrasonic non-spatial composite image data with different deflection angles; interleave processing ultrasonic non-spatial composite image data; and perform ultrasonic non-spatial composite image data after interleaving processing Performing coordinate correction to restore the image data of the angular deflection to the image data of the zero angle deflection; deinterleaving the image data after the coordinate correction; and corresponding to the ultrasonic non-spatial composite image data of the different deflection angles after deinterleaving
- the pixels are superimposed to generate a spatial composite image.
- the application reduces the complexity of the spatial composite algorithm implementation, can effectively improve the speed of the spatial composite operation, reduce the processing time of the image, and improve the imaging frame rate.
- FIG. 1 is a flow chart of a specific implementation manner of an FPGA-based spatial composite imaging method provided by the present invention
- FIG. 2 is a flowchart of another specific implementation manner of an FPGA-based spatial composite imaging method provided by the present invention.
- 3 is a flow chart showing steps of performing interleaving processing on ultrasonic non-spatial composite image data in an embodiment of the present invention
- FIG. 5 is a schematic diagram of still another specific implementation manner of the FPGA-based spatial composite imaging method provided by the present invention.
- FIG. 6 is a schematic diagram of an overall framework of a spatial composite imaging FPGA implementation
- FIG. 7 is a flow chart of data processing implemented by a spatial composite imaging FPGA
- FIG. 8 is a structural block diagram of an FPGA-based spatial composite imaging apparatus according to an embodiment of the present invention.
- FIG. 1 A flowchart of a specific implementation manner of the FPGA-based spatial composite imaging method provided by the present invention is shown in FIG. 1 , and the method includes:
- Step S101 Acquire ultrasonic non-spatial composite image data of different deflection angles.
- the ultrasonic probe emits an ultrasonic signal and receives ultrasonic echo signals of different deflection angles.
- the ultrasonic echo signals include ultrasonic non-spatial composite image data, that is, ultrasonic B image echo signals, and include beam signals and ADC acquisition echoes. Other signals such as signals.
- the ultrasonic echo signals of different deflection angles are stored, and may be stored in an external storage device such as a DDR memory.
- the ultrasonic non-spatial composite image data is read from the external storage device in step S101 of the embodiment.
- Step S102 Perform interleaving processing on the ultrasound non-spatial composite image data.
- Interleaving is essentially the realization of changing the information structure to the maximum without changing the content of the information.
- the complexity of the spatial composite algorithm can be reduced.
- Step S103 Perform coordinate correction on the interlaced ultrasonic non-spatial composite image data to restore the angularly deflected image data to the image data of the zero angle deflection.
- Correcting the coordinates can uniformly restore the images of different deflection angles to the undeflected image for subsequent spatial recombination operations.
- Step S104 Perform deinterleave processing on the coordinate-corrected image data.
- the deinterleaving process is the inverse process of the interleaving process in step S102, and the image data is restored to the original image rule and size by deinterleaving the image data.
- Step S105 superimposing pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving to generate a spatial composite image.
- the spatial composite image can be obtained by fusing the ultrasonic non-spatial composite image data with different deflection angles.
- the FPGA-based spatial composite imaging method acquires ultrasonic non-spatial composite image data with different deflection angles; interleaves the ultrasonic non-spatial composite image data; coordinates the interlaced ultrasonic non-spatial composite image data Correcting to restore image data of the angular deflection to image data of zero angle deflection; deinterleaving the image data after coordinate correction; and pixel points corresponding to ultrasonic non-spatial composite image data of different deflection angles after deinterleaving Perform superposition to generate a spatial composite image.
- the application reduces the complexity of the spatial composite algorithm implementation, can effectively improve the speed of the spatial composite operation, reduce the processing time of the image, and improve the imaging frame rate.
- the spatial composite algorithm is implemented by FPGA. Since the internal buffer space of the FPGA is limited, an external storage device is required for data caching.
- the external storage device may be specifically an off-chip high-speed DDR3 memory (Double Data Rate SDRAM 3rd generation), and other devices capable of realizing storage are not limited to this one.
- the process may specifically include:
- Step S201 Acquire ultrasonic non-spatial composite image data of different deflection angles from an external storage device.
- the ultrasound probe transmits an ultrasound signal and stores the received ultrasound non-spatial composite image data of different deflection angles into an external storage device.
- the ultrasound non-spatial composite image data may be subjected to signal amplification, analog/digital conversion processing, beam synthesis, and digital before storage. Signal processing and other processes.
- Step S202 performing pipelined interleaving processing on the ultrasonic non-spatial composite image data, and storing the interleaved image data to an external storage device;
- the process of performing interleaving processing on the ultrasound non-spatial composite image data may specifically include:
- Step S2021 Segmenting the ultrasound non-spatial composite image data by using a preset window.
- the FPGA When processing large-size image data, the FPGA requires a large amount of storage space if it is processed once. Therefore, in practical applications, due to limitations of storage space and logical resources, it is necessary to divide large-sized image data into small-sized image data for processing. Therefore, in this embodiment, the ultrasound non-spatial composite image data is segmented by a preset window.
- the size of the preset window may be 64 x 64 windows, or 128 x 128 or 32 x 32, which does not affect the implementation of the present invention.
- Step S2022 Perform pipelined interleaving processing on the segmented ultrasound non-spatial composite image data.
- the logic hardware circuit inside the FPGA divides the ultrasonic non-spatial composite image data into a plurality of 64 ⁇ 64 windows, and uses the window as a minimum unit to carry the data of the B image DDR3 buffer space to the logical internal Ram storage space.
- the logical hardware circuit performs pipelined interleaving processing on the logically buffered data in a window unit, and uploads the interleaved window data to an external cache device.
- the ping-pong operation is performed by using the dual-chip Ram, so that the logical internal cache data can be continuously pipelined.
- the ping-pong operation is applied to the processing technique of data flow control, and the data stream is equally allocated to two data buffers, and then switched according to the beat and the mutual cooperation, so that the data stream is in a continuous state to achieve data-free.
- Pipelined processing is a technical means in high-speed design that divides the data processing flow into several steps, the output of the previous step being the input to the next step.
- the biggest feature is that the processing of each step of the data stream is continuous in time, and the efficiency of data processing can be further improved by such setting.
- a window of 64 ⁇ 64 is used as a minimum processing unit, which is a preferred solution of the present invention in combination with the characteristics of the Xilinx tdpram resource, and the interlace can also be used.
- the window size of the control is modified to 128 ⁇ 128 or 32 ⁇ 32, etc., that is, the window size is different from the preferred scheme in the present invention, which does not affect the implementation of the present invention.
- Step S203 Acquiring the data after the interleaving process from the external storage device, performing coordinate correction, and storing the corrected data to the external storage device;
- the step of performing coordinate correction may be specifically as follows:
- Step S2031 Acquire interlaced image data from an external storage device according to a preset direction.
- the preset direction may be specifically: according to the depth direction of the ultrasonic echo signal, that is, the direction in which the ultrasonic probe is measured by the human body.
- the logic hardware circuit reads the interleaved data from an external storage device in units of different lines of the same depth.
- the scheduling method can also be implemented by using queue flow control to improve the efficiency of data processing.
- Step S2032 Read parameter information required for coordinate correction from an external storage device.
- processing may be performed in parallel between step S2031 and step S2032, and there is no order limitation.
- Step S2033 Perform coordinate correction on the interlaced ultrasonic non-spatial composite image data when the interleaved image data and the parameter information satisfy a preset threshold condition.
- a threshold condition is set before the data signal processing, and the preset threshold condition is that both the parameter information required for the coordinate correction and the image data after the interleaving are reached.
- the preset threshold condition may specifically determine whether the interleaved buffer area data and the status of the queue corresponding to the parameter are not empty, that is, whether data and parameters exist in the corresponding queue. The subsequent operations can be performed only when the interleaving buffer data and the parameter information arrive and the determination meets the preset threshold condition. When the preset threshold condition is not met, the process of performing task scheduling is returned until the preset threshold condition is met.
- Coordinate correction includes two processes: registration processing and boundary processing.
- the registration process can use a linear interpolation algorithm to register the image data.
- a Catmull-ROM linear interpolation algorithm is employed.
- the Catmull-ROM interpolation algorithm performs coordinate correction by first performing linear interpolation in the longitudinal direction (depth direction) and then interpolating in the lateral direction (line direction).
- depth direction longitudinal direction
- line direction lateral direction
- other interpolation algorithms can also be used, which will not be repeated here.
- the boundary processing includes: calculating coordinates of pixel points in the image data after registration processing.
- the pixel points are processed according to the preset boundary of the ultrasonic non-spatial composite image, and the pixel points beyond the preset boundary are removed, and the pixel points of the area not covered by the preset boundary are zero-padded.
- the preset boundary is the boundary size range corresponding to the originally acquired ultrasonic non-spatial composite image.
- the parameter information required for coordinate correction specifically includes parameters required for the linear interpolation algorithm in the registration process and size parameters of the preset boundary in the boundary processing process.
- Step S204 The coordinate-corrected image data is read from the external storage device, the coordinate-corrected image data is deinterleaved, and the de-interleaved data is stored to the external storage device.
- Step S205 Acquire deinterleaved data from the external storage device, and superimpose the pixels corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving to generate a spatial composite image.
- the normalization processing is performed, and the superimposed region of the ultrasonic non-spatial composite image is weakened.
- the embodiment of the present application implements a spatial composite algorithm by using an FPGA, and uses an off-chip high-speed DDR memory for data buffering, and the FPGA implements data scheduling and composite operations through a high-speed access interface.
- the application can effectively improve the speed of the spatial composite operation, greatly reduce the processing time of the image, improve the imaging frame rate, and effectively improve the effect of the spatial composite imaging without replacing the CPU.
- the embodiment uses the external storage device buffer interleaving, coordinate correction, and deinterleaving image processing results, and the image data is scheduled to different degrees, thereby effectively improving the scale and efficiency of the FPGA image processing.
- the step of performing the coordinate correction on the ultrasonic non-spatial composite image data after the interleaving process may further include:
- a process of performing line smoothing and/or edge weakening processing on ultrasonic non-spatial composite image data is a process of performing line smoothing and/or edge weakening processing on ultrasonic non-spatial composite image data.
- FIG. 5 is another specific example of the FPGA-based spatial composite imaging method provided by the present invention.
- the embodiment of the present invention uses an off-chip high-speed DDR memory for data caching, and the FPGA uses a high-speed access interface to implement a data scheduling and a compound operation by using a pipelined task control scheduling method.
- the process can specifically include:
- Step S301 After the spatial recombination starts, the ultrasonic non-spatial composite image data in the DDR3 is read for interleaving processing.
- Step S302 Read the interleave buffer data and the required parameters from the DDR3.
- the parameters required for processing include parameters and coefficients required in data processing such as line smoothing, edge weakening, and registration.
- Step S303 determining whether the interleave buffer data and the parameter information in the data processing process meet the preset threshold condition, that is, whether the data and the parameters have all arrived; if yes, proceed to step S304; if not, return to step S302.
- the logic hardware circuit When reading the line data of all lines and the same depth after interleaving, the logic hardware circuit carries data and digital signal processing from the buffer DDR3 space in units of different lines of the same depth.
- the scheduling mode is implemented by using queue flow control, and threshold value judgment is set in the front end of the digital signal processing.
- the depth direction refers to the vertical direction, that is, the direction in which the probe is measured to the human body, and the horizontal direction is the method of line scanning.
- Step S304 When the threshold judgment condition is satisfied, the interlaced buffer area data is subjected to line smoothing processing to reduce the near field jagged shape and the far field mosaic of the ultrasonic image.
- segmentation processing may be specifically adopted, including maintaining the original value and the intermediate line data FIR filtering at the line edge, wherein the number of buses is N, and the current line number is n, Then the line edge is defined as n ⁇ 2, n ⁇ N-2, and the middle section is defined as 3 ⁇ n ⁇ n-3.
- the logic hardware circuit performs special processing on the edge of the line by changing the weighting coefficient at the edge of the line, that is, keeping the original value unchanged; the middle line is smoothed, and the third-order FIR filtering is used, and the result of the smoothed output is smoothed.
- the third-order FIR filtering method is adopted, which satisfies the requirement of smoothing the ultrasound image, and is easier to implement than the fifth-order filtering and the seventh-order filtering. Passing on ultrasound Like line smoothing, it can reduce near-field jagged and far-field mosaic.
- Step S305 Perform edge weakening processing on the data after the interleaving buffer smoothing to reduce artifacts after the ultrasound image is merged.
- Step S306 Perform registration and boundary processing on the data of the interleaved buffer edge weakening process to perform coordinate correction on the images of different deflection angles.
- the embodiment of the present invention adopts the Catmull-ROM interpolation algorithm for coordinate correction, first performing linear interpolation in the longitudinal direction, and then performing interpolation calculation in the lateral direction (line direction).
- Registration processing of line data at the same depth including FPGA implementation of Catmull-ROM linear interpolation algorithm, synchronization operation between interpolation coefficient and echo data required by Catmull-ROM algorithm, and dual parallel registration processing unit.
- the two-way parallel registration processing unit multiplies the data and the interpolation coefficient by a 2 ⁇ 4 operation window in a 4 ⁇ 4 shift window according to the Catmull-ROM interpolation rule, and outputs a registration process every four clock cycles. result.
- the registration processing may specifically adopt a two-way parallel registration processing unit, which is a preferred scheme made by the present invention in combination with the FPGA system logic resources of the ultrasound system, and may have other methods.
- the registration processing unit can be extended to four channels, six channels, and eight channels to increase the rate of registration processing, that is, the number of paths of the processing unit to be registered in the preferred scheme of the present invention is not the same value, which does not affect the present invention. Implementation of the invention.
- the logic hardware circuit calculates the coordinates of the pixel points in the image data after the registration interpolation, and then constrains the registered pixel points according to the boundary of the ultrasonic B echo image. Discard pixels that are beyond the boundary of the B image, and zeros are added to areas that are not registered.
- Step S307 Deinterleave the interleave buffer data to restore the image to the original image rule and size.
- Step S308 Normalize the data after the deinterleaving process, and fuse the ultrasound images with different deflection angles.
- the embodiment of the present invention further includes: performing signal amplification, analog/digital conversion processing, beam synthesis, and the ultrasonic echo signal. Digital signal processing.
- the data stream processing is performed according to interleaving, line smoothing, edge weakening,
- the registration, the boundary processing, the de-interlacing, and the normalization are sequentially performed.
- This method is a preferred solution implemented by the present invention in combination with an algorithm, and of course, there may be other methods.
- the position of the line smoothing and the edge weakening may be interchanged, or the position of the line smoothing and the edge weakening may be completed before the interleaving, that is, the data stream processing sequence of the present invention may be exchanged.
- FIG. 6 is a schematic diagram of the overall framework of the space composite imaging FPGA implementation.
- the upstream of the system is the user configuration interface and the DDR3 high-speed read interface, and the downstream of the system is the DDR3 high-speed storage interface.
- the entire system data stream is from left to right.
- the interleaving processing module acquires ultrasonic echo signals of different deflection angles from DDR3 for interleaving processing; then the signal processing module (SCON_PRO) performs the interleaved ultrasonic echo signals.
- the de-interleaving module restores the signal-processed echo signal to the ultrasound signal in the usual format; finally, the different deflection angles are returned by the fusion (NOR) module. Boss is normalized to obtain a spatial composite image.
- Figure 7 is a flow chart of data processing implemented by a spatial composite imaging FPGA.
- the figure identifies the signal name, symbol, bit width and precision of the input/output of the data processing link.
- DownSizer is mainly used for bit width truncation of digital signals
- UpSizer is used for bit width expansion of digital signals
- Edge_weak is edge weakening processing module
- Scon_Cal is registration coordinate transformation module.
- the interleaving preprocessing of the ultrasonic echo signals can greatly reduce the complexity of the spatial composite algorithm implementation, and thus can pass the field programmable logic when the actual imaging frame rate is low and the CPU performance of the medical ultrasound system is insufficient.
- the array FPGA performs storage scheduling, hardware acceleration, and spatial composite processing on the ultrasonic non-spatial composite image data.
- the FPGA-based spatial composite imaging apparatus provided by the embodiment of the present invention is described below.
- the FPGA-based spatial composite imaging apparatus described below and the FPGA-based spatial composite imaging method described above can be mutually referenced.
- FIG. 8 is a structural block diagram of an FPGA-based spatial composite imaging apparatus according to an embodiment of the present invention.
- the space-based composite imaging apparatus based on the FPGA may include:
- the acquiring module 100 is configured to acquire ultrasonic non-spatial composite image data with different deflection angles
- the interleaving module 200 is configured to perform interleaving processing on the ultrasound non-spatial composite image data
- the coordinate correction module 300 is configured to perform coordinate correction on the interlaced ultrasonic non-spatial composite image data to restore the angularly deflected image data to the image data with zero angle deflection;
- the deinterleaving module 400 is configured to perform deinterleaving processing on the coordinate corrected image data
- the composite module 500 is configured to superimpose pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving to generate a spatial composite image.
- the application reduces the complexity of the spatial composite algorithm implementation, can effectively improve the spatial composite operation speed, reduce the processing time of one frame of image, and improve the imaging frame rate.
- this application replaces the traditional MPU processor with a field programmable logic array FPGA to realize the spatial composite function on the medical ultrasound image service.
- FPGA hardware acceleration can effectively improve the rate and imaging frame rate of spatial composite implementation.
- the present invention first pre-interleaves the ultrasonic non-spatial composite image data, and then performs line smoothing, edge weakening, and registration processing, which greatly reduces the complexity of the spatial composite FPGA implementation and saves the logic of the FPGA. Resources.
- the steps of the method or algorithm described in connection with the embodiments disclosed herein may be directly applied to the hard Implemented by a piece of software, a software module executed by a processor, or a combination of both.
- the software module can be placed in random access memory (RAM), memory, read only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or technical field. Any other form of storage medium known.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
Abstract
An FPGA-based space compound imaging method and device. The method comprises: acquiring ultrasonic non-space compound image data at different deflection angles (S101); carrying out interleaving processing on the ultrasonic non-space compound image data (S102); correcting the coordinates of the ultrasonic non-space compound image data, which has undergone interleaving processing, so as to recover image data subjected to angle deflection to image data with zero angle deflection (S103); carrying out de-interleaving processing on the image data, the coordinates thereof having been corrected (S104); and superposing pixel points corresponding to the de-interleaved ultrasonic non-space compound image data with different deflection angles, and generating a space compound image (S105). The present invention reduces the complexity of implementing a space compound algorithm, can effectively improve a space compound calculation speed, reduces the processing time for an image, and improves an imaging frame rate.
Description
本申请要求于2016年11月21日提交中国专利局、申请号为201611023419.4、发明名称为“基于FPGA的空间复合成像方法以及装置”的中国专利申请的优先权,其全部内容通过引用结合在本申请中。The present application claims priority to Chinese Patent Application No. 201611023419.4, entitled "Synthesis-Based Space Imaging Method and Apparatus Based on FPGA", filed on November 21, 2016, the entire contents of which are incorporated herein by reference. In the application.
本发明涉及医学超声波成像技术领域,特别是涉及一种基于FPGA的空间复合成像方法以及装置。The present invention relates to the field of medical ultrasonic imaging technology, and in particular to an FPGA-based spatial composite imaging method and apparatus.
医学超声成像系统中的空间复合成像是一种沿不同角度对扫描对象进行扫描,然后将这些不同角度的图像对应的像素点进行叠加,形成一幅图像的成像方法。空间复合成像技术可以提高图像的细节表现,同时降低图像中的斑点噪声、杂波以及其他超声伪像对图像质量的影响,增强组织间的分辨能力,可明显地提高超声图像中低对比度组织及微小病变的清晰度,清晰地显示组织间的边界。Spatial composite imaging in medical ultrasound imaging systems is an imaging method that scans scanned objects at different angles and then superimposes the corresponding pixels of these different angles to form an image. Space composite imaging technology can improve the detail performance of the image, reduce the influence of speckle noise, clutter and other ultrasonic artifacts on the image quality, enhance the resolving power between tissues, and obviously improve the low contrast tissue in the ultrasound image. The sharpness of the tiny lesions clearly shows the boundaries between the tissues.
现有的空间复合成像技术是通过扫描线的偏转,从不同角度对目标组织进行扫描,获取这些不同扫描角度的图像后,将这些图像对应的像素点加权求和以合成一幅图像。通常,空间复合算法都是在CPU通过软件编程实现。空间复合成像在应用实现时,往往不会等所有角度都扫描完了之后再进行一次复合图像处理,而是采取一种滚动处理。以三个角度偏转的空间复合成像为例,角度A/B/C。上一个周期扫描获得A0/B0/C0三个角度的图像从而获得了一幅空间复合的图像Y0,当下一个扫描周期的A1采集到后,Y1的输出就可以用B0/C0/A1获得,原理上可以简单的表示为Y1=Y0-A0+A1,以此类推,每采集到一个偏转角度的图像即可获得一帧空间复合的输出。The existing spatial composite imaging technology scans the target tissue from different angles by the deflection of the scan line, and obtains images of these different scan angles, and then weights the pixels corresponding to the images to form an image. Usually, spatial compounding algorithms are implemented in the CPU through software programming. When the spatial composite imaging is implemented, it is not necessary to wait for all angles to be scanned before performing a composite image processing, but a scrolling process is adopted. Taking spatial composite imaging with three angles of deflection as an example, the angle A/B/C. The previous cycle scan obtains the image of three angles of A0/B0/C0 to obtain a spatial composite image Y0. When A1 of the next scan cycle is acquired, the output of Y1 can be obtained by B0/C0/A1. It can be simply represented as Y1=Y0-A0+A1, and so on, a single spatial composite output can be obtained for each image with a deflection angle.
如果空间复合成像的帧率要达到普通B模式成像(非空间复合)的帧率,对CPU的运算能力和数据调度性能有一定要求。如果CPU性能不够,就会导致实现空间复合算法时间较长,导致超声成像的帧率变低。因此,
提供一种提高处理速度以及成像帧率的空间复合成像方法是非常有必要的。If the frame rate of spatial composite imaging is to reach the frame rate of ordinary B-mode imaging (non-spatial composite), there are certain requirements for the computing power and data scheduling performance of the CPU. If the CPU performance is not enough, it will lead to a longer time to implement the spatial composite algorithm, resulting in lower frame rate of ultrasound imaging. Therefore,
It is very necessary to provide a spatial composite imaging method that increases processing speed and imaging frame rate.
发明内容Summary of the invention
本发明的目的是提供一种空间复合成像方法以及装置,以解决现有空间复合的处理速度较慢、成像帧率较低的问题。It is an object of the present invention to provide a spatial composite imaging method and apparatus for solving the problem of slow processing speed and low imaging frame rate of the existing spatial composite.
为解决上述技术问题,本发明提供一种基于FPGA的空间复合成像方法,包括:To solve the above technical problem, the present invention provides an FPGA-based spatial composite imaging method, including:
获取不同偏转角度的超声非空间复合图像数据;Acquiring ultrasonic non-spatial composite image data of different deflection angles;
对所述超声非空间复合图像数据进行交织处理;Interleaving the ultrasound non-spatial composite image data;
对交织处理后的超声非空间复合图像数据进行坐标校正,以将发生角度偏转的图像数据恢复至零角度偏转的图像数据;Performing coordinate correction on the interlaced ultrasonic non-spatial composite image data to restore the angularly deflected image data to the zero angle deflected image data;
对坐标校正后的图像数据进行解交织处理;Performing deinterleaving processing on the coordinate corrected image data;
将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像。The pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving are superimposed to generate a spatial composite image.
可选地,所述对所述超声非空间复合图像数据进行交织处理的步骤包括:Optionally, the step of performing interleaving processing on the ultrasound non-spatial composite image data includes:
采用预设窗口对所述超声非空间复合图像数据进行分割;Separating the ultrasound non-spatial composite image data by using a preset window;
对分割后的所述超声非空间复合图像数据进行流水线式的交织处理;Performing pipelined interleaving processing on the segmented ultrasound non-spatial composite image data;
将交织处理后的图像数据存储至外部存储设备。The interleaved image data is stored to an external storage device.
可选地,所述对交织处理后的超声非空间复合图像数据进行坐标校正的步骤包括:Optionally, the step of performing coordinate correction on the interlaced ultrasonic non-spatial composite image data comprises:
按照预设方向从外部存储设备获取交织处理后的图像数据;Acquiring the interleaved image data from the external storage device according to a preset direction;
从所述外部存储设备读取坐标校正所需的参数信息;Reading parameter information required for coordinate correction from the external storage device;
当交织处理后的图像数据以及所述参数信息满足预设阈值条件时,对交织处理后的超声非空间复合图像数据进行坐标校正;Performing coordinate correction on the interlaced ultrasonic non-spatial composite image data when the interleaved image data and the parameter information satisfy a preset threshold condition;
将坐标校正后的图像数据存储至所述外部存储设备。The coordinate corrected image data is stored to the external storage device.
可选地,所述对交织处理后的超声非空间复合图像数据进行坐标校正的步骤包括:
Optionally, the step of performing coordinate correction on the interlaced ultrasonic non-spatial composite image data comprises:
采用线性插值算法对交织处理后的图像数据进行配准;The interpolated image data is registered by a linear interpolation algorithm;
对配准处理后的图像数据中像素点的坐标进行运算,根据超声非空间复合图像的预设边界对所述像素点进行处理,去除超出所述预设边界的像素点,未被所述预设边界覆盖的区域的像素点补零。Performing operations on the coordinates of the pixel points in the image data after the registration processing, processing the pixel points according to a preset boundary of the ultrasonic non-spatial composite image, and removing pixel points beyond the preset boundary, without being pre-predicted Set the pixel of the area covered by the boundary to zero.
可选地,所述对坐标校正后的图像数据进行解交织处理的步骤包括:Optionally, the step of performing deinterleaving processing on the coordinate-corrected image data includes:
从所述外部存储设备读取坐标校正后的图像数据;Reading coordinate corrected image data from the external storage device;
对坐标校正后的图像数据进行解交织处理;Performing deinterleaving processing on the coordinate corrected image data;
将解交织处理后的图像数据存储至所述外部存储设备。The image data after deinterleaving is stored to the external storage device.
可选地,所述将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像包括:Optionally, the pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving are superimposed, and generating the spatial composite image includes:
从所述外部存储设备获取解交织处理后的图像数据;Acquiring image data after deinterleaving processing from the external storage device;
将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像。The pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving are superimposed to generate a spatial composite image.
可选地,在所述获取不同偏转角度的超声非空间复合图像数据的步骤之后,在所述对交织处理后的超声非空间复合图像数据进行坐标校正的步骤之前还包括:Optionally, after the step of acquiring the ultrasonic non-spatial composite image data of different deflection angles, before the step of performing coordinate correction on the ultrasonic non-spatial composite image data after the interleaving processing, the method further includes:
对所述超声非空间复合图像数据进行线平滑处理和/或边缘弱化处理。The ultrasonic non-spatial composite image data is subjected to line smoothing processing and/or edge weakening processing.
可选地,对所述超声非空间复合图像数据进行线平滑处理的步骤包括:Optionally, the step of performing line smoothing on the ultrasound non-spatial composite image data comprises:
采用分段式进行线平滑预处理,对线边缘处保持原始数值,对中间段线采用三阶FIR滤波的方式进行平滑处理;The segmentation method is used for line smoothing preprocessing, and the original value is maintained at the edge of the line, and the third stage FIR filter is used for smoothing the middle segment line;
其中,设总线数为N,当前线号为n,则线边缘处定义为n≤2,n≥N-2,中间段线定位为3≤n≤n-3。Wherein, if the number of buses is N and the current line number is n, the edge of the line is defined as n≤2, n≥N-2, and the middle segment line is positioned as 3≤n≤n-3.
本发明还提供了一种基于FPGA的空间复合成像装置,包括:The invention also provides an FPGA-based spatial composite imaging device, comprising:
获取模块,用于获取不同偏转角度的超声非空间复合图像数据;Obtaining module for acquiring ultrasonic non-spatial composite image data with different deflection angles;
交织模块,用于对所述超声非空间复合图像数据进行交织处理;An interleaving module, configured to perform interleaving processing on the ultrasound non-spatial composite image data;
坐标校正模块,用于对交织处理后的超声非空间复合图像数据进行坐标校正,以将发生角度偏转的图像数据恢复至零角度偏转的图像数据;a coordinate correction module, configured to perform coordinate correction on the interlaced ultrasonic non-spatial composite image data to restore the angularly deflected image data to the image data with zero angle deflection;
解交织模块,用于对坐标校正后的图像数据进行解交织处理;a deinterleaving module, configured to perform deinterleaving processing on the coordinate corrected image data;
复合模块,用于将解交织后的不同偏转角度的超声非空间复合图像数
据对应的像素点进行叠加,生成空间复合图像。Composite module for ultrasonic non-spatial composite image with different deflection angles after deinterleaving
A spatial composite image is generated by superimposing corresponding pixel points.
可选地,还包括:Optionally, it also includes:
预处理模块,用于在获取不同偏转角度的超声非空间复合图像数据的步骤之后,在对交织处理后的超声非空间复合图像数据进行坐标校正的步骤之前,对所述超声非空间复合图像数据进行线平滑处理和/或边缘弱化处理。a pre-processing module, configured to: after acquiring the ultrasonic non-spatial composite image data of different deflection angles, the ultrasonic non-spatial composite image data before the step of performing coordinate correction on the interlaced ultrasonic non-spatial composite image data Perform line smoothing and/or edge weakening.
本发明所提供的基于FPGA的空间复合成像方法及装置,通过获取不同偏转角度的超声非空间复合图像数据;对超声非空间复合图像数据进行交织处理;对交织处理后的超声非空间复合图像数据进行坐标校正,以将发生角度偏转的图像数据恢复至零角度偏转的图像数据;对坐标校正后的图像数据进行解交织处理;将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像。本申请降低了空间复合算法实现的复杂度,能够有效提高空间复合运算的速度,降低图像的处理时间,提高成像帧率。The FPGA-based spatial composite imaging method and device provided by the invention obtain ultrasonic non-spatial composite image data with different deflection angles; interleave processing ultrasonic non-spatial composite image data; and perform ultrasonic non-spatial composite image data after interleaving processing Performing coordinate correction to restore the image data of the angular deflection to the image data of the zero angle deflection; deinterleaving the image data after the coordinate correction; and corresponding to the ultrasonic non-spatial composite image data of the different deflection angles after deinterleaving The pixels are superimposed to generate a spatial composite image. The application reduces the complexity of the spatial composite algorithm implementation, can effectively improve the speed of the spatial composite operation, reduce the processing time of the image, and improve the imaging frame rate.
为了更清楚的说明本发明实施例或现有技术的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the prior art, the drawings used in the embodiments or the description of the prior art will be briefly described below. Obviously, the drawings in the following description are merely Some embodiments of the present invention may also be used to obtain other drawings based on these drawings without departing from the art.
图1为本发明所提供的基于FPGA的空间复合成像方法的一种具体实施方式的流程图;1 is a flow chart of a specific implementation manner of an FPGA-based spatial composite imaging method provided by the present invention;
图2为本发明所提供的基于FPGA的空间复合成像方法的另一种具体实施方式的流程图;2 is a flowchart of another specific implementation manner of an FPGA-based spatial composite imaging method provided by the present invention;
图3为本发明实施例中对超声非空间复合图像数据进行交织处理步骤的流程图;3 is a flow chart showing steps of performing interleaving processing on ultrasonic non-spatial composite image data in an embodiment of the present invention;
图4为本发明实施例中进行坐标校正步骤的流程图;4 is a flowchart of performing a coordinate correction step in an embodiment of the present invention;
图5本发明所提供的基于FPGA的空间复合成像方法的又一种具体实施方式示意图;
FIG. 5 is a schematic diagram of still another specific implementation manner of the FPGA-based spatial composite imaging method provided by the present invention; FIG.
图6为空间复合成像FPGA实现的整体框架示意图;6 is a schematic diagram of an overall framework of a spatial composite imaging FPGA implementation;
图7为空间复合成像FPGA实现的数据处理流程图;7 is a flow chart of data processing implemented by a spatial composite imaging FPGA;
图8为本发明实施例提供的基于FPGA的空间复合成像装置的结构框图。FIG. 8 is a structural block diagram of an FPGA-based spatial composite imaging apparatus according to an embodiment of the present invention.
为了使本技术领域的人员更好地理解本发明方案,下面结合附图和具体实施方式对本发明作进一步的详细说明。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The present invention will be further described in detail below in conjunction with the drawings and embodiments. It is apparent that the described embodiments are only a part of the embodiments of the invention, and not all of the embodiments. 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.
本发明所提供的基于FPGA的空间复合成像方法的一种具体实施方式的流程图,如图1所示,该方法包括:A flowchart of a specific implementation manner of the FPGA-based spatial composite imaging method provided by the present invention is shown in FIG. 1 , and the method includes:
步骤S101:获取不同偏转角度的超声非空间复合图像数据。Step S101: Acquire ultrasonic non-spatial composite image data of different deflection angles.
超声探头发射超声信号,并接收不同偏转角度的超声回波信号,超声回波信号中除了包含超声非空间复合图像数据,即超声B图像回波信号以外,还包括有波束信号、ADC采集回波信号等其他信号。将不同偏转角度的超声回波信号进行存储,具体可以存储在外部存储设备,如DDR存储器中。在本实施例步骤S101中从该外部存储设备中读取超声非空间复合图像数据。The ultrasonic probe emits an ultrasonic signal and receives ultrasonic echo signals of different deflection angles. The ultrasonic echo signals include ultrasonic non-spatial composite image data, that is, ultrasonic B image echo signals, and include beam signals and ADC acquisition echoes. Other signals such as signals. The ultrasonic echo signals of different deflection angles are stored, and may be stored in an external storage device such as a DDR memory. The ultrasonic non-spatial composite image data is read from the external storage device in step S101 of the embodiment.
步骤S102:对所述超声非空间复合图像数据进行交织处理。Step S102: Perform interleaving processing on the ultrasound non-spatial composite image data.
交织处理从其本质上来说就是实现最大限度的改变信息结构而不改变信息内容。通过对超声非空间复合图像数据进行交织预处理,能够降低空间复合算法实现的复杂度。Interleaving is essentially the realization of changing the information structure to the maximum without changing the content of the information. By interpolating the ultrasound non-spatial composite image data, the complexity of the spatial composite algorithm can be reduced.
步骤S103:对交织处理后的超声非空间复合图像数据进行坐标校正,以将发生角度偏转的图像数据恢复至零角度偏转的图像数据。Step S103: Perform coordinate correction on the interlaced ultrasonic non-spatial composite image data to restore the angularly deflected image data to the image data of the zero angle deflection.
对坐标进行校正能够将不同偏转角度的图像统一恢复至未偏转时的图像,以便后续进行空间复合的操作。
Correcting the coordinates can uniformly restore the images of different deflection angles to the undeflected image for subsequent spatial recombination operations.
步骤S104:对坐标校正后的图像数据进行解交织处理。Step S104: Perform deinterleave processing on the coordinate-corrected image data.
本步骤中解交织处理为步骤S102中的交织处理的逆过程,通过对图像数据进行解交织处理,将图像数据还原为原始图像规则以及尺寸。In this step, the deinterleaving process is the inverse process of the interleaving process in step S102, and the image data is restored to the original image rule and size by deinterleaving the image data.
步骤S105:将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像。Step S105: superimposing pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving to generate a spatial composite image.
将不同偏转角度的超声非空间复合图像数据进行融合,即可得到空间复合图像。The spatial composite image can be obtained by fusing the ultrasonic non-spatial composite image data with different deflection angles.
本发明所提供的基于FPGA的空间复合成像方法,通过获取不同偏转角度的超声非空间复合图像数据;对超声非空间复合图像数据进行交织处理;对交织处理后的超声非空间复合图像数据进行坐标校正,以将发生角度偏转的图像数据恢复至零角度偏转的图像数据;对坐标校正后的图像数据进行解交织处理;将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像。本申请降低了空间复合算法实现的复杂度,能够有效提高空间复合运算的速度,降低图像的处理时间,提高成像帧率。The FPGA-based spatial composite imaging method provided by the invention acquires ultrasonic non-spatial composite image data with different deflection angles; interleaves the ultrasonic non-spatial composite image data; coordinates the interlaced ultrasonic non-spatial composite image data Correcting to restore image data of the angular deflection to image data of zero angle deflection; deinterleaving the image data after coordinate correction; and pixel points corresponding to ultrasonic non-spatial composite image data of different deflection angles after deinterleaving Perform superposition to generate a spatial composite image. The application reduces the complexity of the spatial composite algorithm implementation, can effectively improve the speed of the spatial composite operation, reduce the processing time of the image, and improve the imaging frame rate.
本申请将空间复合算法采用FPGA进行实现,由于FPGA内部的缓存空间有限,因此需要采用外部存储设备进行数据的缓存。外部存储设备可以具体为片外高速DDR3存储器(Double Data Rate SDRAM第3代),其他能够实现存储的设备均可,并不限于这一种。In this application, the spatial composite algorithm is implemented by FPGA. Since the internal buffer space of the FPGA is limited, an external storage device is required for data caching. The external storage device may be specifically an off-chip high-speed DDR3 memory (Double Data Rate SDRAM 3rd generation), and other devices capable of realizing storage are not limited to this one.
下面通过具体实施例对FPGA与外部存储设备进行空间复合任务调度的过程进行进一步详细阐述。The process of spatial composite task scheduling of an FPGA and an external storage device is further elaborated below through a specific embodiment.
如图2本发明所提供的基于FPGA的空间复合成像方法的另一种具体实施方式的流程图所示,该过程具体可以包括:As shown in the flowchart of another embodiment of the FPGA-based spatial composite imaging method provided by the present invention, the process may specifically include:
步骤S201:从外部存储设备中获取不同偏转角度的超声非空间复合图像数据。Step S201: Acquire ultrasonic non-spatial composite image data of different deflection angles from an external storage device.
超声探头发射超声信号,并将接收到的不同偏转角度的超声非空间复合图像数据存储至外部存储设备中。优选地,在存储之前还可以包括对超声非空间复合图像数据进行信号放大、模/数转换处理、波束合成以及数字
信号处理等过程。The ultrasound probe transmits an ultrasound signal and stores the received ultrasound non-spatial composite image data of different deflection angles into an external storage device. Preferably, the ultrasound non-spatial composite image data may be subjected to signal amplification, analog/digital conversion processing, beam synthesis, and digital before storage.
Signal processing and other processes.
步骤S202:对超声非空间复合图像数据进行流水线式的交织处理,将交织处理后的图像数据存储至外部存储设备;Step S202: performing pipelined interleaving processing on the ultrasonic non-spatial composite image data, and storing the interleaved image data to an external storage device;
本实施例中,如图3所示,对超声非空间复合图像数据进行交织处理的过程可以具体包括:In this embodiment, as shown in FIG. 3, the process of performing interleaving processing on the ultrasound non-spatial composite image data may specifically include:
步骤S2021:采用预设窗口对所述超声非空间复合图像数据进行分割。Step S2021: Segmenting the ultrasound non-spatial composite image data by using a preset window.
FPGA在处理大尺寸的图像数据时,如果一次性处理完则需要很大的存储空间。因此,在实际应用中由于存储空间和逻辑资源的限制,需要将大尺寸的图像数据分割成小尺寸的图像数据进行处理。因此,本实施例中以预设窗口对超声非空间复合图像数据进行分割。预设窗口的大小可以为64×64的窗口,也可以为128×128或者32×32,这均不影响本发明的实现。When processing large-size image data, the FPGA requires a large amount of storage space if it is processed once. Therefore, in practical applications, due to limitations of storage space and logical resources, it is necessary to divide large-sized image data into small-sized image data for processing. Therefore, in this embodiment, the ultrasound non-spatial composite image data is segmented by a preset window. The size of the preset window may be 64 x 64 windows, or 128 x 128 or 32 x 32, which does not affect the implementation of the present invention.
步骤S2022:对分割后的所述超声非空间复合图像数据进行流水线式的交织处理。Step S2022: Perform pipelined interleaving processing on the segmented ultrasound non-spatial composite image data.
FPGA内部的逻辑硬件电路将超声非空间复合图像数据分割为多个64×64窗口,并以此窗口为最小单位,将B图像DDR3缓存空间的数据搬运至逻辑内部Ram存储空间。逻辑硬件电路以窗口为单位,对逻辑内部缓存的数据进行流水线式的交织处理,并将交织后的窗口数据上传至外部缓存设备。The logic hardware circuit inside the FPGA divides the ultrasonic non-spatial composite image data into a plurality of 64×64 windows, and uses the window as a minimum unit to carry the data of the B image DDR3 buffer space to the logical internal Ram storage space. The logical hardware circuit performs pipelined interleaving processing on the logically buffered data in a window unit, and uploads the interleaved window data to an external cache device.
其中,在对缓存数据进行交织处理时,采用双片Ram进行乒乓操作,以实现逻辑内部缓存数据能够持续不断地流水线式处理。Among them, in the interleaving process of the buffer data, the ping-pong operation is performed by using the dual-chip Ram, so that the logical internal cache data can be continuously pipelined.
乒乓操作为应用于数据流控制的处理技巧,将数据流等时分配到两个数据缓存区,然后按节拍、相互配合的切换,从而将数据流处于一种连续不断的状态,以达到数据无缝缓冲和处理的效果。The ping-pong operation is applied to the processing technique of data flow control, and the data stream is equally allocated to two data buffers, and then switched according to the beat and the mutual cooperation, so that the data stream is in a continuous state to achieve data-free. The effect of seam buffering and processing.
流水线式处理的方式是高速设计中的一种技术手段,即将数据处理流程分为若干步骤,前一个步骤的输出是下一个步骤的输入。其最大的特点是数据流的每个步骤的处理从时间上来看是连续的,通过这样的设置能够进一步提高数据处理的效率。Pipelined processing is a technical means in high-speed design that divides the data processing flow into several steps, the output of the previous step being the input to the next step. The biggest feature is that the processing of each step of the data stream is continuous in time, and the efficiency of data processing can be further improved by such setting.
需要指出的是,本实施例中以64×64为窗口作为最小处理单元,该方式是本发明结合Xilinx tdpram资源特性的一种优选方案,还可以将交织处
理的窗口大小修改为128×128或者32×32等,即与本发明中优选方案设置窗口大小不同,这均不影响本发明的实现。It should be noted that in this embodiment, a window of 64×64 is used as a minimum processing unit, which is a preferred solution of the present invention in combination with the characteristics of the Xilinx tdpram resource, and the interlace can also be used.
The window size of the control is modified to 128×128 or 32×32, etc., that is, the window size is different from the preferred scheme in the present invention, which does not affect the implementation of the present invention.
步骤S203:从外部存储设备中获取交织处理后的数据,进行坐标校正,将校正后的数据存储至外部存储设备;Step S203: Acquiring the data after the interleaving process from the external storage device, performing coordinate correction, and storing the corrected data to the external storage device;
本实施例中,如图4所示,进行坐标校正的步骤可以具体为:In this embodiment, as shown in FIG. 4, the step of performing coordinate correction may be specifically as follows:
步骤S2031:按照预设方向从外部存储设备获取交织处理后的图像数据。Step S2031: Acquire interlaced image data from an external storage device according to a preset direction.
预设方向可以具体为:按照超声回波信号的深度方向,即,超声探头(probe)往人体测量的方向。逻辑硬件电路以同一深度不同线为单位,从外部存储设备中读取交织后的数据。同样地,调度的方式也可以采用队列流控的方式来实现,以提高数据处理的效率。The preset direction may be specifically: according to the depth direction of the ultrasonic echo signal, that is, the direction in which the ultrasonic probe is measured by the human body. The logic hardware circuit reads the interleaved data from an external storage device in units of different lines of the same depth. Similarly, the scheduling method can also be implemented by using queue flow control to improve the efficiency of data processing.
步骤S2032:从外部存储设备读取坐标校正所需的参数信息。Step S2032: Read parameter information required for coordinate correction from an external storage device.
需要指出的是,步骤S2031以及步骤S2032之间可以并行进行处理,并没有顺序的限制。It should be noted that the processing may be performed in parallel between step S2031 and step S2032, and there is no order limitation.
步骤S2033:当交织后的图像数据以及参数信息满足预设阈值条件时,对交织处理后的超声非空间复合图像数据进行坐标校正。Step S2033: Perform coordinate correction on the interlaced ultrasonic non-spatial composite image data when the interleaved image data and the parameter information satisfy a preset threshold condition.
本实施例在数据信号处理之前设置阈值条件,预设阈值条件为坐标校正所需的参数信息和交织处理后的图像数据二者均到达。In this embodiment, a threshold condition is set before the data signal processing, and the preset threshold condition is that both the parameter information required for the coordinate correction and the image data after the interleaving are reached.
预设阈值条件具体可以为判断交织缓存区数据以及参数对应的队列的状态是否均不为空,即相应队列中是否存在数据和参数。只有在交织缓存区数据以及参数信息均到达,判定满足预设阈值条件时,才能够执行后续操作。在不满足预设阈值条件时,则返回执行任务调度的过程,直至满足预设阈值条件为止。The preset threshold condition may specifically determine whether the interleaved buffer area data and the status of the queue corresponding to the parameter are not empty, that is, whether data and parameters exist in the corresponding queue. The subsequent operations can be performed only when the interleaving buffer data and the parameter information arrive and the determination meets the preset threshold condition. When the preset threshold condition is not met, the process of performing task scheduling is returned until the preset threshold condition is met.
坐标校正包括:配准处理以及边界处理两个过程。Coordinate correction includes two processes: registration processing and boundary processing.
其中,配准处理可以采用线性插值算法对图像数据进行配准。在本实施例中,采用Catmull-ROM线性插值算法。Catmull-ROM插值算法进行坐标矫正,首先在纵向方向(深度方向)进行线性插值,然后再横向(线方向)进行插值计算。当然还可以采用其他插值算法,在此不一一赘述。Among them, the registration process can use a linear interpolation algorithm to register the image data. In this embodiment, a Catmull-ROM linear interpolation algorithm is employed. The Catmull-ROM interpolation algorithm performs coordinate correction by first performing linear interpolation in the longitudinal direction (depth direction) and then interpolating in the lateral direction (line direction). Of course, other interpolation algorithms can also be used, which will not be repeated here.
边界处理包括:对配准处理后的图像数据中像素点的坐标进行运算。
根据超声非空间复合图像的预设边界对像素点进行处理,去除超出预设边界的像素点,未被预设边界覆盖的区域的像素点补零。预设边界为原始采集到的超声非空间复合图像对应的边界尺寸范围。The boundary processing includes: calculating coordinates of pixel points in the image data after registration processing.
The pixel points are processed according to the preset boundary of the ultrasonic non-spatial composite image, and the pixel points beyond the preset boundary are removed, and the pixel points of the area not covered by the preset boundary are zero-padded. The preset boundary is the boundary size range corresponding to the originally acquired ultrasonic non-spatial composite image.
本实施例中,坐标校正所需的参数信息具体包括配准过程中线性插值算法所需的参数以及边界处理过程中的预设边界的尺寸参数。In this embodiment, the parameter information required for coordinate correction specifically includes parameters required for the linear interpolation algorithm in the registration process and size parameters of the preset boundary in the boundary processing process.
步骤S204:从外部存储设备读取坐标校正后的图像数据,对坐标校正后的图像数据进行解交织处理,并将解交织处理后的数据存储至外部存储设备。Step S204: The coordinate-corrected image data is read from the external storage device, the coordinate-corrected image data is deinterleaved, and the de-interleaved data is stored to the external storage device.
步骤S205:从外部存储设备获取解交织处理后的数据,将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像。Step S205: Acquire deinterleaved data from the external storage device, and superimpose the pixels corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving to generate a spatial composite image.
通过对不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,即进行归一化处理,实现了对超声非空间复合图像叠加的区域进行弱化的目的。By superimposing the pixels corresponding to the ultrasonic non-spatial composite image data with different deflection angles, the normalization processing is performed, and the superimposed region of the ultrasonic non-spatial composite image is weakened.
本申请实施例采用FPGA实现空间复合算法,采用片外高速DDR存储器进行数据缓存,FPGA通过高速访问接口实现数据调度和复合运算。本申请能够有效提高空间复合运算的速度,大大降低了图像的处理时间,提高了成像帧率,在不更换CPU的情况下有效提升了空间复合成像的效果。另外,本实施例采用外部存储设备缓存交织、坐标校正以及解交织图像处理结果,将图像数据进行不同程度上的调度,有效地提高了FPGA图像处理的规模和效率。The embodiment of the present application implements a spatial composite algorithm by using an FPGA, and uses an off-chip high-speed DDR memory for data buffering, and the FPGA implements data scheduling and composite operations through a high-speed access interface. The application can effectively improve the speed of the spatial composite operation, greatly reduce the processing time of the image, improve the imaging frame rate, and effectively improve the effect of the spatial composite imaging without replacing the CPU. In addition, the embodiment uses the external storage device buffer interleaving, coordinate correction, and deinterleaving image processing results, and the image data is scheduled to different degrees, thereby effectively improving the scale and efficiency of the FPGA image processing.
在上述任一实施例的基础上,在获取不同偏转角度的超声非空间复合图像数据的步骤之后,在对交织处理后的所述超声非空间复合图像数据进行坐标校正的步骤之前还可以包括:On the basis of any of the above embodiments, after the step of acquiring the ultrasonic non-spatial composite image data of the different deflection angles, the step of performing the coordinate correction on the ultrasonic non-spatial composite image data after the interleaving process may further include:
对超声非空间复合图像数据进行线平滑处理和/或边缘弱化处理的过程。A process of performing line smoothing and/or edge weakening processing on ultrasonic non-spatial composite image data.
如图5本发明所提供的基于FPGA的空间复合成像方法的又一种具体
实施方式示意图所示,本发明实施例采用片外高速DDR存储器进行数据缓存,FPGA通过高速访问接口,采用流水的任务控制调度方式来实现数据调度和复合运算。该过程可以具体包括:FIG. 5 is another specific example of the FPGA-based spatial composite imaging method provided by the present invention.
As shown in the schematic diagram of the embodiment, the embodiment of the present invention uses an off-chip high-speed DDR memory for data caching, and the FPGA uses a high-speed access interface to implement a data scheduling and a compound operation by using a pipelined task control scheduling method. The process can specifically include:
步骤S301:空间复合开始后,读取DDR3中超声非空间复合图像数据进行交织处理。Step S301: After the spatial recombination starts, the ultrasonic non-spatial composite image data in the DDR3 is read for interleaving processing.
步骤S302:从DDR3中读取交织缓存区数据和处理所需参数。Step S302: Read the interleave buffer data and the required parameters from the DDR3.
本实施例中,处理所需的参数包括线平滑、边缘弱化和配准等数据处理过程中所需要的参数以及系数。In this embodiment, the parameters required for processing include parameters and coefficients required in data processing such as line smoothing, edge weakening, and registration.
步骤S303:判断交织缓存区数据以及数据处理过程中的参数信息是否满足预设阈值条件,即判断数据和参数是否已经全部到达;如果是,则执行步骤S304;如果否,则返回步骤S302。Step S303: determining whether the interleave buffer data and the parameter information in the data processing process meet the preset threshold condition, that is, whether the data and the parameters have all arrived; if yes, proceed to step S304; if not, return to step S302.
读取交织后所有线且同一深度的线数据时,逻辑硬件电路以同一深度不同线为单位,从缓存DDR3空间搬运数据、数字信号处理。其调度方式采用队列流控的方式来实现,并在数字信号处理前端设置阈值判断,当缓存数据事件队列和缓存参数队列的状态同时不为空时,即在缓存数据和缓存参数的队列容器中同时存在数据和参数时,则认为满足阈值判断条件,才开始线平滑、边缘弱化和配准处理。When reading the line data of all lines and the same depth after interleaving, the logic hardware circuit carries data and digital signal processing from the buffer DDR3 space in units of different lines of the same depth. The scheduling mode is implemented by using queue flow control, and threshold value judgment is set in the front end of the digital signal processing. When the state of the cached data event queue and the cache parameter queue are not empty at the same time, that is, in the queue container of the cached data and the cache parameter. When there are data and parameters at the same time, it is considered that the threshold judgment condition is satisfied, and line smoothing, edge weakening, and registration processing are started.
本实施例中深度方向指垂直方向,即探头往人体测量的方向,水平方向为线扫查的方法。In this embodiment, the depth direction refers to the vertical direction, that is, the direction in which the probe is measured to the human body, and the horizontal direction is the method of line scanning.
步骤S304:当满足阈值判断条件,则对交织缓存区数据进行线平滑处理,以降低超声图像近场锯齿状以及远场马赛克。Step S304: When the threshold judgment condition is satisfied, the interlaced buffer area data is subjected to line smoothing processing to reduce the near field jagged shape and the far field mosaic of the ultrasonic image.
对同一深度不同线的线数据进行线平滑处理,可以具体采用分段式处理,包括线边缘处保持原值、中间段线数据FIR滤波,其中,设总线数为N,当前线号为n,则线边缘处定义为n≤2,n≥N-2,中间段定义为3≤n≤n-3。逻辑硬件电路通过改变线边缘处加权系数,对线边缘处进行特殊处理,即保持原始数值不变;中间段线平滑处理,采用三阶FIR滤波的方式,流水式输出平滑后的结果。For line smoothing of line data of different lines at the same depth, segmentation processing may be specifically adopted, including maintaining the original value and the intermediate line data FIR filtering at the line edge, wherein the number of buses is N, and the current line number is n, Then the line edge is defined as n ≤ 2, n ≥ N-2, and the middle section is defined as 3 ≤ n ≤ n-3. The logic hardware circuit performs special processing on the edge of the line by changing the weighting coefficient at the edge of the line, that is, keeping the original value unchanged; the middle line is smoothed, and the third-order FIR filtering is used, and the result of the smoothed output is smoothed.
本实施例采用三阶FIR滤波的方式,满足了对超声图像进行平滑处理的需求,与五阶滤波、七阶滤波等方式相比,更容易实现。通过对超声图
像进行线平滑处理,能够降低近场锯齿状以及远场马赛克现象。In this embodiment, the third-order FIR filtering method is adopted, which satisfies the requirement of smoothing the ultrasound image, and is easier to implement than the fifth-order filtering and the seventh-order filtering. Passing on ultrasound
Like line smoothing, it can reduce near-field jagged and far-field mosaic.
步骤S305:对交织缓存区平滑处理后的数据进行边缘弱化处理,以降低超声图像融合后的伪影。Step S305: Perform edge weakening processing on the data after the interleaving buffer smoothing to reduce artifacts after the ultrasound image is merged.
步骤S306:对交织缓存区边缘弱化处理后的数据进行配准和边界处理,以对不同偏转角度的图像进行坐标校正。Step S306: Perform registration and boundary processing on the data of the interleaved buffer edge weakening process to perform coordinate correction on the images of different deflection angles.
本发明实施例采用Catmull-ROM插值算法进行坐标矫正,首先在纵向方向进行线性插值,然后在横向(线方向)进行插值计算。The embodiment of the present invention adopts the Catmull-ROM interpolation algorithm for coordinate correction, first performing linear interpolation in the longitudinal direction, and then performing interpolation calculation in the lateral direction (line direction).
对同一深度上的线数据进行配准处理,包括Catmull-ROM线性插值算法的FPGA实现、Catmull-ROM算法所需的插值系数与回波数据之间的同步运算以及双路并行配准处理单元。其中,双路并行配准处理单元按照Catmull-ROM插值规律,采用2×4运算窗口以4bit移位寄存器方式,将数据与插值系数进行乘法累加,每4个是时钟周期输出一个配准处理的结果。Registration processing of line data at the same depth, including FPGA implementation of Catmull-ROM linear interpolation algorithm, synchronization operation between interpolation coefficient and echo data required by Catmull-ROM algorithm, and dual parallel registration processing unit. The two-way parallel registration processing unit multiplies the data and the interpolation coefficient by a 2×4 operation window in a 4×4 shift window according to the Catmull-ROM interpolation rule, and outputs a registration process every four clock cycles. result.
配准处理可具体采用双路并行配准处理单元,该方式是本发明结合超声系统FPGA逻辑资源做出的优选方案,还可以有其他方法。例如:可以将配准处理单元拓展到四通道、六通道和八通道,以提高配准处理的速率,即与本发明中优选方案配准处理单元通路个数不是同一值,这均不影响本发明的实现。The registration processing may specifically adopt a two-way parallel registration processing unit, which is a preferred scheme made by the present invention in combination with the FPGA system logic resources of the ultrasound system, and may have other methods. For example, the registration processing unit can be extended to four channels, six channels, and eight channels to increase the rate of registration processing, that is, the number of paths of the processing unit to be registered in the preferred scheme of the present invention is not the same value, which does not affect the present invention. Implementation of the invention.
对同一深度上的线数据进行边界处理时,逻辑硬件电路对配准插值后的图像数据中像素点的坐标进行运算,然后根据超声B回波图像的边界对配准后的像素点进行约束,将超出B图像边界的像素点抛弃,未被配准覆盖的区域补充零。When performing boundary processing on the line data at the same depth, the logic hardware circuit calculates the coordinates of the pixel points in the image data after the registration interpolation, and then constrains the registered pixel points according to the boundary of the ultrasonic B echo image. Discard pixels that are beyond the boundary of the B image, and zeros are added to areas that are not registered.
步骤S307:对交织缓存区数据进行解交织处理,将图像还原为原始图像规则和尺寸。Step S307: Deinterleave the interleave buffer data to restore the image to the original image rule and size.
步骤S308:对解交织处理后的数据进行归一化处理,将不同偏转角度的超声图像进行融合。Step S308: Normalize the data after the deinterleaving process, and fuse the ultrasound images with different deflection angles.
在上述任一实施例的基础上,本发明实施例在对不同偏转角度的超声回波信号数据进行存储之前,还包括:对超声回波信号进行信号放大、模/数转换处理、波束合成以及数字信号处理。On the basis of any of the above embodiments, before storing the ultrasonic echo signal data of different deflection angles, the embodiment of the present invention further includes: performing signal amplification, analog/digital conversion processing, beam synthesis, and the ultrasonic echo signal. Digital signal processing.
需要指出的是,本实施例中数据流处理按照交织、线平滑、边缘弱化、
配准、边界处理、解交织和归一化先后顺序进行,该方式是本发明结合算法进行实现的优选方案,当然还可以有其他方法。例如:可以将线平滑与边缘弱化的位置进行互换,或则将线平滑和边缘弱化的位置放在交织前完成,即本发明的数据流处理顺序可以进行调换。It should be noted that, in this embodiment, the data stream processing is performed according to interleaving, line smoothing, edge weakening,
The registration, the boundary processing, the de-interlacing, and the normalization are sequentially performed. This method is a preferred solution implemented by the present invention in combination with an algorithm, and of course, there may be other methods. For example, the position of the line smoothing and the edge weakening may be interchanged, or the position of the line smoothing and the edge weakening may be completed before the interleaving, that is, the data stream processing sequence of the present invention may be exchanged.
图6为空间复合成像FPGA实现的整体框架示意图,该系统上游为用户配置接口以及DDR3高速读取接口,系统下游则为DDR3高速存储接口。而整个系统数据流从左至右,首先交织处理模块(INLE)从DDR3中获取不同偏转角度的超声回波信号,进行交织处理;然后信号处理模块(SCON_PRO)对交织后的超声回波信号进行边缘弱化、线平滑、配准和边界等处理;而后解交织模块(INV_INLE)将信号处理后的回波信号还原为通常格式上的超声信号;最后通过融合(NOR)模块将不同偏转角度的回波信进行归一化,得到空间复合图像。Figure 6 is a schematic diagram of the overall framework of the space composite imaging FPGA implementation. The upstream of the system is the user configuration interface and the DDR3 high-speed read interface, and the downstream of the system is the DDR3 high-speed storage interface. The entire system data stream is from left to right. First, the interleaving processing module (INLE) acquires ultrasonic echo signals of different deflection angles from DDR3 for interleaving processing; then the signal processing module (SCON_PRO) performs the interleaved ultrasonic echo signals. Edge weakening, line smoothing, registration, and boundary processing; then the de-interleaving module (INV_INLE) restores the signal-processed echo signal to the ultrasound signal in the usual format; finally, the different deflection angles are returned by the fusion (NOR) module. Boss is normalized to obtain a spatial composite image.
图7为空间复合成像FPGA实现的数据处理流程图,图中标识了数据处理环节输入/出的信号名称、符号、位宽和精度,以sample_shift s15.5为例,其中sampe_shift为信号名称、s代表signed有符号数、15代表整数部分位宽为15,.5则代表小数部分位宽为5,而信号整体位宽则为20。图中所示,DownSizer主要用于数字信号的位宽截断,UpSizer则用于数字信号的位宽拓展,Edge_weak为边缘弱化处理模块,Scon_Cal则为配准坐标变换模块。Figure 7 is a flow chart of data processing implemented by a spatial composite imaging FPGA. The figure identifies the signal name, symbol, bit width and precision of the input/output of the data processing link. Take sample_shift s15.5 as an example, where sampe_shift is the signal name and s. Represents a signed signed number, 15 represents an integer part with a bit width of 15, and .5 represents a fractional part with a bit width of 5 and a signal overall bit width of 20. As shown in the figure, DownSizer is mainly used for bit width truncation of digital signals, UpSizer is used for bit width expansion of digital signals, Edge_weak is edge weakening processing module, and Scon_Cal is registration coordinate transformation module.
本实施例通过对超声回波信号进行交织预处理,能够大大降低空间复合算法实现的复杂度,因此在实际成像帧率较低以及医学超声系统CPU性能不足的情况下,能够通过现场可编程逻辑阵列FPGA对超声非空间复合图像数据进行存储调度、硬件加速以及空间复合处理。In this embodiment, the interleaving preprocessing of the ultrasonic echo signals can greatly reduce the complexity of the spatial composite algorithm implementation, and thus can pass the field programmable logic when the actual imaging frame rate is low and the CPU performance of the medical ultrasound system is insufficient. The array FPGA performs storage scheduling, hardware acceleration, and spatial composite processing on the ultrasonic non-spatial composite image data.
下面对本发明实施例提供的基于FPGA的空间复合成像装置进行介绍,下文描述的基于FPGA的空间复合成像装置与上文描述的基于FPGA的空间复合成像方法可相互对应参照。The FPGA-based spatial composite imaging apparatus provided by the embodiment of the present invention is described below. The FPGA-based spatial composite imaging apparatus described below and the FPGA-based spatial composite imaging method described above can be mutually referenced.
图8为本发明实施例提供的基于FPGA的空间复合成像装置的结构框图,参照图8基于FPGA的空间复合成像装置可以包括:
FIG. 8 is a structural block diagram of an FPGA-based spatial composite imaging apparatus according to an embodiment of the present invention. The space-based composite imaging apparatus based on the FPGA may include:
获取模块100,用于获取不同偏转角度的超声非空间复合图像数据;The acquiring module 100 is configured to acquire ultrasonic non-spatial composite image data with different deflection angles;
交织模块200,用于对所述超声非空间复合图像数据进行交织处理;The interleaving module 200 is configured to perform interleaving processing on the ultrasound non-spatial composite image data;
坐标校正模块300,用于对交织处理后的超声非空间复合图像数据进行坐标校正,以将发生角度偏转的图像数据恢复至零角度偏转的图像数据;The coordinate correction module 300 is configured to perform coordinate correction on the interlaced ultrasonic non-spatial composite image data to restore the angularly deflected image data to the image data with zero angle deflection;
解交织模块400,用于对坐标校正后的图像数据进行解交织处理;The deinterleaving module 400 is configured to perform deinterleaving processing on the coordinate corrected image data;
复合模块500,用于将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像。The composite module 500 is configured to superimpose pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving to generate a spatial composite image.
本申请降低了空间复合算法实现的复杂度,能够有效提高空间复合运算速度,降低一帧图像的处理时间,提高成像帧率。The application reduces the complexity of the spatial composite algorithm implementation, can effectively improve the spatial composite operation speed, reduce the processing time of one frame of image, and improve the imaging frame rate.
综上,本申请采用现场可编程逻辑阵列FPGA取代传统MPU处理器,来实现医学超声图像业务上的空间复合功能。一方面给超声图像处理——空间复合,提供了多元化的实现方式;另一方面采用FPGA硬件加速的方式,能够有效地提高空间复合实现的速率、成像帧率。此外,本发明首先将超声非空间复合图像数据进行预交织处理,然后再进行线平滑、边缘弱化和配准等处理,很大程度上降低了空间复合FPGA实现的复杂度,节省了FPGA的逻辑资源。In summary, this application replaces the traditional MPU processor with a field programmable logic array FPGA to realize the spatial composite function on the medical ultrasound image service. On the one hand, it provides a variety of implementation methods for ultrasound image processing-space composite; on the other hand, FPGA hardware acceleration can effectively improve the rate and imaging frame rate of spatial composite implementation. In addition, the present invention first pre-interleaves the ultrasonic non-spatial composite image data, and then performs line smoothing, edge weakening, and registration processing, which greatly reduces the complexity of the spatial composite FPGA implementation and saves the logic of the FPGA. Resources.
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其它实施例的不同之处,各个实施例之间相同或相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。The various embodiments in the specification are described in a progressive manner, and each embodiment focuses on differences from other embodiments, and the same or similar parts of the respective embodiments may be referred to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.
专业人员还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。A person skilled in the art will further appreciate that the elements and algorithm steps of the various examples described in connection with the embodiments disclosed herein can be implemented in electronic hardware, computer software or a combination of both, in order to clearly illustrate the hardware and software. Interchangeability, the composition and steps of the various examples have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the solution. A person skilled in the art can use different methods for implementing the described functions for each particular application, but such implementation should not be considered to be beyond the scope of the present invention.
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬
件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。The steps of the method or algorithm described in connection with the embodiments disclosed herein may be directly applied to the hard
Implemented by a piece of software, a software module executed by a processor, or a combination of both. The software module can be placed in random access memory (RAM), memory, read only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or technical field. Any other form of storage medium known.
以上对本发明所提供的基于FPGA的空间复合成像方法以及装置进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以对本发明进行若干改进和修饰,这些改进和修饰也落入本发明权利要求的保护范围内。
The FPGA-based spatial composite imaging method and apparatus provided by the present invention are described in detail above. The principles and embodiments of the present invention have been described herein with reference to specific examples, and the description of the above embodiments is only to assist in understanding the method of the present invention and its core idea. It should be noted that those skilled in the art can make various modifications and changes to the present invention without departing from the spirit and scope of the invention.
Claims (10)
- 一种基于FPGA的空间复合成像方法,其特征在于,包括:An FPGA-based spatial composite imaging method, comprising:获取不同偏转角度的超声非空间复合图像数据;Acquiring ultrasonic non-spatial composite image data of different deflection angles;对所述超声非空间复合图像数据进行交织处理;Interleaving the ultrasound non-spatial composite image data;对交织处理后的超声非空间复合图像数据进行坐标校正,以将发生角度偏转的图像数据恢复至零角度偏转的图像数据;Performing coordinate correction on the interlaced ultrasonic non-spatial composite image data to restore the angularly deflected image data to the zero angle deflected image data;对坐标校正后的图像数据进行解交织处理;Performing deinterleaving processing on the coordinate corrected image data;将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像。The pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving are superimposed to generate a spatial composite image.
- 如权利要求1所述的基于FPGA的空间复合成像方法,其特征在于,所述对所述超声非空间复合图像数据进行交织处理的步骤包括:The FPGA-based spatial composite imaging method according to claim 1, wherein the step of performing interleaving processing on the ultrasound non-spatial composite image data comprises:采用预设窗口对所述超声非空间复合图像数据进行分割;Separating the ultrasound non-spatial composite image data by using a preset window;对分割后的所述超声非空间复合图像数据进行流水线式的交织处理;Performing pipelined interleaving processing on the segmented ultrasound non-spatial composite image data;将交织处理后的图像数据存储至外部存储设备。The interleaved image data is stored to an external storage device.
- 如权利要求1所述的基于FPGA的空间复合成像方法,其特征在于,所述对交织处理后的超声非空间复合图像数据进行坐标校正的步骤包括:The FPGA-based spatial composite imaging method according to claim 1, wherein the step of performing coordinate correction on the interlaced ultrasonic non-spatial composite image data comprises:按照预设方向从外部存储设备获取交织处理后的图像数据;Acquiring the interleaved image data from the external storage device according to a preset direction;从所述外部存储设备读取坐标校正所需的参数信息;Reading parameter information required for coordinate correction from the external storage device;当交织处理后的图像数据以及所述参数信息满足预设阈值条件时,对交织处理后的超声非空间复合图像数据进行坐标校正;Performing coordinate correction on the interlaced ultrasonic non-spatial composite image data when the interleaved image data and the parameter information satisfy a preset threshold condition;将坐标校正后的图像数据存储至所述外部存储设备。The coordinate corrected image data is stored to the external storage device.
- 如权利要求1所述的基于FPGA的空间复合成像方法,其特征在于,所述对交织处理后的超声非空间复合图像数据进行坐标校正的步骤包括:The FPGA-based spatial composite imaging method according to claim 1, wherein the step of performing coordinate correction on the interlaced ultrasonic non-spatial composite image data comprises:采用线性插值算法对交织处理后的图像数据进行配准;The interpolated image data is registered by a linear interpolation algorithm;对配准处理后的图像数据中像素点的坐标进行运算,根据超声非空间复合图像的预设边界对所述像素点进行处理,去除超出所述预设边界的像素点,未被所述预设边界覆盖的区域的像素点补零。Performing operations on the coordinates of the pixel points in the image data after the registration processing, processing the pixel points according to a preset boundary of the ultrasonic non-spatial composite image, and removing pixel points beyond the preset boundary, without being pre-predicted Set the pixel of the area covered by the boundary to zero.
- 如权利要求3所述的基于FPGA的空间复合成像方法,其特征在于,所述对坐标校正后的图像数据进行解交织处理的步骤包括: The FPGA-based spatial composite imaging method according to claim 3, wherein the step of deinterleaving the coordinate-corrected image data comprises:从所述外部存储设备读取坐标校正后的图像数据;Reading coordinate corrected image data from the external storage device;对坐标校正后的图像数据进行解交织处理;Performing deinterleaving processing on the coordinate corrected image data;将解交织处理后的图像数据存储至所述外部存储设备。The image data after deinterleaving is stored to the external storage device.
- 如权利要求5所述的基于FPGA的空间复合成像方法,其特征在于,所述将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像包括:The FPGA-based spatial composite imaging method according to claim 5, wherein the superimposing pixels corresponding to different non-spatial composite image data of different deflection angles are superimposed to generate a spatial composite image comprising:从所述外部存储设备获取解交织处理后的图像数据;Acquiring image data after deinterleaving processing from the external storage device;将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像。The pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving are superimposed to generate a spatial composite image.
- 如权利要求1至6任一项所述的基于FPGA的空间复合成像方法,其特征在于,在所述获取不同偏转角度的超声非空间复合图像数据的步骤之后,在所述对交织处理后的超声非空间复合图像数据进行坐标校正的步骤之前还包括:The FPGA-based spatial composite imaging method according to any one of claims 1 to 6, wherein after the step of acquiring ultrasonic non-spatial composite image data of different deflection angles, after the pair of interleaving processes Before the step of performing coordinate correction on the ultrasonic non-spatial composite image data, the method further includes:对所述超声非空间复合图像数据进行线平滑处理和/或边缘弱化处理。The ultrasonic non-spatial composite image data is subjected to line smoothing processing and/or edge weakening processing.
- 如权利要求7所述的基于FPGA的空间复合成像方法,其特征在于,对所述超声非空间复合图像数据进行线平滑处理的步骤包括:The FPGA-based spatial composite imaging method according to claim 7, wherein the step of performing line smoothing on the ultrasonic non-spatial composite image data comprises:采用分段式进行线平滑预处理,对线边缘处保持原始数值,对中间段线采用三阶FIR滤波的方式进行平滑处理;The segmentation method is used for line smoothing preprocessing, and the original value is maintained at the edge of the line, and the third stage FIR filter is used for smoothing the middle segment line;其中,设总线数为N,当前线号为n,则线边缘处定义为n≤2,n≥N-2,中间段线定位为3≤n≤n-3。Wherein, if the number of buses is N and the current line number is n, the edge of the line is defined as n≤2, n≥N-2, and the middle segment line is positioned as 3≤n≤n-3.
- 一种基于FPGA的空间复合成像装置,其特征在于,包括:An FPGA-based spatial composite imaging device, comprising:获取模块,用于获取不同偏转角度的超声非空间复合图像数据;Obtaining module for acquiring ultrasonic non-spatial composite image data with different deflection angles;交织模块,用于对所述超声非空间复合图像数据进行交织处理;An interleaving module, configured to perform interleaving processing on the ultrasound non-spatial composite image data;坐标校正模块,用于对交织处理后的超声非空间复合图像数据进行坐标校正,以将发生角度偏转的图像数据恢复至零角度偏转的图像数据;a coordinate correction module, configured to perform coordinate correction on the interlaced ultrasonic non-spatial composite image data to restore the angularly deflected image data to the image data with zero angle deflection;解交织模块,用于对坐标校正后的图像数据进行解交织处理;a deinterleaving module, configured to perform deinterleaving processing on the coordinate corrected image data;复合模块,用于将解交织后的不同偏转角度的超声非空间复合图像数据对应的像素点进行叠加,生成空间复合图像。The composite module is configured to superimpose pixel points corresponding to the ultrasonic non-spatial composite image data of different deflection angles after deinterleaving to generate a spatial composite image.
- 如权利要求9所述的基于FPGA的空间复合成像装置,其特征在 于,还包括:The FPGA-based spatial composite imaging apparatus according to claim 9, wherein It also includes:预处理模块,用于在获取不同偏转角度的超声非空间复合图像数据的步骤之后,在对交织处理后的超声非空间复合图像数据进行坐标校正的步骤之前,对所述超声非空间复合图像数据进行线平滑处理和/或边缘弱化处理。 a pre-processing module, configured to: after acquiring the ultrasonic non-spatial composite image data of different deflection angles, the ultrasonic non-spatial composite image data before the step of performing coordinate correction on the interlaced ultrasonic non-spatial composite image data Perform line smoothing and/or edge weakening.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611023419.4A CN106780349B (en) | 2016-11-21 | 2016-11-21 | Space composite imaging method and device based on FPGA |
CN201611023419.4 | 2016-11-21 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2018090654A1 true WO2018090654A1 (en) | 2018-05-24 |
Family
ID=58969829
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2017/092936 WO2018090654A1 (en) | 2016-11-21 | 2017-07-14 | Fpga-based space compound imaging method and device |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN106780349B (en) |
WO (1) | WO2018090654A1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113934393A (en) * | 2021-10-13 | 2022-01-14 | 北京青云航空仪表有限公司 | Image transformation method based on FPGA and aircraft cabin display system |
CN114414662A (en) * | 2021-12-21 | 2022-04-29 | 上海船舶工程质量检测有限公司 | System and method for ultrasonic plane scanning imaging |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106780349B (en) * | 2016-11-21 | 2020-02-14 | 深圳开立生物医疗科技股份有限公司 | Space composite imaging method and device based on FPGA |
CN107970043B (en) * | 2017-12-28 | 2021-01-19 | 深圳开立生物医疗科技股份有限公司 | Shear wave detection method and device |
CN110101411B (en) * | 2019-05-28 | 2020-11-10 | 飞依诺科技(苏州)有限公司 | Ultrasonic imaging space compounding method and system |
CN111110273A (en) * | 2019-12-31 | 2020-05-08 | 无锡祥生医疗科技股份有限公司 | Ultrasonic image processing method and device and storage medium |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101617946A (en) * | 2008-07-04 | 2010-01-06 | 深圳迈瑞生物医疗电子股份有限公司 | Spatial fusion method and spatial fusion device for ultrasonic imaging and ultrasonic imaging system |
CN102038519A (en) * | 2010-12-31 | 2011-05-04 | 飞依诺科技(苏州)有限公司 | Space compound imaging method in ultrasonic diagnosis |
CN102727255A (en) * | 2012-07-13 | 2012-10-17 | 深圳市理邦精密仪器股份有限公司 | Method and device for ultrasonic image space compound imaging |
US20130253325A1 (en) * | 2010-04-14 | 2013-09-26 | Josef R. Call | Systems and methods for improving ultrasound image quality by applying weighting factors |
CN106780349A (en) * | 2016-11-21 | 2017-05-31 | 深圳开立生物医疗科技股份有限公司 | Spatial compound imaging method and device based on FPGA |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101199430B (en) * | 2006-12-15 | 2011-12-28 | 深圳迈瑞生物医疗电子股份有限公司 | Spatial compound imaging method and equipment and ultrasonic imaging system thereof |
KR101374452B1 (en) * | 2008-10-30 | 2014-03-17 | 노키아 코포레이션 | Method and apparatus for interleaving a data block |
-
2016
- 2016-11-21 CN CN201611023419.4A patent/CN106780349B/en active Active
-
2017
- 2017-07-14 WO PCT/CN2017/092936 patent/WO2018090654A1/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101617946A (en) * | 2008-07-04 | 2010-01-06 | 深圳迈瑞生物医疗电子股份有限公司 | Spatial fusion method and spatial fusion device for ultrasonic imaging and ultrasonic imaging system |
US20130253325A1 (en) * | 2010-04-14 | 2013-09-26 | Josef R. Call | Systems and methods for improving ultrasound image quality by applying weighting factors |
CN102038519A (en) * | 2010-12-31 | 2011-05-04 | 飞依诺科技(苏州)有限公司 | Space compound imaging method in ultrasonic diagnosis |
CN102727255A (en) * | 2012-07-13 | 2012-10-17 | 深圳市理邦精密仪器股份有限公司 | Method and device for ultrasonic image space compound imaging |
CN106780349A (en) * | 2016-11-21 | 2017-05-31 | 深圳开立生物医疗科技股份有限公司 | Spatial compound imaging method and device based on FPGA |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113934393A (en) * | 2021-10-13 | 2022-01-14 | 北京青云航空仪表有限公司 | Image transformation method based on FPGA and aircraft cabin display system |
CN114414662A (en) * | 2021-12-21 | 2022-04-29 | 上海船舶工程质量检测有限公司 | System and method for ultrasonic plane scanning imaging |
Also Published As
Publication number | Publication date |
---|---|
CN106780349A (en) | 2017-05-31 |
CN106780349B (en) | 2020-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2018090654A1 (en) | Fpga-based space compound imaging method and device | |
US10729407B2 (en) | Ultrasonic diagnostic apparatus, image processing apparatus, and image processing method | |
CN102727255B (en) | Method and device for ultrasonic image space compound imaging | |
JP5002397B2 (en) | Ultrasonic diagnostic apparatus and program | |
KR101911734B1 (en) | Software-based ultrasound imaging system | |
JP2006325686A (en) | Ultrasonic diagnostic apparatus and its image processing method | |
JP5918200B2 (en) | Ultrasonic diagnostic equipment | |
KR101978728B1 (en) | Method and device for classifying medical ultrasound image based on deep learning using smart device | |
CN111839588A (en) | Ultrasonic image generation method, device, equipment and storage medium | |
KR20110054924A (en) | Ultrasound system and method for providing ultrasound spatial compound image | |
US8811774B1 (en) | Super resolution using an interpretive scaler | |
JPH09220228A (en) | Blood flow information image display method and ultrasonic diagnostic system | |
CN104586433B (en) | First-harmonic based on frequency conversion/harmonic wave merges the formation method combined with space compound | |
JP2008229226A (en) | Image processor | |
Biswas et al. | An FPGA-based architecture of DSC–SRI units specially for motion blind ultrasound systems | |
WO2013150917A1 (en) | Ultrasound diagnostic device and ultrasound image super-resolution image generation method | |
CN104546008B (en) | A kind of first-harmonic/harmonic wave merges the formation method combined with space compound | |
CN108389174B (en) | Ultrasonic imaging system and ultrasonic imaging method | |
JP2018023798A (en) | Ultrasound diagnostic apparatus, image processing apparatus, and image processing method | |
JP4879084B2 (en) | Image processing device | |
JP6045866B2 (en) | Ultrasonic image processing device | |
CN104921756A (en) | Improved digital scanning conversion method | |
JP4616711B2 (en) | Ultrasonic diagnostic equipment | |
JP5072585B2 (en) | Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program | |
JP2002113001A (en) | Ultrasonic diagnostic device, and ultrasonic signal processing 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: 17872469 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 14.10.2019) |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 17872469 Country of ref document: EP Kind code of ref document: A1 |