CN102949201A - Method for phase correlated tomographic image reconstruction from data measured by medical imaging, involves performing phase-correlated iterative image reconstruction which is initialized with non-phase related image volume - Google Patents

Method for phase correlated tomographic image reconstruction from data measured by medical imaging, involves performing phase-correlated iterative image reconstruction which is initialized with non-phase related image volume Download PDF

Info

Publication number
CN102949201A
CN102949201A CN2012102846487A CN201210284648A CN102949201A CN 102949201 A CN102949201 A CN 102949201A CN 2012102846487 A CN2012102846487 A CN 2012102846487A CN 201210284648 A CN201210284648 A CN 201210284648A CN 102949201 A CN102949201 A CN 102949201A
Authority
CN
China
Prior art keywords
measurement data
image reconstruction
phase
reconstruction
time
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.)
Pending
Application number
CN2012102846487A
Other languages
Chinese (zh)
Inventor
H.布鲁德
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens AG
Original Assignee
Siemens AG
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Siemens AG filed Critical Siemens AG
Publication of CN102949201A publication Critical patent/CN102949201A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • 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
    • G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/412—Dynamic
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/424—Iterative

Abstract

The method involves reconstructing (6) a non-phase-related image volume using a portion or all data. The phase-correlated iterative image reconstruction which is initialized with the non-phase related image volume is performed (8). The reconstructed image volume is calculated by using specific equation relating the phase-correlated measurement data, settling term for stabilization of iteration and the variable corresponding to the projector. An independent claim is included for a device for phase correlated tomographic image reconstruction from data measured by medical imaging.

Description

Be used for medical imaging with the time method and apparatus of image reconstruction relatively mutually
Technical field
The present invention relates to a kind of for from the measurement data of the tomography of medical imaging with the time method of (phasenkorrelierten) image reconstruction relatively mutually, this measurement data is comprising the time period record of a plurality of cardiac cycles.The invention still further relates to a kind of device for carrying out the method for advising.The method and device are suitable for CT (computer tomography) (CT), particularly heart CT(Kardio-CT) in image reconstruction.
Background technology
In heart CT, use at present following acquisition scheme and method for reconstructing: perhaps by the triggering that can expect or by retrieve data select (gate) realize about the time demonstration heart mutually.Only recording measurement data during mutually during the specific heart at each cardiac cycle in the situation of the triggering that can expect.In all the other time periods, do not move or only with the Power operation X-ray tube of obvious reduction.Produce trigger impulse by the EKG signal.Can in the image of rebuilding, avoid because the motion artifacts that heart movement causes thus.In principle can be about whole cardiac cycle record measurement data in the situation that retrieve data is selected, so as to be chosen in afterwards each of cardiac cycle identical the time mutually in the measurement data that recorded of institute and be the image reconstruction introducing.
When fast frame rotation and/or when the employing two-pipe system, reach high temporal resolution.But in the situation of higher heart rate in the triggering that can expect the high galvanic areas of X-ray tube must expand the relevant time range of diagnosis because phase must adjust the heart of the best herein the time with recalling.For functional check, also need irradiation at heart contracting end-systole (Endsystole) in addition.Select the same so-called overshoot that needs with retrieve data in both cases, this overshoot is not exclusively by clinical use.But use the radiation that all applies in order to improve quantum statistics expectation in the image of rebuilding, in order to use best the dosage that applies.If can use the radiation that applies fully, then can reduce patient dose in constant contrast in to the situation of noise (Kontrast-zu-Rauschen).
In order to improve the utilization of the x-ray dose that applies, can registration when different the measurement data of phase.During for this reason at first for different heart mutually respectively about the time reconstructed image and then determine sports ground by the 3D-3D registration mutually.Then this sports ground is found out by measurement data motion compensation ground and is carried out beginning in the image reconstruction.But by about the time image rebuild mutually determine that sports ground must be high-precision, to guarantee the accuracy registration of heart coronary artery.
Summary of the invention
The technical problem to be solved in the present invention is, be provided for from the measurement data of the tomography of medical imaging with the time mutually relatively a kind of method and a kind of device of image reconstruction, the method need not to determine that sports ground just can use the dosage that applies fully.
Above-mentioned technical problem is by solving according to method and apparatus of the present invention.Preferred embodiment being the content of dependent claims or can from subsequently description and embodiment, drawing of the method and device.
In the method for advising, according to about the measurement data of the tomography of the time period record that comprises a plurality of cardiac cycles at first by whole or at least a portion measurement data about a plurality of cardiac cycles rebuild not with the time mutually relevant 3D rendering or image volume f StdImage reconstruction is undertaken by suitable back projection with known method.Then rebuild iteratively with the time mutually relevant image, this image pass through not with the time mutually relevant image volume f StdInitialize.This iteration with the time in k iterative step, recomputate the image volume f of reconstruction in the mutually relevant image reconstruction according to following equation (renewal equation) k:
f k+1=f k+X -1 pc*(p-X pc*f k)+R(f k)
X wherein -1 PcCorresponding to the back projection to image reconstruction, X PcCorresponding to the projection of correspondence, p when being used for mutually correct back projection with the time mutually relevant measurement data, and R(f k) corresponding to the regularization term that is used for stablizing iteration.When being used for mutually correct back projection with the time relevant measurement data p be respectively in measurement data during mutually during same cardiac at each cardiac cycle of identical time period.This measurement data is thus also corresponding to the measurement data that is used in the situation about selecting in retrieve data rebuilding.Phase or identical time period in the time of can deriving identical heart by the EKG triggering signal that during gathering measurement data, records simultaneously in common mode equally.Then draw and data that further processed before reconstructed image in case of necessity rebuild 3D rendering or image volume according to these measurement data during this measurement data is understood to be in the execution imaging.
In the method for advising, after the iterative step of fixing predetermined quantity or meeting or exceeding can finishing iteration after the consistent degree between the image volume of the reconstruction of the image volume of the reconstruction of current iteration step and previous iteration step.After iteration finishes the display demonstration with the time mutually relevant reconstruction image volume.
The method of passing through to advise based on the time mutually relevant image reconstruction reduce because the motion artifacts that heart movement causes, but use simultaneously the measurement data of whole records or therefrom selection portion assign to reduce picture noise.Realize thus with the time mutually relevant image reconstruction, enough and however this image reconstruction need not accurately to determine that sports ground is with regard to the but best radiation that applies of use.
The method of advising preferred embodiment in, the high pass filter that is used for stablizing iteration is used as regularization term.In the embodiment of replacing, also can with TV(Total Variation, change fully) the formal representation regularization of Norm minimum, this regularization also use not with the time mutually relevant image volume f StdIn each iterative step, make cost function at this
Figure BDA00001996188500031
0＜α＜1 minimizes.At this, parameter alpha represents weighting factor.L1 represents the L1 norm of vector.In this embodiment will thus R ( f k ) = α · | ▿ f std | L 1 + ( 1 - α ) · | ▿ f k | L 1 Be set to above-mentioned renewal equation.
The device of advising comprises be used to the storage device that is stored in the measurement data that records on a plurality of cardiac cycles and computation processor, utilizes this computation processor to be implemented as the required calculation procedure for image reconstruction of method.This device also comprises reconstruction unit, and this reconstruction unit is constructed like this so that its at first from the measurement data of storage or be enough to rebuild most of measurement data of noise decrease subsequently not with the time mutually relevant image volume f StdReconstruction unit is also constructed like this so that its carry out iteratively subsequently with the time mutually relevant reconstructed image, this image reconstruction by not with the time mutually relevant image volume f StdInitialize and wherein in k iterative step, recomputate according to above-mentioned renewal equation the image volume f of reconstruction kPreferably, this device is connected with display, on this display after iteration finishes demonstration with the time mutually relevant reconstruction image volume.
The method both can be used for only having the CT system an of camera system (i.e. an X-ray tube and an X-ray detector), also can be for so-called multi radiator system.The method demonstration also is applicable to other CT system or C C-arm system C.
Description of drawings
The below is described further with relevant device the method for advising by reference to the accompanying drawings according to embodiment.In the accompanying drawing:
Fig. 1 shows in heart CT the example as the tube current of the function of time;
Fig. 2 shows the sketch map of method flow in the method for advising; With
Fig. 3 shows the sketch map of device.
The specific embodiment
In heart CT, recalling in the situation of gate record about the measurement data of a plurality of cardiac cycles.Fig. 1 there is shown the current I about the X-ray tube of complete cardiac cycle on a left side.But for the time mutually relevant image reconstruction only use in many cases measurement data near time point 1 the zone 2, should CT image reconstruction at this time point 1 place.All the other zones 3 of measurement data are not used in image reconstruction.High galvanic areas can closely be limited in for the important window (Fenster) of diagnosis in the situation of the triggering that can expect.But must further expand high galvanic areas for higher heart rate or in the situation that needs additional function information.This point causes again having larger time range in the measurement data that is not used.
In the method for advising, in the process of reconstruction of iteration, optimally or best use the x-ray dose that applies for gathering measurement data.This point is delineated in the right figure of Fig. 1, wherein for the image reconstruction of time point 1, uses total zone 4 of measurement data in cardiac cycle.This point obviously is applicable to whole cardiac cycles, about this cardiac cycle measurement data is carried out record.
In the method for advising, after gathering 5 measurement data at first not with the time mutually relevant reconstruction 6 in generation not with the time image volume f that is correlated with Std, the below is also referred to as formerly image volume (Prior-Bildvolumen).This is image volume f formerly StdAccumulated the dosage that all applies, namely considered the time range of measurement data 180 ° (with respect to the frame rotation).Phase when this scope of data comprises whole high electric current for the X-ray tube of each cardiac cycle under extreme case.Then image volume f formerly StdInitialization about the time phase iterative approximation 8, select 7 measurement data p to be incorporated into this iterative approximation.By correction term p-X Pc* f kWith the time mutually relevant reconstruction reduce motion artifacts in the image formerly, but simultaneously formerly the little noise transmission of image in the update image of iteration.In k iterative step, set up following renewal equation:
f k+1=f k+X -1 pc*(p–X pc*f k)+R(f k)
f kThe image volume of rebuilding in step of updating k has been described.Projection X PcOr the X of back projection -1 PcGuaranteed forward projection or the reconstruction of phase when correct, and p represents the scope of measurement data, this scope when being mutually correct cardiac reconstruction required.By image volume f formerly StdInitialize iteration.R(f k) represent for stablizing the high pass filter of iteration (regularization).Regularization is associated with the gray value of pixel adjacent.
After iteration finishes, finishing iteration and in image volume or the CT image (step 9) among Fig. 2 of the current reconstruction of display demonstration the iterative step of fixed qty after or in the situation about relatively measuring of the image volume of the image volume that surpasses each reconstruction and previous iteration step for example.Schematically method flow can be found out by Fig. 2.
Fig. 3 has schematically shown the exemplary structure of this device or data processing equipment strongly, has for the storage device 10 of storing the measurement data that records, utilizes for example software module of its computation processor 11 that is implemented as the required calculation procedure of image reconstruction and reconstruction unit 12().The image reconstruction that reconstruction unit 12 is carried out according to the method for advising in the situation of using memory element and computation processor.Then the image of rebuilding can show at the display screen 13 that connects.
Although the present invention is further described and description by preferred embodiment, the present invention is not subjected to the restriction of disclosed example and can derives other by the professional and be out of shape, and does not break away from protection scope of the present invention.

Claims (5)

  1. One kind be used for from the measurement data of the tomography of medical imaging with the time method of image reconstruction relatively mutually, this measurement data is comprising the time period record of a plurality of cardiac cycles, wherein
    At first under using the situation of part or all of measurement data, rebuild not with the time mutually relevant image volume f Std(6), and
    Then carry out iteratively with the time mutually relevant image reconstruction (8), by described not with the time mutually relevant image volume f StdInitialize this image reconstruction, and wherein in k iterative step, recomputate the image volume f of reconstruction according to following equation k:
    f k+1=f k+X -1 pc*(p–X pc*f k)+R(f k)
    X wherein -1 PcCorresponding to the back projection to image reconstruction, X PcCorresponding to the projection of correspondence, p when being used for mutually correct back projection with the time mutually relevant measurement data, and R(f k) corresponding to the regularization term that is used for stablizing iteration.
  2. 2. method according to claim 1 is characterized in that, high pass filter is used as regularization term.
  3. 3. method according to claim 1 is characterized in that, for TV Norm minimum technology is used in regularization.
  4. 4. each described method in 3 according to claim 1 is characterized in that, minimizes cost function in each iterative step
    Figure FDA00001996188400011
    0＜α≤1 wherein.
  5. One kind be used for from the measurement data of the tomography of medical imaging with the time device of image reconstruction relatively mutually, comprising having time period this measurement data of record of a plurality of cardiac cycles
    -storage device (10) is used for storage of measurement data,
    -computation processor (11) utilizes this computation processor to be implemented as the required calculation procedure of image reconstruction, and
    -reconstruction unit (12), this reconstruction unit is constructed like this so that its at first under using the situation of the part or all of measurement data of storage, rebuild not with the time mutually relevant image volume f Std(6) and subsequently according to the described method of the claims carry out iteratively with the time mutually relevant image reconstruction (8).
CN2012102846487A 2011-08-10 2012-08-10 Method for phase correlated tomographic image reconstruction from data measured by medical imaging, involves performing phase-correlated iterative image reconstruction which is initialized with non-phase related image volume Pending CN102949201A (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102011080733.0 2011-08-10
DE201110080733 DE102011080733A1 (en) 2011-08-10 2011-08-10 Method for phase correlated tomographic image reconstruction from data measured by medical imaging, involves performing phase-correlated iterative image reconstruction which is initialized with non-phase related image volume

Publications (1)

Publication Number Publication Date
CN102949201A true CN102949201A (en) 2013-03-06

Family

ID=47595501

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012102846487A Pending CN102949201A (en) 2011-08-10 2012-08-10 Method for phase correlated tomographic image reconstruction from data measured by medical imaging, involves performing phase-correlated iterative image reconstruction which is initialized with non-phase related image volume

Country Status (2)

Country Link
CN (1) CN102949201A (en)
DE (1) DE102011080733A1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103310472A (en) * 2013-06-21 2013-09-18 中国科学院自动化研究所 Limited angle photoacoustic imaging reconstruction method and device on basis of regularization iteration
CN107862722A (en) * 2016-09-21 2018-03-30 株式会社岛津制作所 Successive approximation image rebuilding method, successive approximation image reconstruction program and laminagraph device
CN109523602A (en) * 2018-10-22 2019-03-26 上海联影医疗科技有限公司 Iterative approximation optimization method, device and computer-readable medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101305921A (en) * 2007-05-17 2008-11-19 上海西门子医疗器械有限公司 Method and device for implement CT respiration gate control
US20100128958A1 (en) * 2008-11-26 2010-05-27 Guang-Hong Chen Method for prior image constrained image reconstruction in cardiac cone beam computed tomography
CN102100565A (en) * 2009-12-17 2011-06-22 通用电气公司 System and method to correct motion in gated-pet images using non-rigid registration

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101305921A (en) * 2007-05-17 2008-11-19 上海西门子医疗器械有限公司 Method and device for implement CT respiration gate control
US20100128958A1 (en) * 2008-11-26 2010-05-27 Guang-Hong Chen Method for prior image constrained image reconstruction in cardiac cone beam computed tomography
CN102100565A (en) * 2009-12-17 2011-06-22 通用电气公司 System and method to correct motion in gated-pet images using non-rigid registration

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
GUANG-HONG CHEN 等: "High temporal resolution cardiac cone-beam CT using a slowly rotating C-arm gantry", 《MEDICAL IMAGING 2009： PHYSICS OF MEDICAL IMAGING, PROC OF SPIE》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103310472A (en) * 2013-06-21 2013-09-18 中国科学院自动化研究所 Limited angle photoacoustic imaging reconstruction method and device on basis of regularization iteration
CN103310472B (en) * 2013-06-21 2016-01-06 中国科学院自动化研究所 Based on limited angle photoacoustic imaging method for reconstructing and the device of regularization iteration
CN107862722A (en) * 2016-09-21 2018-03-30 株式会社岛津制作所 Successive approximation image rebuilding method, successive approximation image reconstruction program and laminagraph device
CN109523602A (en) * 2018-10-22 2019-03-26 上海联影医疗科技有限公司 Iterative approximation optimization method, device and computer-readable medium
CN109523602B (en) * 2018-10-22 2023-02-17 上海联影医疗科技股份有限公司 Iterative reconstruction optimization method, device and computer readable medium

Also Published As

Publication number Publication date
DE102011080733A1 (en) 2013-02-14

Similar Documents

Publication Publication Date Title
US8421021B2 (en) Motion correction of SPECT images
EP2174294B1 (en) Motion correction in nuclear imaging
US7558439B2 (en) Motion artifact correction of tomographical images
JP5734664B2 (en) Image Restoration Method Using Dilution Constraint Correction
JP5172347B2 (en) Reconstruction of 2D planar images of nuclear medicine by iterative constraint deconvolution
KR101636041B1 (en) Tomography apparatus and method for reconstructing a tomography image thereof
US10213179B2 (en) Tomography apparatus and method of reconstructing tomography image
CN110168612B (en) Normalized uptake value directed reconstruction control for improved robustness of results in positron emission tomography imaging
CN111540025B (en) Predicting images for image processing
JP2011503570A (en) Apparatus and method for forming an attenuation map
CN104424659A (en) Image-based motion compensation of image data
CN103156629A (en) Image processing apparatus, image processing method, and program
CN102376084B (en) Use anisotropic noise model to the iterative image filtering of CT image
WO2017132465A1 (en) System and method for denoising medical images by enforcing low rank spatial-temporal or spatial-spectral image matrices
CN106462987B (en) The improvement image reconstruction for volume based on data for projection collection
US9953440B2 (en) Method for tomographic reconstruction
CN102949201A (en) Method for phase correlated tomographic image reconstruction from data measured by medical imaging, involves performing phase-correlated iterative image reconstruction which is initialized with non-phase related image volume
JP6143533B2 (en) Nuclear medicine image reconstruction device, nuclear medicine image reconstruction method, and program
JP7382496B2 (en) Continuous bed movement acquisition using short-axis phantom for PET imaging system setup and quality control
JP4360601B2 (en) Image reconstruction processing program and recording medium
US11663758B2 (en) Systems and methods for motion estimation in PET imaging using AI image reconstructions
CN102599927B (en) The X-ray tomographic of Moving Objects is rebuild
Taubmann et al. Background Matters: A Correction Scheme for Dynamic Iterative CBCT with Limited Grid Size
CN115038382A (en) PET imaging using multi-organ specific short CT scans
Verhaeghe et al. PET projection data supersets for reconstruction with acquisition motion

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130306