CN103006218A - Data processing method and system for cardiac magnetic resonance real-time film imaging - Google Patents
Data processing method and system for cardiac magnetic resonance real-time film imaging Download PDFInfo
- Publication number
- CN103006218A CN103006218A CN2012105630855A CN201210563085A CN103006218A CN 103006218 A CN103006218 A CN 103006218A CN 2012105630855 A CN2012105630855 A CN 2012105630855A CN 201210563085 A CN201210563085 A CN 201210563085A CN 103006218 A CN103006218 A CN 103006218A
- Authority
- CN
- China
- Prior art keywords
- cardiac
- left ventricle
- magnetic resonance
- imaging
- 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.)
- Granted
Links
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention provides a data processing method and a data processing system for cardiac magnetic resonance real-time film imaging. The method comprises the following steps that a selected imaging layer is obtained, and an image of the imaging layer is read; a reference line set on the image is obtained; a dynamic variation chart of a signal contour of the reference line along with acquisition time is obtained; a time interval from the beginning to the end of the end exhalation period is obtained from the dynamic variation chart; the diameter of the left ventricle within the time interval is calculated through the dynamic variation chart, and the heart phases of the end diastolic and the end systolic are determined in accordance with the diameter of the left ventricle; and images corresponding to the heart phases of the end diastolic and the end systolic are extracted and saved. By adopting the data processing method and the data processing system, the data processing speed can be increased.
Description
Technical field
The present invention relates to the nuclear magnetic resonance treatment technology, particularly relate to data processing method and the system of the real-time film imaging of cardiac magnetic resonance.
Background technology
The imaging of cardiac magnetic resonance film is the formation method of the measurement cardiac function comparatively commonly used clinically, generally collecting part K spatial data in a cardiac cycle.Because complete K spatial data is finished in a plurality of cardiac cycles, thereby the method is comparatively responsive for respiratory movement.In order to eliminate respiratory movement to the impact of image quality, the experimenter needs repeated multiple times holding one's breath in scanning process.This requirement is difficult to realize suffering from serious heart disease, respiratory disorder and children's with it, so the method has larger clinical limitation.
The real-time film imaging of cardiac magnetic resonance is by the suitable reduced spatial resolution of technology, employing parallel imaging technique and special algorithm for reconstructing (sharing method etc. such as sliding window method, echo), in a cardiac cycle, can finish the collection of whole K spatial datas, therefore picture quality be subjected to respirometric affect less, can under the free breathing state of experimenter, realize data acquisition, greatly improve the clinical application range of cine cardiac imaging.
Yet the respiratory movement meeting causes the locus of same aspect between different cardiac phases and different aspects not mated in the in real time film imaging, measures for follow-up cardiac function and introduces error.In order to solve in the real-time film imaging of cardiac magnetic resonance because the aspect problem of misalignment that respiratory movement causes, usually the method for doing that adopts is to browse the image of all collections, by range estimation, diastole latter stage and end-systole picture when manually being chosen at EEP (being that respiratory movement is minimum).But this processing method is comparatively loaded down with trivial details time-consuming, greatly reduces its clinical service efficiency and using value.
Summary of the invention
Based on this, being necessary provides the data processing method of the faster real-time film imaging of cardiac magnetic resonance of a kind of processing speed for the slower problem of date processing.
Simultaneously, also be necessary to provide the data handling system of the faster real-time film imaging of cardiac magnetic resonance of a kind of processing speed.
The data processing method of the real-time film imaging of a kind of cardiac magnetic resonance comprises:
Obtain selected imaging aspect, read the image of described imaging aspect;
Obtain the reference line that arranges on the described image;
Obtain the signal profile of described reference line with the dynamic change figure of acquisition time;
From described dynamic change figure, obtain the time interval that EEP begins and finishes;
Calculate left ventricle diameter in the described time interval by described dynamic change figure, and determine cardiac phase and end systole cardiac phase place diastasis according to described left ventricle diameter;
Extract described diastasis of cardiac phase and image corresponding to described end systole cardiac phase place and preserving.
The data handling system of the real-time film imaging of a kind of cardiac magnetic resonance comprises:
Imaging level selection module is used for obtaining selected imaging aspect, reads the image of imaging aspect;
The Guide Settings module is used for obtaining the reference line that arranges on the image;
Dynamic change figure extraction module is used for obtaining the signal profile of reference line with the dynamic change figure of acquisition time;
The time interval determination module is used for obtaining the time interval that EEP begins and finishes from dynamic variation diagram;
Left ventricle diameter computing module is used for by the left ventricle diameter in dynamic change figure interval computation time, and determines cardiac phase and end systole cardiac phase place diastasis according to the left ventricle diameter;
The Image Saving module is used for image corresponding to extraction cardiac phase diastasis and end systole cardiac phase place and preservation.
Adopt the method and system of this programme, read the image on the imaging aspect, and at this image reference line is set, draw the signal profile of reference line with the dynamic change figure of acquisition time, determine cardiac phase and end systole cardiac phase place diastasis by calculating left ventricle diameter in time interval, process fast the data of the real-time film imaging of cardiac magnetic resonance, need not to choose manually, process more easy and efficient is higher.
Description of drawings
Fig. 1 is the flow chart of the data processing method of the real-time film imaging of cardiac magnetic resonance;
Fig. 2 is the image of imaging aspect;
Fig. 3 is that the signal profile of reference line is with the dynamic change figure of acquisition time;
Fig. 4 is left ventricle diameter dynamic change figure;
Fig. 5 is the structured flowchart of the data handling system of the real-time film imaging of cardiac magnetic resonance.
The specific embodiment
In cardiac magnetic resonance film imaging process, because the respiratory movement meeting causes the locus of same aspect between different cardiac phases and different aspects not mated, diastole latter stage and heart contraction picture in latter stage in the time of need to being chosen at EEP (being that respiratory movement is minimum).
The data processing method of the real-time film imaging of cardiac magnetic resonance, consult accompanying drawing 1 ~ 4 and comprise:
S10: obtain selected imaging aspect, read the image of imaging aspect.The imaging of cardiac magnetic resonance film at first need to gather the K spatial data in cardiac cycle, selected some imaging aspects, according to sweep time order automatically read all images of this imaging aspect.In this programme, do not do restriction for the selection of imaging aspect, the selection of for example imaging aspect can be that the user sets.In other embodiments, this image is the T2 weighted graph.
S20: obtain the reference line that arranges on the image.In the image that generates, be provided with reference line.By reference to the accompanying drawings 2, reference line is through thoracic wall, left ventricle center and diaphragm; Be appreciated that reference line can only can reflect the organizational structure of movable information through thoracic wall, left ventricle center or diaphragm etc.In this programme, do not do restriction for the selection of reference line, for example the selection of reference line can be that the user sets.
S30: obtain the signal profile of reference line with the dynamic change figure of acquisition time.By reference to the accompanying drawings 3, be pcolor, the power of color representation signal, blue representation signal intensity a little less than, red representation signal intensity is stronger.Cine cardiac imaging is the T2 weighted graph, so chambers of the heart inner blood signal is high signal (redness), and pulmonary is low signal (blueness), and other organs (such as liver) signal is placed in the middle.The signal intensity that reference line on all images on this aspect is corresponding launches (that is, launching by the imaging time with reference to the line signal intensity), gets final product the picked up signal profile with the dynamic change figure of acquisition time.The y direction of this dynamic change figure represents the signal profile of reference line, and X direction is acquisition time, and comprises a plurality of cardiac cycles.
S40: from dynamic variation diagram, obtain the time interval that EEP begins and finishes.EEP is that corresponding time point was interval when thoracic wall extreme lower position or diaphragm moved to the extreme higher position among the dynamic change figure.Because respiration motion cycle is longer than cardiac cycle, therefore, EEP generally comprises a complete cardiac cycle at least in time interval.
S50: by the left ventricle diameter in dynamic change figure interval computation time, and determine cardiac phase and end systole cardiac phase place diastasis according to the left ventricle diameter.Dynamic change figure in time interval calculates left ventricle diameter among the dynamic change figure.The left ventricle diameter value of all acquisitions is compared, and wherein left ventricle diameter maximum is cardiac phase diastasis, left ventricle diameter minimum be the end systole cardiac phase place.By reference to the accompanying drawings 4, the computational methods of left ventricle diameter are: obtain the pixel number of left ventricle inner blood signal on the reference line, calculating pixel is counted and is obtained the left ventricle diameter.This shows that these left ventricle diameter computational methods are not only convenient but also accurate.
In other embodiments, determining of diastasis cardiac phase and end systole cardiac phase place, can pass through this scheme: setting threshold, the left ventricle diameter is cardiac phase diastasis greater than threshold value, the left ventricle diameter is the end systole cardiac phase place less than threshold value.The selection that is appreciated that threshold value can be a value range, also can be a value; Determining of this threshold value can be obtaining of experience accumulation, also can be that the user sets.In this application, the size of threshold value is not done restriction.
S60: the image that extraction cardiac phase diastasis and end systole cardiac phase place are corresponding is also preserved.Extract the image of corresponding phase and be saved to specific file, serve follow-up graphical analysis.
Adopt this programme, read the image on the imaging aspect, and at this image reference line is set, draw the signal profile of reference line with the dynamic change figure of acquisition time, determine cardiac phase and end systole cardiac phase place diastasis by calculating left ventricle diameter in time interval, process fast the data of the real-time film imaging of cardiac magnetic resonance, need not to choose manually, process more easy and efficient is higher.
Further, the data processing method of the real-time film imaging of above-mentioned cardiac magnetic resonance also comprises: choose other imaging aspect, and return the step of the image that reads the imaging aspect.Finish the date processing of this aspect, according to this scheme other imaging aspect is being disposed.
Based on the data processing method of the real-time film imaging of above-mentioned cardiac magnetic resonance, the data handling system of the real-time film imaging of cardiac magnetic resonance also is provided, consult accompanying drawing 5, comprising:
Imaging level selection module 10 is used for obtaining selected imaging aspect, reads the image of imaging aspect.The imaging of cardiac magnetic resonance film at first need to gather the K spatial data in cardiac cycle, selected some imaging aspects, according to sweep time order automatically read all images of this imaging aspect.
Dynamic change figure extraction module 30 is used for obtaining the signal profile of reference line with the dynamic change figure of acquisition time.
Time interval determination module 40 is used for obtaining the time interval that EEP begins and finishes from dynamic variation diagram.EEP is that corresponding time point was interval when thoracic wall extreme lower position or diaphragm moved to the extreme higher position among the dynamic change figure.Because respiration motion cycle is longer than cardiac cycle, therefore, the EEP time interval generally comprises a complete cardiac cycle at least.
The direct computing module 50 of left ventricle is used for by the left ventricle diameter in dynamic change figure interval computation time, and determines cardiac phase and end systole cardiac phase place diastasis according to the left ventricle diameter.Dynamic change figure in time interval calculates left ventricle diameter among the dynamic change figure.The left ventricle diameter value of all acquisitions is compared, and wherein left ventricle diameter maximum is cardiac phase diastasis, left ventricle diameter minimum be the end systole cardiac phase place.
In other embodiments, determining of diastasis cardiac phase and end systole cardiac phase place, can pass through setting threshold determining unit (not shown): setting threshold, the left ventricle diameter is cardiac phase diastasis greater than threshold value, and the left ventricle diameter is the end systole cardiac phase place less than threshold value.The selection that is appreciated that threshold value can be a value range, also can be a value; Determining of this threshold value can be obtaining of experience accumulation, also can be that the user sets.In this application, the size of threshold value is not done restriction.
Adopt native system, read the image on the imaging aspect, and at this image reference line is set, draw the signal profile of reference line with the dynamic change figure of acquisition time, determine cardiac phase and end systole cardiac phase place diastasis by calculating left ventricle diameter in time interval, process fast the data of the real-time film imaging of cardiac magnetic resonance, need not to choose manually, process more easy and efficient is higher.
In other embodiments, the data handling system of the real-time film imaging of above-mentioned cardiac magnetic resonance also comprises: other imaging aspect is chosen the module (not shown), is used for choosing other imaging aspect, and is back into picture level selection module 10.Finish the date processing of this aspect, according to this scheme other imaging aspect is being disposed.
The above embodiment has only expressed several embodiment of the present invention, and it describes comparatively concrete and detailed, but can not therefore be interpreted as the restriction to claim of the present invention.Should be pointed out that for the person of ordinary skill of the art without departing from the inventive concept of the premise, can also make some distortion and improvement, these all belong to protection scope of the present invention.Therefore, the protection domain of patent of the present invention should be as the criterion with claims.
Claims (10)
1. the data processing method of the real-time film imaging of cardiac magnetic resonance comprises:
Obtain selected imaging aspect, read the image of described imaging aspect;
Obtain the reference line that arranges on the described image;
Obtain the signal profile of described reference line with the dynamic change figure of acquisition time;
From described dynamic change figure, obtain the time interval that EEP begins and finishes;
Calculate left ventricle diameter in the described time interval by described dynamic change figure, and determine cardiac phase and end systole cardiac phase place diastasis according to described left ventricle diameter;
Extract described diastasis of cardiac phase and image corresponding to described end systole cardiac phase place and preserving.
2. the data processing method of the real-time film imaging of cardiac magnetic resonance according to claim 1 is characterized in that, described reference line is through thoracic wall, left ventricle center and diaphragm.
3. the data processing method of the real-time film imaging of cardiac magnetic resonance according to claim 2 is characterized in that, described EEP is: corresponding time point was interval when thoracic wall extreme lower position or diaphragm moved to the extreme higher position.
4. the data processing method of the real-time film imaging of cardiac magnetic resonance according to claim 3 is characterized in that, described time interval comprises a cardiac cycle at least.
5. the data processing method of the real-time film imaging of cardiac magnetic resonance according to claim 4, it is characterized in that, describedly determine cardiac phase and end systole cardiac phase place diastasis according to described left ventricle diameter, be specially: described left ventricle diameter maximum be cardiac phase diastasis; Described left ventricle diameter minimum be the end systole cardiac phase place.
6. the data handling system of the real-time film imaging of cardiac magnetic resonance is characterized in that, comprising:
Imaging level selection module is used for obtaining selected imaging aspect, reads the image of imaging aspect;
The Guide Settings module is used for obtaining the reference line that arranges on the image;
Dynamic change figure extraction module is used for obtaining the signal profile of reference line with the dynamic change figure of acquisition time;
The time interval determination module is used for obtaining the time interval that EEP begins and finishes from dynamic variation diagram;
Left ventricle diameter computing module is used for by the left ventricle diameter in dynamic change figure interval computation time, and determines cardiac phase and end systole cardiac phase place diastasis according to the left ventricle diameter;
The Image Saving module is used for image corresponding to extraction cardiac phase diastasis and end systole cardiac phase place and preservation.
7. the data handling system of the real-time film imaging of cardiac magnetic resonance according to claim 6 is characterized in that, described reference line is through thoracic wall, left ventricle center and diaphragm.
8. the data handling system of the real-time film imaging of cardiac magnetic resonance according to claim 7 is characterized in that, described EEP is: corresponding time point was interval when thoracic wall extreme lower position or diaphragm moved to the extreme higher position.
9. the data handling system of the real-time film imaging of cardiac magnetic resonance according to claim 8 is characterized in that, described time interval comprises a cardiac cycle at least.
10. the data handling system of the real-time film imaging of cardiac magnetic resonance according to claim 9, it is characterized in that, describedly determine cardiac phase and end systole cardiac phase place diastasis according to described left ventricle diameter, be specially: described left ventricle diameter maximum be cardiac phase diastasis; Described left ventricle diameter minimum be the end systole cardiac phase place.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210563085.5A CN103006218B (en) | 2012-12-22 | 2012-12-22 | The data processing method of cardiac magnetic resonance real-time film imaging and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210563085.5A CN103006218B (en) | 2012-12-22 | 2012-12-22 | The data processing method of cardiac magnetic resonance real-time film imaging and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103006218A true CN103006218A (en) | 2013-04-03 |
| CN103006218B CN103006218B (en) | 2015-12-23 |
Family
ID=47955714
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210563085.5A Active CN103006218B (en) | 2012-12-22 | 2012-12-22 | The data processing method of cardiac magnetic resonance real-time film imaging and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103006218B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103340628A (en) * | 2013-06-28 | 2013-10-09 | 中国科学院深圳先进技术研究院 | Method and system for processing heart real-time film imaged picture |
| CN109001660A (en) * | 2018-06-12 | 2018-12-14 | 上海联影医疗科技有限公司 | Film imaging method and magnetic resonance imaging system |
| CN113133756A (en) * | 2021-04-23 | 2021-07-20 | 上海联影医疗科技股份有限公司 | Three-dimensional cardiac cine imaging method, magnetic resonance imaging system and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007160108A (en) * | 2005-12-12 | 2007-06-28 | Siemens Medical Solutions Usa Inc | System and method for image based physiological monitoring of cardiovascular function |
| US20100264922A1 (en) * | 2009-04-15 | 2010-10-21 | Siemens Medical Solutions Usa, Inc. | System for Multi Nucleus Cardiac MR Imaging and Spectroscopy |
-
2012
- 2012-12-22 CN CN201210563085.5A patent/CN103006218B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007160108A (en) * | 2005-12-12 | 2007-06-28 | Siemens Medical Solutions Usa Inc | System and method for image based physiological monitoring of cardiovascular function |
| US20100264922A1 (en) * | 2009-04-15 | 2010-10-21 | Siemens Medical Solutions Usa, Inc. | System for Multi Nucleus Cardiac MR Imaging and Spectroscopy |
Non-Patent Citations (1)
| Title |
|---|
| MEINRAD BEER ET AL.: "Free breathing cardiac real-time cine MR without ECG triggering", 《INTERNATIONAL JOURNAL OF CARDIOLOGY》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103340628A (en) * | 2013-06-28 | 2013-10-09 | 中国科学院深圳先进技术研究院 | Method and system for processing heart real-time film imaged picture |
| CN103340628B (en) * | 2013-06-28 | 2015-07-29 | 中国科学院深圳先进技术研究院 | Heart movie real-time image processing method and system |
| CN109001660A (en) * | 2018-06-12 | 2018-12-14 | 上海联影医疗科技有限公司 | Film imaging method and magnetic resonance imaging system |
| CN109001660B (en) * | 2018-06-12 | 2020-07-28 | 上海联影医疗科技有限公司 | Film imaging method and magnetic resonance imaging system |
| CN113133756A (en) * | 2021-04-23 | 2021-07-20 | 上海联影医疗科技股份有限公司 | Three-dimensional cardiac cine imaging method, magnetic resonance imaging system and storage medium |
| CN113133756B (en) * | 2021-04-23 | 2023-08-15 | 上海联影医疗科技股份有限公司 | Three-dimensional cardiac cine imaging method, magnetic resonance imaging system and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103006218B (en) | 2015-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101084528B (en) | A method and a system for integration of medical diagnostic information and a geometric model of a movable body | |
| US11026636B2 (en) | Method for generating a medical data set of a moving body part | |
| CN108364290B (en) | Method, medium, and system for analyzing a sequence of images of periodic physiological activity | |
| CN103792502B (en) | Automatic system for sequential in imaging | |
| EP3413075A1 (en) | Mri involving the generation of a physiological motion signal using a pilot tone navigator | |
| US9804246B2 (en) | Magnetic resonance imaging apparatus | |
| US9462952B2 (en) | System and method for estimating artery compliance and resistance from 4D cardiac images and pressure measurements | |
| US9603554B2 (en) | Systems, methods, and media for monitoring the condition of a patient's heart | |
| CN110246137A (en) | A kind of imaging method, device and storage medium | |
| CN102008307A (en) | Magnetic resonance diffusion tensor imaging method and system | |
| JP2008529611A (en) | Method and system for specimen motion artifact correction | |
| Veronesi et al. | Tracking of left ventricular long axis from real-time three-dimensional echocardiography using optical flow techniques | |
| KR20130079694A (en) | Apparatus and method for measureing velocity vector imaging of blood vessel | |
| CN103340628B (en) | Heart movie real-time image processing method and system | |
| WO2016069602A1 (en) | A method and device for magnetic resonance imaging data acquisition guided by physiologic feedback | |
| US10548494B2 (en) | Method for determining a personalized cardiac model using a magnetic resonance imaging sequence | |
| CN103006218A (en) | Data processing method and system for cardiac magnetic resonance real-time film imaging | |
| Marin et al. | Numerical surrogates for human observers in myocardial motion evaluation from SPECT images | |
| CN116849691B (en) | Method, equipment and storage medium for automatically identifying global optimal phase of cardiac CT imaging | |
| Hennersperger et al. | Vascular 3D+ T freehand ultrasound using correlation of doppler and pulse-oximetry data | |
| Christodoulou et al. | Four-dimensional MR cardiovascular imaging: Method and applications | |
| JP2013040829A (en) | Volume data processor and method | |
| CN109118526A (en) | A kind of senile dementia image analysis system and analysis method based on virtual reality | |
| CN104361606B (en) | Cardiac cycle restoration methods in x-ray imaging image sequence | |
| CN102764126B (en) | Method and device for generating T1 weighted image |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |