CN111445412A - Two-dimensional geometric correction method for magnetic resonance image - Google Patents

Two-dimensional geometric correction method for magnetic resonance image Download PDF

Info

Publication number
CN111445412A
CN111445412A CN202010225659.2A CN202010225659A CN111445412A CN 111445412 A CN111445412 A CN 111445412A CN 202010225659 A CN202010225659 A CN 202010225659A CN 111445412 A CN111445412 A CN 111445412A
Authority
CN
China
Prior art keywords
image
magnetic resonance
deformation
points
control points
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.)
Pending
Application number
CN202010225659.2A
Other languages
Chinese (zh)
Inventor
袁双虎
李玮
李莉
韩毅
刘宁
张云
胡建建
袁朔
吕慧颖
于金明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jinan Bishan Network Technology Co ltd
Shandong Cancer Hospital & Institute (shandong Cancer Hospital)
Beijing Yikang Medical Technology Co ltd
Shandong University
Original Assignee
Jinan Bishan Network Technology Co ltd
Shandong Cancer Hospital & Institute (shandong Cancer Hospital)
Beijing Yikang Medical Technology Co ltd
Shandong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jinan Bishan Network Technology Co ltd, Shandong Cancer Hospital & Institute (shandong Cancer Hospital), Beijing Yikang Medical Technology Co ltd, Shandong University filed Critical Jinan Bishan Network Technology Co ltd
Priority to CN202010225659.2A priority Critical patent/CN111445412A/en
Publication of CN111445412A publication Critical patent/CN111445412A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/80Geometric correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

A two-dimensional geometric correction method for magnetic resonance images is characterized in that when geometric distortion occurs in a magnetic resonance water model image, an image obtained by using an improved automatic window adjusting algorithm is used, the water model is used for displaying the image, basic parameters of longitudinal and transverse stretching states are further obtained by combining the displayed definition and contrast with the automatic window adjusting algorithm, and the automatic window adjusting algorithm is used for carrying out audit correction so that the shape deviation caused by organ self-weight is improved. Moreover, the problems can be effectively solved, and clear and rich-level magnetic resonance images can be quickly and conveniently obtained.

Description

Two-dimensional geometric correction method for magnetic resonance image
Technical Field
The invention relates to the field of magnetic resonance imaging, in particular to an algorithm for two-dimensional geometric correction of a magnetic resonance image.
Background
An algorithm for two-dimensional geometric correction of magnetic resonance images is a digital image processing algorithm that corrects images that are geometrically distorted in magnetic resonance. The geometric distortion exists widely at the edge of a television screen, the edge part of a camera lens and a magnetic resonance image, and the invention particularly refers to the geometric distortion of the magnetic resonance image. The geometric distortion of the magnetic resonance image is mainly caused by the nonlinearity of a gradient field generated on the hardware of a magnetic resonance imaging instrument, the distortion is more serious in a region which is farther away from the center of a magnet, and the difference is also large according to the personal size. Geometric correction algorithms for magnetic resonance images are of great help for medical diagnosis.
The general algorithm for solving the problem of geometric distortion of the magnetic resonance image is to utilize magnetic resonance to perform imaging experiment, estimate the deformation quantity of a control point by adopting a proper algorithm and further perform geometric correction by utilizing geometric transformation of digital image processing. The algorithm is characterized in that the algorithm can be predicted according to the position information of the control points of the magnetic resonance image, the actual positions of the control points are used for determining the primary positioning of the positions of the control points in the image, and the positions of the control points in the image are accurately searched near the positions. However, when the visual field of the image is large, the image deformation is also large, the control points may deform from the line where the control points are located to another line, and at this time, the above algorithm cannot effectively extract the corresponding control points, and may use the horizontal and vertical control points of another line as the control points of the current line, which is the first defect, the algorithm can perform initial positioning on the control points with small deformation amount, and when the initial positioning is performed on the points with large deformation amount, the obtained points are the points of the adjacent lines or the points of the adjacent columns; this geometry correction algorithm cannot estimate the amount of deformation outside the control point coverage using a simple interpolation algorithm, which is a second drawback.
Disclosure of Invention
Aiming at the defects that the prior art can not effectively extract the control points of the magnetic resonance image and can not estimate the deformation field outside the coverage range of the control points of the magnetic resonance image, the invention provides an algorithm for geometric correction of the magnetic resonance image, which comprises the following steps:
the method comprises the following steps: acquiring experimental data and recording the number and position information of grid points of a water model by using a magnetic resonance imaging image and the water model as the number and position information of initial transverse and longitudinal control points of the image;
step two: acquiring a plurality of water model images by utilizing different areas of the water model effect of the magnetic resonance image, and extracting the number and position information of control points in each image, wherein the position information of the transverse control points and the longitudinal control points is the position of the control points in the image after the image is deformed; the control points of the image correspond to the stretching grid points of the transverse direction and the longitudinal direction of the water model image;
step three: combining the plurality of images according to different regions of the magnetic resonance image, calculating position information and deformation quantity of the transverse and longitudinal stretching grid points, namely the control points, and estimating a deformation field according to the deformation quantity;
step four: and correcting the magnetic resonance image with geometric deformation by using an improved automatic window adjusting algorithm according to the deformation field.
Wherein, the step of using the number and position information of the distortion control points of the magnetic resonance image in the magnetic resonance water model image in the step two comprises the following steps:
step A1: performing convolution processing on the magnetic resonance water model image by adopting a cross template, wherein the convolution processing is used for highlighting signals at the centers of the grid cross points, so that the accurate position of the control point of the magnetic resonance image can be conveniently determined in the subsequent steps;
step A2: setting a threshold value to filter the background of the water model image, and setting the initial value of a control point of the background smaller than the threshold value as 0;
step A3: recording the number and position information of control points according to the maximum brightness value in a rectangular window range with a certain size in the magnetic resonance water model image as the control points;
step A4: if the number of the magnetic resonance image control points is larger than the number of the grid points, increasing the threshold value and re-executing the step A2; if the number of control points is less than the number of grid points, decreasing the threshold and re-executing step A2; and if the number of the control points is equal to the number of the grid points, finishing extracting the number and the position information of the control points.
In the third step, the deformation quantity of each control point within the range of the control point of the magnetic resonance water model image is estimated by adopting a bilinear interpolation algorithm to estimate the transverse and longitudinal variables of each control point, so that a deformation field is formed; and fitting and estimating the distortion degree of each control point by adopting a least square method for each point outside the range of the transverse and longitudinal control points, and finally forming a deformation field containing distortion values.
The step of correcting the magnetic resonance image in the fourth step is a process of calculating the intensity of each pixel point on the final image (i.e., the corrected image) by using the two-dimensional geometric correction algorithm of the magnetic resonance image with geometric distortion and combining an automatic window adjustment algorithm, and specifically includes:
step C1: for a control point on the final image, calculating the position of the control point in the magnetic resonance image corresponding to the geometric deformation according to the deformation field;
step C2: assigning the last control point of the final magnetic resonance water model image as a specific value of the deformation position in the magnetic resonance image with the geometric deformation;
step C3: and repeating the steps C1 to C2 until all control points on the final image are processed and then the process is terminated.
And when the real position of the control point exceeds the magnetic resonance image, the variable value of the pixel point is assigned to be 0.
When the true positions of the transverse control point and the longitudinal control point are within the image range and exceed the deformation field and the distortion description range, at least two distortion data points are selected in the adjacent area of the pixel point, and the deformation position of the pixel point is obtained by adopting least square normal fitting. This is the second innovation of the present invention.
Wherein the two-dimensional geometric correction algorithm is expressed by the following formula of a least square method:
f(x,y)=bx+y+a
Figure BDA0002427546360000041
Figure BDA0002427546360000042
in the formula, x represents the abscissa of the deformation position, y represents the ordinate of the deformation position, and xi, yi represent the fitted numerical values.
The invention adopts the iterative search algorithm of the local maximum value on the extraction of the control point of the magnetic resonance water model image, and overcomes the defect that the conventional algorithm is difficult to extract the control point with larger deformation; according to the method, the deformation field is estimated for the region outside the control point range by adopting the least square normal fitting algorithm, so that the defect of region condition limitation in the conventional algorithm is overcome effectively, and a good correction effect is obtained.
Drawings
FIG. 1 is a flow chart of the two-dimensional geometry correction algorithm of the present invention.
Fig. 2 is a schematic view of a positioning plate used in the water model data acquisition process.
FIG. 3 is a flowchart illustrating the processing of the water model image according to the present embodiment.
Detailed Description
The invention is further described with reference to the following figures and examples.
The present invention will be described in further detail with reference to the following specific examples and the accompanying drawings. The procedures, conditions, experimental methods and the like for carrying out the present invention are general knowledge and common general knowledge in the art except for the contents specifically mentioned below, and the present invention is not particularly limited.
The algorithm flow of the two-dimensional geometric correction of the magnetic resonance image of the present invention, as shown in fig. 1, includes:
the method comprises the following steps: and forming a water model according to the magnetic resonance image, and recording the number and the position information of grid points of the magnetic resonance water model image as the number and the position information of the initial magnetic resonance image simulation control points, wherein the grid points are selected as the control points.
Step two: acquiring a plurality of magnetic resonance water model images according to different magnetic resonance imaging areas by using the magnetic resonance water model images, and extracting the number and position information of control points in each magnetic resonance water model image; and the control points of the magnetic resonance water model image correspond to grid points of each water model.
Step three: and combining a plurality of magnetic resonance water model images according to different areas of magnetic resonance imaging, calculating the deformation quantity between the position information of the grid points and the position information of the control points, namely the distortion effect, and estimating the deformation field according to the deformation quantity.
Step four: and correcting the magnetic resonance water model image according to the deformation field.
In this embodiment, the internal size of the magnetic resonance water-mode image designed based on the magnetic resonance image is 190mm × 190mm, the external size is 210mm × 210mm, the water-mode is filled with the aqueous solution to form 18 rows and 18 columns, the grid points are selected as control points, and the number of the grid points is 324 after the magnetic resonance water-mode images are spliced together.
In this embodiment, it is preferable that the direction of the main magnetic field is a depth direction, the left-right direction of the human body is an x-axis direction, and the front-back direction of the human body is a y-axis direction, as shown in fig. 2. The experiment is divided into four times, and the water model is placed at four positions, namely, the upper left position, the lower left position, the upper right position and the lower right position of the position O where the cross point of the bed surface of the instrument and the central positioning line of the magnet is positioned, and the water model can be spliced together without gaps. The four magnetic resonance water model image experiments are completed on a Siemens 3.0T magnetic resonance scanner by adopting the same imaging protocol. Four water model images with y being 0 layers are respectively selected from four experimental data. 324 x 4 control points of the magnetic resonance water model image cover the range of 380mm x 380mm on the plane with y being 0 and centered at the center of the magnet. The acquired water model image is then processed by the following steps, including:
step A1: and (3) convolving the magnetic resonance water mode image by adopting the normalized mean template of the cross-shaped templates of 3 rows and 9 columns and 9 rows and 3 columns to emphasize the variable values on the grid points.
Step A2: and setting a threshold T to filter the background of the magnetic resonance water model image, and setting the variable value of the background smaller than the threshold T as 0.
Step A3: and recording the number and position information of the control points by combining the maximum brightness value in a rectangular window range with a certain size on the magnetic resonance water model image as the control points. The rectangular range is a rectangular range of the automatic window adjusting window W × Wmm, and the interval between the control points in this embodiment is theoretically 10mm × 10mm, so that W is set to 9. Therefore, only one control point exists in each rectangular range, and the highest brightness value in the rectangular range is obtained as the control point. And sequencing the control points from top to bottom and from left to right, and storing the positions of the transverse control points and the longitudinal control points.
Step A4: if the number of the control points of the magnetic resonance water model image is larger than that of the grid points, increasing the threshold value T and executing the step A2 again; if the number of control points is less than the number of grid points, decreasing the threshold and re-executing step A2; and if the number of the control points of the magnetic resonance image is equal to the number of the control points of the magnetic resonance image, finishing extracting the number and the position information of the control points.
The step of estimating the deformation field in the present invention comprises the step B: estimating a deformation field by combining deformation quantities of all points within the coverage range of the control points by adopting a double-line interpolation method; and fitting and estimating the deformation field by adopting a least square method for each point outside the control point range. The control point range is an area surrounded by the outermost control points.
The step of correcting the magnetic resonance image in the fourth step of the present invention is a process of calculating each control point on the final image (i.e., the corrected image) using the magnetic resonance image with geometric deformation, and specifically includes:
step C1: calculating the position of a pixel point on the final image, which corresponds to the magnetic resonance image with geometric deformation, according to the deformation field;
step C2: assigning the variable value of the pixel point on the final magnetic resonance image to the variable value of the deformation position in the magnetic resonance image with geometric deformation;
step C3: and C1 to C2 are repeatedly executed until all the pixels on the final image are processed and then the process is terminated.
And when the deformation positions of the transverse control points and the longitudinal control points exceed the magnetic resonance water model image, the variable value of the control points is assigned to be 0.
When the real positions of the transverse control points and the longitudinal control points exceed the range of the deformation field, data points of at least two deformation fields are selected in the adjacent areas of different control points. In this embodiment, 10 data points of the deformation field are selected, and the deformation position of the control point is obtained by least square normal linear fitting. And correcting the points in the control point range by adopting an algorithm of automatic two-dimensional geometric correction of the magnetic resonance image, and correcting the deformation equivalent extension of the control point boundary by adopting the points outside the control point range.
And the result obtained by the steps is the result corrected by the two-dimensional geometric correction square algorithm of the magnetic resonance image. Compared with the acquired corrected result image, the acquired corrected result image has better correction effect at a position far away from the center of the magnet, the geometric deformation can cause variable unevenness, the gray correction, the image enhancement and other operations need to be carried out on the image, the variable correction is introduced into the geometric correction by the current geometric correction algorithm, and the invention aims to improve the geometric correction algorithm.
Example (b):
in this embodiment, the magnetic resonance water model image is placed at the upper left position of the positioning plate shown in fig. 2, the outer contour of the water model is attached to the cross mark position line, the scanning sequence adopted is tr _ mpr _ tra _ p2_ iso, the voxel is 1.0mm × 1.0mm, and the image size is: 256mm by 256mm, image center shifted from magnet center r101.7mm, p18.8mm, and h101.0mm, and the number of collected layers was 144, resulting in Data 1.
The magnetic resonance water model image is placed at the upper right position of the positioning plate shown in fig. 2, the outer contour of the magnetic resonance water model image is attached to the cross mark bit line, the adopted scanning sequence is tr _ mpr _ tra _ p2_ iso, the voxel is 1.0mm, the image size is 256mm, the image center is deviated from the magnet center by L120.7.7 mm, P18.8mm and H109.0mm, the number of collected layers is 144, and the Data set Data2 is obtained.
The magnetic resonance water model image is placed at the lower right position of the positioning plate shown in fig. 2, the outer contour of the magnetic resonance water model image is attached to the cross mark bit line, the adopted scanning sequence is tr _ mpr _ tra _ p2_ iso, the voxel is 1.0mm, the image size is 256mm, the image center is deviated from the magnet center by L120.7.7 mm, P18.8mm and F118.0mm, and the number of collected layers is 144, so that a Data set Data3 is obtained.
The magnetic resonance water model image is placed at the lower left position of the positioning plate image shown in fig. 2, the outer contour of the magnetic resonance water model image is attached to the cross mark position line, the adopted scanning sequence is tr _ mpr _ tra _ p2_ iso, the voxel is 1.0mm, and the image size is as follows: 256mm, image center shifted from the magnet center r123.7mm, a12.3mm, and f123.7mm, and the number of collected layers was 144, yielding Data set Data 4.
And selecting a Data layer where y is 0 in the Data1 according to the parameter information Slice location of the image, extracting 324 control points according to the second step of the method, and obtaining the positions P2 after the deformation of the 324 control points according to the pixel spacing and the image orientation of the parameter image. The exact position of the control point P1 is measurable with the magnet center as the origin. The deformation amount of the upper left 324 control points can be calculated according to the formula D-P1-P2.
The deformation quantities of the four acquired magnetic resonance water model images are stored in a matrix map, and the storage mode is as follows: assuming that the interval between two adjacent control points in two same images is 9 data points, the distance between two adjacent control points in two same images is 1mm, the two adjacent images are the upper left and the upper right, the interval between the rightmost control point extracted by the upper left water model image and the leftmost control point extracted by the upper right magnetic resonance water model image is 40mm, and the interval is 39 data points. The values of the data points between the control points are filled in with neighboring data points by means of bilinear interpolation. The formula is shown below
Figure BDA0002427546360000081
Wherein, Q11(x1, y1), Q21(x2, y1), Q12(x1, y2) and Q22(x2, y2) are four corner points of a rectangular frame, I (x, y) refers to the variable value of the position (x, y), and x and y are the horizontal and vertical coordinates of the internal points of the rectangle. For example, I (Q11(1,2)) -1, I (Q21(2,2)) -2, I (Q12(1, 1)) -3, I (Q22(2, 1)) -4, can be obtained by calculation using the formula
Figure BDA0002427546360000082
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.

