WO2014173068A1 - 一种超声影像分析系统及其分析方法 - Google Patents
一种超声影像分析系统及其分析方法 Download PDFInfo
- Publication number
- WO2014173068A1 WO2014173068A1 PCT/CN2013/084044 CN2013084044W WO2014173068A1 WO 2014173068 A1 WO2014173068 A1 WO 2014173068A1 CN 2013084044 W CN2013084044 W CN 2013084044W WO 2014173068 A1 WO2014173068 A1 WO 2014173068A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- image data
- dimensional
- dimensional ultrasound
- ultrasound image
- ultrasonic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Ceased
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/466—Displaying means of special interest adapted to display 3D data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/481—Diagnostic techniques involving the use of contrast agents, e.g. microbubbles introduced into the bloodstream
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/483—Diagnostic techniques involving the acquisition of a 3D volume of data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/56—Details of data transmission or power supply
Definitions
- the invention relates to the field of medical ultrasound image analysis, and in particular to an ultrasound image analysis system and an analysis method thereof.
- liver cancer is one of the most common malignant tumors.
- the annual death rate of liver cancer is 300,000, and the mortality rate of liver cancer accounts for the second place in malignant tumors.
- liver cancer Among the many therapies for treating liver cancer, surgical resection is still the first. 80% of patients with liver cancer have poor liver function and coagulation mechanism due to cirrhosis, or cannot be surgically removed due to constraints such as lumps and poor heart and kidney function. The remaining 20% of patients with liver cancer undergo resection without recurrence. The rate is still as high as 70%. Liver cancer cells are usually not sensitive to chemotherapy and radiation therapy. Therefore, various non-surgical interventional treatments have become an important means for clinical treatment of liver cancer.
- contrast-enhanced ultrasound is used to evaluate interventional outcomes.
- the main evaluation method is to measure the long diameter of the tumor area by pre- and post-operative ultrasound two-dimensional angiography. This measurement can not guarantee the consistency of the ultrasound section and position during the two measurements, especially in the multi-needle ablation of large tumors. When a simple long diameter measurement does not represent the entire ablation area.
- Most of the current three-dimensional ultrasound imaging techniques are used to display the individual states of the target area separately, or to use 4D imaging to display the dynamic blood flow perfusion pattern of the target area.
- An additional interventional outcome assessment program is based on CT/MRI equipment for tomographic imaging, but this protocol does not meet the requirements for immediate assessment after clinical treatment, and there is also a lack of preoperative and postoperative three-dimensional volumetric quantitative contrast for tumor interventional therapy. Moreover, most of the clinical observations of the three-dimensional angiography of a certain phase. Obviously, the accuracy of such observations is not high.
- Embodiments of the present invention provide an ultrasound image analysis system and an analysis method, which can directly and clearly evaluate a patient's treatment effect.
- an ultrasound image analysis system including:
- An image collection module configured to transmit an ultrasonic wave and receive an ultrasonic echo, and obtain a first set of three-dimensional ultrasonic image data according to the ultrasonic echo, and obtain the first set of three-dimensional ultrasonic image numbers After that, the ultrasonic wave is transmitted and the ultrasonic echo is received, and the second set of three-dimensional ultrasonic image data is obtained according to the ultrasonic echo;
- An image analysis module configured to perform an analytical registration on the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data, to obtain the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound images Registration mapping relationship between data;
- an image display module configured to display the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data according to the registration mapping relationship.
- the ultrasound image analysis system further includes an image storage module for storing the first set of three-dimensional ultrasound image data, the second set of three-dimensional ultrasound image data, and/or the registration mapping relationship.
- first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data acquired by the image collection module are three-dimensional ultrasound tissue image data and/or three-dimensional ultrasound contrast image data.
- the image analysis module includes: an automatic analysis sub-module, configured to perform automatic analysis and registration on the second set of three-dimensional ultrasonic image data and the first set of three-dimensional ultrasonic image data, and obtain the registration map relationship.
- the image analysis module includes: a semi-automatic analysis sub-module, configured to perform semi-automatic processing on the second set of three-dimensional ultrasonic image data and the first set of three-dimensional ultrasonic image data in combination with the human-machine interaction module and the image display module; The registration is analyzed and the registration mapping relationship is obtained.
- a semi-automatic analysis sub-module configured to perform semi-automatic processing on the second set of three-dimensional ultrasonic image data and the first set of three-dimensional ultrasonic image data in combination with the human-machine interaction module and the image display module; The registration is analyzed and the registration mapping relationship is obtained.
- the image display module further includes:
- An image reconstruction sub-module configured to receive a selection signal of the selected target area, and reconstruct the target area by using the first set of three-dimensional ultrasound image data and the set of three-dimensional ultrasound image data of the second set of three-dimensional ultrasound image data a first image; and reconstructing the first image with another set of three-dimensional ultrasound image data of the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data according to the registration mapping relationship Corresponding second image of the target area;
- the image display module displays the first image in a first area of the display and the second image in a second area of the display.
- first image and the second image are three-dimensional images and/or at least one slice image of the target area.
- An embodiment of the present invention further provides an ultrasound image analysis method, comprising: collecting a first set of three-dimensional ultrasound image data;
- the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data are displayed according to the registration mapping relationship.
- the ultrasound image analysis method further includes: storing the first set of three-dimensional ultrasound image data, the second set of three-dimensional ultrasound image data, and/or the registration mapping relationship.
- first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data are three-dimensional ultrasound tissue image data and/or three-dimensional ultrasound contrast image data.
- performing analytical registration on the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data comprises: performing the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data The registration is automatically analyzed and the registration mapping relationship is obtained.
- performing analytical registration on the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data comprises: performing the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data The registration is semi-automatically analyzed and the registration mapping relationship is obtained.
- the displaying the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data according to the registration mapping relationship comprises:
- the second image is displayed in a second area of the display.
- first image and the second image are three-dimensional images and/or at least one slice image of the target area.
- Embodiments of the present invention also provide an ultrasound image analysis method, including:
- the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data are displayed according to the registration mapping relationship.
- the displaying the first set of three-dimensional ultrasound image data and guiding the second set of three-dimensional ultrasound image data according to the first set of three-dimensional ultrasound image data comprises:
- the ultrasound is transmitted and received to obtain a sliced image, and the position of the ultrasound probe is changed until the same slice image as the at least one slice image is obtained, and the second set of three-dimensional ultrasound image data is collected.
- the displaying the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data according to the registration mapping relationship comprises:
- the second image is displayed in a second area of the display.
- first image and the second image are three-dimensional images and/or at least one slice image of the target area.
- Embodiments of the present invention provide an ultrasound image analysis system and an analysis method thereof. At least two imaging of the target area is performed using three-dimensional ultrasound imaging techniques. Get the target area Two sets of data for different phases. Based on the image registration technology, the correspondence between the two sets of data is established. The state of the target area at different phases is displayed based on the comparison of the registration results. At the same time, based on the mapping relationship between the two sets of data established by the registration, the operations of measurement, tracing, and segmentation performed on a certain set of data can be mapped to the corresponding positions of another set of data, thereby visually comparing the two sets of data. , evaluate the treatment effect.
- FIG. 1 is a schematic diagram of an embodiment of an ultrasound image analysis system according to an embodiment of the present invention.
- FIG 2 is an interface distribution diagram of an image display module of an embodiment of an ultrasonic image analysis system according to an embodiment of the present invention
- FIG. 3 is a first schematic flow chart of an ultrasonic image analysis method according to an embodiment of the present invention.
- FIG. 4 is a second simplified flow chart of an ultrasonic image analysis method according to an embodiment of the present invention.
- FIG. 5 is a third schematic flowchart of an ultrasonic image analysis method according to an embodiment of the invention.
- Embodiments of the present invention provide an ultrasound image analysis system and an analysis method thereof. At least two imaging of the target area is performed using three-dimensional ultrasound imaging techniques. Obtain three-dimensional ultrasound image data of different time phases in the target area. Based on image registration technology, the correspondence between two sets of three-dimensional ultrasound image data is established. The state of the target area at different times is displayed based on the comparison of the registration results. At the same time, based on the mapping relationship between the two sets of data established by the registration, the operations of measurement, tracing, and segmentation performed on a certain set of data can be mapped to the corresponding positions of another set of data, thereby visually comparing the two sets of data. , evaluate the treatment effect.
- an ultrasound image analysis system includes:
- the image collection module 101 is configured to receive and transmit an ultrasonic shape, and process the same to obtain three-dimensional ultrasound image data.
- the three-dimensional ultrasound image data acquired by the image collection module 101 may be three-dimensional ultrasound tissue image data or
- the three-dimensional ultrasound contrast image data may also include three-dimensional ultrasound tissue image data and three-dimensional ultrasound contrast image data.
- the three-dimensional ultrasound image data acquisition step at least two sets of three-dimensional ultrasound image data (i.e., the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data) of the scanned object are collected by the image collection module 101.
- the two sets of three-dimensional ultrasound image data sets can be timed at intervals, and other scans or tumor ablation operations can be performed during the time interval of the two sets. Try to ensure that the patient's position is consistent when collecting the two times, and try to ensure that the position of the probe is consistent with the direction of the probe.
- the image analysis module 102 is configured to perform analysis and registration on the first set of three-dimensional ultrasonic image data and the second set of three-dimensional ultrasonic image data, and obtain registration between the first set of three-dimensional ultrasonic image data and the second set of three-dimensional ultrasonic image data. Mapping relations.
- the analytical registration calculation can be based on three-dimensional ultrasound tissue image data, three-dimensional ultrasound contrast image data, or a combination of both.
- said image analysis module 102 comprises:
- the automatic analysis sub-module is configured to perform automatic analysis and registration on the first set of three-dimensional ultrasonic image data and the second set of three-dimensional ultrasonic image data, and obtain a registration mapping relationship.
- said image analysis module 102 comprises:
- the human-computer interaction module 105 is configured to process, by the user, the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data by using an external hardware device;
- the image display module 103 is configured to display the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data and the result processed by the human-machine interaction module 105 according to the registration mapping relationship.
- the ultrasound image analysis system further comprises:
- the image storage module 104 is configured to store a first set of three-dimensional ultrasound image data, a second set of three-dimensional ultrasound image data, and/or a registration mapping relationship obtained by the image analysis module 102.
- the image display module 103 comprises:
- An image reconstruction sub-module configured to receive a selection signal of the selected target area, and reconstruct a first image of the target area by using the first set of three-dimensional ultrasound image data and the set of three-dimensional ultrasound image data of the second set of three-dimensional ultrasound image data;
- the registration mapping relationship reconstructs a second image of the target area corresponding to the first image by using the first set of three-dimensional ultrasound image data and another set of three-dimensional ultrasound image data of the second set of three-dimensional ultrasound image data.
- the image display module 103 displays the first image in the first area of the display and the second image in the second area of the display.
- the image display module 103 may further include:
- An image measurement sub-module is configured to measure and edit the registration relationship generated by the first set of three-dimensional ultrasound image data, the second set of three-dimensional ultrasound image data, and the improved image analysis module 102.
- the image automatic registration algorithm mainly includes two parts, namely, the similarity measurement method between images and the mapping method between images. According to different mapping methods, it can be divided into rigid body transformation (rotation and translation, ie rigid body registration), radiation transformation (scaling, rotation and translation, ie radiation registration) and nonlinear transformation (to establish different mappings for different local images) , ie nonlinear registration).
- rigid body transformation rotation and translation, ie rigid body registration
- radiation transformation scaling, rotation and translation, ie radiation registration
- nonlinear transformation to establish different mappings for different local images
- the registration algorithm can be defined as rigid body registration, that is, including rotation and translation. If the depth settings are inconsistent when the two sets are set, that is, the size of the image pixels collected by the set is inconsistent, then the bilinear difference algorithm is used to scale the two sets of data to the same scale, and then the registration is performed according to the rigid body registration. .
- the image brightness of a group of three-dimensional ultrasound contrast image data is f (X, and the image brightness of another set of three-dimensional ultrasound contrast image data pixel ⁇ is g (1 ⁇ 4), and the mapping between two sets of three-dimensional ultrasound contrast image data can be Expressed as: Can define two at the same time
- the above similarity measure is a typical "minimum absolute difference sum", and similar definitions are: “minimum square error sum” and “maximum mutual ⁇ ..) off, and minimum absolute improvement based on ultrasonic noise Rayleigh distribution characteristics Difference algorithm.
- f (X and g in the measurement function Y can also be defined as the local gradient size of the corresponding data, local gray quotient, etc.
- the specific automatic registration algorithm has a lot of literature in the field of image processing, only given here. An easy to implement implementation.
- a semi-automatic registration (ie, interactive registration) implementation that is, an interactive interaction between the human-computer interaction module 105 and the image display module 103, may be employed.
- Registration Under the assumption of rigid body registration, the upper left corner 3 x 3 matrix of the mapping matrix A is a symmetric orthogonal matrix.
- the most direct solution is to use interactive operation to select 4 or more pairs of corresponding pairs on the two sets of three-dimensional ultrasound data, and to solve the problem by least squares.
- the method can be used to solve the optimal mapping matrix A; or one of the two sets of three-dimensional ultrasonic data can be selected separately, and the corresponding relationship of the two aspects can be established by means of fusion display, and then one corresponding point is selected in the two sets of data. Yes, the corresponding point pair is outside the two known planes.
- the above two interactive registration methods are two common solutions, and only an easy-to-implement embodiment is given here. It should be noted that the image registration operation can also be a combination of interactive registration and automatic registration algorithms.
- the registration result can be displayed based on the image display module 103 after the registration operation is completed. It should be noted that, as shown in FIG. 2, the screen area of the image display module 103 includes at least two display sub-windows for comparing two sets of data.
- the registration result displayed may be three-dimensional ultrasound tissue image data, three-dimensional ultrasound contrast image data, or a combination of both.
- the image display module 103 simultaneously includes an image multi-planar reconstruction related function, and combines the functions of the human-machine interaction module 105 to select a target area in the image and perform multi-planar reconstruction on the target.
- the reconstructed section may be a three-dimensional ultrasound tissue image or a three-dimensional ultrasound contrast image. Displaying two sets of reconstructed images in the display sub-window, for example: displaying a reconstructed image from the first set of three-dimensional ultrasound contrast image data in one of the display sub-windows For example, the image reconstructed from the second set of three-dimensional ultrasound contrast image data corresponding to the image is displayed in another display sub-window, and the pixels in the image displayed in the two display sub-windows are in one-to-one correspondence. It should be noted that the three-dimensional ultrasound tissue image data may be displayed in the two display sub-windows, or may be three-dimensional ultrasound contrast image data, and may also support switching between the two data.
- the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data displayed by the image display module 103 include, but are not limited to: a first set of three-dimensional ultrasound image data and a second set of three-dimensional ultrasound image data. A portion, a section, a plurality of sections, and a projection of the first set of three-dimensional ultrasound image data or the second set of three-dimensional ultrasound image data in a certain direction.
- the two display sub-windows can also support measurement functions, such as distance measurement, area measurement, hooking, tracing, etc.; two display sub-windows can also support measurement function linkage display, ie in which One display sub-window performs the measurement operation, and the same operation is displayed in the same position of the other display sub-window; at the same time, the de-linking operation function can be supported, and the measurement operation of one display sub-window is performed separately.
- the above display sub-window can display an image cut surface, or can display multiple cut surfaces at the same time, that is, images arranged in the same order in the two display sub-windows.
- the screen area of the image display module 103 allows the user to select different cut surfaces (reconstructed) for display and measurement operations.
- This embodiment also preferably provides the following ultrasonic image analysis methods. Further, the ultrasonic image analysis methods will be further described below.
- an ultrasonic image analysis method includes:
- the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data are collected.
- the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data are collected in the image collection module 101, and the specific collection method may be consistent with the above-mentioned three-dimensional ultrasound image data acquisition link, Let me repeat.
- the first set of three-dimensional ultrasonic image data and the second set of three-dimensional ultrasonic image data are analyzed and registered, and a registration mapping relationship between the first set of three-dimensional ultrasonic image data and the second set of three-dimensional ultrasonic image data is obtained.
- the specific three-dimensional ultrasonic image data analysis registration method may adopt the above-mentioned automatic registration method or semi-automatic registration method, or may simultaneously apply the combination of the automatic registration method and the semi-automatic registration method to the two groups collected.
- the three-dimensional ultrasound image data is analyzed and registered, and the specific registration is not repeated here.
- the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data are displayed in the image display module 103 described above in accordance with the registration mapping relationship.
- displaying the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data according to the registration mapping relationship may include the following steps:
- a selection signal for selecting a target area is received.
- the selection signal may be input by the user through the human-machine interaction module 105, or may be system default.
- Through the selection signal Determine the target area that the user is interested in.
- the target area may be part of a first set of three-dimensional ultrasound image data or a second set of three-dimensional ultrasound image data. That is, the target area may be selected in the first set of three-dimensional ultrasound image data, or may be selected in the second set of three-dimensional ultrasound image data.
- the second image is then displayed in the second area of the display.
- the first image and the second image herein may be a three-dimensional image of the target area and/or at least one cut surface image.
- an ultrasonic image analysis method includes:
- the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data are displayed according to the registration mapping relationship.
- displaying the first set of three-dimensional ultrasound image data and guiding the second set of three-dimensional ultrasound image data according to the first set of three-dimensional ultrasound image data may include the following steps: displaying at least one of the first set of three-dimensional ultrasound image data Cut surface image
- the ultrasound is transmitted and received to obtain a sliced image, and the position of the ultrasound probe is changed until the same slice image as the at least one slice image is obtained, and the second set of three-dimensional ultrasound image data is collected.
- the first set of three-dimensional ultrasound image data when collecting the first set of three-dimensional ultrasound image data, it generally includes the target area to be evaluated. Importing the collected first set of three-dimensional ultrasound image data to the image display module 103, displaying a certain display sub-window in the screen area of the image display module 103, and then guiding The second set of three-dimensional ultrasound image data is collected. For example, in one of the display sub-windows, one of the first set of three-dimensional ultrasound image data is displayed (the slice in a certain direction, including but not limited to the axial intermediate slice in the three-dimensional ultrasound image data), and the second three-dimensional ultrasound image Before the data is collected, the moving probe finds the same slice and then collects the second set of 3D ultrasound image data.
- the image analysis registration method in the image analysis module 102 is used to register two sets of three-dimensional data.
- the specific registration scheme may be automatic registration based on image processing algorithms or semi-automatically combined with manual interactive operations. Registration (ie interactive registration) mode.
- the method steps of displaying the first set of three-dimensional ultrasound image data and the second set of three-dimensional ultrasound image data according to the registration mapping relationship may be the same as in the foregoing second method.
- the third type is the third type.
- an ultrasonic image analysis method includes:
- Two sets of three-dimensional ultrasound image data are collected; two sets of three-dimensional ultrasound image data are collected in the image collection module 101 described above.
- the specific how to collect the three-dimensional ultrasound image data is consistent with the data acquisition of the three-dimensional ultrasound image described above, and details are not described herein again.
- the specific three-dimensional ultrasound image data analysis registration method can use the above automatic registration method or
- the semi-automatic registration method can also simultaneously analyze and register the two sets of three-dimensional ultrasonic image data collected by the combination of the automatic registration method and the semi-automatic registration method. The specific registration is not repeated here.
- a set of three-dimensional ultrasound image data reconstructs a section corresponding to a previously reconstructed set of three-dimensional ultrasound image data, and displays the section.
- the image reconstruction sub-module is used to reconstruct image data obtained by the target region at two sets of time, select a target region in the three-dimensional ultrasound image data, and perform multi-plane reconstruction on the target, and the reconstructed image data is corresponding of.
- the reconstructed section may be three-dimensional ultrasound tissue image data or three-dimensional ultrasound contrast image data. For example: displaying an image reconstructed from the first set of three-dimensional ultrasound image data in one of the display sub-windows, and displaying the reconstructed image in the second set of three-dimensional ultrasound image data corresponding thereto in another display sub-window, The pixels in the image displayed in the two windows are - corresponding.
- any of the above ultrasonic image analysis methods further includes:
- the three-dimensional ultrasound image data may be three-dimensional ultrasound tissue image data, or three-dimensional ultrasound contrast image data, or To simultaneously include three-dimensional ultrasound tissue image data and three-dimensional ultrasound contrast image data.
- performing image registration on the image comprises: automatically analyzing and registering the three-dimensional ultrasound image data, and generating a registration result.
- the specific automatic analysis registration method is consistent with the automatic registration method described above, and will not be described herein.
- analyzing the image for registration comprises: performing semi-automatic analysis registration on the three-dimensional ultrasound image data, and generating a registration result.
- the specific automatic analysis registration method is consistent with the semi-automatic registration method described above, and will not be described again.
- any one of the above ultrasonic image analysis methods performs multi-planar reconstruction on the three-dimensional ultrasound image data and the registration result.
- the specific multi-planar reconstruction function is consistent with the multi-planar reconstruction described above, and is not described here.
- any one of the above ultrasonic image analysis methods measures and edits the three-dimensional ultrasonic image data and the registration result.
- the specific measurement and editing functions are consistent with the measurement and editing described above, and are not described here.
- the three-dimensional ultrasound contrast imaging apparatus is used to image the clinical tumor intervention treatment before and after surgery, and the target area is imaged at least twice. Obtain two sets of data for different phases in the target area. ⁇ Use image processing technology to register the two sets of data and establish the correspondence between the two sets of data. The state of the target area at different times is displayed based on the comparison of the registration results. At the same time, the mapping relationship between the two sets of data based on the registration is mapped to the corresponding position display of another set of data, thereby visually comparing the two sets of data, and the method of length measurement, tumor segmentation, etc. can be used to evaluate whether the ablation effect area covers the whole The tumor area, which visually compares the two sets of data, evaluates the treatment effect. The ablation effect is assessed on-site and if the entire tumor area is not covered by ablation, needle ablation can be performed on site. Avoid based on
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Hematology (AREA)
- Computer Graphics (AREA)
- General Engineering & Computer Science (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
提供一种超声影像分析系统及其分析方法。该分析系统包括：图像采集模块（101），用于接收和发射超声波形，并对其进行处理获取两组三维超声图像数据；图像分析模块（102），用于对两组三维超声图像数据进行分析配准，获得两组三维超声图像数据之间的配准映射关系；人机交互模块（105），用于响应使用者通过外部硬件设备对两组三维超声图像数据进行的处理；图像显示模块（103），用于根据配准映射关系显示两组三维超声图像数据，以及人机交互模块（105）所处理的结果；图像存储模块（104），用于存储两组三维超声图像数据和配准映射关系。利用该超声影像分析系统，可以有效地对病人的治疗效果进行评估。
Description
一种超声影像分析系统及其分析方法 技术领域
本发明涉及医用超声图像分析领域， 尤其涉及一种超声影像分析系 统及其分析方法。
背景技术
联合国世界卫生组织下属的国际癌症研究机构（IARC )报告 2008 年 确诊癌症 1270 万， 760 万人死于癌症， 2030 年全世界确诊癌症病人将 有大约 2140 万， 1320 万人将死于癌症。 当前中国肿瘤病发率和死亡率 都呈加速发展态势。 肿瘤介入治疗已作为临床最主要治疗手段之一， 超 声介入治疗又是其中发展最块， 对肝癌、 肺癌、 曱状腺癌有着明确治疗 效果的手段。 在中国， 肝癌是最常见的恶性肿瘤之一， 每年肝癌死亡达 30 万人， 肝癌死亡率占恶性肿瘤的第二位。 在治肝癌的诸多疗法中， 目 前仍首推手术切除。 8 0 %的肝癌患者因为肝硬化导致肝功能和凝血机 制差， 或者由于肿块部位及心肾功能差等制约因素， 不能做手术切除， 剩下的 2 0 %的肝癌患者即使接受手术切除， 复发率仍高达 7 0 %。 而 肝癌细胞通常对化疗和放疗不敏感。 因此各种非手术的介入性治疗， 成 为了临床治疗肝癌的重要手段。
在超声引导肿瘤消融治疗过程中， 超声造影被用来对介入治疗结果 进行评估。 目前主要的评估手段是釆用对术前、 术后的超声二维造影， 测量肿瘤区域的长径， 该测量不能保证两次测量时超声切面和位置的一 致性， 尤其在大肿瘤多针消融时， 简单的长径测量不能代表整个消融区 域。 目前的三维超声造影技术， 大都用来单独显示目标区域的单独状态， 或则釆用 4D 成像， 显示目标区域的动态血流灌注形式。 另外的介入效 果评估方案是基于 CT/MRI 设备进行断层成像， 但是该方案不能满足临 床手术治疗后即时评估的要求，而且同样缺少针对肿瘤介入治疗的术前、 术后的三维体积定量对比功能。 而且， 临床大都单独观察某一时相的三 维造影结果。 显然， 这样的观察结果准确率不高。
不难看出， 现有技术还存在一定缺陷。
发明内容
本发明实施例提供一种超声影像分析系统及分析方法， 可以直接明 确地评估病人的治疗效果。
为实现上述发明目的， 本发明的实施例中提供了一种超声影像分析 系统， 包括：
图像釆集模块， 用于发射超声波并接收超声回波， 并根据超声回波 获得第一组三维超声图像数据， 以及在获得所述第一组三维超声图像数
据之后， 发射超声波并接收超声回波， 并根据超声回波获得第二组三维 超声图像数据；
图像分析模块， 用于对所述第一组三维超声图像数据和所述第二组 三维超声图像数据进行分析配准， 获得所述第一组三维超声图像数据和 所述第二组三维超声图像数据之间的配准映射关系；
图像显示模块， 用于根据所述配准映射关系显示所述第一组三维超 声图像数据和所述第二组三维超声图像数据。
进一步地， 该超声影像分析系统还包括图像存储模块， 用于存储所 述第一组三维超声图像数据、第二组三维超声图像数据和 /或所述配准映 射关系。
进一步地， 所述图像釆集模块所获取的所述第一组三维超声图像数 据和所述第二组三维超声图像数据为三维超声组织图像数据和 /或三维 超声造影图像数据。
进一步地， 所述图像分析模块包括： 自动分析子模块， 用于对所述 第二组三维超声图像数据和所述第一组三维超声图像数据进行自动分析 配准， 并获得所述配准映射关系。
进一步地， 所述图像分析模块包括： 半自动分析子模块， 用于结合 人机交互模块和所述图像显示模块对所述第二组三维超声图像数据与所 述第一组三维超声图像数据进行半自动分析配准， 并获得所述配准映射 关系。
进一步地， 所述图像显示模块还包括：
图像重建子模块， 用于接收选择目标区域的选择信号， 并用所述第 一组三维超声图像数据和所述第二组三维超声图像数据中的一组三维超 声图像数据重建出所述目标区域的第一图像； 以及根据所述配准映射关 系， 用所述第一组三维超声图像数据和所述第二组三维超声图像数据中 的另一组三维超声图像数据重建出与所述第一图像对应的所述目标区域 的第二图像；
所述图像显示模块在显示器的第一区域中显示所述第一图像， 以及 在显示器的第二区域中显示所述第二图像。
进一步地， 所述第一图像和所述第二图像是所述目标区域的三维图 像和 /或至少一个切面图像。
本发明的实施例中还提供了一种超声影像分析方法， 包括： 釆集第一组三维超声图像数据；
在釆集所述第一组三维超声图像数据之后， 釆集第二组三维超声图 像数据；
对所述第二组三维超声图像数据和所述第一组三维数据进行分析配 准， 获得所述第一组三维超声图像数据和所述第二组三维超声图像数据
之间的配准映射关系；
根据所述配准映射关系显示所述第一组三维超声图像数据和所述第 二组三维超声图像数据。
进一步地， 该超声影像分析方法， 还包括： 存储所述第一组三维超 声图像数据、 所述第二组三维超声图像数据和 /或所述配准映射关系。
进一步地， 所述第一组三维超声图像数据和所述第二组三维超声图 像数据为三维超声组织图像数据和 /或三维超声造影图像数据。
进一步地， 对所述第一组三维超声图像数据和所述第二组三维超声 图像数据进行分析配准包括： 对所述第一组三维超声图像数据和所述第 二组三维超声图像数据进行自动分析配准， 并获得所述配准映射关系。
进一步地， 对所述第一组三维超声图像数据和所述第二组三维超声 图像数据进行分析配准包括： 对所述第一组三维超声图像数据和所述第 二组三维超声图像数据进行半自动分析配准，并获得所述配准映射关系。
进一步地， 所述根据所述配准映射关系显示所述第一组三维超声图 像数据和所述第二组三维超声图像数据包括：
接收选择目标区域的选择信号；
用所述第一组三维超声图像数据和所述第二组三维超声图像数据中 的一组三维超声图像数据重建出所述目标区域的第一图像；
在显示器的第一区域中显示所述第一图像；
根据所述配准映射关系， 用所述第一组三维超声图像数据和所述第 二组三维超声图像数据中的另一组三维超声图像数据重建出与所述第一 图像对应的所述目标区域的第二图像；
在显示器的第二区域中显示所述第二图像。
进一步地， 所述第一图像和所述第二图像是所述目标区域的三维图 像和 /或至少一个切面图像。
本发明的实施例还提供了一种超声影像分析方法， 包括：
釆集第一组三维超声图像数据；
显示所述第一组三维超声图像数据并依据所述第一组三维超声图像 数据引导釆集第二组三维超声图像数据；
对所述第一组三维超声图像数据和所述第二组三维超声图像数据进 行分析配准， 获得所述第一组三维超声图像数据和所述第二组三维超声 图像数据之间的配准映射关系；
根据所述配准映射关系显示所述第一组三维超声图像数据和所述第 二组三维超声图像数据。
进一步地， 所述显示所述第一组三维超声图像数据并依据所述第一 组三维超声图像数据引导釆集第二组三维超声图像数据包括：
显示所述第一组三维超声图像数据的至少一个切面图像；
发射超声波并接收超声波以获得切面图像，并改变超声探头的位置， 直到获得与所述至少一个切面图像相同的切面图像时， 开始釆集所述第 二组三维超声图像数据。
进一步地， 所述根据所述配准映射关系显示所述第一组三维超声图 像数据和所述第二组三维超声图像数据包括：
接收选择目标区域的选择信号；
用所述第一组三维超声图像数据和所述第二组三维超声图像数据中 的一组三维超声图像数据重建出所述目标区域的第一图像；
在显示器的第一区域中显示所述第一图像；
根据所述配准映射关系， 用所述第一组三维超声图像数据和所述第 二组三维超声图像数据中的另一组三维超声图像数据重建出与所述第一 图像对应的所述目标区域的第二图像；
在显示器的第二区域中显示所述第二图像。
进一步地， 所述第一图像和所述第二图像是所述目标区域的三维图 像和 /或至少一个切面图像。
本发明的实施例提出了一种超声影像分析系统及其分析方法。 利用 三维超声造影成像技术， 对目标区域进行至少两次成像。 获得目标区域 不同时相两组数据。基于图像配准技术建立两组数据间的——对应关系。 基于配准结果对比显示目标区域在不同时相的状态。 同时基于配准建立 的两组数据间的映射关系， 可以把在某一组数据上进行的测量、 描迹、 分割等操作， 映射到另外一组数据的对应位置显示， 从而直观对比两组 数据， 对治疗效果进行评估。
附图说明
图 1 是本发明实施例提供的一种超声影像分析系统的实施例的模 拟图；
图 2 是本发明实施例提供的一种超声影像分析系统的实施例的图 像显示模块的界面分布图；
图 3 是发明实施例提供的一种超声影像分析方法的第一种简要流 程图；
图 4 是发明实施例提供的一种超声影像分析方法的第二种简要流 程图；
图 5 是发明实施例提供的一种超声影像分析方法的第三种简要流 程图。
具体实施方式
为使本发明实施例的目的、 技术方案和优点更加清楚， 下面将结合 本发明实施例， 对本发明实施例中的技术方案进行清楚、 完整地描述。 需要说明的是， 在附图或说明书中， 相似或相同的元件使用相同的附图
标记。
在临床介入治疗过程中， 超声被广泛的用来进行介入引导和术中观 察。 热消融是目前开展最广泛的肿瘤介入治疗技术， 而临床缺少有效的 对消融治疗效果进行即时评估的功能。
本发明的实施例提出了一种超声影像分析系统及其分析方法。 利用 三维超声造影成像技术， 对目标区域进行至少两次成像。 获得目标区域 不同时相两组三维超声图像数据。 基于图像配准技术建立两组三维超声 图像数据间的——对应关系。 基于配准结果对比显示目标区域在不同时 相的状态。 同时基于配准建立的两组数据间的映射关系， 可以把在某一 组数据上进行的测量、 描迹、 分割等操作， 映射到另外一组数据的对应 位置显示， 从而直观对比两组数据， 对治疗效果进行评估。
下面对本发明的实施例中提出的一种超声影像分析系统做进一步说 明：
如图 1 所示， 一种超声影像分析系统， 包括：
图像釆集模块 101 , 用于接收和发射超声波形， 并对其进行处理获 取三维超声图像数据； 作为优选， 所述图像釆集模块 101所获取的三维 超声图像数据可以是三维超声组织图像数据或者三维超声造影图像数 据， 也可以同时包含三维超声组织图像数据和三维超声造影图像数据。
在三维超声图像数据获取环节， 利用图像釆集模块 101釆集扫描对 象的至少两组三维超声图像数据 （即第一组三维超声图像数据和第二组 三维超声图像数据）。在釆集的过程中尽量利用较大的釆集角度和釆集范 围， 为后续的配准流程操作获取尽量丰富的数据信息。 两组三维超声图 像数据釆集可以有一定的时间间隔， 在两次釆集的时间间隔内可以进行 其它的扫描或者肿瘤消融手术等操作。 两次釆集时尽量保证病人体位一 致， 尽量保证探头位置和探头方向一致。
需要说明的是， 具体的造影成像以及三维超声造影成像原理是本技 术领域内熟知的技术， 在此不做进一步阐述说明。
图像分析模块 102 , 用于对第一组三维超声图像数据和第二组三维 超声图像数据进行分析配准， 并获得第一组三维超声图像数据和第二组 三维超声图像数据之间的配准映射关系。 分析配准计算可以是基于三维 超声组织图像数据、 三维超声造影图像数据或者同时基于两者的组合。
作为优选， 所述图像分析模块 102包括：
自动分析子模块， 用于对第一组三维超声图像数据和第二组三维超 声图像数据进行自动分析配准， 并获得配准映射关系。
作为优选， 所述图像分析模块 102包括：
半自动分析子模块，用于结合人机交互模块 105和图像显示模块 103 对第一组三维超声图像数据和第二组三维超声图像数据进行半自动分析
配准， 并获得配准映射关系。
人机交互模块 105 , 用于响应使用者通过外部硬件设备对第一组三 维超声图像数据和第二组三维超声图像数据进行的处理；
图像显示模块 103 , 用于根据配准映射关系显示第一组三维超声图 像数据和第二组三维超声图像数据以及所述人机交互模块 105所处理的 结果。
作为优选， 所述超声影像分析系统还包括：
图像存储模块 104, 用于存储第一组三维超声图像数据、 第二组三 维超声图像数据和 /或所述图像分析模块 102所获得的配准映射关系。
作为优选， 所述图像显示模块 103包括：
图像重建子模块， 用于接收选择目标区域的选择信号， 并用第一组 三维超声图像数据和第二组三维超声图像数据中的一组三维超声图像数 据重建出目标区域的第一图像； 以及根据配准映射关系， 用第一组三维 超声图像数据和第二组三维超声图像数据中的另一组三维超声图像数据 重建出与该第一图像对应的该目标区域的第二图像。
图像显示模块 103在显示器的第一区域中显示第一图像， 以及在显 示器的第二区域中显示第二图像。
作为优选， 所述图像显示模块 103还可以包括：
图像测量子模块， 用于对所述第一组三维超声图像数据、 第二组三 维超声图像数据和改进图像分析模块 102所生成的配准映射关系进行测 量和编辑。
图像自动配准算法主要包括两部分， 即图像间的相似性度量方法和 图像间的映射方法。根据映射方法不同，可以分为刚体变换 (旋转和平移， 即刚体配准）、 放射变换 (放缩、 旋转和平移， 即放射配准)和非线性变换 (对不同局部的图像建立不同的映射， 即非线性配准）。 结合超声三维造 影成像特点， 下面给出刚体配准和放射配准的两种实施方案， 需要说明 的是， 本实施例不限定于某一种配准算法。 如果两组造影数据是在同一 深度下釆集， 即釆集的像素是同一尺度， 配准算法可以限定为刚体配准， 即包括旋转和平移。 如果两次釆集时深度设定不一致， 即釆集到的图像 像素大小尺度不一致， 则釆用双线性差值算法把两组数据放缩到同一尺 度， 然后按照刚体配准进行配准计算。
假设某组三维超声造影图像数据中像素点为 的图像亮度为 f (X , 另一组三维超声造影图像数据像素点 ^ 的图像亮度为 g (¼) , 两组三维超声造影图像数据间映射可表示为：
同时可定义两
E二 /m - g
ί=ί…- Ν
上述相似性度量为典型的 "最小绝对差值和" ， 类似的定义方式还 有： "最小平方误差和"和"最大互相 \..)关"，以及基于超声波噪声 Rayleigh 分布特性改进的最小绝对差值算法。 同时度量函数中的 f (X 和 g (Y 还可以定义为对应数据的局部梯度大小、 局部灰度商等等形式。 具体的 自动配准算法在图像处理领域有很多文献资料， 在此只给出易于实现的 实施方案。
进行图像配准操作时， 除了可以釆用上述自动配准算法外还可以釆 用半自动配准 （即交互式配准） 实施方案， 即结合人机交互模块 105和 图像显示模块 103 , 进行交互式配准。 在假设刚体配准要求下， 映射矩 阵 A 的左上角 3 x 3 矩阵为对称正交矩阵。 求解两组三维超声数据间的 刚体配准方法有多种， 最直接方案是利用交互式操作， 分别在两组三维 超声数据上选择 4 组或者 4 组以上的对应点对，通过最小二乘拟合方法 求解最优映射矩阵 A; 也可以是在两组三维超声数据上分别选择一个切 面， 用融合显示的方式， 建立两个切面的——对应关系， 然后在两组数 据中选择一个对应点对， 该对应点对在已知的两个切面外。 上述两种交 互式配准方法是常见的两种方案， 在此只给出易于实现的实施方案。 需 要说明的是， 进行图像配准操作时还可以是交互式配准与自动配准算法 的组合应用。
配准操作完成后可以基于图像显示模块 103显示配准结果。 需要说 明的是， 如图 2 所示， 所述图像显示模块 103的屏幕区域至少包括两个 显示子窗口， 用于对比显示两组数据。 显示的配准结果可以是三维超声 组织图像数据、 三维超声造影图像数据或者是两者的组合。
图像显示模块 103同时包含图像多平面重建相关功能， 结合人机交 互模块 105功能， 在图像中选择一个目标区域， 并对目标进行多平面重 建。 重建出来的切面可以是三维超声组织图像也可以是三维超声造影图 像。 在所述显示子窗口中分别显示重建后的两组图像， 例如： 在其中一 个显示子窗口中显示从第一组三维超声造影图像数据中重建出来的图
像， 在另一个显示子窗口中显示与之对应的从第二组三维超声造影图像 数据中重建出来的图像， 这两个显示子窗口中显示的图像中的像素是一 一对应的。 需要说明的是， 所述两个显示子窗口中显示的可以是三维超 声组织图像数据也可以是三维超声造影图像数据， 也可以支持两种数据 间的切换。
需要说明的是， 所述图像显示模块 103所显示的第一组三维超声图 像数据和第二组三维超声图像数据包括但不限于： 第一组三维超声图像 数据和第二组三维超声图像数据的一部分、 一个切面、 多个切面和第一 组三维超声图像数据或第二组三维超声图像数据在某个方向的投影。
与此同时， 作为优选， 所述两个显示子窗口还可以支持测量功能， 例如距离测量、 面积测量、 勾边、 描迹等； 两个显示子窗口还可以支持 测量功能联动显示， 即在其中一个显示子窗口进行测量操作， 在另外一 个显示子窗口的相同位置显示同样操作； 同时也可以支持解除联动操作 功能， 单独对某一个显示子窗口进行测量操作。 上述显示子窗口可以显 示一个图像切面， 也可以同时显示多个切面， 即两个显示子窗口中按照 相同顺序排列的图像。 所述图像显示模块 103的屏幕区域以支持使用者 选择不同的切面 （重建后的） 进行显示和测量操作。
本实施例还优选提供如下几种超声影像分析方法， 下面将对这几种 超声影像分析方法做进一步阐述。
第一种：
如图 3 所示， 一种超声影像分析方法， 包括：
釆集第一组三维超声图像数据和第二组三维超声图像数据。 在上述 的图像釆集模块 101中釆集第一组三维超声图像数据和第二组三维超声 图像数据， 具体的釆集方法可以与上文所述的三维超声图像数据获取环 节一致， 在此不再赘述。
对所釆集的第一组三维超声图像数据和第二组三维超声图像数据进 行分析配准， 获得第一组三维超声图像数据和第二组三维超声图像数据 之间的配准映射关系。 具体的三维超声图像数据分析配准方法可以釆用 上述的自动配准方法或者是半自动配准方法， 也可以同时釆用自动配准 方法和半自动配准方法的组合应用对所釆集的两组三维超声图像数据进 行分析配准， 具体的如何配准， 在此不再赘述。
根据配准映射关系在上述的图像显示模块 103中显示第一组三维超 声图像数据和第二组三维超声图像数据。
本实施例中， 根据配准映射关系显示第一组三维超声图像数据和第 二组三维超声图像数据可以包括下列步骤：
首先， 接收选择目标区域的选择信号。 该选择信号可以是由用户通 过人机交互模块 105输入的， 也可以是系统默认的。 通过该选择信号可
以确定用户感兴趣的目标区域。 该目标区域可以是第一组三维超声图像 数据或者第二组三维超声图像数据中的一部分。 也就是说， 目标区域可 以在第一组三维超声图像数据中选择， 也可以在第二组三维超声图像数 据中选择。
然后， 用第一组三维超声图像数据和所述第二组三维超声图像数据 中的其中一组三维超声图像数据重建出该目标区域的第一图像， 并在显 示器的第一区域中显示该第一图像。
然后， 根据前述的配准映射关系， 用第一组三维超声图像数据和所 述第二组三维超声图像数据中的另一组三维超声图像数据重建出与该第 一图像对应的目标区域的第二图像。 容易理解， 由前文所述， 已经获得 了两组三维超声图像数据的配准映射关系， 则， 在其中一组三维超声图 像数据中选择了目标区域后， 根据该配准映射关系， 即可获得该目标区 域在另一组三维超声图像数据中的对应区域， 用另一组三维超声图像数 据重建出的该对应区域的图像， 即为前文所说的与第一图像对应的目标 区域的第二图像。
然后， 在显示器的第二区域中显示该第二图像。
容易理解， 本实施例中， 这里的第一图像和第二图像可以是目标区 域的三维图像和 /或至少一个切面图像。
第二种：
如图 4 所示， 一种超声影像分析方法， 包括：
釆集第一组三维超声图像数据；
显示第一组三维超声图像数据并依据第一组三维超声图像数据引导 釆集第二组三维超声图像数据；
对第一组三维超声图像数据和第二组三维超声图像数据进行分析配 准， 获得第一组三维超声图像数据和第二组三维超声图像数据之间的配 准映射关系；
根据配准映射关系显示第一组三维超声图像数据和第二组三维超声 图像数据。
本发明的实施例中， 显示第一组三维超声图像数据并依据第一组三 维超声图像数据引导釆集第二组三维超声图像数据可以包括下列步骤： 显示第一组三维超声图像数据的至少一个切面图像；
发射超声波并接收超声波以获得切面图像，并改变超声探头的位置， 直到获得与该至少一个切面图像相同的切面图像时， 开始釆集第二组三 维超声图像数据。
例如， 釆集第一组三维超声图像数据时， 一般要包含所要评估的目 标区域。 导入所釆集的第一组三维超声图像数据到所述图像显示模块 103 ,显示在图像显示模块 103屏幕区域的某一个显示子窗口，进而引导
第二组三维超声图像数据的釆集。 例如在其中一个显示子窗口显示第一 组三维超声图像数据的某一个切面（在一定方向上的切面，包括但不限于 是三维超声图像数据中轴向中间切面），在第二次三维超声图像数据釆集 之前， 移动探头找到相同的切面， 然后釆集第二组三维超声图像数据。 保证两次数据釆集时探头的位置和探头方向一致， 从而保证釆集的数据 基本代表人体的相同部位。 利用图像分析模块 102中图像分析配准的方 法， 对两组三维数据进行配准， 具体的配准方案， 可以是基于图像处理 算法的全自动配准， 也可以是结合手动交互式操作的半自动配准 （即交 互式配准）模式。
本实施例中， 根据配准映射关系显示第一组三维超声图像数据和第 二组三维超声图像数据的方法步骤可以与前述第二种方法中的相同。
第三种：
如图 5 所示， 一种超声影像分析方法， 包括：
釆集两组三维超声图像数据； 在上述的图像釆集模块 101中釆集两 组三维超声图像数据。 具体的如何釆集三维超声图像数据与上文所述的 三维超声图像数据获取环节一致， 在此不再赘述。
对所釆集的两组三维超声图像数据进行分析配准， 确定两组三维超 声图像数据的配准映射关系； 具体的三维超声图像数据分析配准方法可 以釆用上述的自动配准方法或者是半自动配准方法， 也可以同时釆用自 动配准方法和半自动配准方法的组合应用对所釆集的两组三维超声图像 数据进行分析配准， 具体的如何配准， 在此不再赘述。
在其中一组三维超声图像数据中选择目标区域， 按照选定方向重建 出所述目标区域的切面， 并显示该组三维超声图像数据的切面； 基于两 组图像的配准映射关系， 在另一组三维超声图像数据中重建出与之前重 建出的一组三维超声图像数据切面对应的切面， 并显示该切面。
运用所述图像重建子模块， 重建目标区域在两次釆集时候获得的图 像数据， 在三维超声图像数据中选择一个目标区域， 并对目标进行多平 面重建， 重建后的图像数据是——对应的。 重建出来的切面可以是三维 超声组织图像数据也可以是三维超声造影图像数据。 例如： 在其中一个 显示子窗口中显示从第一组三维超声图像数据中重建出来的图像， 在另 一个显示子窗口中显示与之对应的第二组三维超声图像数据中重建出来 的图像， 这两个窗口中显示的图像中的像素是——对应的。
作为优选， 上述任意一种超声影像分析方法， 均还包括：
存储所述三维超声图像数据； 存储所述配准结果。 在上述的图像存 储模块 104内存储所述三维超声图像数据以及所述配准结果。
作为优选， 上述任意一种超声影像分析方法， 所述三维超声图像数 据可以是三维超声组织图像数据， 或者是三维超声造影图像数据， 也可
以同时包含三维超声组织图像数据和三维超声造影图像数据。 作为优选， 上述任意一种超声影像分析方法， 对所述图像进行分析 配准包括： 对所述三维超声图像数据进行自动分析配准， 并生成配准结 果。 具体的自动分析配准方法与上文所述的自动配准方法一致， 在此不 再赘述。
作为优选， 上述任意一种超声影像分析方法， 对所述图像进行分析 配准包括： 对所述三维超声图像数据进行半自动分析配准， 并生成配准 结果。 具体的自动分析配准方法与上文所述的半自动配准方法一致， 在 此不再赘述。
作为优选， 上述任意一种超声影像分析方法， 对所述三维超声图像 数据和配准结果进行多平面重建。 具体的多平面重建功能与上文所述的 多平面重建一致， 在此不再赘述。
作为优选， 上述任意一种超声影像分析方法， 对所述三维超声图像 数据和所述配准结果进行测量和编辑。 具体的测量和编辑功能如上文所 述的测量和编辑一致， 在此不再赘述。
本发明的实施例中， 利用三维超声造影成像设备对临床肿瘤介入治 疗的术前、 术后分别进行成像， 对目标区域进行至少两次成像。 获得目 标区域不同时相两组数据。 釆用图像处理技术对两组数据进行配准， 建 立两组数据的——对应关系。 基于配准结果对比显示目标区域在不同时 相的状态。 同时基于配准建立的两组数据间的映射关系， 映射到另外一 组数据的对应位置显示， 从而直观对比两组数据， 可以釆用长度测量， 肿瘤分割等方法， 评估消融效果区域是否涵盖整个肿瘤区域， 从而直观 地对比两组数据， 对治疗效果进行评估。 现场对消融效果进行评估， 如 果消融没有涵盖整个肿瘤区域， 则可以现场进行补针消融。 避免基于
CT/MRI 等进行临床评估需要延后若干时间，并且给病人二次介入手术。
最后应说明的是： 以上实施例仅用以说明本发明的技术方案， 而非 对其限制； 尽管参照前述实施例对本发明进行了详细的说明， 本领域的 普通技术人员应当理解： 其依然可以对前述实施例所记载的技术方案进 行修改，或者对其中部分技术特征进行等同替换； 而这些修改或者替换， 并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范 围。
Claims
1、 一种超声影像分析系统， 其特征在于， 包括：
图像釆集模块（ 101 ), 用于发射超声波并接收超声回波， 并根据超 声回波获得第一组三维超声图像数据， 以及在获得所述第一组三维超声 图像数据之后， 发射超声波并接收超声回波， 并根据超声回波获得第二 组三维超声图像数据；
图像分析模块（ 102 ), 用于对所述第一组三维超声图像数据和所述 第二组三维超声图像数据进行分析配准， 获得所述第一组三维超声图像 数据和所述第二组三维超声图像数据之间的配准映射关系；
图像显示模块（ 103 ), 用于根据所述配准映射关系显示所述第一组 三维超声图像数据和所述第二组三维超声图像数据。
2、 根据权利要求 1 所述的超声影像分析系统， 其特征在于， 还包 括：
图像存储模块（ 104 ), 用于存储所述第一组三维超声图像数据、 第 二组三维超声图像数据和 /或所述配准映射关系。
3、 根据权利要求 1 所述的超声影像分析系统， 其特征在于， 所述 图像釆集模块（ 101 )所获取的所述第一组三维超声图像数据和所述第二 组三维超声图像数据为三维超声组织图像数据和 /或三维超声造影图像 数据。
4、 根据权利要求 1 所述的超声影像分析系统， 其特征在于， 所述 图像分析模块 （ 102 ) 包括：
自动分析子模块， 用于对所述第二组三维超声图像数据和所述第一 组三维超声图像数据进行自动分析配准， 并获得所述配准映射关系。
5、 根据权利要求 1 所述的超声影像分析系统， 其特征在于， 所述 图像分析模块 （ 102 ) 包括：
半自动分析子模块， 用于结合人机交互模块（ 105 )和所述图像显示 模块（ 103 )对所述第二组三维超声图像数据与所述第一组三维超声图像 数据进行半自动分析配准， 并获得所述配准映射关系。
6、 根据权利要求 1 所述的超声影像分析系统， 其特征在于， 所述 图像显示模块 （ 103 )还包括：
图像重建子模块， 用于接收选择目标区域的选择信号， 并用所述第 一组三维超声图像数据和所述第二组三维超声图像数据中的一组三维超 声图像数据重建出所述目标区域的第一图像； 以及根据所述配准映射关 系， 用所述第一组三维超声图像数据和所述第二组三维超声图像数据中 的另一组三维超声图像数据重建出与所述第一图像对应的所述目标区域 的第二图像；
所述图像显示模块（ 103 )在显示器的第一区域中显示所述第一图像，
以及在显示器的第二区域中显示所述第二图像。
7、 根据权利要求 6 所述的超声影像分析系统， 其特征在于： 所述 第一图像和所述第二图像是所述目标区域的三维图像和 /或至少一个切 面图像。
8、 一种超声影像分析方法， 其特征在于， 包括：
釆集第一组三维超声图像数据；
在釆集所述第一组三维超声图像数据之后， 釆集第二组三维超声图 像数据；
对所述第二组三维超声图像数据和所述第一组三维数据进行分析配 准， 获得所述第一组三维超声图像数据和所述第二组三维超声图像数据 之间的配准映射关系；
根据所述配准映射关系显示所述第一组三维超声图像数据和所述第 二组三维超声图像数据。
9、 根据权利要求 8 所述的超声影像分析方法， 其特征在于， 还包 括：
存储所述第一组三维超声图像数据、 所述第二组三维超声图像数据 和 /或所述配准映射关系。
10、 根据权利要求 8 所述的超声影像分析方法， 其特征在于， 所述 第一组三维超声图像数据和所述第二组三维超声图像数据为三维超声组 织图像数据和 /或三维超声造影图像数据。
11、 根据权利要求 8 所述的超声影像分析方法， 其特征在于， 对所 述第一组三维超声图像数据和所述第二组三维超声图像数据进行分析配 准包括：
对所述第一组三维超声图像数据和所述第二组三维超声图像数据进 行自动分析配准， 并获得所述配准映射关系。
12、 根据权利要求 8 所述的超声影像分析方法， 其特征在于， 对所 述第一组三维超声图像数据和所述第二组三维超声图像数据进行分析配 准包括：
对所述第一组三维超声图像数据和所述第二组三维超声图像数据进 行半自动分析配准， 并获得所述配准映射关系。
13、 根据权利要求 8 所述的超声影像分析方法， 其特征在于， 所述 根据所述配准映射关系显示所述第一组三维超声图像数据和所述第二组 三维超声图像数据包括：
接收选择目标区域的选择信号；
用所述第一组三维超声图像数据和所述第二组三维超声图像数据中 的一组三维超声图像数据重建出所述目标区域的第一图像；
在显示器的第一区域中显示所述第一图像；
根据所述配准映射关系， 用所述第一组三维超声图像数据和所述第 二组三维超声图像数据中的另一组三维超声图像数据重建出与所述第一 图像对应的所述目标区域的第二图像；
在显示器的第二区域中显示所述第二图像。
14、 根据权利要求 13 所述的超声影像分析方法， 其特征在于： 所 述第一图像和所述第二图像是所述目标区域的三维图像和 /或至少一个 切面图像。
15、 一种超声影像分析方法， 其特征在于， 包括：
釆集第一组三维超声图像数据；
显示所述第一组三维超声图像数据并依据所述第一组三维超声图像 数据引导釆集第二组三维超声图像数据；
对所述第一组三维超声图像数据和所述第二组三维超声图像数据进 行分析配准， 获得所述第一组三维超声图像数据和所述第二组三维超声 图像数据之间的配准映射关系；
根据所述配准映射关系显示所述第一组三维超声图像数据和所述第 二组三维超声图像数据。
16、 如权利要求 15所述的超声影像分析方法， 其特征在于， 所述显 示所述第一组三维超声图像数据并依据所述第一组三维超声图像数据引 导釆集第二组三维超声图像数据包括：
显示所述第一组三维超声图像数据的至少一个切面图像；
发射超声波并接收超声波以获得切面图像，并改变超声探头的位置， 直到获得与所述至少一个切面图像相同的切面图像时， 开始釆集所述第 二组三维超声图像数据。
17、 根据权利要求 15 所述的超声影像分析方法， 其特征在于， 所 述根据所述配准映射关系显示所述第一组三维超声图像数据和所述第二 组三维超声图像数据包括：
接收选择目标区域的选择信号；
用所述第一组三维超声图像数据和所述第二组三维超声图像数据中 的一组三维超声图像数据重建出所述目标区域的第一图像；
在显示器的第一区域中显示所述第一图像；
根据所述配准映射关系， 用所述第一组三维超声图像数据和所述第 二组三维超声图像数据中的另一组三维超声图像数据重建出与所述第一 图像对应的所述目标区域的第二图像；
在显示器的第二区域中显示所述第二图像。
18、 根据权利要求 17 所述的超声影像分析方法， 其特征在于： 所 述第一图像和所述第二图像是所述目标区域的三维图像和 /或至少一个 切面图像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US14/921,915 US11083436B2 (en) | 2013-04-25 | 2015-10-23 | Ultrasonic image analysis systems and analysis methods thereof |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310147583.6A CN104116523B (zh) | 2013-04-25 | 2013-04-25 | 一种超声影像分析系统及其分析方法 |
| CN201310147583.6 | 2013-04-25 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| US14/921,915 Continuation US11083436B2 (en) | 2013-04-25 | 2015-10-23 | Ultrasonic image analysis systems and analysis methods thereof |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2014173068A1 true WO2014173068A1 (zh) | 2014-10-30 |
Family
ID=51762310
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/CN2013/084044 Ceased WO2014173068A1 (zh) | 2013-04-25 | 2013-09-24 | 一种超声影像分析系统及其分析方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US11083436B2 (zh) |
| CN (1) | CN104116523B (zh) |
| WO (1) | WO2014173068A1 (zh) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11109843B2 (en) * | 2015-02-26 | 2021-09-07 | Hitachi, Ltd. | Ultrasonic image pickup device and image processing device |
| CN106157282A (zh) * | 2015-03-31 | 2016-11-23 | 深圳迈瑞生物医疗电子股份有限公司 | 图像处理系统及方法 |
| WO2016175758A2 (en) * | 2015-04-28 | 2016-11-03 | Analogic Corporation | Image guided steering of a transducer array and/or an instrument |
| WO2016192114A1 (zh) * | 2015-06-05 | 2016-12-08 | 深圳迈瑞生物医疗电子股份有限公司 | 超声流体成像方法及超声流体成像系统 |
| CN106934807B (zh) * | 2015-12-31 | 2022-03-01 | 深圳迈瑞生物医疗电子股份有限公司 | 一种医学影像分析方法、系统及医疗设备 |
| CN106971055B (zh) * | 2016-01-12 | 2019-10-25 | 深圳迈瑞生物医疗电子股份有限公司 | 超声造影成像方法和系统 |
| CN105894442A (zh) * | 2016-03-30 | 2016-08-24 | 武汉声锐达科技有限责任公司 | 应激超声图像自动配准方法 |
| CN109310399B (zh) * | 2016-06-06 | 2022-12-06 | 皇家飞利浦有限公司 | 医学超声图像处理设备 |
| CN108053887A (zh) * | 2017-11-02 | 2018-05-18 | 杭州英放生物科技有限公司 | 一种实体瘤疗效评估方法、电子设备、存储介质、系统 |
| CN114007513B (zh) * | 2019-07-10 | 2025-04-08 | 深圳迈瑞生物医疗电子股份有限公司 | 超声成像设备及检测b线的方法、装置、存储介质 |
| CN110428407A (zh) * | 2019-07-31 | 2019-11-08 | 深圳市人民医院 | 一种超声影像分析系统及其分析方法 |
| CN112634191A (zh) * | 2019-09-24 | 2021-04-09 | 深圳迈瑞生物医疗电子股份有限公司 | 医学影像分析方法、超声成像设备及计算机存储介质 |
| US11980501B2 (en) * | 2020-02-06 | 2024-05-14 | GE Precision Healthcare LLC | Method and system for providing enhanced ultrasound images simulating acquisition at high acoustic power by processing ultrasound images acquired at low acoustic power |
| CN112345646B (zh) * | 2020-10-30 | 2022-03-08 | 安徽理工大学 | 一种细杆中微缺陷重构方法 |
| CN114881974A (zh) * | 2022-05-11 | 2022-08-09 | 上海联影医疗科技股份有限公司 | 一种医学图像分析方法和系统 |
| CN116687353B (zh) * | 2023-08-01 | 2023-12-19 | 宁波杜比医疗科技有限公司 | 新辅助化疗疗效评估系统、设备及介质 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060020204A1 (en) * | 2004-07-01 | 2006-01-26 | Bracco Imaging, S.P.A. | System and method for three-dimensional space management and visualization of ultrasound data ("SonoDEX") |
| CN1961841A (zh) * | 2006-11-30 | 2007-05-16 | 上海交通大学 | 图像引导的水冷式射频消融肿瘤治疗一体机 |
| CN101252886A (zh) * | 2005-06-29 | 2008-08-27 | 艾可瑞公司 | 不用基准标记而用超声波图像对软组织目标的动态跟踪 |
| CN101568942A (zh) * | 2006-12-29 | 2009-10-28 | 皇家飞利浦电子股份有限公司 | 图像引导介入过程中改进的用于补偿手术中发生的运动的图像配准和方法 |
| US20090304252A1 (en) * | 2008-06-05 | 2009-12-10 | Dong Gyu Hyun | Non-Rigid Registration Between CT Images And Ultrasound Images |
| CN101681504A (zh) * | 2006-11-27 | 2010-03-24 | 皇家飞利浦电子股份有限公司 | 用于将实时超声图像与预先获取的医学图像进行融合的系统和方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2460474B1 (en) * | 2003-05-08 | 2015-12-16 | Hitachi Medical Corporation | Reference image display method for ultrasonography and ultrasonic diagnosis apparatus |
| US20090156933A1 (en) * | 2005-09-07 | 2009-06-18 | Koninklijke Philips Electronics, N.V. | Ultrasound system for reliable 3d assessment of right ventricle of the heart and method of doing the same |
| US8771188B2 (en) * | 2007-06-20 | 2014-07-08 | Perception Raisonnement Action En Medecine | Ultrasonic bone motion tracking system |
| US9144461B2 (en) * | 2008-12-03 | 2015-09-29 | Koninklijke Philips N.V. | Feedback system for integrating interventional planning and navigation |
| KR101194282B1 (ko) * | 2010-05-17 | 2012-10-24 | 삼성메디슨 주식회사 | 3차원 초음파 검사기의 디스플레이 시스템 및 방법 |
| GB201019694D0 (en) * | 2010-11-19 | 2011-01-05 | Setred As | Camera and visualisation of images and video for autostereoscopic display |
| BR112014029565B1 (pt) * | 2012-05-31 | 2021-12-28 | Koninklijke Philips N.V. | Sistema de formação de imagens de ultrassom e método de fornecimento de imagem de ultrassom com informações de vasos |
-
2013
- 2013-04-25 CN CN201310147583.6A patent/CN104116523B/zh active Active
- 2013-09-24 WO PCT/CN2013/084044 patent/WO2014173068A1/zh not_active Ceased
-
2015
- 2015-10-23 US US14/921,915 patent/US11083436B2/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060020204A1 (en) * | 2004-07-01 | 2006-01-26 | Bracco Imaging, S.P.A. | System and method for three-dimensional space management and visualization of ultrasound data ("SonoDEX") |
| CN101252886A (zh) * | 2005-06-29 | 2008-08-27 | 艾可瑞公司 | 不用基准标记而用超声波图像对软组织目标的动态跟踪 |
| CN101681504A (zh) * | 2006-11-27 | 2010-03-24 | 皇家飞利浦电子股份有限公司 | 用于将实时超声图像与预先获取的医学图像进行融合的系统和方法 |
| CN1961841A (zh) * | 2006-11-30 | 2007-05-16 | 上海交通大学 | 图像引导的水冷式射频消融肿瘤治疗一体机 |
| CN101568942A (zh) * | 2006-12-29 | 2009-10-28 | 皇家飞利浦电子股份有限公司 | 图像引导介入过程中改进的用于补偿手术中发生的运动的图像配准和方法 |
| US20090304252A1 (en) * | 2008-06-05 | 2009-12-10 | Dong Gyu Hyun | Non-Rigid Registration Between CT Images And Ultrasound Images |
Also Published As
| Publication number | Publication date |
|---|---|
| US11083436B2 (en) | 2021-08-10 |
| US20160045186A1 (en) | 2016-02-18 |
| CN104116523B (zh) | 2016-08-03 |
| CN104116523A (zh) | 2014-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104116523B (zh) | 一种超声影像分析系统及其分析方法 | |
| CN106934807B (zh) | 一种医学影像分析方法、系统及医疗设备 | |
| JP2021041268A (ja) | Ｘ線誘導胸部生検をナビゲートするためのシステムおよび方法 | |
| US20200085412A1 (en) | System and method for using medical image fusion | |
| US8102392B2 (en) | Image processing/displaying apparatus having free moving control unit and limited moving control unit and method of controlling the same | |
| CN107835661B (zh) | 超声图像处理系统和方法及其装置、超声诊断装置 | |
| JP4571187B2 (ja) | 病気の進行または治療効果を解析するために関心領域を複数の時点にわたってリンクさせるシステムおよび方法 | |
| JP5513790B2 (ja) | 超音波診断装置 | |
| CN115811961A (zh) | 三维显示方法和超声成像系统 | |
| US20130257910A1 (en) | Apparatus and method for lesion diagnosis | |
| CN101336844A (zh) | 医用图像处理装置以及医用图像诊断装置 | |
| JP2008086767A (ja) | ３次元及び４次元コントラスト撮像のためのシステム及び方法 | |
| US10916010B2 (en) | Learning data creation support apparatus, learning data creation support method, and learning data creation support program | |
| JP2016049327A (ja) | 医用画像計測装置および方法並びにプログラム | |
| CN108937831A (zh) | 用于淋巴液样本跟踪、引流确定、可视化和处理的系统、装置和方法 | |
| CN116263948A (zh) | 用于图像融合的系统和方法 | |
| JP2007144162A (ja) | 改良型腫瘍アブレーションのシステム及び方法 | |
| US20120230575A1 (en) | Quantification results in multiplane imaging | |
| JP2000350722A (ja) | 器官の注目する要素の配置および三次元表現の方法 | |
| CN110458948A (zh) | 基于智能化3d重建系统中图像的处理方法 | |
| CN102646266A (zh) | 一种处理图像的方法 | |
| Walk et al. | High-frequency ultrasound imaging of mouse cervical lymph nodes | |
| CN108877922A (zh) | 病变程度判断系统及其方法 | |
| CN115205459B (zh) | 一种医学影像处理显示方法及系统 | |
| CN115998423A (zh) | 模拟消融的显示方法和超声成像系统 |
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: 13882727 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 1205 DATED 22.03.2016) |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 13882727 Country of ref document: EP Kind code of ref document: A1 |