CN1663526A - Method for removing Gibbs ring pseudo-image generated in MRI reconstruction process - Google Patents
Method for removing Gibbs ring pseudo-image generated in MRI reconstruction process Download PDFInfo
- Publication number
- CN1663526A CN1663526A CN 200510033773 CN200510033773A CN1663526A CN 1663526 A CN1663526 A CN 1663526A CN 200510033773 CN200510033773 CN 200510033773 CN 200510033773 A CN200510033773 A CN 200510033773A CN 1663526 A CN1663526 A CN 1663526A
- Authority
- CN
- China
- Prior art keywords
- matrix
- data
- image
- frequency domain
- data matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Analysis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
Claims (5)
- The removing method of the Gibbs ring artifact that produces when 1, a kind of magnetic resonance image (MRI) is rebuild is characterized in that may further comprise the steps:(1) obtains the original frequency domain vectors matrix data of the part that is used for nuclear magnetic resonance, i.e. Partial K spatial data;(2) be the Kx direction with the frequency coding direction, phase-encoding direction is the Ky direction, is decided to be Ky when phase encoding gradient Gy is zero and equals zero and determine the K space coordinates; The cyclophysis of utilization Fourier transform, with each original frequency domain vectors matrix conversion of obtain become in the K space coordinates with phase-encoding direction Ky equal zero be axis of symmetry grip symmetric matrix data 1 up and down altogether;(3) frequency domain vectors matrix data 1 is carried out a low-pass filtering and obtain filtered frequency domain vectors matrix data 2;(4) utilize 2 pairs of images of frequency domain vectors matrix data to carry out rim detection;(5) be traditional Fourier imaging equation that the image function equation is rewritten into the second class Chebyshev multinomial with fourier progression expanding method be that basic function carries out unfolded image imaging equation;(6) marginal information of obtaining by step (4) in each continuum, is brought frequency domain vectors matrix data 1 into try to achieve in the step (5) image imaging equation, obtains the image area data matrix asked with this;(7) will obtain the image area data matrix at last, and vector will be converted into vector, the image that obtains revising.
- The removing method of the Gibbs ring artifact that produces when 2, magnetic resonance image (MRI) according to claim 1 is rebuild, it is characterized in that: the concrete steps of described frequency domain filtering are as follows:(1) selection index low pass filter in the filtering, wherein the attenuation rate coefficient is selected between 4~8, and the filtering strength factor alpha is that empirical value is generally selected between 25~50;(2) according to above-mentioned rule structure filtering data matrix, matrix size is consistent with frequency domain vectors matrix data 1;(3) wait until filtered frequency domain vectors matrix data 2 with multiply each other corresponding of frequency domain vectors data matrix 1 with filtering matrix.
- The removing method of the Gibbs ring artifact that produces when 3, magnetic resonance image (MRI) according to claim 1 is rebuild, it is characterized in that: described rim detection is to carry out rim detection respectively to obtain a marginal information on the direction on frequency coding Kx direction and phase code Ky direction, and concrete steps are as follows:(1) on frequency coding Kx direction, selected concentration factor τ, the concentration factor data matrix of formation;(2) utilize concentration factor structure to concentrate kernel function K, make and concentrate kernel function K to have following characteristic: carry out the jump that its value behind the convolution algorithm converges on function to be detected with image function;(3) to concentrating the nuclear matrix data to carry out Fourier transform, with data matrix after the conversion and frequency domain data matrix 2 corresponding multiplying each other, the data matrix that obtains carries out inverse fourier transform again;(4) setting threshold is got rid of the false discontinuous point in the detection, has only test value just to be judged to be real testing result greater than the result of threshold value; The size of threshold value is 9% of the maximum jump value of function;(5) matrix data that carries out obtaining after the threshold test is the final result of frequency coding Kx direction edge detection;(6) according to above-mentioned steps phase code Ky direction is detected matrix data after obtaining detecting;
- The removing method of the Gibbs ring artifact that produces when 4, magnetic resonance image (MRI) according to claim 1 is rebuild, it is characterized in that: the concrete steps of construct image imaging equation are as follows:(1) respectively in frequency coding Kx direction, the selected second class Chebyshev multinomial U (x) on the phase code Ky direction, U (y);(2) data matrix that utilizes Chebyshev multinomial U (x) to constitute carries out Fourier and changes and obtain transformation matrix w (x), and data matrix U (y) carries out Fourier transform and obtains transformation matrix w (y);(3) on frequency coding direction and phase-encoding direction, the transform data matrix is done inverse transformation and obtain w -1(x), w -1(y);(4) the inverse transformation data matrix w that obtains -1(x), w -1(y) obtain image imaging equation expansion coefficient matrix g with frequency domain vectors data matrix 1 corresponding multiplying each other respectively l(x), g l(y);(5), use function expansion coefficient data matrix g on the phase code Ky direction respectively in frequency coding Kx direction l(x), g l(y) with Chebyshev polynomial data matrix U (x), U (y) correspondence multiplies each other and can obtain image imaging equation G on the frequency coding direction l(x) and the image imaging equation G on the phase-encoding direction l(y).
- The removing method of the Gibbs ring artifact that produces when 5, magnetic resonance image (MRI) according to claim 1 is rebuild, it is characterized in that: the process that obtains the image area data matrix is: respectively in frequency coding Kx direction, on the phase code Ky direction, by the marginal information of obtaining, calling graph is as imaging equation in the continuum, obtain the image area data matrix F (x) on frequency coding direction and the phase-encoding direction, F (y); With the image area data matrix F (x) on the both direction, F (y) does weighted average and obtains final image numeric field data matrix F.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100337730A CN1327810C (en) | 2005-03-28 | 2005-03-28 | Method for removing Gibbs ring pseudo-image generated in MRI reconstruction process |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100337730A CN1327810C (en) | 2005-03-28 | 2005-03-28 | Method for removing Gibbs ring pseudo-image generated in MRI reconstruction process |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1663526A true CN1663526A (en) | 2005-09-07 |
CN1327810C CN1327810C (en) | 2007-07-25 |
Family
ID=35034965
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100337730A Expired - Fee Related CN1327810C (en) | 2005-03-28 | 2005-03-28 | Method for removing Gibbs ring pseudo-image generated in MRI reconstruction process |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1327810C (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100409807C (en) * | 2005-11-29 | 2008-08-13 | 东南大学 | Artifact correction method based on even marker in plane echo imaging technique |
CN101908204A (en) * | 2010-05-25 | 2010-12-08 | 南方医科大学 | Inverse diffusion method for eliminating Gibbs annular artifact of magnetic resonance image |
CN101995561A (en) * | 2010-10-22 | 2011-03-30 | 南方医科大学 | Reconstruction method of periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) magnetic resonance data based on image domain overlapping |
CN102096055A (en) * | 2010-12-14 | 2011-06-15 | 南方医科大学 | Rapid and accurate reconstructing method for non-uniform sampling data of magnetic resonance imaging |
CN102324089A (en) * | 2011-07-13 | 2012-01-18 | 南方医科大学 | Maximum posteriori reconstruction method of PET (positron emission tomography) image based on generalized entropy and MR (magnetic resonance) prior |
CN103942763A (en) * | 2014-05-03 | 2014-07-23 | 南方医科大学 | Voxel level PET (positron emission tomography) image partial volume correction method based on MR (magnetic resonance) information guide |
CN104020430A (en) * | 2014-04-15 | 2014-09-03 | 清华大学 | Correction method and system for movement artifacts of magnetic resonance imaging |
CN104635187A (en) * | 2015-02-28 | 2015-05-20 | 深圳先进技术研究院 | Adaptive active automatic shimming method and system of magnetic resonance system |
CN110133556A (en) * | 2019-05-29 | 2019-08-16 | 上海联影医疗科技有限公司 | A kind of magnetic resonance image processing method, device, equipment and storage medium |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IL86570A (en) * | 1988-05-31 | 1991-07-18 | Elscint Ltd | Reduction of truncation caused artifacts |
US5285157A (en) * | 1991-05-03 | 1994-02-08 | Elscint Ltd. | Reduction of truncation caused artifacts |
CN1050427C (en) * | 1995-06-17 | 2000-03-15 | 浙江大学 | Fuzzy mutticriteria image reconstruction method |
-
2005
- 2005-03-28 CN CNB2005100337730A patent/CN1327810C/en not_active Expired - Fee Related
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100409807C (en) * | 2005-11-29 | 2008-08-13 | 东南大学 | Artifact correction method based on even marker in plane echo imaging technique |
CN101908204B (en) * | 2010-05-25 | 2012-01-04 | 南方医科大学 | Inverse diffusion method for eliminating Gibbs annular artifact of magnetic resonance image |
CN101908204A (en) * | 2010-05-25 | 2010-12-08 | 南方医科大学 | Inverse diffusion method for eliminating Gibbs annular artifact of magnetic resonance image |
CN101995561A (en) * | 2010-10-22 | 2011-03-30 | 南方医科大学 | Reconstruction method of periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) magnetic resonance data based on image domain overlapping |
CN102096055B (en) * | 2010-12-14 | 2013-01-23 | 南方医科大学 | Rapid and accurate reconstructing method for non-uniform sampling data of magnetic resonance imaging |
CN102096055A (en) * | 2010-12-14 | 2011-06-15 | 南方医科大学 | Rapid and accurate reconstructing method for non-uniform sampling data of magnetic resonance imaging |
CN102324089A (en) * | 2011-07-13 | 2012-01-18 | 南方医科大学 | Maximum posteriori reconstruction method of PET (positron emission tomography) image based on generalized entropy and MR (magnetic resonance) prior |
CN104020430A (en) * | 2014-04-15 | 2014-09-03 | 清华大学 | Correction method and system for movement artifacts of magnetic resonance imaging |
CN104020430B (en) * | 2014-04-15 | 2017-01-25 | 清华大学 | Correction method and system for movement artifacts of magnetic resonance imaging |
CN103942763A (en) * | 2014-05-03 | 2014-07-23 | 南方医科大学 | Voxel level PET (positron emission tomography) image partial volume correction method based on MR (magnetic resonance) information guide |
CN104635187A (en) * | 2015-02-28 | 2015-05-20 | 深圳先进技术研究院 | Adaptive active automatic shimming method and system of magnetic resonance system |
CN104635187B (en) * | 2015-02-28 | 2017-07-07 | 深圳先进技术研究院 | The adaptive active automatic shimming method and system of magnetic resonance system |
CN110133556A (en) * | 2019-05-29 | 2019-08-16 | 上海联影医疗科技有限公司 | A kind of magnetic resonance image processing method, device, equipment and storage medium |
CN110133556B (en) * | 2019-05-29 | 2021-01-19 | 上海联影医疗科技股份有限公司 | Magnetic resonance image processing method, device, equipment and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN1327810C (en) | 2007-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1663526A (en) | Method for removing Gibbs ring pseudo-image generated in MRI reconstruction process | |
CN110490832A (en) | A kind of MR image reconstruction method based on regularization depth image transcendental method | |
Cheon et al. | Generative adversarial network-based image super-resolution using perceptual content losses | |
CN105631807B (en) | The single-frame image super-resolution reconstruction method chosen based on sparse domain | |
CN111932460B (en) | MR image super-resolution reconstruction method, device, computer equipment and storage medium | |
CN107576925B (en) | Magnetic resonance multi-contrast image reconstruction method and device | |
CN112396672B (en) | Sparse angle cone-beam CT image reconstruction method based on deep learning | |
CN107274462A (en) | The many dictionary learning MR image reconstruction methods of classification based on entropy and geometric direction | |
CN106485764A (en) | The quick exact reconstruction methods of MRI image | |
CN105118053B (en) | A kind of full reference picture assessment method for encoding quality based on compressed sensing | |
CN104700440B (en) | Magnetic resonant part K spatial image reconstruction method | |
CN113160380B (en) | Three-dimensional magnetic resonance image super-resolution reconstruction method, electronic equipment and storage medium | |
Wang et al. | High‐field mr diffusion‐weighted image denoising using a joint denoising convolutional neural network | |
US8634671B2 (en) | Methods and apparatus to perform multi-focal plane image acquisition and compression | |
Brun et al. | Effective implementation of ring artifacts removal filters for synchrotron radiation microtomographic images | |
CN110942496A (en) | Propeller sampling and neural network-based magnetic resonance image reconstruction method and system | |
CN117523095A (en) | Sparse angle THz-CT image reconstruction method based on deep learning | |
CN112150360A (en) | IVUS image super-resolution reconstruction method based on dense residual error network | |
Schrank et al. | Deep learning MPI super-resolution by implicit representation of the system matrix | |
CN114897691A (en) | Image compressed sensing reconstruction method and system based on standardized stream | |
CN115239836A (en) | Extreme sparse view angle CT reconstruction method based on end-to-end neural network | |
Abbas et al. | PET image reconstruction and denoising on hexagonal lattices | |
CN116934618B (en) | Image halftone method, system and medium based on improved residual error network | |
CN112801886B (en) | Dynamic PET image denoising method and system based on image wavelet transformation | |
CN117036533B (en) | Sparse angle CT reconstruction method and system for coaxial phase contrast imaging |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract |
Assignee: Ningbo Xingaoyi Magnetism Co., Ltd. Assignor: Southern Medical University Contract fulfillment period: 2007.11.15 to 2013.11.14 contract change Contract record no.: 2008330000178 Denomination of invention: Method for removing Gibbs ring pseudo-image generated in MRI reconstruction process Granted publication date: 20070725 License type: Exclusive license Record date: 2008.8.15 |
|
LIC | Patent licence contract for exploitation submitted for record |
Free format text: EXCLUSIVE LICENCE; TIME LIMIT OF IMPLEMENTING CONTACT: 2007.11.15 TO 2013.11.14 Name of requester: NINGBO XIN GAO YI MAGNETIC CO., LTD. Effective date: 20080815 |
|
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070725 Termination date: 20150328 |
|
EXPY | Termination of patent right or utility model |