Claims (7)

1. A two-dimensional geometric correction method for magnetic resonance images is characterized by comprising the following steps:
step 1: acquiring experimental data and recording the number and position information of grid points of a water model by using a magnetic resonance imaging image and the water model as the number and position information of initial transverse and longitudinal control points of the image;
step 2: acquiring a plurality of water model images by utilizing different areas of the water model effect of the magnetic resonance image, and extracting the number and the position information of control points in each image, wherein the position information of the transverse control points and the longitudinal control points is the position of the control points in the image after the image is deformed; the control points of the image correspond to the stretching grid points of the horizontal direction and the longitudinal direction of the water model image;
and step 3: combining the plurality of images according to different regions of the magnetic resonance image, calculating position information and deformation quantities of the transverse and longitudinal stretching and shrinking grid points, namely the control points, and estimating a deformation field according to the deformation quantities;
and 4, step 4: and correcting the magnetic resonance image with geometric deformation by using an improved automatic window adjusting algorithm according to the deformation field.
2. The method as claimed in claim 1, wherein the water model has a square shape, the grid of the water model image is a grid-like glass filler, the filler in the water model image is water or solution, the inner dimension of the water model has a side length of 190mm, the outer dimension has a side length of 210mm, the outer wall has a thickness of 10mm, and the number of the grid points of the magnetic resonance water model image is 81.
3. The method of claim 1, wherein the step of extracting the number and position information of the control points in each image comprises:
step A1: performing water modeling image convolution processing on the magnetic resonance image by adopting a cross template;
step A2: setting a threshold value to filter out the background of the magnetic resonance water model image;
step A3: traversing pixel points with the maximum brightness value in a rectangular window range with a certain size in the magnetic resonance water model image as the transverse and longitudinal initial control points, and recording the quantity and position information of the control points;
step A4: if the number of the transverse and longitudinal control points is larger than the number of the grid points, increasing the threshold value and executing the step A2 again; if the number of control points is less than the number of grid points, decreasing the threshold and re-executing the step a 2; and if the number of the control points is equal to the number of the grid points, finishing the extraction of the number and the position information of the transverse and longitudinal control points.
4. The method of claim 1, wherein the step of estimating the deformation field comprises: estimating the deformation quantity of each transverse and longitudinal control point by adopting an automatic window adjusting algorithm for the deformation quantity of each control point within the range of the transverse and longitudinal control points, calculating the stretching and shrinking quantities and forming a deformation field; and fitting and estimating the deformation quantity of each control point by adopting a least square method for each point outside the range of the control points, and calculating the stretching and shrinking quantities to form a deformation field.
5. The method of claim 1, wherein the step of modifying the geometrically deformed magnetic resonance image comprises:
step C1: calculating the position of a transverse and longitudinal control point on the final image, which corresponds to the magnetic resonance image with the geometric deformation, according to the deformation field;
step C2: assigning the stretching and shrinking variable quantity of the control point on the final magnetic resonance image as the variable quantity of the control point of the deformation position in the magnetic resonance image with the geometric deformation;
step C3: and repeating the steps C1 to C2 until all the image control points on the final image are processed and then the process is terminated.
6. The method of claim 5, wherein when the deformation position of the pixel point exceeds the magnetic resonance image, the initial value of the image control point is 0; and when the real position of each control point is in the image range and exceeds the deformation field description range, selecting at least two deformation field data points in the adjacent area of the control point, and obtaining the deformation positions of the transverse and longitudinal control points by adopting least square normal fitting.
7. The method of claim 6, wherein the least squares method is formulated as follows: f (x, y) ═ bx + y + a
Figure FDA0002427546350000021
Figure FDA0002427546350000031
In the formula, x represents the abscissa of the deformation position, y represents the ordinate of the deformation position, and xi, yi represent the fitted numerical values.
CN202010225659.2A 2020-03-26 2020-03-26 Two-dimensional geometric correction method for magnetic resonance image Pending CN111445412A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010225659.2A CN111445412A (en) 2020-03-26 2020-03-26 Two-dimensional geometric correction method for magnetic resonance image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010225659.2A CN111445412A (en) 2020-03-26 2020-03-26 Two-dimensional geometric correction method for magnetic resonance image

Publications (1)

Publication Number Publication Date
CN111445412A true CN111445412A (en) 2020-07-24

Family

ID=71648236

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010225659.2A Pending CN111445412A (en) 2020-03-26 2020-03-26 Two-dimensional geometric correction method for magnetic resonance image

Country Status (1)

Country Link
CN (1) CN111445412A (en)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070018645A1 (en) * 2003-05-29 2007-01-25 Deming Wang Method and apparatus for mapping and correcting geometric distortion in mri
CN101046506A (en) * 2007-03-28 2007-10-03 新奥博为技术有限公司 Method of correcting gradient deformation of magnetic resonant image based on spherical harmonic function
CN101216541A (en) * 2008-01-15 2008-07-09 新奥博为技术有限公司 Magnetic resonance image-forming system gradient field spherical harmonic coefficient obtaining method
CN102283651A (en) * 2011-07-11 2011-12-21 苏州安科医疗***有限公司 Magnetic resonance imaging gradient field correction method
CN102521809A (en) * 2011-12-08 2012-06-27 沈阳工业大学 Regularization correction method of magnetic resonance phase array coil image uniformity
CN103020944A (en) * 2013-01-08 2013-04-03 华东师范大学 Method for two-dimensional geometric correction of magnetic resonance image
US20140177933A1 (en) * 2012-12-21 2014-06-26 Industry-University Cooperation Foundation Hanyang University Method and apparatus for mapping image information of target object onto interface of the target object in magnetic resonance imaging
WO2016084964A1 (en) * 2014-11-27 2016-06-02 学校法人岩手医科大学 Phantom for magnetic resonance imaging device
CN108334733A (en) * 2017-09-01 2018-07-27 上海联影医疗科技有限公司 Medical image display method, display system and computer readable storage medium
US20200020082A1 (en) * 2018-07-10 2020-01-16 The Board Of Trustees Of The Leland Stanford Junior University Un-supervised convolutional neural network for distortion map estimation and correction in MRI
US10557911B2 (en) * 2014-04-24 2020-02-11 David W. Holdsworth Method and apparatus for measuring 3D geometric distortion in MRI and CT images with a 3D physical phantom
CN110807492A (en) * 2019-11-06 2020-02-18 厦门大学 Magnetic resonance multi-parameter simultaneous quantitative imaging method and system

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070018645A1 (en) * 2003-05-29 2007-01-25 Deming Wang Method and apparatus for mapping and correcting geometric distortion in mri
CN101046506A (en) * 2007-03-28 2007-10-03 新奥博为技术有限公司 Method of correcting gradient deformation of magnetic resonant image based on spherical harmonic function
CN101216541A (en) * 2008-01-15 2008-07-09 新奥博为技术有限公司 Magnetic resonance image-forming system gradient field spherical harmonic coefficient obtaining method
CN102283651A (en) * 2011-07-11 2011-12-21 苏州安科医疗***有限公司 Magnetic resonance imaging gradient field correction method
CN102521809A (en) * 2011-12-08 2012-06-27 沈阳工业大学 Regularization correction method of magnetic resonance phase array coil image uniformity
US20140177933A1 (en) * 2012-12-21 2014-06-26 Industry-University Cooperation Foundation Hanyang University Method and apparatus for mapping image information of target object onto interface of the target object in magnetic resonance imaging
CN103020944A (en) * 2013-01-08 2013-04-03 华东师范大学 Method for two-dimensional geometric correction of magnetic resonance image
US10557911B2 (en) * 2014-04-24 2020-02-11 David W. Holdsworth Method and apparatus for measuring 3D geometric distortion in MRI and CT images with a 3D physical phantom
WO2016084964A1 (en) * 2014-11-27 2016-06-02 学校法人岩手医科大学 Phantom for magnetic resonance imaging device
CN108334733A (en) * 2017-09-01 2018-07-27 上海联影医疗科技有限公司 Medical image display method, display system and computer readable storage medium
US20200020082A1 (en) * 2018-07-10 2020-01-16 The Board Of Trustees Of The Leland Stanford Junior University Un-supervised convolutional neural network for distortion map estimation and correction in MRI
CN110807492A (en) * 2019-11-06 2020-02-18 厦门大学 Magnetic resonance multi-parameter simultaneous quantitative imaging method and system

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
DEMING WANG: "A novel phantom and method for comprehensive 3-dimensional measurement and correction of geometric distortion in magnetic resonance imaging", 《MAGNETIC RESONANCE IMAGING》 *
吕磊 等: "一种16位灰度图像自动调窗算法", 《光电技术应用》 *
康晓东: "《现代医学影像技术》", 31 January 2000, 天津:天津科技翻译出版公司 *
张彬: "最小二乘法的未知控制点检测", 《中国组织工程研究》 *
沈琴 等: "一种改进的磁共振图像自动调窗算法", 《中国医疗器械杂志》 *
陈华莉: "基于背景信息的磁共振成像尖峰噪声消除", 《波谱学杂志》 *

Similar Documents

Publication Publication Date Title
CN109146930B (en) Infrared and visible light image registration method for electric power machine room equipment
US10224955B2 (en) Data compression and decompression method of demura table, and mura compensation method
CN102790841B (en) Method of detecting and correcting digital images of books in the book spine area
CN102883175B (en) Methods for extracting depth map, judging video scene change and optimizing edge of depth map
CN106960449B (en) Heterogeneous registration method based on multi-feature constraint
CN106611416B (en) Method and device for segmenting lung in medical image
CN105005981B (en) Based on multiple dimensioned Light stripes center extraction method and device
CN102156972B (en) Image tilting correcting method and system
CN113284189B (en) Distortion parameter calibration method, device, equipment and storage medium
CN108895959A (en) A kind of camera calibration plate angle point calculating method based on sub-pix
CN110414385A (en) A kind of method for detecting lane lines and system based on homography conversion and characteristic window
CN112087970A (en) Information processing apparatus, information processing method, and computer program
JP6830712B1 (en) Random sampling Consistency-based effective area extraction method for fisheye images
CN111489383B (en) Depth image up-sampling method and system based on depth marginal point and color image
CN114092468A (en) Standard target counting method based on machine vision
CN117576219A (en) Camera calibration equipment and calibration method for single shot image of large wide-angle fish-eye lens
CN103020944A (en) Method for two-dimensional geometric correction of magnetic resonance image
CN112037284A (en) Checkerboard grid center point detection method, camera distortion calibration method and system
CN113723399A (en) License plate image correction method, license plate image correction device and storage medium
CN110223356A (en) A kind of monocular camera full automatic calibration method based on energy growth
CN114612412A (en) Processing method of three-dimensional point cloud data, application of processing method, electronic device and storage medium
CN108564622A (en) The method for realizing sub-pixel Corner character in positioning plate image
CN111445412A (en) Two-dimensional geometric correction method for magnetic resonance image
CN112285876A (en) Camera automatic focusing method based on image processing and bubble detection
CN116012322A (en) Camera dirt detection method, device, equipment and medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20200724

RJ01 Rejection of invention patent application after publication