CN101238979B - MRI involving a cine prescan for motion analysis - Google Patents
MRI involving a cine prescan for motion analysis Download PDFInfo
- Publication number
- CN101238979B CN101238979B CN2008100045914A CN200810004591A CN101238979B CN 101238979 B CN101238979 B CN 101238979B CN 2008100045914 A CN2008100045914 A CN 2008100045914A CN 200810004591 A CN200810004591 A CN 200810004591A CN 101238979 B CN101238979 B CN 101238979B
- Authority
- CN
- China
- Prior art keywords
- mentioned
- image
- imaging apparatus
- time series
- magnetic field
- 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.)
- Expired - Fee Related
Links
- 230000033001 locomotion Effects 0.000 title claims description 23
- 238000004458 analytical method Methods 0.000 title claims description 10
- 238000003384 imaging method Methods 0.000 claims abstract description 29
- 239000000284 extract Substances 0.000 claims abstract description 17
- 210000004351 coronary vessel Anatomy 0.000 claims description 42
- 238000000605 extraction Methods 0.000 claims description 22
- 238000012360 testing method Methods 0.000 claims description 17
- 238000003703 image analysis method Methods 0.000 claims description 12
- 238000001514 detection method Methods 0.000 claims description 5
- 230000002269 spontaneous effect Effects 0.000 claims description 4
- 238000002595 magnetic resonance imaging Methods 0.000 abstract 1
- 238000000034 method Methods 0.000 description 17
- 230000009471 action Effects 0.000 description 9
- 238000012545 processing Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000003068 static effect Effects 0.000 description 7
- 238000005481 NMR spectroscopy Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 238000013480 data collection Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000010247 heart contraction Effects 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000005311 nuclear magnetism Effects 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000010626 work up procedure Methods 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
A magnetic resonance imaging apparatus comprises a cine-imaging unit (12). The cine-imaging unit (12) collects the time-series images of a region of interest in a subject (10) and reconstructs an image. Next, the characterizing-part detecting unit (14) detects the characteristics of the time-series images. The motion-analyzing unit (16) analyzes motion-characteristic values of the characteristics extracted by the characterizing-part detecting unit (14). The image-extracting unit (18) extracts a specified time-series image in accordance with the motion-characteristic values.
Description
The cross-reference of related application
The application is based on the previous Japanese patent application of submitting on February 6th, 2007 2007-27050 number and require its priority, introduces its full content herein as a reference.
Technical field
MR imaging apparatus that the present invention relates to use in medical diagnostic imaging apparatus and image analysis method thereof in more detail, relate to the MR imaging apparatus and the image analysis method thereof of the heart coronary artery shooting that is used for being undertaken by MR imaging apparatus.
Background technology
Image diagnosing system has in recent years developed into can be at clinical field conduct cardiac work up.In nuclear magnetic resonance (MRI), motion picture pickup, delay radiography, the shooting of load heart muscle perfusion etc. can have been used.But, although it gets along with, coronary artery shooting aspect picture quality or stable aspect still have the shooting of higher technology barrier.This is owing to have that coronary artery is thin, it move towards the complicated movement of complexity, heart and according to the patient and not equal many difficult points, the cause of high spatial resolution, high time resolution must be arranged.Therefore, from the adjustment of patient's the setting that is arranged to photography conditions, magnetic field bump etc. etc., require the height technology respectively.
Recently, for example in Magnetic Resonance In Resonance magazine, 50,1223-1228,2003 or by assistant, in image information MEDICAL magazine, Vol.36, No.10, proposed to use SSFP (the Steady-state FreePrecession steady state free precession) sequence that was balanced between for a long time etc. and merged to use and breathe navigator and follow the trail of the method that covers heart integral body by Weber etc.
In such method, because because of the per 1 heart beating data in synchronization acquisition time of sectional cause is as short as about 100ms, so be necessary to make during this data collection just in time consistent with the time phase of immobilized heart coronarius.Therefore, in the prior art,, before the coronary artery shooting, implement the heart motion picture pickup in advance, and very attentional view is examined this view data in order to obtain the time phase of immobilized heart coronarius.At this moment, be necessary to observe simultaneously right coronary artery and left coronary artery.
But under most situation, right coronary artery is different with phase place quiescent time of left coronary artery, and in addition, time width is also different.Therefore, the operator is necessary the observation of exactissima diligentia more film image.Therefore, carrying out needing the more time before the actual coronary artery shooting, become and check the throughput main reasons for decrease.In addition, because of the review time lengthening, patient's feeling of fatigue increases, and the result is because electrocardiogram or breathing pattern become unstable, so become the main cause that causes image quality decrease.
Summary of the invention
Thereby the object of the present invention is to provide can be easily and to obtain coronary artery exactly be static or move MR imaging apparatus and image analysis method thereof during few.
Promptly, the present invention is a kind of MR imaging apparatus, magnetostatic field takes place in this device, for the overlapping leaning magnetic field that utilizes gradient magnetic field coil to take place of this magnetostatic field, subject in the above-mentioned magnetostatic field is detected the magnetic resonance signal that the high frequency magnetic field be used to spontaneous sending coil produces in above-mentioned subject, obtain image from this magnetic resonance signal, it is characterized in that, possess: collect the film image image pickup part that the time series image at the object position of above-mentioned subject is reconstructed; The feature portion test section that the characteristic of above-mentioned time series image is detected; To the motion analysis portion that analyzes with the kinetic characteristic value of the detected characteristic of above-mentioned feature portion's test section; And according to the extraction unit of above-mentioned kinetic characteristic value extraction about the information of specific time phase.
In addition, the present invention is a kind of MR imaging apparatus, it is characterized in that, possesses: the applying unit that applies leaning magnetic field and high-frequency impulse in magnetostatic field for the subject of placing; According to corresponding to the applying of above-mentioned leaning magnetic field and above-mentioned high-frequency impulse, collect the film image image pickup part that the time series image at the object position of above-mentioned subject is reconstructed from the magnetic resonance signal of above-mentioned subject radiation; The feature portion test section that the characteristic of above-mentioned time series image is detected; To the motion analysis portion that analyzes with the kinetic characteristic value of the detected characteristic of above-mentioned feature portion's test section; And according to the extraction unit of above-mentioned kinetic characteristic value extraction about the information of specific time phase.
Moreover, the present invention is a kind of image analysis method of MR imaging apparatus, this MR imaging apparatus generation magnetostatic field, for the overlapping leaning magnetic field that utilizes gradient magnetic field coil to take place of this magnetostatic field, subject in the above-mentioned magnetostatic field is detected the magnetic resonance signal that the high frequency magnetic field be used to spontaneous sending coil produces in above-mentioned subject, obtain image from this magnetic resonance signal, this image analysis method is characterised in that: the time series image of collecting the object position of above-mentioned subject is reconstructed, characteristic to above-mentioned time series image detects, kinetic characteristic value to above-mentioned detected characteristic is analyzed, according to the information of above-mentioned kinetic characteristic value extraction about specific time phase.
To state additional purpose of the present invention and advantage in the following description, the part of these purposes and advantage is conspicuous according to description, or can recognize from the practice of the present invention.Utilize particularly in following means of pointing out and combination, can realize and obtain objects and advantages of the present invention.
Description of drawings
The accompanying drawing that is incorporated in the description and constitutes the part of description illustrates the preferred embodiments of the present invention now, with the above overall description that provides and below the concrete description of the preferred embodiment that provides, be used for illustrating principle of the present invention.
Figure 1A is the block diagram of the basic structure of the expression MR imaging apparatus relevant with the 1st embodiment of the present invention, and Figure 1B is the block diagram of the basic structure of the expression MR imaging apparatus relevant with the 2nd embodiment of the present invention.
Fig. 2 is the block diagram of the structure of the expression MRI device relevant with the 1st embodiment of the present invention.
Fig. 3 is the flow chart of work that is used for illustrating the MRI device of the 1st embodiment of the present invention.
Fig. 4 A and 4B are the figure of example of the film image of expression heart, and Fig. 4 A is the figure of expression original image, and Fig. 4 B is the figure of expression binary image.
Fig. 5 is that expression is calculated kinetic characteristic value P to right coronary artery
iThe curve chart of data.
Fig. 6 is the flow chart of action that is used for illustrating the MRI device of the 2nd embodiment of the present invention.
The specific embodiment
Below, with reference to the description of drawings embodiments of the present invention.
Figure 1A is the block diagram of the basic structure of the expression MR imaging apparatus (be designated hereinafter simply as MRI device) relevant with the 1st embodiment of the present invention.
In Figure 1A, this MRI device is made of film image image pickup part 12, feature portion test section 14, motion analysis portion 16 and extraction unit 18.Above-mentioned film image image pickup part 12 is used for the time series image at intelligence-collecting object position (subject 10) and is reconstructed.Above-mentioned feature portion test section 14 detects the characteristic of above-mentioned time series image, and above-mentioned motion analysis portion 16 is used to analyze the kinetic characteristic value of this characteristic.Have, extraction unit 18 is extracted specific time series image according to above-mentioned kinetic characteristic value again.
In such structure, in film image image pickup part 12 from as the subject 10 acquisition time image series at object position and be reconstructed.Then, in feature portion test section 14, detect the specific part of heart of characteristic, for example subject 10 of the time series image of collecting with above-mentioned film image image pickup part 12.Obtain at the displacement of the specific part of this detected heart and characteristic quantity thereof as the kinetic characteristic value with motion analysis portion 16.Then, in extraction unit 18, extract specific time series image according to above-mentioned kinetic characteristic value.
Fig. 2 is the block diagram of the structure of the expression MRI device relevant with the 1st embodiment of the present invention.
This MRI device 30 possesses: the bed portion of placing subject 10; The magnetostatic field generating unit of magnetostatic field takes place; Be used for leaning magnetic field generating unit to the magnetostatic field additional location information; Send the transmission acceptance division that receives high-frequency signal; Bear the control of entire system and control, the operational part of image reconstruction; Detection is as the electrocardio test section of the ECG signal of the signal of the time phase of the heart of expression subject P; And be used for the instruction department of holding one's breath that instruction is held one's breath to subject 10.
Above-mentioned magnetostatic field generating unit comprises the Magnet (static field magnet) 32 of for example superconduction formula and this Magnet 32 is supplied with the magnetostatic field power supply 34 of electric currents, goes up generation magnetostatic field H at the direction of principal axis (Z-direction) of the peristome cylindraceous that inserts subject 10 (space use in diagnosis)
0In addition, in this Magnet portion, be provided with unifluxor circle 58.Under the control of master computer 42 described later, supply with the electric current that is used to make the magnetostatic field homogenization from 60 pairs of these unifluxor circles of unifluxor circle power supply 58.For bed portion, the top board of having placed subject 10 can be inserted in the peristome of Magnet 32 in the mode that can keep out of the way.
The leaning magnetic field generating unit comprises gradient magnetic field coil unit 36.This gradient magnetic field coil unit 36 possesses 3 groups of coil 36x, coil 36y, the coil 36z of the leaning magnetic field of the X, the Y that are used to take place mutually orthogonal and Z-direction.The leaning magnetic field generating unit also comprises the leaning magnetic field power supply 38 of above-mentioned coil 36x～36z being supplied with electric current.This leaning magnetic field power supply 38 is supplied with the pulse current that is used for run-off the straight magnetic field to coil 36x～36z under the control of sequencer 40 described later.
By adjusting the pulse current of supplying with from 38 couples of coil 36x of leaning magnetic field power supply～36z, each the leaning magnetic field that can synthesize X, Y as physical axis and Z direction is at random set slice direction leaning magnetic field Gs, the phase-encoding direction leaning magnetic field Ge of mutually orthogonal and is read each logic direction of principal axis of direction (frequency coding direction) leaning magnetic field Gr.Slice direction, phase-encoding direction and each leaning magnetic field of reading direction are overlapped magnetostatic field H
0On.
Sending acceptance division is included in the photographing space in the Magnet 32 near the RF of the configuration subject 10 coil 44 and is connected to transmitter 46T and receptor 46R on this RF coil 44.Transmitter 46T and receptor 46R move under the control of sequencer 40 described later.Transmitter 46T supplies with the RF pulse of the Larmor frequency that is used to cause nuclear magnetic resonance, NMR (NMR) to RF coil 44.Receptor 46R is taken into the echo-signal (high-frequency signal) that RF coil 44 receives, behind the various signal processing that it carried out preposition amplification, medium-frequency transformation, phase detection, low frequency amplification, filtering etc., carry out the A/D conversion, generate the echo data (former data) of the digital quantity corresponding with echo-signal.
Control, operational part comprise sequencer (being also referred to as sequence controller) 40, master computer 42, arithmetic element 50, memory element 52, display 54, loader 56 and sound generator 62.Wherein, master computer 42 has the function of utilizing institute's saved software time ordered pair sequencer 40 command pulse sequence informations and synthetically controlling the whole action of this device.
Master computer 42 is in the scanning of imaging after the location is with the preparation work that scans etc.Image scanning is the scanning of collecting the group of echo data necessary in the image reconstruction, is set at two-dimensional scan at this.Can merge and use the ECG door method that depends on the ECG signal to carry out image scanning.In addition, according to circumstances, but also this ECG door method is used in nonjoinder.
Though it is not shown, but sequencer 40 has CPU and memorizer, the pulse train information that storage is sent here from master computer 42, action separately according to this information control leaning magnetic field power supply 38, transmitter 46T and receptor 46R, in case and the echo data of input sink 46R output, just send it in the arithmetic element 50.Pulse train information is in order to make leaning magnetic field power supply 38, transmitter 46T and receptor 46R move necessary full detail according to the series of pulses sequence, the information that for example comprise intensity, application time about the pulse current that coil 36x～36z is applied, applies regularly etc.
Arithmetic element 50 is through the echo data of sequencer 40 input sink 46R output.Arithmetic element 50 is configuration echo data in the fourier space (being also referred to as k space or frequency space) on portion's memorizer within it, on each group in each group this echo data is implemented two dimension or three-dimensional Fourier transform, be reconstructed into the view data of the real space.
Memory element 52 not only can be certainly by the view data of reconstruct, and can carry out the view data of above-mentioned synthetic processing or difference processing certainly.Display 54 is display image under the control of master computer 42.Can be input in the master computer 42 through photography conditions, pulse train, the information synthetic about image or calculus of differences that loader 56 is wished the operator.
In addition, this MRI device 30 possesses the key element of sound generator 62 as the instruction department of holding one's breath.This sound generator 62 can send the hold one's breath beginning and the message that finishes of holding one's breath under the instruction from master computer 42 as sound.
The electrocardio test section comprises the body surface detection ECG signal that is attached to subject and comprises the various processing of digitized processing and export to master computer 42 and the ECG unit 66 of sequencer 40 as the ECG pick off 66 of the signal of telecommunication with to this sensor signal.When carrying out image scanning, use the measuring-signal that obtains by the electrocardio test section by sequencer 40.Thus, synchronization timing can be suitably set, image scanning according to the ECG door method of this synchronization timing can be carried out to carry out data collection based on ECG door method (the same footwork of electrocardio).
Secondly, the action of the MRI device 30 that constitutes like that more than illustrating with reference to the flow chart of Fig. 3.
This sequence is then at first carried out the film image shooting that utilizes the time series image at film image image pickup part 12 intelligence-collecting object positions and be reconstructed at the beginning in step S1.For example, in MRI device 30, realize as described below.That is, utilize the nuclear-magnetism spin of static field magnet 32 magnetization shooting objects, send coil by RF and apply RF magnetic field, nmr phenomena takes place.Detect the resonance signal that takes place this moment with the RF receiving coil.
In order to make a video recording, promptly, send reception with RF and regularly synchronously drive gradient magnetic field coil unit 36 with pulse type in order to give spatial information to received signal, in static field magnet 32, generate leaning magnetic field thus.Then, make the control information of having put down in writing with the situation that applies and receive of predetermined timing execution RF magnetic field or leaning magnetic field, send it to sequencer 40 with master computer 42.
The received signal that has been endowed above-mentioned spatial information is carried out image reconstruction in reconfiguration device (arithmetic element 50).In order to carry out motion picture pickup, repeat the shooting in same cross section by time series.With high time resolution the big position of motions such as heart is being carried out under the situation of motion picture pickup, the data collection of 1 width of cloth size is being carried out segmentation, the data collection that repeats respectively to be partitioned into by time series.The result who does like this, the film image that can obtain having time series information.
For example, be under the situation of purpose with the heart coronary artery shooting, the coronary artery about the motion picture pickup cross section is mostly is discernible cross section (for example cross section, four chambeies) on piece image.As the pulse train of in motion picture pickup, using, often use the SSFP method that was balanced, can describe in following mode, that is, the signal intensity height of coronary artery and fat is described in the coronary artery and the boundary signal intensity lowland of fat.
Because to run through situation out of shape on the direction in cross section many for coronary artery in cross section, four chambeies, so right coronary artery has the toroidal that has the black surround edge on the time series image.Left coronary artery also is depicted as substantially by the black surround edge and surrounds, but also departs from circle sometimes.No matter under any situation, coronary artery is all described in the high signal area at the edge with low-signal areas.
Secondly, in step S2, detected characteristics part in each image of the time series image that obtains with above-mentioned steps S1.Under the situation of heart coronary artery shooting, utilization coronary artery in the time series image is the regional this point with edge, and characteristic is equivalent to coronary artery itself.In ECG unit 66, sequencer 40 etc., in order to detect coronary artery, make the signal intensity binaryzation of each image of time series image, make binary image.In this binary image, extract small closed-loop structure.The time series image of heart by binaryzation situation under because a plurality of closed-loop structure take place, derive from closed-loop structure coronarius so be necessary to detect effectively.
The method that is used for this has 2.The 1st method is the upper limit of the size of supposition closed-loop configuration, and will be decided to be the method coronarius that derives from smaller or equal to the closed loop of this upper limit.This be according to as if the situation of heart under little closed-loop configuration can only be that the blood vessel this point detects.The 2nd method is to limit in advance to be envisioned for to have zone coronarius, and will this region memory closed-loop configuration be decided to be the method coronarius that derives from.About the zone that is limited, available apparatus is specified, and also can be specified by the operator.In addition, also these 2 methods can be combined, in this case, can further improve and detect precision coronarius.
But, as above-mentioned characteristic, can not necessarily detect part coronarius, so long as the expression specific part just can not limited by this.
Fig. 4 A and 4B are the figure of example of the film image of expression heart, and Fig. 4 A is the figure of expression original image, and Fig. 4 B is the figure of expression binary image.In addition, the part of the white circle of representing with arrow among the figure is a coronary artery.
Secondly, in step S3, analyze the kinetic characteristic value of the feature portion of being detected.Execution analysis method as described below.At first, obtain position of centre of gravity on the image of detected feature portion.To be decided to be about whole position of centre of gravitys of obtaining of time series image:
(X
i, Y
i) (i=1,2 ..., N, N is the number of time series image)
In addition, define kinetic characteristic value P as described below
i
P
i=
{(X
i-X
i-1)
2+(Y
i-Y
i-1)
2}
(i＝1、2、...、N-1)
P
iIt is the amount what degree the interframe position of centre of gravity that is illustrated in the time series image has been offset.If the side-play amount in interframe feature portion is big, then P
iBecome big value, if side-play amount is little, P then
iBecome little value.
Fig. 5 is that expression is calculated kinetic characteristic value P to right coronary artery
iThe curve chart of data.
The motion of latter half heart diminishes in expansionary phase, but P
iBe illustrated in the time phase that the end stops quantitatively.As the kinetic characteristic value, also may be defined as other value.For example, such definition of representing of also useful following formula.
Q
i,k＝
{(X
i-X
k)
2+(Y
i-Y
k)
2}
(i＝1、2、...、N-1，1≤k≤N)
Q
I, kBe equivalent to k image of time series image is decided to be the side-play amount of benchmark.
The kinetic characteristic value also can be considered definition that has utilized the phase cross-correlation etc. except above-mentioned two kinds of definition.
At last, in step S4, extract specific time series image according to above-mentioned kinetic characteristic value.Under the situation of heart coronary artery shooting, know that be important the time phase that coronary artery moves hardly.Thereby, utilize above-mentioned kinetic characteristic value P
i, select P
iValue smaller or equal to the P that approaches zero
LimitI image get final product.P
LimitThe feasible value that is equivalent to move can also can be specified the inner setting value that should be worth as extraction unit 18 by the operator.Display 54 with the time image series shows the image that is extracted selectively, or adds identifying information in incidental information.By doing like this, the operator can easily judge the extraction result.
Secondly, the variation of the 1st embodiment of the present invention is described.
Extracting under the situation of image, also be useful as the variation of following narration according to above-mentioned kinetic characteristic value.
With the extraction method of the 1st variation, can tackle and use with respect to P
LimitSize judged P
iThe time a plurality of images of being extracted situations about on the time series direction, scattering.In this case, also will be sandwiched in the image that is extracted diffusedly and the image that is not extracted is added on the object that is extracted.Perhaps, P in being sandwiched in the image that is extracted diffusedly and in the image that is not extracted
iSatisfy under the situation of following formula,
P
i＜(1+ α) P
Limit(α be little on the occasion of)
This image is added on the object that is extracted.By doing like this, the image that is not extracted because of the error of calculation can be added to and extract on the object, about the proximate consecutive image of extraction property stably that moves.
Secondly, the 2nd variation of the 1st embodiment of the present invention is described.
In this 2nd variation, utilize a plurality of feature of feature portion test section 14 detections portion.In the heart coronary artery shooting, for example right coronary artery and left coronary artery are detected as feature portion.Moreover, provide kinetic characteristic value R with following formula
i
R
i＝(∏i
{(Xji-Xji-1))2
+(Yji-Yji-1)2})1/M
(i=1,2 ..., N-1, N is the number of time series image)
(j=1,2 ..., M, M is the number of the feature portion detected)
(∏ j represent j=1,2 ..., M long-pending)
By using this characteristic value R
iCarry out extraction unit 18, can extract all static such image of a plurality of feature portions.Certainly, at this R
iResearch R when all static
iDefinition, can realize that the corresponding extraction of motor pattern with feature portion is judged.
As mentioned above, according to present embodiment, in the time series image, can extract about image like its property class of moving.Thereby, when heart coronary artery is made a video recording desired, can be promptly and catch motion coronarius few time phase reliably, can improve the picture quality of coronary artery shooting.
In addition, in the above-described embodiment, under the situation of heart coronary artery shooting, be made as coronary artery almost motionless time phase, but be not limited to 1 this time phase.For example also can be a certain regular between in can obtain the situations of 2 time phases.
(the 2nd embodiment)
Secondly, the 2nd embodiment of the present invention is described.
Figure 1B is the block diagram of the basic structure of the expression MR imaging apparatus (MRI device) relevant with the 2nd embodiment of the present invention.
In addition, in this 2nd embodiment because identical with action about the structure of the MRI device of the 1st embodiment represented among the basic structure of MRI device and action and Figure 1A, so to together a part of attached with same with reference to label, omit its diagram and explanation, different parts only is described.
In Figure 1B, this MRI device is made of following part: film image image pickup part 12; Feature portion test section 14; Motion analysis portion 16; Extraction unit 18; About the specific time series image that extracts with this extraction unit 18 carry out coronary artery photography synchronously or the condition enactment portion 20 of the setting of photography conditions; And the image pickup part 22 of making a video recording with the photography conditions that this condition enactment portion 20 sets.
In such structure, in film image image pickup part 12, be reconstructed from subject 10 acquisition time image series as the object position.Then, in feature portion test section 14, detect the specific part that with the characteristic of the time series image of above-mentioned film image image pickup part 12 collections, for example detects the heart of subject 10.Obtain at the displacement of the specific part of this detected heart and characteristic quantity thereof as the kinetic characteristic value with motion analysis portion 16.Then, in extraction unit 18, extract specific time series image according to above-mentioned kinetic characteristic value.
Then, according to specific time series image or the incidental information that extract with this extraction unit 18, in condition enactment portion 20, carry out the setting of the synchronous or photography conditions of coronary artery photography.According to the photography conditions of setting with this condition enactment portion 20 etc., carry out the shooting of the time phase of necessary heart by image pickup part 22.
Secondly, with reference to the flow chart of Fig. 6, the action of the MRI device of the 2nd embodiment is described.
In addition, because the step S1～S4 in the flow chart of the processing action of the step S11～S14 in this flow chart and above-mentioned Fig. 3 is identical, as moving, in this description will be omitted with reference to the processing of corresponding processing action.
In step S14,, then in the step S15 that follows, carry out the setting of the synchronous or photography conditions of coronary artery photography if extract specific time series image according to above-mentioned kinetic characteristic value.At this, by the operator import as for example in Fig. 4 B, represent like that the time delay of the incidental information 742 shown in the part of image 74 (in this case, being 800ms) as identifying information.Perhaps, also can use at the number of the white circle (coronary artery) shown in the image 74 etc. and set photography conditions as identifying information.
If in step S16, determined photography conditions like this, then in step S17, use based on making a video recording the time phase of the heart of the photography conditions that sets.About this shooting, the camera method of above-mentioned SSFP sequence or the shooting of being undertaken by any camera methods such as SE (spin echo) method, EPI (echo-planar imaging) methods for example can be to use.
In addition, when the setting of the photography conditions that carries out above-mentioned steps S15, can set photography conditions according to such performance plot of representing among Fig. 5.In addition, when showing, not only displayable image 74, and the performance plot of expression in the displayed map 5 simultaneously.
In addition, in the above-described embodiment, suppose that the characteristic of the time series image at the object position of detecting subject illustrates, but be not limited thereto,, then also can come the detected characteristics part from basic image if data volume is no more than allowance.
Moreover, in the above-described embodiment, film image having been carried out shooting back detected characteristics portion, extracted image, but also can use sample image to detect similar degree (characteristic) based on the kinetic characteristic value.Embodiments of the present invention more than have been described, but the present invention also can implement various distortion except above-mentioned embodiment in the scope that does not break away from main idea of the present invention.
Moreover, comprised the invention in each stage in the above-described embodiment, utilize the suitable combination of disclosed a plurality of constitutive requirements can extract various inventions.For example, even the several constitutive requirements of deletion also can solve the problem narrated in the part of narrating the problem that invention plan to solve, also can obtain narrating under the situation of the effect of narrating in the part of effect of invention in whole constitutive requirements of representing from embodiment, the structure of having deleted these constitutive requirements also can be used as invention and extracts.
According to the present invention, can provide can be easily and to obtain coronary artery exactly be static or move MR imaging apparatus and image analysis method and image analysis program during few.
For those skilled in the art, can easily realize additional advantage of the present invention and modification.Thereby, the present invention be not limited to here represent aspect wideer at it with the specific details and the representational embodiment that describe.Therefore, under situation about not departing from, can do various corrections by the spirit and scope of the inventive concept of accompanying Claim and the universality of the present invention that equivalent limited thereof.
Claims (10)
1. MR imaging apparatus, magnetostatic field takes place in this device, for the overlapping leaning magnetic field that utilizes gradient magnetic field coil to take place of this magnetostatic field, subject in the above-mentioned magnetostatic field is detected the magnetic resonance signal that the high frequency magnetic field be used to spontaneous sending coil produces in above-mentioned subject, obtain image from this magnetic resonance signal, it is characterized in that possessing:
Collect the film image image pickup part (1) that the time series image at the object position of above-mentioned subject is reconstructed;
To the feature portion test section (2) that detects as one or more closed-loop configuration on the motion picture pickup cross section of the characteristic of above-mentioned time series image;
To the motion analysis portion (3) that analyzes with the kinetic characteristic value of the detected characteristic of above-mentioned feature portion's test section (2); And
Extract the extraction unit (4) about the information of specific time phase of above-mentioned kinetic characteristic value smaller or equal to the image of threshold value, this threshold value is the permissible value of motion.
2. MR imaging apparatus as claimed in claim 1 is characterized in that:
The information of extracting with above-mentioned extraction unit (4) about above-mentioned specific time phase is specific time series image.
3. MR imaging apparatus as claimed in claim 2 is characterized in that:
Also possesses the display part (22) that shows the above-mentioned time series image that extracts with above-mentioned extraction unit (4).
4. MR imaging apparatus as claimed in claim 1 is characterized in that:
Above-mentioned feature portion's test section (2) detects this closed-loop configuration under the big slight situation of above-mentioned closed-loop configuration ratio regulation.
5. MR imaging apparatus as claimed in claim 1 is characterized in that:
Also possess be used for according to the information that extracts with above-mentioned extraction unit (4) about above-mentioned specific time phase carry out the coronary artery photography synchronously and the condition enactment portion (5) of the setting of photography conditions and use the image pickup part (6) of making a video recording by the photography conditions of above-mentioned condition enactment portion (5) setting.
6. MR imaging apparatus as claimed in claim 1 is characterized in that:
Above-mentioned feature portion's test section (2) detects right coronary artery and left coronary artery as characteristic.
7. the image analysis method of a MR imaging apparatus, this MR imaging apparatus generation magnetostatic field, for the overlapping leaning magnetic field that utilizes gradient magnetic field coil to take place of this magnetostatic field, subject in the above-mentioned magnetostatic field is detected the magnetic resonance signal that the high frequency magnetic field be used to spontaneous sending coil produces in above-mentioned subject, obtain image from this magnetic resonance signal, this image analysis method is characterised in that:
The time series image of collecting the object position of above-mentioned subject is reconstructed,
To detecting as one or more closed-loop configuration on the motion picture pickup cross section of the characteristic of above-mentioned time series image,
Kinetic characteristic value to above-mentioned detected characteristic is analyzed,
Extract the information about specific time phase of kinetic characteristic value smaller or equal to the image of threshold value, this threshold value is the permissible value of motion.
8. the image analysis method of MR imaging apparatus as claimed in claim 7 is characterized in that:
The above-mentioned information about above-mentioned specific time phase that is extracted is specific time series image.
9. the image analysis method of MR imaging apparatus as claimed in claim 8 is characterized in that:
And then show the above-mentioned above-mentioned time series image that is extracted out.
10. the image analysis method of MR imaging apparatus as claimed in claim 7 is characterized in that:
In above-mentioned detection, under the big slight situation of above-mentioned closed-loop configuration ratio regulation, detect this closed-loop configuration.
Applications Claiming Priority (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007027050 | 2007-02-06 | ||
| JP2007027050 | 2007-02-06 | ||
| JP2007-027050 | 2007-02-06 | ||
| JP2007327683 | 2007-12-19 | ||
| JP2007327683A JP2008212634A (en) | 2007-02-06 | 2007-12-19 | Magnetic resonance imaging apparatus and image analysis method and image analysis program therefor |
| JP2007-327683 | 2007-12-19 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101238979A CN101238979A (en) | 2008-08-13 |
| CN101238979B true CN101238979B (en) | 2011-08-10 |
Family
ID=39833386
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2008100045914A Expired - Fee Related CN101238979B (en) | 2007-02-06 | 2008-01-25 | MRI involving a cine prescan for motion analysis |
Country Status (3)
| Country | Link |
|---|---|
| JP (2) | JP2008212634A (en) |
| CN (1) | CN101238979B (en) |
| DE (1) | DE602008003041D1 (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9921285B2 (en) | 2012-04-19 | 2018-03-20 | New York University | System, method and computer-accessible medium for highly-accelerated dynamic magnetic resonance imaging using golden-angle radial sampling and compressed sensing |
| US9671482B2 (en) | 2012-10-18 | 2017-06-06 | Samsung Electronics Co., Ltd. | Method of obtaining image and providing information on screen of magnetic resonance imaging apparatus, and apparatus thereof |
| IN2013CH05587A (en) * | 2013-12-04 | 2015-06-12 | Gen Electric | |
| DE102014206724B4 (en) | 2014-04-08 | 2015-11-12 | Siemens Aktiengesellschaft | Dynamic imaging with variable contrast |
| US11284811B2 (en) * | 2016-06-22 | 2022-03-29 | Viewray Technologies, Inc. | Magnetic resonance volumetric imaging |
| JP7123767B2 (en) * | 2018-11-20 | 2022-08-23 | キヤノンメディカルシステムズ株式会社 | Magnetic resonance imaging device |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6292683B1 (en) * | 1999-05-18 | 2001-09-18 | General Electric Company | Method and apparatus for tracking motion in MR images |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU8523291A (en) * | 1990-08-17 | 1992-03-17 | Trustees Of The University Of Pennsylvania, The | Analyzing heart wall motion using spatial modulation of magnetization |
| US7031504B1 (en) * | 2000-09-26 | 2006-04-18 | Vital Images, Inc. | Image data based retrospective temporal selection of medical images |
| AU2003279484A1 (en) * | 2002-12-13 | 2004-07-09 | Koninklijke Philips Electronics N.V. | System and method for processing a series of image frames representing a cardiac cycle |
| JP2005278919A (en) * | 2004-03-30 | 2005-10-13 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
| JP4507081B2 (en) * | 2004-08-02 | 2010-07-21 | 横河電機株式会社 | Image processing method, image processing apparatus, and program for image processing |
| US8352013B2 (en) * | 2005-01-18 | 2013-01-08 | Siemens Medical Solutions Usa, Inc. | Method and system for motion compensation in magnetic resonance (MR) imaging |
-
2007
- 2007-12-19 JP JP2007327683A patent/JP2008212634A/en active Pending
-
2008
- 2008-01-23 DE DE200860003041 patent/DE602008003041D1/en active Active
- 2008-01-25 CN CN2008100045914A patent/CN101238979B/en not_active Expired - Fee Related
-
2013
- 2013-11-06 JP JP2013230598A patent/JP5684363B2/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6292683B1 (en) * | 1999-05-18 | 2001-09-18 | General Electric Company | Method and apparatus for tracking motion in MR images |
Non-Patent Citations (1)
| Title |
|---|
| K Nehrke, P Bornert, T Netsch.Automatic Selection of Cardiac Acquisition Window Using anImage-Based Global Cross-Correlation of Multi Heart PhaseCine Scans.Proc Int Soc Magn Reson Med.2003,1623. * |
Also Published As
| Publication number | Publication date |
|---|---|
| DE602008003041D1 (en) | 2010-12-02 |
| JP2008212634A (en) | 2008-09-18 |
| CN101238979A (en) | 2008-08-13 |
| JP2014023963A (en) | 2014-02-06 |
| JP5684363B2 (en) | 2015-03-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9804245B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| EP1825808B1 (en) | Magnetic resonance imaging system | |
| US20170127972A1 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| US6144201A (en) | MR imaging utilizing ECG gating technique | |
| JP4912154B2 (en) | Nuclear magnetic resonance imaging system | |
| JP5854575B2 (en) | Magnetic resonance imaging system | |
| JP5100181B2 (en) | Magnetic resonance imaging system | |
| US9684979B2 (en) | MRI 3D cine imaging based on intersecting source and anchor slice data | |
| US8417007B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| US6505064B1 (en) | Diagnostic imaging systems and methods employing temporally resolved intensity tracing | |
| EP1615044A1 (en) | Continuous table motion MRI involving variable fields of view during the scan | |
| CN101238979B (en) | MRI involving a cine prescan for motion analysis | |
| EP2350698B1 (en) | System and method for moving table mri | |
| US8488856B2 (en) | Magnetic resonance imaging apparatus and method of analyzing images provided thereby | |
| CN101138497B (en) | Magnetic resonance imaging apparatus | |
| JP3895972B2 (en) | Magnetic resonance imaging device | |
| CN101558995B (en) | Medical imaging apparatus and medical display image generation method | |
| JP2000350714A (en) | Mri inspection device | |
| JP4349647B2 (en) | Magnetic resonance imaging system | |
| KR102386797B1 (en) | Apparatus and method for magnet resonance imaging using phase contrast flow | |
| JP5289292B2 (en) | Magnetic resonance imaging device | |
| JP4519827B2 (en) | Magnetic resonance imaging device | |
| WO2021183989A1 (en) | Free-breathing abdominal magnetic resonance fingerprinting using a pilot tone navigator | |
| JPWO2003101294A1 (en) | Magnetic resonance imaging apparatus and multi-station CE-MRA method |
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 | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20160720 Address after: Japan Tochigi Patentee after: TOSHIBA MEDICAL SYSTEMS Corp. Address before: Tokyo, Japan Patentee before: Toshiba Corp. Patentee before: TOSHIBA MEDICAL SYSTEMS Corp. |
|
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110810 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |