CN115137383B - Method and device for selecting and retrospective reconstruction data of heart spiral retrospective reconstruction - Google Patents
Method and device for selecting and retrospective reconstruction data of heart spiral retrospective reconstruction Download PDFInfo
- Publication number
- CN115137383B CN115137383B CN202210667055.2A CN202210667055A CN115137383B CN 115137383 B CN115137383 B CN 115137383B CN 202210667055 A CN202210667055 A CN 202210667055A CN 115137383 B CN115137383 B CN 115137383B
- Authority
- CN
- China
- Prior art keywords
- phase
- reconstruction
- cardiac cycle
- image
- cardiac
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- 230000000747 cardiac effect Effects 0.000 claims abstract description 191
- 230000008569 process Effects 0.000 claims description 18
- 238000012552 review Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000010009 beating Methods 0.000 description 8
- 238000002565 electrocardiography Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 5
- 230000018199 S phase Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 210000004351 coronary vessel Anatomy 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 206010011224 Cough Diseases 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Surgery (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Cardiology (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Signal Processing (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- General Physics & Mathematics (AREA)
- Pulmonology (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The embodiment of the invention discloses a method and a device for selecting and retrospectively reconstructing data of a heart spiral retrospective reconstruction, which determine a focus position needing image reconstruction; determining a reconstruction phase of the location of interest; determining a first cardiac cycle in which the reconstruction phase is positioned, and reconstructing a first image sequence of the optimal reconstruction phase in the first cardiac cycle; reconstructing images in a preset phase range in a second cardiac cycle to obtain a second image sequence; comparing the overlapped images of the first image sequence and the second image sequence, and taking the reconstruction phase corresponding to the image in the second image sequence which meets the requirement as the optimal reconstruction phase in the second cardiac cycle; and taking the second cardiac cycle as the first cardiac cycle, and repeating the steps until the optimal reconstruction phase in each cardiac cycle is determined. The invention can improve the motion artifact generated in the heart spiral retrospective reconstruction.
Description
Technical Field
The application relates to the technical field of imaging, in particular to a cardiac spiral review reconstruction data selecting and reviewing reconstruction method and device.
Background
CT provides a clinical method for non-invasive examination of the anatomy of a tissue or organ in a patient, and CT scan reconstruction requires that the scanned object remain stationary during the scan. In clinical practice, due to the physiological characteristics of the heart that does not move continuously, if a conventional scanning reconstruction protocol is adopted, serious motion artifacts exist in the finally obtained image, so that coronary vessels cannot be observed to influence the final clinical diagnosis.
Scanning of the heart is well established in the art by rapidly scanning the heart of a patient while simultaneously acquiring Electrocardiography (ECG) data from the patient. During image reconstruction, according to the acquired ECG gating curve information, data segments of appointed phases in each cardiac cycle are selected, wherein the data segments meet the reconstruction requirements. Phase here refers to the percentage of the designated locations within adjacent R-waves in the ECG signal that are within the RR interval, which is 0% -100%. In this process, we consider that the motion states of the heart are consistent at the same phase time in different heartbeat cycles. The reconstructed images are consistent and have the minimum motion artifact by combining the data of the same phase range in different cardiac cycles.
In practice, however, the trajectory of the human heart motion cannot be perfectly repeated in each cardiac cycle, and the patient's heart rate may vary over several cardiac cycles during the scan, as well as other situations such as patient respiration. This results in that if the same phase data of each cardiac cycle is simply used for reconstruction, the corresponding true heart morphology may also change, eventually leading to additional motion artifacts in the reconstructed image.
No effective solution has been proposed to solve the problem of reconstruction image artifacts caused by the inconsistency between the beating state of the patient's heart and the ECG signal cycle or the aperiodic variation of the patient's heart during the beating in the related art.
Disclosure of Invention
The embodiment of the invention provides a method and a device for selecting and retrospective reconstructing data of a heart spiral retrospective reconstruction, a computer device and a storage medium, which are used for solving the problem of reconstructed image artifacts caused by inconsistency between the beating state of the heart of a patient and the period of an ECG signal or aperiodic change of the heart of the patient in the beating process in the related art.
In order to achieve the above object, in a first aspect of the embodiments of the present invention, there is provided a cardiac spiral review reconstruction data sorting method, including:
step 1: determining a focus position needing image reconstruction, wherein the focus position is a position where an artifact is likely to occur;
step 2: determining an optimal reconstruction phase of the position of interest, wherein the optimal reconstruction phase refers to a cardiac motion phase with the smallest cardiac motion amount among the plurality of cardiac motion phases;
and step 3: determining a first cardiac cycle in which the optimal reconstruction phase is located, and reconstructing a first image sequence of the optimal reconstruction phase in the first cardiac cycle;
and 4, step 4: reconstructing images in a second cardiac cycle within a preset phase range to obtain a second image sequence, wherein the second cardiac cycle comprises a previous cardiac cycle or a next cardiac cycle adjacent to the first cardiac cycle;
and 5: comparing the overlapped images of the first image sequence and the second image sequence, and taking the reconstruction phase corresponding to the image in the second image sequence meeting the requirement as the optimal reconstruction phase in the second cardiac cycle;
step 6: repeating steps 2-5 with the second cardiac cycle as the first cardiac cycle until an optimal reconstruction phase is determined for each cardiac cycle.
Optionally, in a possible implementation manner of the first aspect, comparing the first image sequence with the overlapping images of the second image sequence, and taking a reconstruction phase corresponding to an image in the second image sequence that meets the requirement as an optimal reconstruction phase in the second cardiac cycle includes:
comparing the overlapped images in the first image sequence and the second image sequence, and selecting at least one image of which the difference degree is smaller than a preset threshold value;
and judging whether the phase distance between the corresponding phase of the at least one image and the optimal reconstruction phase of the first cardiac cycle is within a preset phase threshold range, and if so, taking the corresponding phase as the optimal reconstruction phase in the second cardiac cycle.
Optionally, in a possible implementation manner of the first aspect, the determining the optimal reconstruction phase of the location of interest includes:
reconstructing images corresponding to all phases within the range of 0% -99% of the concerned position at intervals according to a preset phase;
and calculating the image difference of all the images at the same attention position in all the phase ranges, and determining the phase corresponding to the image with the minimum image difference as the optimal reconstruction phase.
Optionally, in a possible implementation manner of the first aspect, the determining the optimal reconstruction phase of the location of interest includes:
calculating a target phase range to be reconstructed according to a rule between the heart rate and the reconstructed target phase;
and evaluating the definition of all images in the target phase range by adopting an image information entropy calculation method, and taking the phase of the image with the highest definition as the optimal reconstruction phase.
In a second aspect of the embodiments of the present invention, there is provided a cardiac spiral retrospective reconstruction method, including:
selecting reconstructed data based on the cardiac spiral review reconstructed data selecting method of any one of claims 1 to 5;
and reconstructing the reconstruction data to obtain a heart reconstruction image.
Optionally, in a possible implementation manner of the second aspect, the method further includes:
in the process of carrying out data selection reconstruction according to the optimal reconstruction phase in each cardiac cycle to obtain a cardiac reconstruction image sequence, if the reconstruction images obtained from the data selected in any two adjacent cardiac cycles have overlapping areas with different degrees, calculating difference points between two groups of reconstruction images, and processing the difference points according to a motion deformation field to eliminate motion artifacts generated in the helical review reconstruction of the heart.
In a third aspect of the embodiments of the present invention, there is provided a cardiac spiral review reconstruction data sorting apparatus, including:
the device comprises a focus position determining module, a focus position determining module and a focusing module, wherein the focus position determining module is used for determining a focus position needing image reconstruction, and the focus position is a position where artifacts possibly appear;
a first optimal reconstruction phase determining module, configured to determine an optimal reconstruction phase of the location of interest, where the optimal reconstruction phase is a cardiac motion phase with a smallest cardiac motion amount among the cardiac motion phases;
the first image sequence reconstruction module is used for determining a first cardiac cycle in which the optimal reconstruction phase is positioned and reconstructing a first image sequence of the optimal reconstruction phase in the first cardiac cycle;
the second image sequence reconstruction module is used for reconstructing images in a second cardiac cycle within a preset phase range to obtain a second image sequence, wherein the second cardiac cycle comprises a previous cardiac cycle or a next cardiac cycle which is adjacent to the first cardiac cycle;
a second optimal reconstruction phase determining module, configured to compare the first image sequence with an overlapped image of the second image sequence, and use a reconstruction phase corresponding to an image in the second image sequence that meets the requirement as an optimal reconstruction phase in a second cardiac cycle;
and the expansion module is used for taking the second cardiac cycle as the first cardiac cycle and repeatedly executing the first optimal reconstruction phase determining module, the first image sequence reconstructing module, the second image sequence reconstructing module and the second optimal reconstruction phase determining module until the optimal reconstruction phase in each cardiac cycle is determined.
Optionally, in a possible implementation manner of the third aspect, the first optimal rebuilding period determining module includes:
the image reconstruction unit is used for reconstructing images corresponding to all phases within the range of 0% -99% of the concerned position according to a preset phase interval;
and the optimal reconstruction phase determining unit is used for calculating the image difference of all the images at the same attention position in all the phase ranges and determining the phase corresponding to the image with the minimum image difference as the optimal reconstruction phase.
In a fourth aspect of the embodiments of the present invention, a computer device is provided, which includes a memory and a processor, where the memory stores a computer program that can be run on the processor, and the processor executes the computer program to implement the steps in the above-mentioned method embodiments.
A fifth aspect of the embodiments of the present invention provides a readable storage medium, in which a computer program is stored, which computer program, when being executed by a processor, is adapted to carry out the steps of the method according to the first aspect of the present invention and any of the various possible designs of the first aspect.
According to the method and the device for selecting and reviewing reconstruction data of cardiac spiral review, computer equipment and storage media, the position of interest to be reconstructed by the image is determined, and the position of interest is a position where artifacts may appear; determining an optimal reconstruction phase of the position of interest, wherein the optimal reconstruction phase refers to a cardiac motion phase with the smallest cardiac motion amount among the plurality of cardiac motion phases; determining a first cardiac cycle in which the optimal reconstruction phase is located, and reconstructing a first image sequence of the optimal reconstruction phase in the first cardiac cycle; reconstructing images in a second cardiac cycle within a preset phase range to obtain a second image sequence, wherein the second cardiac cycle comprises a previous cardiac cycle or a next cardiac cycle adjacent to the first cardiac cycle; comparing the overlapped images of the first image sequence and the second image sequence, and taking the reconstruction phase corresponding to the image in the second image sequence meeting the requirement as the optimal reconstruction phase in the second cardiac cycle; and taking the second cardiac cycle as the first cardiac cycle, and repeating the steps until the optimal reconstruction phase in each cardiac cycle is determined. The invention can reduce the probability of inconsistent data motion states selected by the heart spiral retrospective scanning reconstruction of different cardiac cycles, and correct the data inconsistent artifact caused by the selection of the same phase in the adjacent cardiac cycles due to the inconsistent beating states of the heart in different cardiac cycles or other reasons, thereby realizing the technical effect of improving the motion artifact generated in the heart spiral retrospective reconstruction.
Drawings
FIG. 1 is a schematic flow chart diagram illustrating an alternative method for selecting cardiac spiral review reconstruction data according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of data selection from ECG versus image reconstruction location in a helical retrospective cardiac reconstruction;
fig. 3 is a block diagram of an alternative cardiac spiral review reconstruction data sorting apparatus according to an embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
The terms "first," "second," "third," "fourth," and the like in the description and in the claims, as well as in the drawings, if any, are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used is interchangeable under appropriate circumstances such that the embodiments of the invention described herein are capable of operation in other sequences than those illustrated or described herein.
It should be understood that, in various embodiments of the present invention, the sequence numbers of the processes do not mean the execution sequence, and the execution sequence of the processes should be determined by the functions and the internal logic of the processes, and should not constitute any limitation on the implementation process of the embodiments of the present invention.
It should be understood that in the present application, "comprising" and "having" and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
It should be understood that, in the present invention, "a plurality" means two or more. "and/or" is merely an association relationship describing an associated object, meaning that there may be three relationships, for example, and/or B, may mean: a exists alone, A and B exist simultaneously, and B exists alone. The character "/" generally indicates that the former and latter associated objects are in an "or" relationship. "comprises A, B and C" and "comprises A, B, C" means that all three of A, B, C comprise, "comprises A, B or C" means that one of three of A, B, C is comprised, "comprises A, B and/or C" means that any 1 or any 2 or 3 of the three of A, B, C is comprised.
It should be understood that in the present invention, "B corresponding to a", "a corresponds to B", or "B corresponds to a" means that B is associated with a, and B can be determined from a. Determining B from a does not mean determining B from a alone, but may also be determined from a and/or other information. And the matching of A and B means that the similarity of A and B is greater than or equal to a preset threshold value.
As used herein, if can be interpreted as "at … …" or "at … …" or "in response to a determination" or "in response to a detection", depending on the context.
The technical solution of the present invention will be described in detail below with specific examples. The following several specific embodiments may be combined with each other, and details of the same or similar concepts or processes may not be repeated in some embodiments.
As described in the background art, retrospective reconstruction of the heart in the prior art generally employs combining data of the same phase range in different cardiac cycles to reconstruct an image; however, the actual situation is that the motion trajectory of the heart is not perfectly repeated in each cardiac cycle, and during the scanning process, the heart rate of the patient may change, and the patient may cough due to breathing, which may cause the real heart morphology to also change, thereby causing the problem of reconstructed image artifact due to inconsistency between the beating state of the heart of the patient and the ECG signal cycle or aperiodic change of the heart of the patient during the beating process.
Therefore, the invention provides a method for adaptively selecting data ranges in different cardiac cycles near a target reconstruction period selected by a user, so that the selected data ranges in different cardiac cycles correspond to the same heart motion state, which comprises the following specific steps:
example 1:
the invention provides a method for selecting heart spiral retrospective reconstruction data, which is shown as a flow chart in figure 1 and comprises the following steps:
when the helical retrospective scan reconstruction of the heart is performed, firstly, scan data synchronized with ECG information of a patient is obtained, and then when a reconstruction task is performed, the scan data needs to be selected, wherein the specific selection steps are as follows:
step 1, determining a concerned position needing image reconstruction.
In this step, for cardiac reconstruction, the degree of motion artifact at different positions is very different, usually the motion amplitude of the right coronary vessel of the heart is much larger, and the right coronal middle segment is the position where the motion artifact occurs most frequently in cardiac reconstruction, so a position of interest S (the position where the artifact may occur) needs to be determined in advance, as shown by S in fig. 2. The focus position can be obtained by manual marking on the positioning sheet by a user; the position range of the middle section of the right crown can be simply deduced by calculating the heart reconstruction range; it can also be obtained by applying more sophisticated recognition algorithms based on pre-reconstructed whole images of the heart.
And 2, determining the optimal reconstruction phase of the concerned position.
And 3, determining a first cardiac cycle in which the optimal reconstruction phase is located, and reconstructing a first image sequence of the optimal reconstruction phase in the first cardiac cycle.
In step 2-3, after obtaining the image position needing attention, it is necessary to determine an optimal reconstruction phase at the attention position, where the optimal reconstruction phase refers to a motion phase with the smallest amount of cardiac motion in the plurality of cardiac motion phases, that is, the motion artifact generated by the cardiac review reconstruction image under the phase is the smallest.
Methods of determining the optimal reconstruction phase at the location of interest include, but are not limited to: 1. reconstructing images corresponding to all phases within the range of 0% -99% according to the interval of a preset phase (for example, 1% phase) at the concerned position, and determining the phase corresponding to the minimum difference as the reconstruction optimal phase P by calculating the image difference at the same position within the range of all phases S (ii) a 2. Calculating according to the actual application and the heart rate sampled by the patient and the target phase reconstruction rule to obtain a target phase range to be reconstructed, evaluating the definition of all images in the target phase range by adopting an image information entropy calculation method, and taking the phase of the image with the highest definition as the optimal reconstruction phase; 3. providing a corresponding interface for a user in a software interface to obtain relevant range and phase interval parameters, displaying all phase range images to the user, and obtaining a reconstruction phase P which is considered by the user as the best image through the interface S 。
The optimal reconstruction phase at the determined location of interest S is P S The data range used (the data required to reconstruct the image for each phase is within a range of projection angles), and the optimal reconstruction phase is determined as P S The corresponding cardiac cycle P S So as to reconstruct the whole image that can be reconstructed by the data range, i.e. the first image sequence.
And 4, reconstructing the image in the preset phase range in the second cardiac cycle to obtain a second image sequence.
In step 4, the second cardiac cycle includes a previous cardiac cycle or a subsequent cardiac cycle adjacent to the first cardiac cycle. The phase range of the preset period can be specified according to the experience of technical personnel, and the range needs to be ensured to be between 0% and 99%. And reconstructing the images in the preset phase range to obtain a second image sequence.
And 5, comparing the overlapped images of the first image sequence and the second image sequence, and taking the reconstruction phase corresponding to the image in the second image sequence meeting the requirement as the optimal reconstruction phase in the second cardiac cycle.
In step 5, in the process of comparing the overlapped images in the first image sequence and the second image sequence, at least one image of which the image difference degree is smaller than a preset threshold value needs to be selected; then judging whether the phase distance between the corresponding phase of the at least one image and the optimal reconstruction phase of the first cardiac cycle is within a preset phase threshold range, and if so, taking the corresponding phase as a second cardiac cycleOptimal reconstruction phase within the cycle. For example: five images (a first image, a second image, a third image, a fourth image and a fifth image) in the second image sequence are overlapped with the first image sequence, and the five images are compared with the first image sequence, wherein the first image, the second image and the third image are similar to the images in the first image sequence, and the difference degrees are smaller than a certain threshold value; respectively comparing the phase phases corresponding to the first image, the second image and the third image with the optimal reconstruction phase P of the first image sequence S Phase distance judgment is carried out, and the phase of the third image is the phase P of the optimal reconstruction phase S Phase distance of 5%, phase of the second image and optimal reconstruction phase P S Phase distance of 3%, phase of the first image and optimal reconstruction phase P S The phase distance of (1%) and the preset phase threshold range of [0%,2%]Therefore, only the phase corresponding to the first image meets the requirement, i.e. the phase corresponding to the first image is used as the optimal reconstruction phase in the second cardiac cycle.
And 6, taking the second cardiac cycle as the first cardiac cycle, and repeating the steps 2-5 until the optimal reconstruction phase in each cardiac cycle is determined.
In step 6, the optimal reconstruction phase in each cardiac cycle required for the whole heart is determined in turn, extending from the first cardiac cycle to both sides according to the method described above, for example: all cardiac cycles include: r1, R2.. Rs-1, rs +1 and Rs + 2.; after the optimal reconstruction phase of the first cardiac cycle Rs is determined, the optimal reconstruction phase of a second cardiac cycle Rs-1 or Rs +1 adjacent to Rs is determined according to the method; then taking the Rs-1 or the Rs +1 as a first cardiac cycle, and determining the optimal reconstruction phase of a second cardiac cycle Rs-2 adjacent to the Rs-1 and the optimal reconstruction phase of a second cardiac cycle Rs +2 adjacent to the Rs +1 according to the method; finally, the optimal reconstruction phase of each cardiac cycle is determined in turn by the expansion mode. As shown in fig. 2, the heart can be reconstructed with data exposed in 4 cardiac cycles: AA '+ BB' + CC '+ DD'. Wherein the time range available for reconstruction in the range of each cardiac cycle is longer than light grayThe column indicates that the phase at the center of the time range is P 1 、P 2 、P 3 、P 4 The data used for reconstruction are indicated by dark grey short columns.
In a typical implementation, the phase of the selected data segment in each cardiac cycle is the same, and in the present invention, in order to make the cardiac images reconstructed from the data used in different cardiac cycles have more consistent morphology, the reconstruction phases in the finally selected cardiac cycles are not necessarily the same, i.e. P1, P2, P3, P4 are not necessarily the same value.
According to the heart spiral retrospective reconstruction data selection method provided by the invention, the concerned position needing image reconstruction is determined, and the concerned position is a position where an artifact is likely to occur; determining an optimal reconstruction phase of the position of interest, wherein the optimal reconstruction phase refers to a cardiac motion phase with the smallest cardiac motion amount among the plurality of cardiac motion phases; determining a first cardiac cycle in which the optimal reconstruction phase is located, and reconstructing a first image sequence of the optimal reconstruction phase in the first cardiac cycle; reconstructing images in a second cardiac cycle within a preset phase range to obtain a second image sequence, wherein the second cardiac cycle comprises a previous cardiac cycle or a next cardiac cycle adjacent to the first cardiac cycle; comparing the overlapped images of the first image sequence and the second image sequence, and taking the reconstruction phase corresponding to the image in the second image sequence meeting the requirement as the optimal reconstruction phase in the second cardiac cycle; and taking the second cardiac cycle as the first cardiac cycle, and repeating the steps until the optimal reconstruction phase in each cardiac cycle is determined. The invention can reduce the probability of inconsistent data motion states selected by the heart spiral retrospective scanning reconstruction of different cardiac cycles, and correct the data inconsistent artifact caused by the selection of the same phase in the adjacent cardiac cycles due to the inconsistent beating states of the heart in different cardiac cycles or other reasons, thereby realizing the technical effect of improving the motion artifact generated in the heart spiral retrospective reconstruction.
Example 2:
embodiments of the present invention further provide a cardiac spiral retrospective reconstruction method, including:
selecting reconstructed data based on the cardiac spiral review reconstructed data selecting method of any one of claims 1 to 5;
and reconstructing the reconstruction data to obtain a heart reconstruction image.
In one embodiment, the method further comprises:
in the process of carrying out data selection reconstruction according to the optimal reconstruction phase in each cardiac cycle to obtain a heart reconstruction image sequence, if reconstruction images obtained from data selected in any two adjacent cardiac cycles have overlapping areas with different degrees, difference points between two groups of reconstruction images are calculated, and the difference points are processed according to a motion deformation field to eliminate motion artifacts generated in the process of heart spiral review reconstruction.
In this embodiment, the entire cardiac image sequence is reconstructed by selecting data according to the optimal reconstruction phase obtained in each cardiac cycle. During the reconstruction, any two adjacent cardiac cycles R S-1 、R S+1 The reconstructed images obtained from the selected data have different overlapping areas, as shown in the reconstructed data (indicated by dark gray short bars) in fig. 2, and the overlapping area a ' B, BB ' of AA ' and BB ' and the overlapping area B ' C, CC ' and DD ' of CC ' have an overlapping area C ' D. According to the overlapping area, the slight difference between two adjacent groups of reconstructed images can be calculated through an algorithm, and the relative motion field MVF in the adjacent cardiac cycle can be further calculated, so that the motion deformation field MVF is applied in the image reconstruction process to further eliminate the data inconsistency motion artifact caused by the reconstruction of the adjacent data segments.
Example 3:
an embodiment of the present invention further provides a cardiac spiral review reconstruction data selecting apparatus, as shown in fig. 3, which includes:
the device comprises a focus position determining module, a focus position determining module and a focus position determining module, wherein the focus position determining module is used for determining a focus position needing image reconstruction, and the focus position is a position where an artifact is likely to occur;
a first optimal reconstruction phase determining module, configured to determine an optimal reconstruction phase of the location of interest, where the optimal reconstruction phase is a cardiac motion phase with a smallest cardiac motion amount among the cardiac motion phases;
the first image sequence reconstruction module is used for determining a first cardiac cycle in which the optimal reconstruction phase is positioned and reconstructing a first image sequence of the optimal reconstruction phase in the first cardiac cycle;
the second image sequence reconstruction module is used for reconstructing images in a second cardiac cycle within a preset phase range to obtain a second image sequence, wherein the second cardiac cycle comprises a previous cardiac cycle or a next cardiac cycle which is adjacent to the first cardiac cycle;
a second optimal reconstruction phase determining module, configured to compare the first image sequence with the overlapping image of the second image sequence, and use a reconstruction phase corresponding to an image in the second image sequence that meets the requirement as an optimal reconstruction phase in a second cardiac cycle;
and the expansion module is used for taking the second cardiac cycle as the first cardiac cycle and repeatedly executing the first optimal reconstruction phase determining module, the first image sequence reconstructing module, the second image sequence reconstructing module and the second optimal reconstruction phase determining module until the optimal reconstruction phase in each cardiac cycle is determined.
In one embodiment, the first optimal reconstruction period phase determination module includes:
the image reconstruction unit is used for reconstructing images corresponding to all phases within the range of 0% -99% of the concerned position according to a preset phase interval;
and the optimal reconstruction phase determining unit is used for calculating the image difference of all the images at the same attention position in all the phase ranges and determining the phase corresponding to the image with the minimum image difference as the optimal reconstruction phase.
The readable storage medium may be a computer storage medium or a communication medium. Communication media includes any medium that facilitates transfer of a computer program from one place to another. Computer storage media may be any available media that can be accessed by a general purpose or special purpose computer. For example, a readable storage medium is coupled to the processor such that the processor can read information from, and write information to, the readable storage medium. Of course, the readable storage medium may also be an integral part of the processor. The processor and the readable storage medium may reside in an Application Specific Integrated Circuits (ASIC). Additionally, the ASIC may reside in user equipment. Of course, the processor and the readable storage medium may also reside as discrete components in a communication device. The readable storage medium may be read-only memory (ROM), random-access memory (RAM), CD-ROMs, magnetic tapes, floppy disks, optical data storage devices, and the like.
The present invention also provides a program product comprising execution instructions stored in a readable storage medium. The at least one processor of the device may read the execution instructions from the readable storage medium, and the execution of the execution instructions by the at least one processor causes the device to implement the methods provided by the various embodiments described above.
In the above embodiments of the terminal or the server, it should be understood that the Processor may be a Central Processing Unit (CPU), other general-purpose processors, a Digital Signal Processor (DSP), etc. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of a method disclosed in connection with the present invention may be embodied directly in a hardware processor, or in a combination of the hardware and software modules within the processor.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and these modifications or substitutions do not depart from the spirit of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A method for selecting data for cardiac spiral retrospective reconstruction, comprising:
step 1: determining a focus position needing image reconstruction, wherein the focus position is a position where an artifact is likely to occur;
step 2: determining an optimal reconstruction phase of the position of interest, wherein the optimal reconstruction phase refers to a cardiac motion phase with the smallest amount of cardiac motion among the plurality of cardiac motion phases;
and step 3: determining a first cardiac cycle in which the optimal reconstruction phase is located, and reconstructing a first image sequence of the optimal reconstruction phase in the first cardiac cycle;
and 4, step 4: reconstructing images in a second cardiac cycle within a preset phase range to obtain a second image sequence, wherein the second cardiac cycle comprises a previous cardiac cycle or a next cardiac cycle adjacent to the first cardiac cycle;
and 5: comparing the overlapped images of the first image sequence and the second image sequence, and taking the reconstruction phase corresponding to the image in the second image sequence meeting the requirement as the optimal reconstruction phase in the second cardiac cycle;
and 6: repeating steps 2-5 with the second cardiac cycle as the first cardiac cycle until an optimal reconstruction phase is determined for each cardiac cycle.
2. The method for selecting spiral retrospective reconstruction data according to claim 1, wherein comparing the overlapping images of the first image sequence and the second image sequence, and regarding a reconstruction phase corresponding to an image in the second image sequence that satisfies a requirement as an optimal reconstruction phase in the second cardiac cycle comprises:
comparing the overlapped images in the first image sequence and the second image sequence, and selecting at least one image of which the difference degree is smaller than a preset threshold value;
and judging whether the phase distance between the corresponding phase of the at least one image and the optimal reconstruction phase of the first cardiac cycle is within a preset phase threshold range, and if so, taking the corresponding phase as the optimal reconstruction phase in the second cardiac cycle.
3. The method for selecting cardiac helical review reconstruction data according to claim 1, wherein the determining the optimal reconstruction phase for the location of interest includes:
reconstructing images corresponding to all phases within the range of 0% -99% of the concerned position at intervals according to a preset phase;
and calculating the image difference of all the images at the same attention position in all the phase ranges, and determining the phase corresponding to the image with the minimum image difference as the optimal reconstruction phase.
4. The method of claim 1, wherein the determining the optimal reconstruction phase for the location of interest comprises:
calculating a target phase range to be reconstructed according to a rule between the heart rate and the reconstructed target phase;
and evaluating the definition of all images in the target phase range by adopting an image information entropy calculation method, and taking the phase of the image with the highest definition as the optimal reconstruction phase.
5. A cardiac spiral retrospective reconstruction method is characterized in that,
selecting reconstructed data based on the cardiac spiral review reconstructed data selecting method of any one of claims 1 to 4;
and reconstructing the reconstruction data to obtain a heart reconstruction image.
6. A cardiac spiral retrospective reconstruction method as set forth in claim 5, further including:
in the process of carrying out data selection reconstruction according to the optimal reconstruction phase in each cardiac cycle to obtain a heart reconstruction image sequence, if reconstruction images obtained from data selected in any two adjacent cardiac cycles have overlapping areas with different degrees, difference points between two groups of reconstruction images are calculated, and the difference points are processed according to a motion deformation field to eliminate motion artifacts generated in the process of heart spiral review reconstruction.
7. A cardiac spiral review reconstruction data culling apparatus, comprising:
the device comprises a focus position determining module, a focus position determining module and a focus position determining module, wherein the focus position determining module is used for determining a focus position needing image reconstruction, and the focus position is a position where an artifact is likely to occur;
a first optimal reconstruction phase determining module, configured to determine an optimal reconstruction phase of the location of interest, where the optimal reconstruction phase is a cardiac motion phase with a smallest cardiac motion amount among the cardiac motion phases;
the first image sequence reconstruction module is used for determining a first cardiac cycle in which the optimal reconstruction phase is positioned and reconstructing a first image sequence of the optimal reconstruction phase in the first cardiac cycle;
the second image sequence reconstruction module is used for reconstructing images in a second cardiac cycle within a preset phase range to obtain a second image sequence, wherein the second cardiac cycle comprises a previous cardiac cycle or a next cardiac cycle which is adjacent to the first cardiac cycle;
a second optimal reconstruction phase determining module, configured to compare the first image sequence with an overlapped image of the second image sequence, and use a reconstruction phase corresponding to an image in the second image sequence that meets the requirement as an optimal reconstruction phase in a second cardiac cycle;
and the expansion module is used for taking the second cardiac cycle as the first cardiac cycle and repeatedly executing the first optimal reconstruction phase determining module, the first image sequence reconstructing module, the second image sequence reconstructing module and the second optimal reconstruction phase determining module until the optimal reconstruction phase in each cardiac cycle is determined.
8. The cardiac spiral review reconstruction data culling apparatus as set forth in claim 7, wherein the first optimal reconstruction phase determining module includes:
the image reconstruction unit is used for reconstructing images corresponding to all periods within the range of 0% -99% of the concerned position at intervals according to a preset period;
and the optimal reconstruction phase determining unit is used for calculating the image difference of all the images at the same attention position in all the phase ranges and determining the phase corresponding to the image with the minimum image difference as the optimal reconstruction phase.
9. A computer device comprising a memory and a processor, the memory storing a computer program operable on the processor, wherein the processor when executing the computer program implements the method of any one of claims 1 to 4 or the method of claim 5.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the method of any one of claims 1 to 4 or the method of claim 5.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210667055.2A CN115137383B (en) | 2022-06-13 | 2022-06-13 | Method and device for selecting and retrospective reconstruction data of heart spiral retrospective reconstruction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210667055.2A CN115137383B (en) | 2022-06-13 | 2022-06-13 | Method and device for selecting and retrospective reconstruction data of heart spiral retrospective reconstruction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115137383A CN115137383A (en) | 2022-10-04 |
CN115137383B true CN115137383B (en) | 2022-12-23 |
Family
ID=83408033
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210667055.2A Active CN115137383B (en) | 2022-06-13 | 2022-06-13 | Method and device for selecting and retrospective reconstruction data of heart spiral retrospective reconstruction |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115137383B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104545973A (en) * | 2014-11-26 | 2015-04-29 | 沈阳东软医疗***有限公司 | CT (Computed Tomography) scanning method and device for heart |
CN110728730A (en) * | 2019-09-30 | 2020-01-24 | 东软医疗***股份有限公司 | Image reconstruction method and device, CT (computed tomography) equipment and CT system |
CN112603338A (en) * | 2020-12-02 | 2021-04-06 | 赛诺威盛科技(北京)有限公司 | Method and device for selecting and retrospective reconstruction data of heart spiral retrospective reconstruction |
CN113269845A (en) * | 2021-04-26 | 2021-08-17 | 上海东软医疗科技有限公司 | Image reconstruction method and device, storage medium and electronic equipment |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8788020B2 (en) * | 1998-10-23 | 2014-07-22 | Varian Medical Systems, Inc. | Method and system for radiation application |
US20040019275A1 (en) * | 2002-07-23 | 2004-01-29 | Maria Iatrou | Method and apparatus for positioning a CT reconstruction window |
JP4731476B2 (en) * | 2003-07-16 | 2011-07-27 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Reconstruction of moving object images from volumetric data |
DE102010024684B8 (en) * | 2010-06-23 | 2019-08-22 | Siemens Healthcare Gmbh | Method for reconstructing image data of a moving examination object, control and processing unit, computer tomography system and computer program |
CN103596503B (en) * | 2011-06-20 | 2015-10-14 | 株式会社日立医疗器械 | X ray CT device and best heart contracting interval determining method |
US10475216B2 (en) * | 2018-02-26 | 2019-11-12 | General Electric Company | Imaging system and method using learned phase acquisition to acquire images |
US10726588B2 (en) * | 2018-10-31 | 2020-07-28 | Biosense Webster (Israel) Ltd. | Magnetic resonance imaging (MRI) image filtration according to different cardiac rhythms |
-
2022
- 2022-06-13 CN CN202210667055.2A patent/CN115137383B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104545973A (en) * | 2014-11-26 | 2015-04-29 | 沈阳东软医疗***有限公司 | CT (Computed Tomography) scanning method and device for heart |
CN110728730A (en) * | 2019-09-30 | 2020-01-24 | 东软医疗***股份有限公司 | Image reconstruction method and device, CT (computed tomography) equipment and CT system |
CN112603338A (en) * | 2020-12-02 | 2021-04-06 | 赛诺威盛科技(北京)有限公司 | Method and device for selecting and retrospective reconstruction data of heart spiral retrospective reconstruction |
CN113269845A (en) * | 2021-04-26 | 2021-08-17 | 上海东软医疗科技有限公司 | Image reconstruction method and device, storage medium and electronic equipment |
Also Published As
Publication number | Publication date |
---|---|
CN115137383A (en) | 2022-10-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20210390696A1 (en) | Medical image processing apparatus, medical image processing method and computer-readable storage medium | |
US7277565B2 (en) | Method of reconstructing a high-resolution 3D image | |
CN110246137B (en) | Imaging method, imaging device and storage medium | |
NL1024314C2 (en) | Integrated image registration for cardiac magnetic resonance blood flow data. | |
US11922601B2 (en) | Medical image processing apparatus, medical image processing method and computer-readable medium | |
IL158566A (en) | Method and apparatus for medical intervention procedure planning | |
JP7000079B2 (en) | Medical image processing equipment, medical image processing program and X-ray computer tomographic imaging equipment | |
US20100128941A1 (en) | Method and System of Motion Artefact Compensation in a Subject | |
CN110969633B (en) | Automatic optimal phase identification method for cardiac CT imaging | |
Xia et al. | Recovering from missing data in population imaging–Cardiac MR image imputation via conditional generative adversarial nets | |
US10854329B2 (en) | Motion estimation method and apparatus | |
Lee et al. | Synthesis of electrocardiogram V-lead signals from limb-lead measurement using R-peak aligned generative adversarial network | |
US11270434B2 (en) | Motion correction for medical image data | |
CN110728730A (en) | Image reconstruction method and device, CT (computed tomography) equipment and CT system | |
CN115137383B (en) | Method and device for selecting and retrospective reconstruction data of heart spiral retrospective reconstruction | |
EP1421556A2 (en) | Method of and system for the automatic registration of anatomically corresponding positions for perfusion measurements | |
EP3658031B1 (en) | Motion compensated cardiac valve reconstruction | |
JP3647970B2 (en) | Region extraction device | |
US20220292641A1 (en) | Dynamic imaging and motion artifact reduction through deep learning | |
CN116168099A (en) | Medical image reconstruction method and device and nonvolatile storage medium | |
CN115861172A (en) | Wall motion estimation method and device based on self-adaptive regularized optical flow model | |
EP4084010A1 (en) | Method for operating an evaluation system for medical image data sets, evaluation system, computer program and electronically readable storage medium | |
WO2022096867A1 (en) | Image processing of intravascular ultrasound images | |
Chittajallu et al. | Image-based methods for phase estimation, gating, and temporal superresolution of cardiac ultrasound | |
CN112258596A (en) | Image generation method and device, console equipment and CT system |
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 |