Nothing Special   »   [go: up one dir, main page]

CN105662413B - A kind of method and apparatus that cardiac muscle is quantitative T1 - Google Patents

A kind of method and apparatus that cardiac muscle is quantitative T1 Download PDF

Info

Publication number
CN105662413B
CN105662413B CN201511033083.5A CN201511033083A CN105662413B CN 105662413 B CN105662413 B CN 105662413B CN 201511033083 A CN201511033083 A CN 201511033083A CN 105662413 B CN105662413 B CN 105662413B
Authority
CN
China
Prior art keywords
sampling line
image
line
rebuild
sampling
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
Application number
CN201511033083.5A
Other languages
Chinese (zh)
Other versions
CN105662413A (en
Inventor
刘新
朱燕杰
钟耀祖
郑海荣
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.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
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 Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201511033083.5A priority Critical patent/CN105662413B/en
Publication of CN105662413A publication Critical patent/CN105662413A/en
Application granted granted Critical
Publication of CN105662413B publication Critical patent/CN105662413B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A61B5/7207Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Signal Processing (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Radiology & Medical Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

This application discloses a kind of methods that cardiac muscle is quantitative T1, including:After ecg-gating trigger delay, apply non-layer choosing inversion pulse;At least two layers or more of image staggeredly acquire in real time using the Fast spoiled gradient echo low-angle excitation sequence of radial sample track, captures the inversion recovery process of signal;Select the sampling line in diastole as K space center line;Centered on the sampling line selected, sampling line is symmetrically selected to carry out image reconstruction according to window size is rebuild;Quantitatively schemed with the image fitting T1 reconstructed.Disclosed herein as well is a kind of devices based on this method.The application may be implemented in acquisition multilayer T1 quantitative images in single-breath hold, and the discomfort for whole-heartedly reducing the hold one's breath time waste and patient that bring is covered in can holding one's breath at 2-3 times.

Description

A kind of method and apparatus that cardiac muscle is quantitative T1
Technical field
A kind of the application angiocarpy imaging field, and in particular to method and apparatus that cardiac muscle is quantitative T1.
Background technology
Cardiovascular magnetic resonance imaging (CMR) is a kind of imaging technique of Noninvasive, function, shape that can be to cardiovascular system State and structure etc. are assessed.The quantization of myocardium T1 values can be used in the diagnosis of a variety of diseases such as heart infarction, myocardial fibrosis.Due to The limitation of heartbeat and respiratory movement to acquisition time, in body-centered, the T1 of flesh is quantitative still has prodigious challenge.
The common method for measuring cardiac muscle T1 values has following several:
(1)MOLLI(Modified Look-Locker inversion recovery):After inversion pulse (IR), not Same delay (TI) time bssfp (balance steady state free precession) imaging sequences, and will be more The data of secondary acquisition merge into one group, to be fitted T1 values.
(2)SASHA(Saturation recovery single-shot acquisition):It is replaced with saturation pulse Inversion pulse is completed to adopt with bssfp imaging sequences through different delay times after saturation pulse in 10 heart beat cycles Collection.
In the above method, single layer image can only be acquired every time by holding one's breath, and multi-layer image could be acquired by needing repeatedly to hold one's breath, and be held one's breath Often, acquisition time is long, and patient comfort is low.
Invention content
The application provides a kind of method of cardiac muscle T1 value quantitative measurments, including:
After ecg-gating trigger delay, apply non-layer choosing inversion pulse, mutually ladder is quickly disturbed using radial sample track The image for spending echo low-angle excitation sequence at least two layers or more staggeredly acquire in real time, captures the inversion recovery mistake of signal Journey;Select the sampling line in diastole as K space center line;Centered on selected sampling line, according to Rebuilding window size symmetrically selects sampling line to carry out image reconstruction;Quantitatively schemed with the image fitting T1 reconstructed.
It is above-mentioned to apply non-layer choosing inversion pulse after ecg-gating trigger delay, using quickly disturbing for radial sample track Phase gradient echo low-angle excitation sequence at least two layers or more of image staggeredly acquire in real time, and the reversion for capturing signal is extensive Multiple process, including:The inversion pulse for applying non-layer choosing excites sequence with the Fast spoiled gradient echo low-angle of radial sample track Row are acquired in real time, are sampled to the inversion recovery process of signal, and the phase gradient of quickly disturbing of the radial direction sample track returns Wave low-angle excitation sequence uses multi-layer intercrossed acquisition mode, and at least two layers or more of image can be acquired after an inversion pulse.
Above-mentioned radial direction sample track uses small gold angle sampling configuration, and the azimuth of i-th sampling line is (i-1) ψN The azimuth of degree, i.e., every sampling line increases a fixed angle ψN, whereinN is determined ψNSize, it is empirically determined.Sampling line number mesh must be more than N.
The above-mentioned Fast spoiled gradient echo low-angle excitation sequence with radial sample track acquire in real time and is being breathed It is carried out during holding one's breath, the most short TE and TR of use.
It is above-mentioned to select the sampling line in diastole as K space center line, including:According to sliding window size, Acquisition time is selected retrospectively is located at the sampling line of diastole as the K space center line rebuild.
Centered on selected sampling line, sampling line is symmetrically selected to carry out image weight according to window size is rebuild Build including:Sampling line using the acquisition time positioned at diastole is rebuild as the K space center line rebuild to determine The position of window is weighted K space data using KWIC methods the reconstruction window, and quick algorithm for reconstructing pair can be used Data are rebuild.
According to the second aspect of the application, a kind of device that cardiac muscle is quantitative T1 is provided, including:Data acquisition module is used for After ecg-gating trigger delay, apply non-layer choosing inversion pulse, the Fast spoiled gradient echo using radial sample track is small Angle excitation sequence at least two layers or more of image staggeredly acquire in real time, captures the inversion recovery process of signal;Selection Module, for selecting the sampling line in diastole as K space center line;Image reconstruction module, for institute Centered on stating the sampling line selected, sampling line is symmetrically selected to carry out image reconstruction according to window size is rebuild;Fitting module, For using the image fitting T1 reconstructed quantitatively to scheme.
Above-mentioned data acquisition module is additionally operable to apply the inversion pulse of non-layer choosing, and mutually ladder is quickly disturbed with radial sample track Degree echo low-angle excitation sequence is acquired in real time, is sampled to the inversion recovery process of signal, the radial sampling rail The Fast spoiled gradient echo low-angle excitation sequence of mark uses multi-layer intercrossed acquisition mode, can acquire after an inversion pulse to Few two layers or more of image.
Above-mentioned radial direction sample track uses small gold angle sampling configuration, and the azimuth of i-th sampling line is (i-1) ψN The azimuth of degree, i.e., every sampling line increases a fixed angle ψN, whereinN is determined ψNSize, for N by empirically determined, sampling line number mesh must be more than N.
Above-mentioned data acquisition module is additionally operable to carry out data acquisition, the most short TE and TR of use during breathing is held one's breath.
Above-mentioned selecting module is additionally operable to, according to sliding window size, select acquisition time retrospectively and be located at diastole The sampling line of phase is as the K space center line rebuild.
Above-mentioned image reconstruction module is additionally operable to be located at the sampling line of diastole as reconstruction using the acquisition time K space center line determine the position of reconstruction window, to the reconstruction window, K space data is added using KWIC methods Power, rebuilds data using quick algorithm for reconstructing.
As a result of above technical scheme, the advantageous effect that the application has is made to be:
In the specific implementation mode of the application, due to the use of small using the Fast spoiled gradient echo of radial sample track Angle excitation sequence and small gold angle sampling configuration carry out at least two layers or more of image after non-layer choosing inversion pulse In real time staggeredly acquisition, it can be achieved that in single-breath hold acquire multilayer T1 quantitative images, cover whole-heartedly, subtract in can holding one's breath at 2-3 times The discomfort of the hold one's breath time waste and patient that bring is lacked;Using KWIC methods be weighted K space data, Ke Yiming Really the time of acquisition K space center point, the deviation of measured value can be effectively prevent.
Description of the drawings
Fig. 1 is the flow chart according to the application method one embodiment;
Fig. 2 is the acquisition schematic diagram according to the application method one embodiment;
Fig. 3 is the radial sampled gradients oscillogram according to the application method one embodiment;
Fig. 4 is the sample track figure according to the application method one embodiment;
Fig. 5 is to be chosen to scheme according to the diastole of the application method one embodiment;
Fig. 6 is to weight weight map according to the KWIC of the application method one embodiment;
Fig. 7 is the structural schematic diagram according to the application device one embodiment.
Specific implementation mode
Below by specific implementation mode combination attached drawing, invention is further described in detail.
Embodiment one:
Fig. 1 shows the flow chart according to the application method one embodiment, including:
Step 102:After ecg-gating trigger delay, apply non-layer choosing inversion pulse, using the fast of radial sample track Speed disturbs phase gradient echo low-angle excitation sequence and at least two layers or more of image staggeredly acquire in real time, captures the anti-of signal Turn recovery process.
Apply the inversion pulse (IR) of non-layer choosing, is then swashed with the Fast spoiled gradient echo low-angle of radial sample track Hair sequence (turbo-flash) is acquired in real time, is sampled to the inversion recovery process of signal, and turbo-flash is used Multi-layer intercrossed acquisition mode (interleaved acquisition) can acquire multi-layer image, such as Fig. 2 after an inversion pulse It is shown.
Radial sampled gradients oscillogram is as shown in Figure 3.In radial sampling, every sampling line can all pass through K space center.It is small The radial sample track of two dimension of gold angle is that the deflection of every sampling line increases a fixed angle ψN,It will not be all overlapped with every sampling line of this sampling configuration, and as sampling line number Mu > When N, sampling line distribution is similar to be uniformly distributed, and can rebuild arbitrary sampling line as initial position, also in Chuan Kou great little > Arbitrarily window size, a kind of embodiment N=5 are rebuild in setting in the range of=N, as shown in figure 4, radial direction when it is N=5 is small The sample track of gold angle.
The trigger delay of ecg-gating is set as 60% (about 500~650ms) of patient's heart beat cycle (RR), ecg-gating Delay time be to ensure that initial inversion recovery signal acquires within heartbeat diastole, therefore delay time is by being imaged What the heart beat cycle of object determined.Entire gatherer process is completed during breathing is held one's breath.To improve the temporal resolution of imaging, Using most short TE (time of repetition) and TR (time of echo), i.e. TR/TE is by the hard of magnetic resonance scanner Part limits the shortest time that (greatest gradient amplitude, gradient rate of climb etc.) determines, TR increases with the increase of the acquisition number of plies, has Body is the number of plies that TR=acquires one layer most short TR × acquisition, and flip angle is 5 degree.
Step 104:Select the sampling line in diastole as K space center line.
Data acquisition carries out in entire heart beat cycle, but due to heart movement, and directly carrying out image reconstruction can not be direct Image is fitted point by point.Heart is run slowly in diastole, can be selected acquisition time retrospectively and is located at diastole The sampling line of phase is fitted for subsequent image reconstruction and data.Choosing method is as follows:When gathered data, initial data head text As starting point at the time of being triggered using current electrocardio in part, current acquisition time is recorded, it, will when next electrocardio is triggered to next The time sets to 0, then maximum time is the time of heart beat cycle, is diastole between heart beat cycle 65%-95%, when selection acquires Gathering line within that range is carved, as effective gathering line, as shown in Figure 5.The effectively acquisition of first of each heart beat cycle Line determines the position of subsequent K space center line, i.e., the adjacent spaces K according to sliding window size as first K space center line The spacing of center line is sliding window size.
Step 106:Centered on selected sampling line, according to rebuild window size symmetrically select sample line into Row image reconstruction.
The sampling line positioned at center line-window size/2~center line+window size/2 is selected to rebuild present image;One Item sampling line can be used for rebuilding the image of different centers;For centreline space away from for sliding step, those skilled in the art should It is clear, it rebuilds window size and sliding step determines based on experience value.It is located at diastole from what is selected in step 104 Selection K space center line in effective gathering line, centreline space is using sliding step as spacing.To each reconstruction window, using KWIC (k-space weighted image contrast) method is weighted K space data., it is empty that KWIC can clearly acquire K Between center line time, Fig. 6 is KWIC weight maps when rebuilding window size=53, black=0 in figure, white=1.With reversion When measuring T1 values, if not using KWIC methods of weighting, measured value can be caused to go out with radial sampling acquisition in real time after pulse Existing deviation.
Step 108:Quantitatively schemed with the image fitting T1 reconstructed.
Image after reconstruction obtains T1 value images through fitting.Fitting formula isT is IR modules Terminate the time to acquisition each image k-space center line, i.e. t=TR × position of center line, obtains three parametersThen T1 value calculation formula areIt carries out curve fitting to the every bit of image spatial coordinates, obtains T1 value figures Picture.Due in the realization of actual sequence, needing that damage gradient (Spoiler) is added after inversion pulse to eliminate due to endless The remaining transverse magnetization vector that beautiful inversion pulse introduces causes have certain delay between inversion pulse and data acquisition, The T1 values measured are influenced, therefore error caused by this delay need to be corrected, formula is as follows:T1true=T1+ 2* Δ t, Δ t are Delay between IR and data acquisition, obtains the T1 value images finally measured.
Heart T1 quantitative measurments at present are all made of single-breath hold, acquire the mode of a tomographic image.Hold one's breath need operator and The interactive process of subject, causes waste of time, and repeatedly holding one's breath also can enable subject/patient feel exhausted.The application proposes Method single breath-hold can acquire the acquisition sequence and method of multilayer cardiac muscle T1 quantitative images, covered in can holding one's breath at 2~3 times complete The heart reduces repeatedly hold one's breath the time waste brought and patient's discomfort.
Embodiment two:
Fig. 7 is the structural schematic diagram according to the application device one embodiment, including:Data acquisition module, selecting module, Image reconstruction module and fitting module.
Data acquisition module, for after ecg-gating trigger delay, applying non-layer choosing inversion pulse, using radial sampling The Fast spoiled gradient echo low-angle excitation sequence of track at least two layers or more of image staggeredly acquire in real time, capture The inversion recovery process of signal.A kind of embodiment is additionally operable to apply the inversion pulse of non-layer choosing, fast with radial sample track Speed is disturbed phase gradient echo low-angle excitation sequence and is acquired in real time, is sampled to the inversion recovery process of signal.Radial direction is adopted The Fast spoiled gradient echo low-angle excitation sequence of sample track uses multi-layer intercrossed acquisition mode, can be adopted after an inversion pulse Collect at least two layers or more of image.A kind of embodiment, radial sample track are small gold angle sampling configuration, i-th sampling The azimuth of line is (i-1) ψNThe azimuth of degree, i.e., every sampling line increases a fixed angle ψN, whereinN determines ψNSize, for N by empirically determined, sampling line number mesh must be more than N.A kind of reality Mode is applied, data acquisition module is additionally operable to carry out data acquisition, the most short TE and TR of use during breathing is held one's breath.
Selecting module, for selecting the sampling line in diastole as K space center line.A kind of embodiment party Formula selects acquisition time and is located at the sampling line of diastole as the K space center line rebuild retrospectively.
Image reconstruction module, for centered on the sampling line selected, symmetrically selection to be adopted according to window size is rebuild Line-transect carries out image reconstruction.A kind of embodiment, the sampling line that diastole is located at using acquisition time are empty as the K rebuild Between center line determine the position of reconstruction window, to the reconstruction window, K space data is weighted using KWIC methods, is made Data are rebuild with quick algorithm for reconstructing.In one embodiment, the application can be used cg-SENSE algorithms and carry out weight It builds, the application can also be used other methods and rebuild.
Fitting module, for using the image reconstructed fitting T1 quantitatively to scheme.Those skilled in the art can use some maturations Method be fitted.
The above content is combining, specific embodiment is made for the present invention to be further described, and it cannot be said that this hair Bright specific implementation is confined to these explanations.For those of ordinary skill in the art to which the present invention belongs, it is not taking off Under the premise of from present inventive concept, a number of simple deductions or replacements can also be made.

Claims (8)

1. a kind of method that cardiac muscle is quantitative T1, which is characterized in that including:
After ecg-gating trigger delay, apply non-layer choosing inversion pulse, is returned using the phase gradient of quickly disturbing of radial sample track Wave low-angle excitation sequence at least two layers or more of image staggeredly acquire in real time, captures the inversion recovery process of signal;
Select the sampling line in diastole as K space center line;
Centered on selected sampling line, sampling line is symmetrically selected to carry out image reconstruction according to window size is rebuild;
Quantitatively schemed with the image fitting T1 reconstructed;
It is described to apply non-layer choosing inversion pulse after ecg-gating trigger delay, mutually ladder is quickly disturbed using radial sample track The image for spending echo low-angle excitation sequence at least two layers or more staggeredly acquire in real time, captures the inversion recovery mistake of signal Journey, including:
The inversion pulse for applying non-layer choosing is carried out real with the Fast spoiled gradient echo low-angle excitation sequence of radial sample track When acquire, the inversion recovery process of signal is sampled, it is described radial direction sample track Fast spoiled gradient echo low-angle Excitation sequence uses multi-layer intercrossed acquisition mode, and at least two layers or more of image can be acquired after an inversion pulse;
The wherein described radial sample track uses small gold angle sampling configuration, and the azimuth of i-th sampling line is (i-1) ψNDegree, The azimuth of i.e. every sampling line increases a fixed angle ψN, whereinN determines ψN Size, for N by empirically determined, sampling line number mesh must be more than N.
2. the method as described in claim 1, which is characterized in that the small angle of Fast spoiled gradient echo of the radial direction sample track Degree excitation sequence carries out acquisition in real time and is carried out during breathing is held one's breath, the most short TE and TR of use.
3. the method as described in claim 1, which is characterized in that described to select the sampling line in diastole as K skies Between center line include:
According to sliding window size, acquisition time is selected retrospectively and is located at the sampling line of diastole as the K rebuild Space center's line.
4. method as claimed in claim 3, which is characterized in that wherein centered on selected sampling line, according to weight Build window size symmetrically select sampling line carry out image reconstruction include:
Sampling line using the acquisition time positioned at diastole determines reconstruction window as the K space center line rebuild Position is weighted K space data using K- spatial weighting image comparison methods the reconstruction window, uses quick weight Algorithm is built to rebuild data.
5. a kind of device that cardiac muscle is quantitative T1, which is characterized in that including:
Data acquisition module, for after ecg-gating trigger delay, applying non-layer choosing inversion pulse, using radial sample track Fast spoiled gradient echo low-angle excitation sequence at least two layers or more of image staggeredly acquire in real time, capture signal Inversion recovery process;
Selecting module, for selecting the sampling line in diastole as K space center line;
Image reconstruction module, for centered on selected sampling line, symmetrically selection to be adopted according to window size is rebuild Line-transect carries out image reconstruction;
Fitting module, the image fitting T1 for being reconstructed described in use quantitatively scheme;
The data acquisition module is additionally operable to apply the inversion pulse of non-layer choosing, is returned with the phase gradient of quickly disturbing of radial sample track Wave low-angle excitation sequence is acquired in real time, is sampled to the inversion recovery process of signal, the radial direction sample track Fast spoiled gradient echo low-angle excitation sequence uses multi-layer intercrossed acquisition mode, and at least two can be acquired after an inversion pulse Layer or more image;
The wherein described radial sample track uses small gold angle sampling configuration, and the azimuth of i-th sampling line is (i-1) ψNDegree, The azimuth of i.e. every sampling line increases a fixed angle ψN, whereinN determines ψN Size, for N by empirically determined, sampling line number mesh must be more than N.
6. device as claimed in claim 5, which is characterized in that the data acquisition module is additionally operable to during breathing is held one's breath Carry out data acquisition, the most short TE and TR of use.
7. device as claimed in claim 5, which is characterized in that the selecting module is additionally operable to, according to sliding window size, return It selects to Gu property acquisition time and is located at the sampling line of diastole as the K space center line rebuild.
8. device as claimed in claim 7, which is characterized in that described image rebuilds module and is additionally operable to use the acquisition time Sampling line positioned at diastole determines the position of reconstruction window as the K space center line rebuild, to the reconstruction window, K space data is weighted using K- spatial weighting image comparison methods, data are rebuild using quick algorithm for reconstructing.
CN201511033083.5A 2015-12-31 2015-12-31 A kind of method and apparatus that cardiac muscle is quantitative T1 Active CN105662413B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511033083.5A CN105662413B (en) 2015-12-31 2015-12-31 A kind of method and apparatus that cardiac muscle is quantitative T1

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511033083.5A CN105662413B (en) 2015-12-31 2015-12-31 A kind of method and apparatus that cardiac muscle is quantitative T1

Publications (2)

Publication Number Publication Date
CN105662413A CN105662413A (en) 2016-06-15
CN105662413B true CN105662413B (en) 2018-10-26

Family

ID=56189976

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511033083.5A Active CN105662413B (en) 2015-12-31 2015-12-31 A kind of method and apparatus that cardiac muscle is quantitative T1

Country Status (1)

Country Link
CN (1) CN105662413B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108175404B (en) * 2018-01-11 2023-08-01 成都心吉康科技有限公司 Electrocardiosignal acquisition device and method for inhibiting touch interference
CN109171670B (en) * 2018-06-25 2021-02-05 天津海仁医疗技术有限公司 3D blood vessel imaging algorithm based on reverse principal component analysis method
CN109143132B (en) * 2018-09-12 2020-11-27 上海联影医疗科技股份有限公司 Magnetic resonance signal acquisition method, imaging method, system, and medium
CN109507622B (en) * 2018-11-16 2021-08-17 上海联影医疗科技股份有限公司 Magnetic resonance imaging method and system
WO2023122919A1 (en) * 2021-12-27 2023-07-06 中国科学院深圳先进技术研究院 Cardiac magnetic resonance diffusion tensor imaging method, apparatus and device, and storage medium

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5545992A (en) * 1995-08-03 1996-08-13 General Electric Company Fast cardiac gated NMR acquisition with improved T1 contrast
US8446148B2 (en) * 2008-08-29 2013-05-21 Siemens Medical Solutions Usa, Inc. Long T1 artifact suppression techniques for magnetic resonance imaging
DE102010001145A1 (en) * 2010-01-22 2011-07-28 Siemens Aktiengesellschaft, 80333 Method and apparatus for magnetic resonance imaging for creating T1 cards
US8848990B2 (en) * 2010-09-28 2014-09-30 Siemens Aktiengesellschaft Automatic registration of image series with varying contrast based on synthetic images derived from intensity behavior model
US9285446B2 (en) * 2010-09-29 2016-03-15 Isis Innovation Limited Systems and methods for shortened look locker inversion recovery (Sh-MOLLI) cardiac gated mapping of T1
US9081073B2 (en) * 2011-01-27 2015-07-14 Siemens Medical Solutions Usa, Inc. System for suppression of artifacts in MR imaging
CN102908143B (en) * 2011-04-22 2016-04-20 纽约大学 To the method and system of at least one image anatomy
CN102764126B (en) * 2011-12-08 2014-03-26 中国科学院深圳先进技术研究院 Method and device for generating T1 weighted image
CN104013405B (en) * 2014-06-09 2017-01-25 深圳先进技术研究院 Dynamic myocardium activity detection method and system

Also Published As

Publication number Publication date
CN105662413A (en) 2016-06-15

Similar Documents

Publication Publication Date Title
CN105662413B (en) A kind of method and apparatus that cardiac muscle is quantitative T1
Scott et al. Motion in cardiovascular MR imaging
US6683454B2 (en) Shifting of artifacts by reordering of k-space
CN102143707B (en) Magnetic resonance imaging apparatus
US20190154785A1 (en) Cardiac and Respiratory Self-Gated Motion-Corrected Free-Breathing Spiral Cine Imaging
CN110031786B (en) Magnetic resonance image reconstruction method, magnetic resonance imaging apparatus, and medium
CN105793722B (en) Real-time adaptive physiology for stable state MR sequence is synchronous and gates
CN104067137B (en) MR imagings with the suppression to flow artefacts
CN105283774B (en) The System and method for of subject's cardiac imaging is improved under unfavorable cardiac condition
WO2013159044A1 (en) System, method and computer-accessible medium for highly-accelerated dynamic magnetic resonance imaging using golden-angle radial samplng and compressed sensing
EP1451603A1 (en) Multiple preparatory excitations and readouts distributed over the cardiac cycle
CN104013404B (en) Cardiac muscle T1 value measurement method and system
KR101809213B1 (en) Method and apparatus for multi-slice imaging of t2-relaxation time
CN109507622A (en) MR imaging method and system
CN106443533B (en) Based on the navigation magnetic resonance diffusion imaging method repeatedly excited and device
Brinegar et al. Real-time cardiac MRI without triggering, gating, or breath holding
Christodoulou et al. Improved subspace estimation for low-rank model-based accelerated cardiac imaging
KR101787369B1 (en) Medical imaging apparatus, magnetic reasonance imaging apparatus and method of controlling the same
CN111631713B (en) Magnetic resonance imaging method, apparatus and storage medium
Lingala et al. Accelerated first pass cardiac perfusion MRI using improved k− t SLR
US8509874B2 (en) Method for non-contrast enhanced pulmonary vein magnetic resonance imaging
CN104013405B (en) Dynamic myocardium activity detection method and system
US7283862B1 (en) Rapid multi-slice MR perfusion imaging with large dynamic range
CN108742626A (en) T1 Parameter Maps imaging method and magnetic resonance imaging system
CN109247938B (en) Quantitative myocardial magnetic resonance imaging method, apparatus and storage medium

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Liu Xin

Inventor after: Zhu Yanjie

Inventor after: Zhong Yaozu

Inventor after: Zheng Hairong

Inventor before: Zhu Yanjie

Inventor before: Zhong Yaozu

Inventor before: Liu Xin

Inventor before: Zheng Hairong

GR01 Patent grant
GR01 Patent grant