CN112951381A - Registration method of liver CT and MRI multi-stage enhancement examination images - Google Patents
Registration method of liver CT and MRI multi-stage enhancement examination images Download PDFInfo
- Publication number
- CN112951381A CN112951381A CN202110096850.6A CN202110096850A CN112951381A CN 112951381 A CN112951381 A CN 112951381A CN 202110096850 A CN202110096850 A CN 202110096850A CN 112951381 A CN112951381 A CN 112951381A
- Authority
- CN
- China
- Prior art keywords
- delta
- period
- images
- portal vein
- layer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
- G16H10/60—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
Abstract
A registration method suitable for a liver CT and MRI multi-stage enhanced examination image solves the problem that the starting positions of the multi-stage enhanced examination image are inconsistent in the scanning process due to various reasons, namely, pictures with the same sequence number in the four stages are not on the same plane or have certain deviation on the longitudinal axis of a human body, namely, the Xth pictures in a flat scanning period, an artery period, a portal vein period and a delay period are not the closest on the longitudinal axis of the human body, a method for positioning a reference layer in the four stages is adopted, and a method for deleting redundant layers is adopted, so that the images with the fixed sequence numbers in the four stages are the closest images with the minimum deviation on the longitudinal axis of the human body.
Description
Technical Field
The invention relates to the field of medical image picture processing, in particular to a registration method of liver CT and MRI multi-stage enhancement examination images
Background
Medical imaging is a medical discipline for diagnosing and treating human diseases by applying medical imaging technology and applying interventional devices under the guidance of medical imaging technology, and is an important component of clinical medicine. Diagnostic imaging (diagnostic imaging) is a major branch of medical imaging. X-ray Computed Tomography (CT)/Magnetic Resonance Imaging (MRI) is the primary modality of imaging diagnostics. The CT/MRI multi-stage enhancement examination technology is mainly used for qualitatively diagnosing diseases of the pelvis and the abdomen and is also a main means for diagnosing liver diseases.
The liver CT multi-stage enhancement examination is that after vein rapid bolus injection of contrast medium, liver artery phase, portal vein phase and delay phase scanning are respectively carried out at different delay time points, and the CT flat scanning examination is combined, so that the method can be used for analyzing the enhancement mode, the enhancement degree and the change of a focus, evaluating the blood supply conditions of the hepatic artery and the portal vein of the focus, and being beneficial to qualitative diagnosis of lesions. Liver MRI enhancement examination is used for flat scan lesion finding, but diagnosis of the case with difficulty, conventional Gd-DTPA injection is carried out, liver T1WI multi-stage enhancement examination is carried out, and the effect and the meaning of the liver MRI enhancement examination are consistent with those of CT multi-stage enhancement examination. Therefore, the diagnosis of liver diseases is critical to the generation of the arterial phase, portal vein phase and delayed phase images by CT/MRI flat scan images and multi-phase enhancement examination techniques.
In the actual operation process, due to the influence of factors such as the respiratory movement of the patient, the posture change of the patient, the operation of the radiography technician and the like, the pictures with the same sequence number in the four phases may not be in the same plane or have large deviation on the longitudinal axis of the human body, namely, the Xth picture in the flat scanning phase, the artery phase, the portal vein phase and the delay phase is not the closest on the longitudinal axis of the human body, and the distance of the deviation is often larger than the layer thickness of the scanning. Although there is also a deviation in the horizontal axis and in the front-back direction, such a deviation is generally small and substantially negligible with respect to the vertical axis.
Images with such a longitudinal axis deviation have the following disadvantages:
in image-related scientific research, particularly in the field of image recognition, manual or semi-automatic labeling of lesions is usually required, when images for multi-stage enhancement examination are researched, multi-stage images need to be observed simultaneously, if longitudinal axis deviation exists, displayed images may have different sizes and positions, and in order to label multi-stage lesions simultaneously in a labeling process, a large labeling frame is often required to be simultaneously brought into the multi-stage lesions, that is, labeling is inaccurate and irregular.
Because the images are not on the same level or the closest level, the images cannot reflect the real level and are not comparable, the images used for research can cause misleading to a computer, and the images with poor registration are generally directly discarded, so that the waste of image resources is caused.
At present, no effective method for solving the problem exists, so that a registration method of images of liver CT/MRI multi-phase enhanced examination (flat scanning phase, arterial phase, portal venous phase and delay phase) is researched, and the method has wide application prospect and great market value.
Disclosure of Invention
In view of the above, the present invention provides a registration method for liver CT and MRI multi-phase enhanced examination images,
in order to achieve the purpose, the first technical scheme adopted by the invention is as follows: a registration method suitable for liver CT and MRI multi-stage enhancement examination images comprises the following steps:
s1, recording the hospitalization number and the examination time corresponding to the case to be registered;
s2, downloading the CT or MRI image with format dicom in the image storage server according to the hospital number and the examination time, and storing the image in a folder named by the hospital number;
s3, opening a folder of the case by using medical image processing software, importing four sequence images of a sweep period, an artery period, a portal vein period and a delay period, and dividing the images into four display interfaces to be displayed in a display;
s4, respectively finding out the layer of the right portal vein trunk in the four sequence images, respectively measuring the long axis length Lx of the right portal vein trunk layer by layer when the layer of the right portal vein trunk appears in the horizontal scanning period, taking the layer with the largest Lx value as the reference layer of the horizontal scanning period, and repeating the steps to respectively record the several pictures of the arterial period, the portal vein period and the delay period, wherein the pictures are respectively in the phase phases of the reference layer and are marked as a, b, c and d;
s5, taking the minimum Nmin and the maximum Nmax of a, b, c and d, and allowing a plurality of same maximum values or minimum values to exist;
s6, calculating differences delta Xa, delta Xb, delta Xc and delta Xd between a, b, c and d and Nmin respectively, deleting delta Xa, delta Xb, delta Xc and delta Xd pictures before the corresponding period respectively, and if the difference between a certain period and Nmin is 0, deleting is not needed;
s7, calculating differences delta Ya, delta Yb, delta Yc and delta Yd between a, b, c and d and Nmax respectively, deleting pictures delta Ya, delta Yb, delta Yc and delta Yd after corresponding periods respectively, and if the difference between a certain period and Nmax is 0, deleting is not needed;
and S8, obtaining the images of the flat scan stage, the artery stage, the portal vein stage and the delay stage after the registration after the deletion.
Further, the medical image processing software in step S3 is a dicom browser.
Further, the four serial images in step S3 are all 72 images.
Further, the step S4 further includes fine-tuning the selected reference level:
from the adjacent images of the slice with the largest Lx value, the liver shape and the relative position of the liver and other organs are finely adjusted by using the auxiliary means, and the final reference slice of the four sequence images is determined.
The key concept of the first technical scheme is as follows: the reference layer of the four-stage image is accurately positioned through the characteristics of the right branch trunk of the portal vein, and the redundant layers are deleted, so that the images with the fixed sequence numbers of the four stages are the closest images with the minimum deviation on the longitudinal axis of the human body.
The beneficial effect of above-mentioned technical scheme one is: the deviation of the registered four-phase image on the longitudinal axis of the human body is minimum, the problem that the initial positions of the multi-phase enhanced inspection images are inconsistent in the scanning process due to various reasons is solved, and the method can be used for research in the field of image identification.
Drawings
FIG. 1 illustrates a registration flow diagram of an embodiment of the present invention;
Detailed Description
The invention is further described below with reference to the following figures and specific examples:
as shown in fig. 1, the registration method for multi-stage enhancement examination images of liver CT and MRI provided by the present invention includes the following steps:
s1, recording the hospitalization number and the examination time corresponding to the case to be registered;
s2, downloading the CT or MRI image with format dicom in the image storage server according to the hospital number and the examination time, and storing the image in a folder named by the hospital number;
s3, opening a folder of the case by using a dicom browser, importing four sequence images of a sweep period, an artery period, a portal vein period and a delay period, and displaying the images in four display interfaces in the display;
s4, respectively finding out the layer of the right portal vein trunk in the four sequence images, respectively measuring the length Lx of the long axis of the right portal vein trunk layer by layer when the layer of the right portal vein trunk appears in the flat scanning period, respectively taking the layer with the largest Lx value as the reference layer of the flat scanning period, finely adjusting the layer with the shape of the liver and the relative positions of the liver and other organs as auxiliary means from the adjacent images of the layer with the largest Lx value, determining the final reference layer of the four sequence images, and repeating the steps to respectively record the pictures of the arterial period, the portal vein period and the reference layer of the delay period in the respective phase, namely a, b, c and d;
s5, taking the minimum Nmin and the maximum Nmax of a, b, c and d, and allowing a plurality of same maximum values or minimum values to exist;
s6, calculating differences delta Xa, delta Xb, delta Xc and delta Xd between a, b, c and d and Nmin respectively, deleting delta Xa, delta Xb, delta Xc and delta Xd pictures before the corresponding period respectively, and if the difference between a certain period and Nmin is 0, deleting is not needed;
s7, calculating differences delta Ya, delta Yb, delta Yc and delta Yd between a, b, c and d and Nmax respectively, deleting pictures delta Ya, delta Yb, delta Yc and delta Yd after corresponding periods respectively, and if the difference between a certain period and Nmax is 0, deleting is not needed;
and S8, obtaining the images of the flat scan stage, the artery stage, the portal vein stage and the delay stage after the registration after the deletion.
Further, the medical image processing software in step S3 is a dicom browser.
Specifically, the working principle of the invention is as follows: the lengths of the scanned images of different longitudinal axis coordinate levels in the CT and MRT images of the right branch trunk of the portal vein are different, so that the right branch trunk of the portal vein can be used for positioning a four-stage image reference level, the number of layers of the level is taken as a reference after the reference level is determined, the total number and the serial number of the registered four-stage images are the same by calculating a difference value and deleting the difference value, and the purpose of registration is achieved.
The following further lists application examples to help those skilled in the art to better understand the technical content of the present invention and the technical contribution of the present invention with respect to the prior art:
the first application embodiment:
and recording the hospitalization number and the examination time corresponding to the case to be registered.
The image in the format dicom (taking the MR image as an example) is downloaded from the image storage server according to the hospital number and the examination time, and the image is saved in the folder named by the hospital number.
A dicom browser (described by medex webmim 6.00.0208, medex technologies ltd, beijing mactex) is used to open the folder of the case, import the image, and 2 × 2 split-screen display shows four sequences of a pan period, an artery period, a portal vein period, and a delay period, wherein the total number of the four sequences is 72, the layer thickness is 3mm, the upper left corner is the pan period, the upper right corner is the artery period, the lower left corner is the portal vein period, and the lower right corner is the delay period.
The number of layers of the right branch trunk of the portal vein is recorded in the four periods respectively, the flat scanning period is 26 th to 30 th layers, the arterial period is 26 th to 29 th layers, the portal vein period is 27 th to 31 th layers, and the delay period is 29 th to 32 th layers.
And measuring the lengths of the long axes of the trunks of the right branches of the portal vein by layers from the 26 th layer to the 30 th layer in the flat scanning period, wherein the lengths are respectively 15.44mm, 17.00mm, 18.09mm, 16.50mm and 14.30mm, and recording the 28 th layer of the layer where the maximum value is 18.09mm as the reference layer in the flat scanning period.
And measuring the lengths of the long axes of the trunks of the right branches of the portal vein at the 26 th layer to the 29 th layer in the arterial period, wherein the lengths are respectively 14.88mm, 17.65mm, 11.76mm and 10.95mm, and recording the 27 th layer of the layer where the maximum value is 17.65mm as the reference layer of the arterial period.
And measuring the lengths of the long axes of the trunks of the right branches of the portal vein layer by layer from 28 th to 30 th in the portal vein period, wherein the lengths are respectively 13.85mm, 15.93mm and 10.08mm, and recording the 29 th layer of the layer where the maximum value is 15.93mm as the portal vein period reference layer.
And measuring the lengths of the long axes of the trunks of the right branches of the portal vein by layers from 29 th to 31 th in the delay period, wherein the lengths are respectively 12.79mm, 13.90mm and 10.15mm, and recording the layer 30 at which the maximum value is 13.90mm as a delay period reference layer.
The minimum value of 28, 27, 29, 30 is taken as 27 and the maximum value is taken as 30.
28. The difference values between 27, 29 and 30 and the minimum value 27 are respectively 1, 0, 2 and 3, and 1, 0, 2 and 3 pictures before the corresponding period are respectively deleted.
28. The differences between 27, 29, 30 and the maximum value 30 are 2, 3, 1, 0 respectively, and 2, 3, 1, 0 pictures after the corresponding period are deleted respectively.
The total number of the deleted images in the flat scanning period, the artery period, the portal vein period and the delay period is the same as 69, and the sequence numbers of the images of the reference layer in the four periods are the same as 27, so that the aim of registration is fulfilled.
Application example two:
and recording the hospitalization number and the examination time corresponding to the case to be registered.
The image in the format dicom (taking the MR image as an example) is downloaded from the image storage server according to the hospital number and the examination time, and the image is saved in the folder named by the hospital number.
A dicom browser (described by medex webmim 6.00.0208, medex technologies ltd, beijing madex technology) is used to open the folder of the case, import the image, and 2 × 2 split-screen display shows four sequences of a pan period, an artery period, a portal vein period, and a delay period, wherein the total number of the four sequences is 72, the upper left corner is the pan period, the upper right corner is the artery period, the lower left corner is the portal vein period, and the lower right corner is the delay period.
The level of the right portal vein trunk is found in the four display interfaces, the display interfaces stay at the level of the right portal vein trunk, and the images displayed in the fourth stage are 36 th, 38 th, 41 th and 43 th images.
In the four display interfaces, the shape of the right branch of the portal vein trunk is taken as the main part, the shape of the liver and the relative position of the liver and other organs are taken as the auxiliary parts, the four display interface layers are compared with each other for fine adjustment, the shapes of the portal vein right branch trunks of the reference layers found by the four sequences are ensured to be consistent as much as possible, the outlines of the surrounding tissues and the liver are also kept consistent as much as possible, namely the deviation on the longitudinal axis of the human body is minimum, and the reference layers of the horizontal scanning period, the artery period, the portal vein period and the delay period are respectively recorded as 35, 37, 40 and 43.
The minimum value of 35, 37, 40, 43 is 35 and the maximum value is 43.
35. The difference values between 37, 40 and 43 and the minimum value 35 are respectively 0, 2, 5 and 8, and 0, 2, 5 and 8 pictures before the corresponding period are respectively deleted.
35. The difference values between 37, 40 and 43 and the maximum value 43 are respectively 8, 6, 3 and 0, and 8, 6, 3 and 0 pictures after the corresponding period are respectively deleted.
The total number of the deleted images in the flat scanning period, the artery period, the portal vein period and the delay period is the same as 64, and the sequence number of the image of the reference layer in the four periods is the same as 35, so that the aim of registration is fulfilled.
The present invention has been described with reference to the above embodiments and the accompanying drawings, however, the above embodiments are only examples for carrying out the present invention. It should be noted that the disclosed embodiments do not limit the scope of the invention. Rather, modifications and equivalent arrangements included within the spirit and scope of the claims are included within the scope of the invention.
Claims (4)
1. A registration method suitable for liver CT and MRI multi-stage enhancement examination images is characterized by comprising the following steps:
s1, recording the hospitalization number and the examination time corresponding to the case to be registered;
s2, downloading the CT or MRI image with format dicom in the image storage server according to the hospital number and the examination time, and storing the image in a folder named by the hospital number;
s3, opening a folder of the case by using medical image processing software, importing four sequence images of a sweep period, an artery period, a portal vein period and a delay period, and dividing the images into four display interfaces to be displayed in a display;
s4, respectively finding out the layer of the right portal vein trunk in the four sequence images, respectively measuring the long axis length Lx of the right portal vein trunk layer by layer when the layer of the right portal vein trunk appears in the horizontal scanning period, taking the layer with the largest Lx value as the reference layer of the horizontal scanning period, and repeating the steps to respectively record the several pictures of the arterial period, the portal vein period and the delay period, wherein the pictures are respectively in the phase phases of the reference layer and are marked as a, b, c and d;
s5, taking the minimum Nmin and the maximum Nmax of a, b, c and d, and allowing a plurality of same maximum values or minimum values to exist;
s6, calculating differences delta Xa, delta Xb, delta Xc and delta Xd between a, b, c and d and Nmin respectively, deleting delta Xa, delta Xb, delta Xc and delta Xd pictures before the corresponding period respectively, and if the difference between a certain period and Nmin is 0, deleting is not needed;
s7, calculating differences delta Ya, delta Yb, delta Yc and delta Yd between a, b, c and d and Nmax respectively, deleting pictures delta Ya, delta Yb, delta Yc and delta Yd after corresponding periods respectively, and if the difference between a certain period and Nmax is 0, deleting is not needed;
and S8, obtaining the images of the flat scan stage, the artery stage, the portal vein stage and the delay stage after the registration after the deletion.
2. The method as claimed in claim 1, wherein the medical image processing software in step S3 is a dicom browser.
3. The method as claimed in claim 1, wherein the number of the four sequential images in step S3 is 72.
4. The registration method for liver CT, MRI multi-stage enhancement examination images as claimed in claim 1, wherein said step S4 further includes fine-tuning the selected reference level:
from the adjacent images of the slice with the largest Lx value, the liver shape and the relative position of the liver and other organs are finely adjusted by using the auxiliary means, and the final reference slice of the four sequence images is determined.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010465017X | 2020-05-28 | ||
| CN202010465017 | 2020-05-28 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112951381A true CN112951381A (en) | 2021-06-11 |
| CN112951381B CN112951381B (en) | 2023-02-28 |
Family
ID=76236395
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110096850.6A Active CN112951381B (en) | 2020-05-28 | 2021-01-25 | Registration method of liver CT and MRI multi-phase enhancement examination images |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112951381B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116563358A (en) * | 2023-07-07 | 2023-08-08 | 浙江大学 | Data alignment preprocessing method for liver enhancement multi-stage CT data AI training |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108053433A (en) * | 2017-11-28 | 2018-05-18 | 浙江工业大学 | A kind of multi-modal arteria carotis MRI method for registering based on physical alignment and outline |
| US20190029561A1 (en) * | 2017-07-26 | 2019-01-31 | Korea Institute Of Science And Technology | Image registration system and method using subject-specific tracker |
| CN110853082A (en) * | 2019-10-21 | 2020-02-28 | 科大讯飞股份有限公司 | Medical image registration method and device, electronic equipment and computer storage medium |
-
2021
- 2021-01-25 CN CN202110096850.6A patent/CN112951381B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190029561A1 (en) * | 2017-07-26 | 2019-01-31 | Korea Institute Of Science And Technology | Image registration system and method using subject-specific tracker |
| CN108053433A (en) * | 2017-11-28 | 2018-05-18 | 浙江工业大学 | A kind of multi-modal arteria carotis MRI method for registering based on physical alignment and outline |
| CN110853082A (en) * | 2019-10-21 | 2020-02-28 | 科大讯飞股份有限公司 | Medical image registration method and device, electronic equipment and computer storage medium |
Non-Patent Citations (1)
| Title |
|---|
| 汪晓妍等: "基于空间对齐和轮廓匹配的颈动脉多对比MRI三维配准方法", 《计算机科学》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116563358A (en) * | 2023-07-07 | 2023-08-08 | 浙江大学 | Data alignment preprocessing method for liver enhancement multi-stage CT data AI training |
| CN116563358B (en) * | 2023-07-07 | 2023-09-22 | 浙江大学 | Data alignment preprocessing method for liver enhancement multi-stage CT data AI training |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112951381B (en) | 2023-02-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8391643B2 (en) | Image registration method | |
| JP5475347B2 (en) | Identification of white matter fiber tracts using magnetic resonance imaging (MRI) | |
| US20090303252A1 (en) | Registration Of CT Image Onto Ultrasound Images | |
| JP4490442B2 (en) | Method and system for affine superposition of an intraoperative 2D image and a preoperative 3D image | |
| US7346381B2 (en) | Method and apparatus for medical intervention procedure planning | |
| Bidaut et al. | Automated registration of dynamic MR images for the quantification of myocardial perfusion | |
| US20190223750A1 (en) | Medical image diagnostic apparatus and medical image display apparatus for volume image correlations | |
| CN109872312B (en) | Medical image segmentation method, device and system and image segmentation method | |
| US20030233039A1 (en) | Physiological model based non-rigid image registration | |
| CN106469449B (en) | Method and device for displaying focus in medical image | |
| US8559686B2 (en) | Medical image data alignment apparatus, method and program | |
| US20050196028A1 (en) | Method of registering a sequence of 2D image data with 3D image data | |
| US20160066795A1 (en) | Stenosis therapy planning | |
| US20090087047A1 (en) | Image display device and image display program storage medium | |
| US7634301B2 (en) | Repeated examination reporting | |
| JP2000185036A (en) | Medical image display device | |
| US8036730B1 (en) | Temporal magnetic resonance imaging | |
| US7627158B2 (en) | Automatic registration of intra-modality medical volume images using affine transformation | |
| WO2009050676A1 (en) | Pathology-related magnetic resonance imaging | |
| CN112767370A (en) | Target region labeling method based on image registration | |
| CN112951381B (en) | Registration method of liver CT and MRI multi-phase enhancement examination images | |
| US20120306862A1 (en) | Image processing device, method and program | |
| Bizino et al. | High spatial resolution coronary magnetic resonance angiography at 7 T: comparison with low spatial resolution bright blood imaging | |
| US20190266727A1 (en) | System and method for using imaging quality metric ranking | |
| Kim et al. | SLIMMER: SLIce MRI motion estimation and reconstruction tool for studies of fetal anatomy |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |