EP1500368B1 - Magnetic resonance imaging device - Google Patents
Magnetic resonance imaging device Download PDFInfo
- Publication number
- EP1500368B1 EP1500368B1 EP03727983A EP03727983A EP1500368B1 EP 1500368 B1 EP1500368 B1 EP 1500368B1 EP 03727983 A EP03727983 A EP 03727983A EP 03727983 A EP03727983 A EP 03727983A EP 1500368 B1 EP1500368 B1 EP 1500368B1
- Authority
- EP
- European Patent Office
- Prior art keywords
- sensitivity
- image data
- receiving coil
- sensitivity distribution
- receiving
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
Definitions
- the present invention relates to a magnetic resonance imaging apparatus (hereinafter referred to as "MRI apparatus") for imaging a tissue structure inside an object to be examined by utilizing a magnetic resonance phenomenon, more particularly to an MRI apparatus suitable for parallel imaging for rapidly imaging the object using a plurality of receiving coils.
- MRI apparatus magnetic resonance imaging apparatus
- a pulse sequence is sequentially executed while a phase encoding amount given to a nuclear magnetic resonance signal (hereinafter referred to as "NMR signal") is varied, and an echo signal necessary for reconstructing one image is thus acquired. Therefore, imaging time depends on time of repetition (TR) and the repetition number of times of pulse sequence.
- TR time of repetition
- a multi-echo pulse sequence in which plural echo signals can be measured by executing the pulse sequence one time and a pulse sequence in which the TR is shortened to several ms to several ten ms are generally used.
- the parallel imaging method is a method of reducing the repetition number of pulse sequence by executing the pulse sequence using a plurality of receiving coils while phase encoding steps are thinned out at even intervals, and thus shortening the imaging time.
- a measured signal obtained by performing the parallel imaging pulse sequence is simply image-reconstructed, an aliasing artifact is generated in the reconstructed image.
- the aliasing artifact is generated due to a difference between a field of view (FOV) of image to be obtained and that in which signals are actually measured in the parallel imaging, and it can be removed by performing a matrix calculation using a pre-calculated sensitivity distribution of each receiving coil.
- FOV field of view
- the imaging time can be shortened depending on the number of receiving coils. That is, it is possible in principle to thin out the number of phase encoding steps as much as the number of receiving coils from the number of phase encoding steps in the usual imaging. As a result, for example, the imaging time can be shortened 1/N times, where the number of receiving coil is N, in comparison with the imaging performed with the usual phase encoding step number.
- a sensitivity distribution of each receiving coil is needed for removing the aliasing artifact generated in the image.
- the pulse sequence for the sensitivity distribution acquisition on every slice position to be imaged Therefore, there is a problem that the total imaging time in the parallel imaging becomes very long.
- Techniques of shortening the total imaging time in the parallel imaging includes one in which the sensitivity distribution of receiving coil is acquired by performing three-dimensional measurement. However, because the sensitivity distribution data are similarly acquired on every imaging slice position, the effect of shortening the time is small, and it is desired that the time be further shortened.
- JP-A-11 262 479 discloses an MRI system using multiple RF coils with distinct spatial sensitivity distributions.
- an MRI system containing means for acquiring examination image data from a plurality of m slices sequentially adjoining on an object, using a phased-array coil. It also contains sensitivity image data acquiring means arranged for acquiring sensitivity data for n slices, with n ⁇ m.
- An object of the present invention is to provide an MRI apparatus which can shorten the total imaging time when plural slices of the imaging portion of the object are pictured using the parallel imaging method.
- a data component of one portion is wrapped and superposed on that of another portion. Accordingly, by forming a simultaneous equation in accordance with a correlation among the image data of the plural receiving coils, the sensitivity distribution data of each receiving coil, and the data components on plural points of the imaging cross section and solving these simultaneous equations on the basis of matrix calculation, the superposing data component at the wrapped portion is separated. The aliasing artifact on the image is removed by expanding thus separated image data.
- the sensitivity distribution data on each slice position of image data are necessary. These sensitivity distribution data are calculated from the sensitivity image data acquired by executing the first pulse sequence for acquiring measured data of central portion (low-frequency portion) in the k space, whereby the time for acquiring the sensitivity distribution of receiving coil can be shortened.
- the slice number of the sensitivity image data acquired in the first pulse sequence is n in number, which is smaller than m in number of examination image data including the aliasing artifact component acquired in the second pulse sequence.
- sensitivity distribution data on a slice position do not exist in the image data for m slices, they are generated by performing interpolation using the above acquired sensitivity distribution data. Therefore, because the repetition number of the first pulse sequence for acquiring the sensitivity distribution of receiving coil can be reduced, the total time of repeatedly executing the first pulse sequence is shortened, and the total imaging time in the parallel imaging including the first pulse sequence and the second pulse sequence can be thus shortened.
- the MRI apparatus includes means for acquiring sensitivity image data of each receiving coil for plural slices on slice positions different from each other by executing the first sequence and acquiring aliasing image data of each receiving coil for plural images by executing the second sequence, and means for generating the sensitivity distribution data of each receiving coil on the slice positions of the above sensitivity image data on the basis thereof and removing the aliasing artifact generated in the examination image on the basis of the sensitivity distribution data on the slice positions of the aliasing image data.
- the plurality of receiving coils may include a first receiving coil and a second receiving coil.
- the sensitivity distribution of the first receiving coil is more uniform than that of the second receiving coil.
- an image having a more uniform sensitivity distribution than that generated on the basis of the measured data acquired by a single receiving coil is generated by combining the measured data obtained by using a plurality of receiving coils.
- this image having a uniform sensitivity distribution the sensitivity distribution of each receiving coil can be accurately calculated.
- Fig.1 is a diagram showing a procedure of signal processing of a signal processing unit in an MRI apparatus according to a first embodiment of the present invention.
- Fig.2 is a diagram showing a schematic structure of the MRI apparatus according to the first embodiment of the present invention.
- Fig.3 is a diagram illustrating a receiving coil in the MRI apparatus according to the first embodiment of the present invention.
- Fig.4 is a diagram showing a schematic structure of the multiple RF receiving coils and a reception unit in the MRI apparatus according to the first embodiment of the present invention.
- Fig.5 is a chart illustrating a pulse sequence used in the MRI apparatus according to the present invention.
- Fig.6 is a diagram illustrating an aliasing artifact to be removed in the parallel imaging method in the MRI apparatus according to the present invention.
- Fig.7 is a diagram illustrating image data acquired in the pulse sequence and those generated by a slice interpolation processing in the MRI apparatus according to the first embodiment of the present invention.
- Fig.8 is a diagram showing a procedure of signal processing in the signal processing unit in the MRI apparatus according to a second embodiment of the present invention.
- Fig.9 is a diagram illustrating sensitivity distribution data generated by a sensitivity distribution calculation processing and a slice interpolation processing and image data acquired in the pulse sequence in the MRI apparatus according to the second embodiment of the present invention.
- Fig.10 is a diagram illustrating a receiving coil in the MRI apparatus according to a third embodiment of the present invention.
- Fig.11 is a diagram showing a schematic structure of multiple receiving coils and a reception unit in the MRI apparatus according to the third embodiment of the present invention.
- the MRI apparatus includes magnet 402 for generating a uniform static magnetic field having a predetermined strength to a space (hereinafter referred to as "imaging space") having a predetermined size for accommodating object 401 to be examined, gradient magnetic field coil 403 for generating gradient magnetic fields to the imaging space, body coil 404 for generating a radio-frequency magnetic field to the imaging space and detecting an NMR signal generated by object 401, and RF probe 405 for detecting the NMR signal generated by object 401.
- imaging space a space having a predetermined size for accommodating object 401 to be examined
- gradient magnetic field coil 403 for generating gradient magnetic fields to the imaging space
- body coil 404 for generating a radio-frequency magnetic field to the imaging space and detecting an NMR signal generated by object 401
- RF probe 405 for detecting the NMR signal generated by object 401.
- RF is an abbreviation of "Radio Frequency.”
- Fig.2 shows an MRI apparatus of horizontal magnetic field type which generates a static magnetic field in a horizontal direction
- the present invention can also be applied to an MRI apparatus of vertical magnetic field type which generates the static magnetic field in the vertical direction.
- Gradient magnetic field coil 403 includes a plurality of coils for generating gradient magnetic fields in X-, Y-, and Z-directions, each of which generates the gradient magnetic field in response to a signal from gradient magnetic field power source 409. These gradient magnetic fields are superposed on the static magnetic field and used for selecting an imaging slice and encoding the NMR signal.
- Body coil 404 is designed to generate an RF magnetic field (RF pulse) in response to a signal output from RF transmission unit 410. Because conductors of body coil 404 are arranged so that its sensitivity distribution becomes substantially uniform, it is also used to receive the NMR signal used to obtain a reference signal for calculating a sensitivity distribution of receiving coil in this embodiment.
- the RF magnetic field is used to excite nuclear spins in the body of object 401 placed in the imaging space.
- the coil for obtaining the reference signal for calculating the sensitivity distribution of receiving coil may be a dedicated receiving coil having a substantially uniform sensitivity distribution, instead of the coil for both transmission and reception.
- the image data dealt with in the present invention includes examination image data used in diagnosis on the object and sensitivity image data for calculating a sensitivity distribution of coil. Among them, the examination image data used in diagnosis of the object is displayed as an image on display unit 408.
- Gradient magnetic field power source 409, RF transmission unit 410, signal processing unit 407, and display unit 408 are controlled by control unit 411.
- a time chart of the control of control unit 411 is generally referred to as "pulse sequence”.
- Object 401 is laid on bed 412, which is designed to move between the imaging space inside the magnets and the outside of the magnets.
- RF probe 405 includes two figure-of-eight shaped coils 4051 and 4052 arranged in the inside space of body coil 405, as shown in Fig.3 .
- the two receiving coils 4051 and 4052 are arranged opposite at a predetermined interval in the Z direction along an X-Y plane, which receive the NMR signal from the object placed inside the space surrounded by these two receiving coils 4051 and 4052.
- elements and processings given a four-digit number having an identical first three digits represent an identical element and processing, and the number of their last one digit is given to distinguish each of plural identical elements and processings. Accordingly, they are hereinafter equivalently treated in this specification.
- RF probe 405 is constructed with a technique referred in the field of MRI as "multiple receiving coils” or “phased array coil” including a plurality of receiving coils.
- Two receiving coils 4051 and 4052 of multiple receiving coils 301 are connected with pre-amplifiers 3021 and 3022 as shown in Fig.4 .
- the reception unit for measured signals received by multiple receiving coils 301 has signal detection unit 406 and signal processing unit 407.
- signal detection unit 406 two AD conversion/quadrature detection circuits 3031 and 3032 are arranged in parallel, each of which is connected to outputs of pre-amplifiers 3021 and 3022.
- Signal processing unit 407 is designed to generate an MR image from signals received by receiving coils 4051 and 4052 by utilizing a Fourier transformation method, a back projection method, a wavelet transformation method, and the like, which includes signal combination processing unit 304 for performing the combination processing on those image signals.
- the pulse sequence for the imaging may be a general spin-echo pulse sequence or a gradient-echo pulse sequence. Furthermore, it may be a spin-echo type multi-slice pulse sequence or a gradient-echo type multi-slice sequence.
- RF pulse 601 applied here is provided a frequency on the imaging slice position determined depending on a magnetic field gradient of the slice-selecting gradient magnetic field superposed on the static magnetic field, and a frequency band determined depending on a slice thickness.
- phase encoding gradient magnetic field pulse 603. After the nuclear spins inside the slice are excited, the excited nuclear spins are given positional information in a phase direction by the application of phase encoding gradient magnetic field pulse 603. After that, echo signal 606 is measured during the application of readout gradient magnetic field pulse 604, i.e. during sampling window 605. This pulse sequence is repeated plural times at time of repetition (TR) 607.
- TR 607 is, for example, an interval between applications of RF pulse 6011 in a first pulse sequence and of RF pulse 6012 in the next pulse sequence.
- the phase encoding step number determined in accordance with the pixel number in the phase encoding direction of the image to be examined is sequentially measured in sampling window 605.
- a set of echo signals necessary for reconstructing one image is measured. That is, in the conventional imaging method, the number of echo signal is determined as 64, 128, 256, 512, or the like in consideration of spatial resolution of field of view (hereinafter referred to as "FOV") in the phase encoding direction. Further, the echo signals are obtained as time-sequential sampling data usually including 128, 256, 512, 1024, signals or so in consideration of spatial resolution of FOV in the readout direction.
- k-space two-dimensional memory space
- a plurality of adjoining slices of, e.g. m in number are imaged on the examining portion of object 401. Accordingly, the above pulse sequence is sequentially executed on each of those slice positions. Meanwhile, the multi-slice imaging method may be employed here as above.
- the pulse sequence is modified for executing the parallel imaging. That is, according to this embodiment, pre-measurement of MR signal is performed with two receiving coils, i.e. body coil 404 and receiving probe 405 to rapidly perform the parallel imaging.
- a pulse sequence of pre-measurement in which the size of phase encoding matrix, i.e. the number of phase encoding step is reduced to be one of several fractions of that in the usual imaging is executed to acquire sensitivity image data for n slices included in the examining portion of the object.
- the sensitivity image data for n slices the data for at least two slices corresponds to the imaging slice position, where the slice number is determined as n ⁇ m.
- the pulse sequence shown in Fig.5 may be applied to the slice positions where the sensitivity image data is acquired, or the multi-slice pulse sequence may be applied so as to measure signals from the slice positions where the sensitivity image data is acquired.
- a pulse sequence of imaging in which the phase encoding steps in the k space are thinned out is executed to acquire image data for m slices (m>n) of the examining portion of object 401 using receiving coils 4051 and 4052 forming multiple receiving coils 405.
- the image data for m slices are acquired from sequentially adjoining slice positions.
- the multi-slice pulse sequence or a pulse sequence for three-dimensional measurement may be used instead of the pulse sequence shown in Fig.5
- the sensitivity image data acquired by body coil 404 and receiving probe 405 and the image data acquired by receiving coils 4051 and 4052 of receiving probe 405 are respectively subjected to two-dimensional Fourier transformation and converted into a sensitivity image and an examination image.
- An aliasing artifact is generated in the examination image reconstructed by the two-dimensional Fourier transformation, which is removed with the following processings.
- the sensitivity distribution data of receiving coils 4051 and 4052 are calculated from the sensitivity image obtained with the pre-measurement pulse sequence, an interpolation processing is performed using the calculated sensitivity distribution data to calculate the sensitivity distribution data corresponding to each slice position of the examination image, the sensitivity distribution of receiving coil is calculated using the calculated sensitivity distribution data, and a matrix calculation is performed using this sensitivity distribution and the examination image, whereby the aliasing artifact on the examination image is removed.
- the repetition number of pulse sequence is reduced by thinning out the intervals between phase encoding steps at a fixed rate.
- This thinning ratio is called "acceleration factor" of the parallel imaging.
- the acceleration factor is 2. That is, when the phase encoding step number is 64, all data of 1st to 64th phase encodes are usually measured. Meanwhile, in the parallel imaging with the acceleration factor of 2, data of 1st, 3rd, 5th, ..., and 63rd phase encodes are measured. At this time, the thinning out ratio is 1/2.
- Control unit 411 executes the real-measurement pulse sequence in which the phase encoding steps are thus thinned out to acquire the image data for m slices with receiving coils 4051 and 4052.
- Fig.6(a) shows the case where imaging is done with usual phase encoding steps (usual imaging), the number of which is not reduced.
- Signals 2021, 2022, ..., and 2027 obtained with each phase encoding amount are arranged in the k space to form measured data 201 for one image.
- Fig.6(c) an image without aliasing is obtained as shown in Fig.6(c) .
- Fig.6(b) in the parallel imaging measurement is performed while the interval between the phase encoding steps are e.g. doubled in the same FOV as in the usual imaging to thin out the phase encoding steps.
- data are measured on every other lines 2041, 2042, ..., and 2043 in a ky direction of the k space, and data on lines 2051, 2052, ..., 2053, and 2054 are not measured. Accordingly, the amount of measured data becomes half and the measured FOV also becomes 1/2.
- An MR signal from a portion of the object out of the measured FOV becomes an aliasing component.
- image 208 in which an aliasing artifact is generated in the phase encoding direction is obtained as shown in Fig.6(b) . That is, object's image 2061 indicated in upper portion 2071 of image 207 and object's image 2062 in lower portion 2072 as shown in Fig.6(c) become aliasing and appear as aliasing image 208 shown in Fig.6(d) .
- SENSE Sensitivity Encoding for Fast MRI ( Klass P. Pruessmann et. al), Magnetic Resonance in Medicine 42: 952-962(1999 ), etc.
- a calculation method of aliasing removal will be described. If the image matrices in the x and y directions are X and Y respectively, a pixel value on a coordinate (x, y) (1 ⁇ x ⁇ X, 1 ⁇ y ⁇ Y) in the image is si (x, y).
- a subscript "i " is the number of receiving coil, where 2 ⁇ i ⁇ N.
- received signal si(x, y) is a product of sensitivity distribution ci(x, y) of receiving coil and proton density distribution p(x, y) of the object, which can be expressed by formula (2).
- s i x y c i x y ⁇ p x y
- formula (1) can be expressed by formula (3) using formula (2).
- Formula (5) can be expressed by formula (6) as a matrix of Nx2.
- S 1 S 2 ⁇ S N 1 a C 11 C 12 C 21 C 22 ⁇ ⁇ C N ⁇ 1 C N ⁇ 2 ⁇ P 1 P 2
- a nucleate factor in measuring the object e.g. proton density distribution P j can be calculated by calculating the inverse matrix.
- a general formula in performing imaging at M-fold speed using N coils is can be expressed by formula (7), where Y' ⁇ Y/ M and 1 ⁇ y' ⁇ Y'.
- character b is a fixed number.
- S 1 S 2 ⁇ S N 1 b C 11 C 12 ... C 1 ⁇ M C 12 C 22 ... C 2 ⁇ M ⁇ ⁇ ⁇ ⁇ C N ⁇ 1 C N ⁇ 2 ... C NM ⁇ P 1 P 2 ⁇ P M
- sensitivity distribution C ij of each of receiving coils 4051 and 4052 is calculated on the basis of sensitivity image data generated from the measured data obtained in advance of the low-frequency band in the phase encoding direction of the k space in the pre-measurement pulse sequence.
- the sensitivity image data for calculating sensitivity distribution C ij of receiving coils 4051 and 4052 are acquired by executing the pre-measurement sequence.
- the interval between the phase encoding steps are not thinned out but the same as in the usual imaging.
- the repetition number of pulse sequence is reduced. For example, when the phase encoding step number in the k space is 64, in the usual imaging all data of 1 st to 64th phase encode steps are measured.
- the phase encoding step number is set to be the same as in the usual imaging, and the data on 16 steps including 25th to 40th phase encode steps are measured in one portion of the k space, e.g. in the central portion (low-frequency area) of the k space.
- the phase encoding step number becomes one fourth of that in the usual imaging, the time for pre-measurement is shortened.
- Control unit 411 executes the pre-measurement pulse sequence to acquire the data to be measured of the central portion being about one fourth of the k space so as to acquire the sensitivity image data for n slices of body coil 404 and receiving coils 4051 and 4052. In this manner, the sensitivity image data for one slice are measured with the small phase encoding step number and the sensitivity image data are acquired with the smaller number of slice position than that in the examination image, whereby the total imaging time of parallel imaging is shortened.
- signal processing unit 407 performs slice interpolation processings 1011, 1012, and 1013 on sensitivity image 701 of n in number of body coil 404 and on sensitivity images 702 and 703, each of n in number, respectively of receiving coils 4051 and 4052 of multiple receiving coils 301 obtained by executing the pre-measurement pulse sequence for sensitivity image data acquisition, and thus generates an interpolated sensitivity image so that all sensitivity images corresponding to the respective slice positions of examination image complete.
- sensitivity distribution data 707 and 708 respectively of receiving coils 4051 and 4052 are calculated by sensitivity distribution calculation processings 7061 and 7062 using a whole sensitivity image including the sensitivity image generated by slice interpolation processing 1011 of body coil 404 and that including the sensitivity image generated by slice interpolation processings 1012 and 1013 of receiving coils 4051 and 4052 of multiple receiving coils 301.
- Calculated sensitivity distribution data 707 and 708 of receiving coils 4051 and 4052 are stored into a memory (not shown) mounted in control unit 411.
- the above sensitivity distribution calculation processings 7061 and 7062 are the processing of dividing the sensitivity images of receiving coils 4051 and 4052 by the sensitivity image of body coil 404 while slice positions are matched.
- a determinant is formed by matrix forming processing 709 on the basis of obtained examination images 704 and 705, each for m slices, of receiving coils 4051 and 4052 of multiple receiving coils 301 and sensitivity distribution data 707 and 708, inverse matrix calculation processing 710 of the determinant is performed, and diagnostic image 711 without an aliasing is resultingly obtained.
- sensitivity images 701, 702, and 703, each for n slices, obtained by executing the pre-measurement are sensitivity images 7011, 7014, ..., and 701m of body coil 404, sensitivity images 7021, 7024, ..., and 702m of receiving coil 4051 of multiple receiving coil 301, and sensitivity images 7031, 7034, ..., 703m of receiving coil 4052 of multiple receiving coils 301, as shown in Fig.7 .
- examination images 704 and 705 including an aliasing artifact obtained by executing the real-measurement pulse sequence are examination images 7041 to 704m of receiving coil 4051 of multiple receiving coils 301 and examination images 7051 to 705m of receiving coil 4052 of multiple receiving coils 301.
- ch.1 corresponds to receiving coil 4051 of multiple receiving coils 301
- ch.2 corresponds to receiving coil 4052 of multiple receiving coils 301.
- Each of those sensitivity images 701, 702, and 703 and examination images 704 and 705 is given a slice number corresponding to its slice position, and images having an identical slice number are those on an identical slice position.
- the slice number corresponds to the last one digit of image number.
- sensitivity images 7011, 7014, ..., and 701m of body coil 404 are respectively given the slice numbers 1, 4, ..., and m. Further, slice intervals of the images obtained in the pre-measurement sequence and the real-measurement pulse sequence are evenly determined.
- sensitivity image data 701 of body coil 404 and sensitivity image data 702 and 703 of receiving coils 4051 and 4052 are not acquired on slice numbers 2, 3, 5, ..., and m-1. Accordingly, in examination images 7042, 7043, 7045, ..., and 704m-1, and 7052, 7053, 7055, ..., and 705m-1 of slice numbers 2, 3, 5, ..., and m-1, the aliasing artifact cannot be removed without sensitivity distribution information on those slice positions.
- the signal processing unit calculates sensitivity images 7012, 7013, 7015, ..., and 701m-1 by performing slice interpolation processing 1011 from sensitivity images 7011, 7014, ..., and 701m of body coil 404, similarly calculates sensitivity image 7022, 7023, 7025, ..., and 702m-1 by performing slice interpolation processing 1012 from obtained sensitivity images 7021, 7024, ..., and 702m of receiving coil 4051, and calculates sensitivity images 7032, 7033, 7035, ..., and 703m-1 by performing slice interpolation processing 1013 from sensitivity images 7031, 7034, ..., and 703m of receiving coil 4052.
- reference number 501 in Fig.7 represents a set of sensitivity image necessary for removing the aliasing artifact in the examination image including the sensitivity images acquired in the pre-measurement pulse sequence and the sensitivity image data acquired by slice interpolation processings 1011, 1012, and 1013
- reference number 502 represents a set of examination image including the aliasing artifact acquired by the real-measurement pulse sequence.
- white sensitivity images in sensitivity image set 501 represent those actually acquired by in the pre-measurement pulse sequence and black ones represent those generated by slice interpolation processing 101.
- a determinant of formula (7) is formed by matrix forming processing 709. After that, diagnostic image 711 in which an aliasing artifact is removed is obtained with reverse matrix calculation processing 710.
- phase encoding step number is set to be one of several fractions of that in the usual imaging in the pre-measurement pulse sequence, imaging for the sensitivity image data and that for the aliasing image data are performed for the same number of slices and on the same slice positions.
- the time for executing the pre-measurement occupies a large proportion of the total imaging time.
- the execution time of the pre-measurement pulse sequence is a major matter in consideration of the efficiency of examination.
- sensitivity image data 701, 702, and 703 are acquired for every several slices, e.g. for every three slices to reduce the slice number of the imaging sensitivity image data to be m being smaller than the slice number m of examination image. Accordingly, the repetition number of the pre-measurement pulse sequence is reduced, so that the imaging time can be shortened for thinned sensitivity image data 701, 702, and 703. Further, when this embodiment is applied to three-dimensional measurement, the slice encoding step number in the pulse sequence for sensitivity distribution acquisition can be reduced.
- the sensitivity image data generated by slice interpolation processing 101 are accurately calculated with a simple calculation.
- setting and control of pulse sequence in control unit 411 is simplified since the slice number is thinned out at even intervals to acquire the sensitivity image data.
- the sensitivity distributions of receiving coils 4051 and 4052 are stored in the memory of control unit 411. Accordingly, when it is necessary to image again 1st to mth slices while slightly shifting the slice positions thereof, it is possible to calculate the sensitivity distribution of each coil on the slice positions to be re-measured by using the sensitivity distribution stored in the memory, without re-measuring them. Therefore, the time for re-measurement can be shortened.
- slice interpolation processing 801 is performed after sensitivity distribution calculation processing 706 is performed using sensitivity image 701 of body coil 404 and sensitivity images 702 and 703 of receiving coils 7051 and 4052 of multiple receiving coils 301 acquired in the pre-measurement pulse sequence for the sensitivity image data acquisition.
- sensitivity distribution data 7061 and 7062 of receiving coils 4051 and 4052 are calculated from sensitivity image 701 of body coil 404 and sensitivity images 702 and 703 of receiving coils 4051 and 4052, the slice interpolation processing 801 is performed next using those sensitivity distribution data 7061 and 7062 to acquire sensitivity distribution data 707 and 708 of receiving coils 4051 and 4052, matrix forming processing 709 is performed after that from those sensitivity distribution data 707 and 708 and examination images 704 and 705, and the aliasing artifact in examination images 704 and 705 are removed by performing inverse matrix calculation processing 710 on the formed matrix.
- sensitivity distribution data of receiving coils 4051 and 4052 calculated by sensitivity distribution calculation processing 706 and those generated by slice interpolation processing 801 will be described.
- Sensitivity distribution data calculated by sensitivity distribution calculation processing 706 are, as shown in Fig.9 , sensitivity distribution data 7071, 7074, ..., and 707m of receiving coil 4051 of multiple receiving coils 301 and sensitivity distribution data 7081, 7084, ..., and 708m of receiving coil 4052.
- the number of measured slices are reduced by setting measuring slice positions as every three slices as in the first embodiment. That is, the sensitivity image data of slice numbers 1, 4, ..., and m are acquired and those of slice numbers 2, 3, 5, ..., and m-1 are not acquired.
- the aliasing artifact appearing in examination images 7042, 7043, 7045, ..., and 704m-1 and 7052, 7053, 7055, ..., and 705m-1 cannot be removed without the sensitivity distribution data of slice numbers 2, 3, 5, ..., and m-1.
- sensitivity distribution data 7072, 7073, 7075, ..., and 707m-1 are calculated by performing slice interpolation processing 8011 on acquired sensitivity distribution data 7071, 7074, ..., and 707m of receiving coil 4051, and sensitivity distribution data 7082, 7083, 7085, ...., and 708m-1 are calculated by performing slice interpolation processing 8012 on sensitivity distribution data 7081, 7084, ..., and 708m of receiving coil 4052.
- reference number 503 in Fig.9 represents a data set necessary for removing the aliasing artifact in the examination image, including sensitivity distribution data calculated by sensitivity distribution calculation processings 7061 and 7062 and those obtained by slice interpolation processings 8011 and 8012.
- the white sensitivity image in data set 503 represents the sensitivity distribution data calculated by sensitivity distribution calculation processings 7061 and 7062, and the black sensitivity distribution data represent sensitivity distribution data generated with by interpolation processing 801.
- signal processing unit 407 performs slice interpolation processings 8011 and 8012 after performing sensitivity distribution calculation processings 7061 and 7062 on the basis of sensitivity image data 701, 702, and 703 acquired in the pre-measurement pulse sequence, forms sensitivity distribution data 7071 to 707m and 7081 to 708m corresponding to each slice position of examination images 7041 to 704m and 7051 to 705m, and after that removes the aliasing artifact in the image by matrix forming processing 709 and inverse matrix calculation processing 710.
- slice interpolation processing 801 is performed after calculating sensitivity distributions of receiving coils 4051 and 4052 on the basis of sensitivity images 701, 702, and 703 acquired in the first pulse sequence, it is unnecessary to perform the slice interpolation processing on sensitivity image 701 of body coil 404. Therefore, in comparison with First embodiment where slice interpolation processing 1011 is performed on sensitivity image 701 of body coil 404, it is possible to simplify sensitivity distribution calculation processing 706 and shorten the time for image reconstruction.
- multiple receiving coils include four receiving coils. That is, as shown in Fig.10 , multiple receiving coils 301A forming an RF probe in the MRI apparatus according to this embodiment include four figure-of-eight shaped receiving coils 4051 to 4054, which are arranged opposite on the X-Y plane and on the T-Z plane at a predetermined distance.
- each receiving coil 4051 to 4054 of multiple RF receiving coils 301A is respectively connected to preamplifiers 3021 to 3024.
- signal detecting unit 406A four AD conversion/quadrature detection circuits 3031 to 3034 are arranged in parallel, which are respectively connected to outputs of preamplifiers 3021 to 3024.
- the matrix calculation of formula (7) may diverge to generate an artifact or deteriorate image quality.
- optimum combination of receiving coil is a combination with which the matrix calculation does not diverge. For example, there may be such conditions that there is not a portion of low sensitivity distribution in the phase encoding direction, that sensitivity distributions of receiving coil groups are not identical, and the like.
- two receiving coils 4051 and 4052 are arranged opposite in the Z direction on the X-Y plane at a predetermined distance. Further, two receiving coils 4053 and 4054 are oppositely arranged in the X direction on a plane parallel to the Z-Y plane at a predetermined distance. Among those four receiving coils, receiving coils 4051 and 4053 are combined to form receiving coil group 4055, and receiving coils 4052 and 4054 are combined to form receiving coil group 4056.
- the receiving coils When the receiving coils are thus arranged and combined, the combined sensitivity distribution data of receiving coils 4051 and 4053 and that of receiving coils 4052 and 4054 become different from each other and the number of low sensitivity portions decreases, whereby the matrix calculation does not easily diverge.
- the apparatus is desirably constructed so that the optimum combination is automatically selected when those imaging conditions are determined.
- sensitivity image data measured by the two selected combinations of receiving coils are picked up, and those sensitivity image data of each set are combined.
- two sets of examination image data measured by the two selected combinations of receiving coils are picked up and those examination image data of each set are combined.
- the sensitivity image data of receiving coil group 4055 generated by thus combining those data correspond to sensitivity image data 7021, 7024, ..., and 702m in ch.1 of Fig.7
- those of receiving coil group 4056 correspond to sensitivity image data 7031, 7034, ..., and 703m in ch.2.
- the combined examination image data of receiving coil group 4055 correspond to examination image data 7041 to 704m in ch.1 of Fig.7 and those of receiving coil group 4056 correspond to examination image data 7051 to 705m in ch.2.
- the sensitivity image data acquired by receiving coils 4051 to 4054 are combined in signal combination processing unit 304A into sensitivity image data 702 and 703 in ch.1 and ch.2 corresponding to the sensitivity distribution of each predetermined receiving coil group, and also made into reference data of the matrix calculation for removing the aliasing artifact in examination image data 704 and 705 of ch.1 and ch.2 combined in each receiving coil group.
- the combination processing based on is done in signal combination processing unit 304A.
- the combination processing may be done by calculating a root sum square of two sets of sensitivity image data and that of two sets of examination image data. The latter case brings an advantage of higher SN ratio in comparison with the complex addition.
- the signals are combined after the AD conversion.
- combining processing may be done on analog signals before being subjected to the AD conversion, and the AD conversion is performed afterward on them.
- one receiving coil may be overlapped in both the first and the second combinations. That is, the plural combinations of receiving coils may include a common receiving coil. Since the plural combinations of receiving coils may thus share a part of receiving coils, the coil combinations with which an optimum sensitivity distribution is obtainable can be more flexibly selected and a better quality image is obtainable.
- the present invention is not limited to the above disclosed embodiments and may employ various embodiments in consideration of the gist of the invention.
- First to Third embodiments are described with regard to the case where the pre-measurement pulse sequence for the sensitivity distribution acquisition is performed on every three slices, it may be done on every arbitrary number of slices.
- First to Third embodiments are described with regard to the case where the slice intervals in acquiring the image data are equivalent, they may be changed to acquire the image data in the MRI apparatus according to the present invention.
- the number of receiving coils is not limited to two and four. That is, it may be eight or other number. However, to simplify the circuitry of the apparatus and to reduce the object's uncomfortable feeling, it is considered that the number of receiving coils is desirably smaller.
- the step interval may be set as three times, four times, or more.
- G the number of selected combination of receiving coils to be the same number as interval M of the phase encoding steps or more.
- First to Third embodiments are described with regard to the case of employing the gradient echo sequence, the type of the parallel imaging and the measurement method are not limited.
- the SE pulse sequence the FSE pulse sequence, the EPI pulse sequence, the SSFP method, and the spiral measurement method.
- the MRI apparatus is applicable not only to the multi-slice measurement method of these embodiments but also to the three-dimensional measurement method.
- the data may be thinned out not in the phase encoding direction but in the slice encoding direction, or in both the phase encoding direction and the slice encoding direction in combination.
- the combination of receiving coil may be selected so that the combination of sensitivity distributions is optimum in the axis direction in which the encoding steps are thinned out.
- slice interpolation processings 101 and 801 may be, e.g. a spline interpolation processing, an interpolation processing based on a Sinc function, a function-based fitting processing, and the like.
- an approximate sensitivity distribution of receiving coil can be calculated by generating a sensitivity image data having a uniform sensitivity distribution by combining sensitivity image data acquired by plural receiving coils instead of the body coil according to those embodiment, and dividing the image data acquired by each receiving coil by thus combined sensitivity image data.
- the sensitivity distribution of receiving coil can be accurately calculated by dividing the image data acquired by each receiving coil by the image data of body coil having a relatively uniform sensitivity distribution described in the above embodiments and calculating the approximate sensitivity distribution of receiving coil. Therefore, the latter is preferable.
- the multiple receiving coils of the MRI apparatus are not limited to those according to First to Third embodiments, and various multiple receiving coils may be utilized, such as ones mentioned in “Array Head Coil for Improved Functional MRI” (Christoph Leussler), 1996 ISMRM abstract p.249, ones mentioned in “Helmet and Cylindrical Shaped CP Array Coils for Brain Imaging: A Comparison of Signal - to - Noise Characteristics” (H. A. Stark, E. M. Haacke), 1996 ISMRM abstract P.1412 as a head QD multiple RF coils for horizontal magnetic field, and ones mentioned in “Four Channel Wrap-Around Coil with Inductive Decoupler for 1.5T Body Imaging”(T. Takahashi et.
- the multiple receiving coils are designed to enlarge an FOV while maintaining high sensitivity of receiving coil to improve the whole sensitivity by combining signals acquired by each of arranged small-sized receiving coils having a relatively high sensitivity. Accordingly, it is possible to enlarge the FOV and improve the whole sensitivity by combining the reconstructed images of each of those receiving coils into one image.
- First to Third embodiments are described for the case of utilizing a horizontal magnetic field type MRI apparatus in which the static magnetic field direction (Z direction) is horizontal, the present invention is not limited thereto and also applicable to an MRI apparatus in which the static magnetic field direction is vertical.
- the sensitivity distribution of receiving coil for removing an aliasing artifact generated in an examination image obtained in the parallel imaging can be acquired by measuring data on fewer slices than imaging slices, it can be acquired within a short time. Therefore, the total imaging time in the parallel imaging can be shortened. Further, since an unacquired sensitivity distribution of the receiving coil on imaging slices is calculated with the interpolation calculation, the calculation can be easily performed.
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002129152 | 2002-04-30 | ||
JP2002129152 | 2002-04-30 | ||
PCT/JP2003/005103 WO2003092497A1 (fr) | 2002-04-30 | 2003-04-22 | Dispositif d'imagerie par resonance magnetique |
Publications (3)
Publication Number | Publication Date |
---|---|
EP1500368A1 EP1500368A1 (en) | 2005-01-26 |
EP1500368A4 EP1500368A4 (en) | 2006-04-12 |
EP1500368B1 true EP1500368B1 (en) | 2011-09-14 |
Family
ID=29397295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP03727983A Expired - Fee Related EP1500368B1 (en) | 2002-04-30 | 2003-04-22 | Magnetic resonance imaging device |
Country Status (5)
Country | Link |
---|---|
US (1) | US7224163B2 (ja) |
EP (1) | EP1500368B1 (ja) |
JP (1) | JP4152381B2 (ja) |
CN (1) | CN100369580C (ja) |
WO (1) | WO2003092497A1 (ja) |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4443079B2 (ja) | 2001-09-13 | 2010-03-31 | 株式会社日立メディコ | 磁気共鳴イメージング装置及び磁気共鳴イメージング装置用rf受信コイル |
US7330027B2 (en) * | 2003-07-02 | 2008-02-12 | Universitat Zurich | System and method of magnetic resonance imaging for producing successive magnetic resonance images |
CN100396240C (zh) * | 2005-02-28 | 2008-06-25 | 西门子(中国)有限公司 | 频域灵敏度编码磁共振并行成像方法 |
EP1790991B1 (en) * | 2005-11-28 | 2011-11-16 | Kabushiki Kaisha Toshiba | Data processing system and data processing method involving a Wiener filter |
JP4350729B2 (ja) * | 2006-06-27 | 2009-10-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | イメージング装置およびイメージング方法 |
AT504060B1 (de) * | 2006-07-24 | 2010-03-15 | Univ Graz Tech | Vorrichtung zur magnetischen induktionstomografie |
DE102006061177B4 (de) * | 2006-12-22 | 2009-04-02 | Siemens Ag | 3D-MR-Bildgebung mit Fettunterdrückung |
DE102007004620B4 (de) * | 2007-01-30 | 2010-02-04 | Siemens Ag | Verbessertes dreidimensionales schichtselektives Mehrschicht-Anregungsverfahren in der MRT-Bildgebung |
JP4995650B2 (ja) * | 2007-06-22 | 2012-08-08 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
CN101470180B (zh) * | 2007-12-29 | 2016-01-20 | 西门子(中国)有限公司 | 磁共振成像中失真校准的方法和装置 |
CN101482600B (zh) * | 2008-01-08 | 2011-10-19 | 西门子(中国)有限公司 | 一种磁共振成像系统的接收装置 |
EP2283373B1 (en) * | 2008-04-28 | 2021-03-10 | Cornell University | Accurate quantification of magnetic susceptibility in molecular mri |
JP5955891B2 (ja) * | 2008-07-01 | 2016-07-20 | 東芝メディカルシステムズ株式会社 | 磁気共鳴イメージング装置および受信経路切り替え方法 |
EP2233941A1 (en) * | 2009-03-25 | 2010-09-29 | Koninklijke Philips Electronics N.V. | MR imaging using parallel signal acquisition |
DE102009053522B4 (de) * | 2009-11-16 | 2012-04-05 | Siemens Aktiengesellschaft | Phasenabhängige Magnetresonanzbilderzeugung mit mehreren Spulen |
JP5721941B2 (ja) * | 2009-11-20 | 2015-05-20 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴イメージング装置 |
DE102010041191B4 (de) * | 2010-09-22 | 2016-02-18 | Siemens Aktiengesellschaft | Erstellung von MR-Bilddaten mit paralleler Schichtanregung und teilweiser Überlappung der Schichten im Frequenzbereich |
DE102011081411B4 (de) * | 2011-08-23 | 2013-04-11 | Friedrich-Alexander-Universität Erlangen-Nürnberg | Abtastmuster für iterative MR-Rekonstruktionsverfahren |
CN103874457B (zh) * | 2011-09-29 | 2016-02-03 | 株式会社日立医疗器械 | 磁共振摄像装置以及高频磁场照射方法 |
JP6042069B2 (ja) * | 2012-01-25 | 2016-12-14 | 東芝メディカルシステムズ株式会社 | 磁気共鳴イメージング装置 |
WO2013182967A1 (en) * | 2012-06-05 | 2013-12-12 | Koninklijke Philips N.V. | Channel by channel artifact reduction in parallel mri |
EP3250935A1 (en) | 2015-01-27 | 2017-12-06 | Koninklijke Philips N.V. | Parallel multi-slice mr imaging with suppression of side band artefacts |
US10132883B2 (en) * | 2016-05-31 | 2018-11-20 | General Electric Company | Foldable coil array |
CN106772167B (zh) * | 2016-12-01 | 2019-05-07 | 中国科学院深圳先进技术研究院 | 核磁共振成像方法及装置 |
JP6820876B2 (ja) | 2018-03-20 | 2021-01-27 | 株式会社日立製作所 | 磁気共鳴撮影装置、感度分布算出方法、及び、感度分布算出プログラム |
JP2022092150A (ja) * | 2020-12-10 | 2022-06-22 | キヤノンメディカルシステムズ株式会社 | 画像生成装置、画像生成方法、および画像生成プログラム |
JP2022178683A (ja) * | 2021-05-20 | 2022-12-02 | 富士フイルムヘルスケア株式会社 | 磁気共鳴撮影装置及び感度分布算出プロブラム |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4770182A (en) * | 1986-11-26 | 1988-09-13 | Fonar Corporation | NMR screening method |
JP3676853B2 (ja) | 1994-08-03 | 2005-07-27 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 表面コイル配設により核磁化分布を決定するmr法 |
JP3526337B2 (ja) * | 1995-01-20 | 2004-05-10 | 株式会社東芝 | 磁気共鳴イメージング装置 |
JP2952228B1 (ja) * | 1998-03-18 | 1999-09-20 | 技術研究組合医療福祉機器研究所 | 連続mri画像再構成方法及び装置 |
US6289232B1 (en) * | 1998-03-30 | 2001-09-11 | Beth Israel Deaconess Medical Center, Inc. | Coil array autocalibration MR imaging |
GB9926918D0 (en) * | 1999-11-15 | 2000-01-12 | Marconi Electronic Syst Ltd | Magnetic resonance imaging |
US6242916B1 (en) * | 1999-11-17 | 2001-06-05 | General Electric Company | Partial fourier acquistion of MR data over a limited field of view and image reconstruction |
JP3952247B2 (ja) | 1999-12-08 | 2007-08-01 | 株式会社日立メディコ | 核磁気共鳴撮影装置 |
US6777934B2 (en) * | 1999-12-08 | 2004-08-17 | Hitachi Medical Corporation | Magnetic resonance imaging method and apparatus |
US6556009B2 (en) * | 2000-12-11 | 2003-04-29 | The United States Of America As Represented By The Department Of Health And Human Services | Accelerated magnetic resonance imaging using a parallel spatial filter |
DE10119660B4 (de) * | 2001-04-20 | 2006-01-05 | Siemens Ag | Verfahren zur schnellen Gewinnung eines Magnetresonanzbildes |
US6559642B2 (en) * | 2001-05-09 | 2003-05-06 | Ge Medical Systems Global Technology Company, Llc | Calibration method for use with sensitivity encoding MRI acquisition |
-
2003
- 2003-04-22 CN CNB038098334A patent/CN100369580C/zh not_active Expired - Lifetime
- 2003-04-22 US US10/511,964 patent/US7224163B2/en not_active Expired - Lifetime
- 2003-04-22 WO PCT/JP2003/005103 patent/WO2003092497A1/ja active Application Filing
- 2003-04-22 EP EP03727983A patent/EP1500368B1/en not_active Expired - Fee Related
- 2003-04-22 JP JP2004500689A patent/JP4152381B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
US20060049829A1 (en) | 2006-03-09 |
EP1500368A1 (en) | 2005-01-26 |
US7224163B2 (en) | 2007-05-29 |
CN1649540A (zh) | 2005-08-03 |
CN100369580C (zh) | 2008-02-20 |
JPWO2003092497A1 (ja) | 2005-09-08 |
JP4152381B2 (ja) | 2008-09-17 |
EP1500368A4 (en) | 2006-04-12 |
WO2003092497A1 (fr) | 2003-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1500368B1 (en) | Magnetic resonance imaging device | |
US6876201B2 (en) | Magnetic resonance imaging apparatus and method | |
KR100553464B1 (ko) | 자기 공명 화상화 방법 및 장치 | |
EP3004908B1 (en) | Parallel mri with b0 distortion correction and multi-echo dixon water-fat separation using regularised sense reconstruction | |
US6559642B2 (en) | Calibration method for use with sensitivity encoding MRI acquisition | |
EP2615470A1 (en) | MR imaging with B1 mapping | |
US9223001B2 (en) | MR imaging using navigators | |
JP2003079595A (ja) | 磁気共鳴イメージング装置及び磁気共鳴イメージング装置用rf受信コイル | |
JP4047553B2 (ja) | 磁気共鳴イメージング装置 | |
JP2023109791A (ja) | 並列マルチスライスmr撮像 | |
US20030055330A1 (en) | Sensitivity encoding MRI acquisition method | |
US10955506B2 (en) | Parallel MR imaging with spectral fat suppression | |
JP3992934B2 (ja) | 磁気共鳴イメージング装置および方法 | |
JP4980662B2 (ja) | 磁気共鳴イメージング装置 | |
EP2581756A1 (en) | MR imaging using parallel signal acquisition | |
EP3185029A1 (en) | Mr imaging using propeller acquisition with t2 decay correction | |
EP3931588A1 (en) | Epi mr imaging with distortion correction | |
CN113544526B (zh) | 使用波编码的并行mr成像 | |
Banerjee et al. | Fast magnetic resonance spectroscopic imaging at 3 Tesla using autocalibrating parallel technique | |
CORBIN et al. | Fast Imaging and Acceleration Techniques | |
Takizawa et al. | Parallel imaging of head with a dedicated multi-coil on a 0.4 T open MRI | |
Öhman | The influence of the reference measurement in MRI image reconstruction using Sensitivity Encoding (SENSE) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
17P | Request for examination filed |
Effective date: 20041029 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PT RO SE SI SK TR |
|
RIC1 | Information provided on ipc code assigned before grant |
Ipc: 7A 61B 5/055 B Ipc: 7G 01R 33/561 A |
|
A4 | Supplementary search report drawn up and despatched |
Effective date: 20060228 |
|
17Q | First examination report despatched |
Effective date: 20081021 |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): DE FR GB IT NL |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R096 Ref document number: 60338361 Country of ref document: DE Effective date: 20111110 |
|
REG | Reference to a national code |
Ref country code: NL Ref legal event code: T3 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20110914 |
|
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
26N | No opposition filed |
Effective date: 20120615 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R097 Ref document number: 60338361 Country of ref document: DE Effective date: 20120615 |
|
GBPC | Gb: european patent ceased through non-payment of renewal fee |
Effective date: 20120422 |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: ST Effective date: 20121228 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: GB Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20120422 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: FR Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20120430 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: NL Payment date: 20140410 Year of fee payment: 12 Ref country code: DE Payment date: 20140430 Year of fee payment: 12 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R119 Ref document number: 60338361 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: NL Ref legal event code: MM Effective date: 20150501 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: DE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20151103 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: NL Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20150501 |