CN101422364A - Magnetic resonance image seamless splicing method - Google Patents

Magnetic resonance image seamless splicing method Download PDF

Info

Publication number
CN101422364A
CN101422364A CNA2007100477523A CN200710047752A CN101422364A CN 101422364 A CN101422364 A CN 101422364A CN A2007100477523 A CNA2007100477523 A CN A2007100477523A CN 200710047752 A CN200710047752 A CN 200710047752A CN 101422364 A CN101422364 A CN 101422364A
Authority
CN
China
Prior art keywords
image
simg
factor
pixel
signal intensities
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
Application number
CNA2007100477523A
Other languages
Chinese (zh)
Other versions
CN101422364B (en
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.)
KALEFU MAGNETIC RESONANCE TECH Co Ltd SHANGHAI
East China Normal University
Original Assignee
KALEFU MAGNETIC RESONANCE TECH Co Ltd SHANGHAI
East China Normal 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 KALEFU MAGNETIC RESONANCE TECH Co Ltd SHANGHAI, East China Normal University filed Critical KALEFU MAGNETIC RESONANCE TECH Co Ltd SHANGHAI
Priority to CN200710047752A priority Critical patent/CN101422364B/en
Publication of CN101422364A publication Critical patent/CN101422364A/en
Application granted granted Critical
Publication of CN101422364B publication Critical patent/CN101422364B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)

Abstract

The invention discloses a method for assembling magnetic resonance images in a seamless manner, pertaining to the technical field of magnetic resonance imaging and seamless assembling of medical images. The method is required to be implemented in a magnetic resonance imaging system and an image processing section of the system comprises a central processing unit and a memorizer. Programs for instructing and commanding the central processing unit to carry out image processing operations are stored in the memorizer, and the programs comprise two procedures of data acquisition and image processing: when magnetic resonance image data of all sections are collected, a marker is arranged beside a patient and the image of the marker is utilized for determining the seam location of the adjacent magnetic resonance images; before the magnetic resonance images are assembled, image data at both sides of the seam are processed for realizing the seamless assembling of the adjacent magnetic resonance images. The method has the advantages of simple and convenient operation, easy and precise image registration as well as fewer artifacts of image assembling.

Description

The method of magnetic resonance image seamless splicing
Technical field
The present invention relates to a kind of method of magnetic resonance image seamless splicing, belong to nuclear magnetic resonance and the seamless spliced technical field of medical image.
Technical background
In general, the each imaging region of nuclear magnetic resonance (observing wild) can be subjected to certain limitation, and the main field uniformity of imaging region, the linearity of gradient magnetic, the factors such as coverage of radio-frequency coil are depended in these restrictions.Sometimes clinical examination needs imaging region to cover very large scope, need observe whole vertebra when checking metastatic tumor to the infringement of vertebra.Solution commonly used at present is: scan a position of patient body earlier, obtain the image at this position; The motion scan bed scans another position of patient body more then, obtains the image at this position; So analogize,, more all image mosaics are become a complete image up to the image that collects needed body part.The shortcoming of background technology is, the data of this image are not handled before splicing above-mentioned image, the factors such as spatially uniform of anamorphose and view data make the width of cloth spliced image flaw occur because of the image of seam crossing is difficult to registration, the discontinuous and pseudo-shadow of splicing as image.
Summary of the invention
At the deficiency of above-mentioned background technology, the purpose of this invention is to provide a kind of method of magnetic resonance image seamless splicing.
Technical scheme of the present invention is, above-mentioned method comprises data acquisition and two processes of Flame Image Process: when gathering the magnetic resonance image data of all sites, at patient's the label of placement by the side of, utilize the picture of label to determine the seaming position of above-mentioned adjacent magnetic resonance images; Before the above-mentioned magnetic resonance image (MRI) of splicing, the view data of above-mentioned seam both sides is handled, realize the seamless spliced of above-mentioned adjacent magnetic resonance images.
Now be described with reference to the accompanying drawings technical scheme of the present invention.
A kind of method of magnetic resonance image seamless splicing, need in magnetic resonance imaging system, to implement, the image processing section of above-mentioned system contains central processing unit and memorizer, store memory has the program of order and the operation of commander's central processing unit carries out image processing, it is characterized in that, comprise two processes of data acquisition and Flame Image Process, the concrete operations of data acquisition are:
The 1st the step patient lie in scanning bed on, on the first adjacent imaging position and the second imaging position respectively coated with first coil and second coil;
The 2nd step was placed two labels between first coil and second coil: imaging sagittal plane or imaging coronalplane, and the former places forward and backward of patient, and the latter places patient's left and right limit;
The 3rd moved further is scanning bed, and the first imaging position is pushed to magnet center;
The 4th step was scanned the first imaging position, and scanner field covers the label between first coil and second coil, and scanning obtains first image 1;
The 5th moved further is scanning bed, and the second imaging position that will imaging pushes to magnet center;
The 6th step was scanned the second imaging position, and scanner field covers the label between first coil and second coil, and scanning obtains second image 2;
The 7th step, scanning obtained the image at more imagings position by that analogy;
The concrete operations of image processing process are:
The 8th step was determined the center of labelling object image 4 on first image 1 and second image 2 respectively; The principle of determining is, 1, the ratio of the image intensity signal meansigma methods of the image intensity signal meansigma methods of label image 4 and label image 4 peripheries must be greater than a certain threshold level, this threshold value is the value that the twice of the image intensity signal meansigma methods of label image 4 obtains divided by the outer image intensity signal meansigma methodss of label image 4,2, the pericentral picture signal symmetry of label image 4 is the highest;
The 9th step determined that the line at the center of the line at center of two label images 4 on 5: the first images 1 in the seam 3 of first image 1 and second image 2 and overlapping region and two label images 4 on second image 2 was exactly a seam 3, and the zone between two lines of 3 liang of back gauge seam 3 certain distances of seam is defined as overlapping region 5;
The 10th goes on foot the distance between the center of measuring labelling object image 4 on first image 1 and second image 2 respectively;
The mode translation that 2: the second images 2 of the 11st step horizontal direction translation second image overlap with the seam 3 of first image 1 with its seam 3, the mid point of the line at the center of two label images 4 on the mid point of the line at the center of two label images 4 on second image 2 and first image 1 overlaps;
The 12nd step, the topography to overlapping region 5 was out of shape correction: adopt the gradual change principle, if the distance between the center of two label images 4 on first image 1 is greater than the distance between the center of two label images 4 on second image 2, then first topography of image 1 in overlapping region 5 diminishes from top to bottom gradually, it is big that second topography of image 2 in overlapping region 5 becomes from top to bottom gradually, if the distance between the center of two label images 4 on first image 1 is less than the distance between the center of two label images 4 on second image 2, then first topography of image 1 in overlapping region 5 becomes big from top to bottom gradually, second topography of image 2 in overlapping region 5 diminishes from top to bottom gradually, after distortion was proofreaied and correct, two label images 4 on first image 1 just in time overlapped with two label images 4 on second image 2 respectively;
The 13rd step was proofreaied and correct all pixel signal intensities on second image 2: measure on first image 1 and second image 2 the pixel signal intensities meansigma methods SAvg_1 and the SAvg_2 of 5 topography in the overlapping region respectively, pixel signal intensities with each pixel of second topography of image 2 in overlapping region 5 all multiply by SAvg_1/SAvg_2 then, and the pixel signal intensities of second image, 2 integral body and first image 1 pixel signal intensities are on the whole coincide;
Level and smooth splicing of the 14th step: vertically respectively the pixel signal intensities of all pixels of first image 1 handled through the 13rd step and second topography of image 2 in overlapping region 5 is handled with the method for fading in-fade out before the splicing, the pixel signal intensities of all pixels of first topography of image 1 in overlapping region 5 is carried out vertically the processing of fading out, the pixel signal intensities of all pixels of second topography of image 2 in overlapping region 5 is carried out vertically the processing of fading in, if gradual change factor factor, the value of factor is between 0~1, the former pixel signal intensity of the respective pixel Img_2 of first image 1 the pixel Img_1 and second image 2 is respectively SImg_1 and SImg_2 in the topography of overlapping region 5, and its pixel signal intensities became factor after pixel Img_1 handled through fading out *Its pixel signal intensities became SImg_2 after SImg_1, pixel Img_2 handled through fading in *(1-factor), pixel Img_1 is SImg_st, then SImg_st=SImg_1 with pixel signal intensities after pixel Img_2 overlaps *Factor+SImg_2 *(1-factor)=SImg_1 *Factor-SImg_2 *Factor+SImg_2, wherein factor is relevant with the vertical direction overlap distance between the image, change to 0 by 1, the pixel signal intensities of the pixel outside the overlapping region 5 still adopts the pixel signal intensities of the former pixel on first image 1 and second image 2 respectively, will smoothly be stitched together through first image 1 and second image of handling 2 that fade in-fade out;
The 15th step was wiped the label on the stitching image: spliced the noise signal strength replacement of the pixel signal intensities of labelling object image 4 on the image that obtains with label image 4 surrounding pixels the 14th step, obtain seamless spliced magnetic resonance image (MRI).
Technical scheme of the present invention is further characterized in that in the 2nd step, described label is the test tube that fills with copper-bath, and the internal diameter of the test tube of copper-bath is 4.25 millimeters.
Technical scheme of the present invention is further characterized in that in the 9th step, described overlapping region 5 is evenly distributed on the both sides of seam 3, and the width of overlapping region 5 on seam 3 each limit is 10 pixels.
Technical scheme of the present invention is further characterized in that, in the 14th step, the pixel signal intensities of all pixels of first topography of image 1 in overlapping region 5 is carried out vertically the processing of fading in, the pixel signal intensities of all pixels of second topography of image 2 in overlapping region 5 is carried out vertically the processing of fading out, and its pixel signal intensities became SImg_1 after pixel Img_1 handled through fading in *(1-factor), its pixel signal intensities became factor after pixel Img_2 handled through fading out *SImg_2, pixel Img_1 is SImg_st=SImg_1 with pixel signal intensities after pixel Img_2 overlaps *(1-factor)+factor *SImg_2=SImg_1-factor *SImg_1+factor *SImg_2.
The invention has the advantages that:
1, simple, convenient.
2, image registration easily, accurately.
3, the pseudo-shadow of image mosaic is few.
Description of drawings
Fig. 1 is the seamless spliced sketch map of image.Wherein, 1 is first image, and 2 is second images, the 3rd, and seam, the 4th, label image, the 5th, overlapping region, the i.e. lap of first image 1 and second image 2.
Fig. 2 is the sketch map of geometry deformation of proofreading and correct TMR's image of overlapping region 5.Among the A figure, the 5th, overlapping region, TMR's image of first image 1 before promptly geometry deformation is proofreaied and correct and the lap of second image 2, the 51st, TMR's image of first image 1 and second image, 2 laps in first image 1 of geometry deformation correction back.Among the B figure, the 52nd, TMR's image of second image 2 and first image, 1 lap in second image 2 of geometry deformation correction back.
Fig. 3 is the magnetic resonance image (MRI) that the first imaging position collects among the embodiment, i.e. first image 1, and two white points among the figure are label images 4.
Fig. 4 is the magnetic resonance image (MRI) that the second imaging position collects among the embodiment, i.e. second image 2, and two white points among the figure are label images 4.
Fig. 5 is first image 1 and second image 2 magnetic resonance image (MRI) after seamless spliced among the embodiment.
The specific embodiment
Now further specify technical scheme of the present invention by drawings and Examples.Embodiment operates according to two processes of above-mentioned image seamless joint method fully, below only enumerates the key technology data of each process and step.
The method of a kind of magnetic resonance image seamless splicing of embodiment: the coronalplane magnetic resonance image (MRI) of head and neck and chest is together seamless spliced.
The concrete operations of data acquisition are:
In the 1st step, patient head and last cervical region cover with first coil, and cervical region and chest cover with second coil under the patient; In the 2nd step, place two labels at first coil and the second coil intersection, label is the test tube that fills with copper-bath, and the internal diameter of test tube is 4.25 millimeters, present embodiment is the imaging coronalplane, and two labels are symmetrically placed in the left side and the right of patient's neck; In the 3rd step, patient head is pushed to magnet center; In the 4th step, the first coil scanning head and last cervical region, scanner field covers above-mentioned label, and in the present embodiment, scanner field is 360 millimeters, and the image that scanning obtains is designated as first image 1, sees Fig. 3; In the 5th step, patient's chest is pushed to magnet center; In the 6th step, the scanning of second coil is cervical region and chest down, and scanner field covers above-mentioned label, and in the present embodiment, scanner field is 360 millimeters, and the image that scanning obtains is designated as second image 2, sees Fig. 4;
The concrete operations of image processing process are:
In the 8th step, determine test tube, the i.e. center of label image 4 on first image 1 and second image 2 respectively; In the 9th step, the line at the center of small test tube image is a seam 3 on first image 1 and second image 2, and two lines on seam 3 both sides and the distance of seam 3 are 10 pixels, and promptly the width of overlapping region 5 is 20 pixels; In the 10th step, the distance that records between the center of labelling object image 4 on first image 1 and second image 2 is respectively 162 and 165 pixels, is 227.8 millimeters and 232.0 millimeters; In the 11st step, horizontal direction translation second image 2, second image, 2 integral levels to right translation 4 pixels, be 360*4/256 millimeter=5.625 millimeter, after the translation, the mid point of the line at the center of two label images 4 on the mid point of the line at the center of two label images 4 on second image 2 and first image 1 overlaps; In the 12nd step, topography to overlapping region 5 is out of shape correction, in the present embodiment, because the distance between the center of two label images 4 on first image 1 is 162 pixels (227.8 millimeters), less than 165 pixels of the distance between the center of two label images 4 on second image 2 (232.0 millimeters), so it is big that first topography of image 1 in overlapping region 5 becomes from top to bottom gradually, second topography of image 2 in overlapping region 5 diminishes from top to bottom gradually, after distortion was proofreaied and correct, two label images 4 on first image 1 just in time overlapped with two label images 4 on second image 2 respectively; In the 13rd step, proofread and correct all pixel signal intensities on second image 2, measure on first image 1 and second image 2 the pixel signal intensities meansigma methods SAvg_1=195 and the SAvg_2=70 of 5 topography respectively in the overlapping region, pixel signal intensities with each pixel of second topography of image 2 in overlapping region 5 all multiply by SAvg_1/SAvg_2=2.79 then, and the pixel signal intensities of second image, 2 integral body and first image 1 are coincide in the pixel signal intensities of integral body; In the 14th step, level and smooth splicing, vertically respectively the pixel signal intensities of all pixels of first image 1 handled through the 13rd step and second topography of image 2 in overlapping region 5 is handled with the method for fading in-fade out before the splicing, the pixel signal intensities of all pixels of first topography of image 1 in overlapping region 5 is carried out vertically the processing of fading out, the pixel signal intensities of all pixels of second topography of image 2 in overlapping region 5 is carried out vertically the processing of fading in, if gradual change factor factor, the value of factor is between 0~1, the former pixel signal intensity of the respective pixel Img_2 of first image 1 the pixel Img_1 and second image 2 is respectively SImg_1 and SImg_2 in the topography of overlapping region 5, and its pixel signal intensities became factor after pixel Img_1 handled through fading out *Its pixel signal intensities became SImg_2 after SImg_1, pixel Img_2 handled through fading in *(1-factor), pixel Img_1 is SImg_st, then SImg_st=SImg_1 with pixel signal intensities after pixel Img_2 overlaps *Factor+SImg_2 *(1-factor)=SImg_1 *Factor-SImg_2 *Factor+SImg_2, wherein factor is relevant with the vertical direction overlap distance between the image, change to 0 by 1, the pixel signal intensities of the pixel outside the overlapping region 5 still adopts the pixel signal intensities of the former pixel on first image 1 and second image 2 respectively, will smoothly be stitched together through first image 1 and second image of handling 2 that fade in-fade out; In the 15th step, wipe the label on the stitching image, spliced of the noise signal strength replacement of the pixel signal intensities of labelling object image 4 on the image that obtains the 14th step, obtain seamless spliced magnetic resonance image (MRI) with label image 4 surrounding pixels.See Fig. 5.

Claims (8)

1, a kind of method of magnetic resonance image seamless splicing, need in magnetic resonance imaging system, to implement, the image processing section of above-mentioned system contains central processing unit and memorizer, store memory has the program of order and the operation of commander's central processing unit carries out image processing, it is characterized in that, comprise two processes of data acquisition and Flame Image Process, the concrete operations of data acquisition are:
The 1st the step patient lie in scanning bed on, on the first adjacent imaging position and the second imaging position respectively coated with first coil and second coil;
The 2nd step was placed two labels between first coil and second coil: imaging sagittal plane or imaging coronalplane, and the former places forward and backward of patient, and the latter places patient's left and right limit;
The 3rd moved further is scanning bed, and the first imaging position is pushed to magnet center;
The 4th step was scanned the first imaging position, and scanner field covers the label between first coil and second coil, and scanning obtains first image (1);
The 5th moved further is scanning bed, and the second imaging position that will imaging pushes to magnet center;
The 6th step was scanned the second imaging position, and scanner field covers the label between first coil and second coil, and scanning obtains second image (2);
The 7th step, scanning obtained the image at more imagings position by that analogy;
The concrete operations of image processing process are:
The 8th step was determined the upward center of labelling object image (4) of first image (1) and second image (2) respectively: definite principle is, 1, the ratio of the image intensity signal meansigma methods of the image intensity signal meansigma methods of label image (4) and label image (4) periphery must be greater than a certain threshold level, this threshold value is the value that the twice of the image intensity signal meansigma methods of label image (4) obtains divided by the outer image intensity signal meansigma methods of label image (4), 2, the pericentral picture signal symmetry of label image (4) is the highest;
The 9th step was determined seam (3) and overlapping region (5) of first image (1) and second image (2): the line at the line at the center of two the label images (4) on first image (1) and the center of two the label images (4) on second image (2) is exactly seam (3), and the zone between two lines of seam (3) two back gauge seam (3) certain distances is defined as overlapping region (5);
The 10th goes on foot the distance between the center of measuring first image (1) and the last labelling object image (4) of second image (2) respectively;
The 11st step horizontal direction translation second image (2): the mode translation that second image (2) overlaps with the seam (3) of first image (1) with its seam (3), the mid point of the line at the center of two the label images (4) on the mid point of the line at the center of two the label images (4) on second image (2) and first image (1) overlaps;
The 12nd step, the topography to overlapping region (5) was out of shape correction: adopt the gradual change principle, if the distance between the center of two the label images (4) on first image (1) is greater than the distance between the center of two the label images (4) on second image (2), then the topography of first image (1) in overlapping region (5) diminishes from top to bottom gradually, it is big that the topography of second image (2) in overlapping region (5) becomes from top to bottom gradually, if the distance between the center of two the label images (4) on first image (1) is less than the distance between the center of two the label images (4) on second image (2), then the topography of first image (1) in overlapping region (5) becomes big from top to bottom gradually, the topography of second image (2) in overlapping region (5) diminishes from top to bottom gradually, after distortion was proofreaied and correct, two the label images (4) on first image (1) just in time overlapped with two label images (4) on second image (2) respectively;
The 13rd step proofreaied and correct second image (2) and goes up all pixel signal intensities: measure pixel signal intensities meansigma methods SAvg_1 and SAvg_2 that first image (1) and second image (2) are gone up the topography of (5) in the overlapping region respectively, pixel signal intensities with each pixel of the topography of second image (2) in overlapping region (5) all multiply by SAvg_1/SAvg_2 then, makes whole pixel signal intensities of second image (2) and first image (1) pixel signal intensities on the whole identical;
Level and smooth splicing of the 14th step: vertically respectively the pixel signal intensities of all pixels of first image (1) handled through the 13rd step and the topography of second image (2) in overlapping region (5) is handled with the method for fading in-fade out before the splicing, the pixel signal intensities of all pixels of the topography of first image (1) in overlapping region (5) is carried out vertically the processing of fading out, the pixel signal intensities of all pixels of the topography of second image (2) in overlapping region (5) is carried out vertically the processing of fading in, if gradual change factor factor, the value of factor is between 0~1, the former pixel signal intensity of the respective pixel Img_2 of pixel Img_1 of first image (1) and second image (2) is respectively SImg_1 and SImg_2 in the topography of overlapping region 5, its pixel signal intensities became factor*SImg_1 after pixel Img_1 handled through fading out, its pixel signal intensities became SImg_2* (1-factor) after pixel Img_2 handled through fading in, pixel Img_1 is SImg_st with pixel signal intensities after pixel Img_2 overlaps, SImg_st=SImg_1*factor+SImg_2* (1-factor)=SImg_1*factor-SImg_2*factor+SImg_2 then, wherein factor is relevant with the vertical direction overlap distance between the image, change to 0 by 1, the pixel signal intensities of the pixel that overlapping region (5) is outer still adopts the pixel signal intensities of the former pixel on first image (1) and second image (2) respectively, will smoothly be stitched together through first image (1) and second image of handling (2) that fade in-fade out;
The 15th step was wiped the label on the stitching image: spliced the noise signal strength replacement of the pixel signal intensities of labelling object image (4) on the image that obtains with label image (4) surrounding pixel the 14th step, obtain seamless spliced magnetic resonance image (MRI).
2, the method for magnetic resonance image seamless splicing according to claim 1 is characterized in that, in the 2nd step, described label is the test tube that fills with copper-bath, and the internal diameter of copper-bath test tube is 4.25 millimeters.
3, the method for magnetic resonance image seamless splicing according to claim 1 is characterized in that, in the 9th step, overlapping region (5) are evenly distributed on the both sides of seam (3), and overlapping region (5) are 10 pixels at the width on each limit of seam (3).
4, the method for magnetic resonance image seamless splicing according to claim 1, it is characterized in that, in the 14th step, the pixel signal intensities of all pixels of the topography of first image (1) in overlapping region (5) is carried out vertically the processing of fading in, to second image
(2) pixel signal intensities of all pixels of the topography in overlapping region (5) is carried out the processing of fading out vertically, its pixel signal intensities became SImg_1* (1-factor) after pixel Img_1 handled through fading in, its pixel signal intensities became factor*SImg_2 after pixel Img_2 handled through fading out, and pixel Img_1 is SImg_st=SImg_1* (1-factor)+factor*SImg_2=SImg_1-factor*SImg_1+factor*SImg_2 with pixel signal intensities after pixel Img_2 overlaps.
5, the method for magnetic resonance image seamless splicing according to claim 1, it is characterized in that, in the 2nd step, described label is the test tube that fills with copper-bath, the internal diameter of copper-bath test tube is 4.25 millimeters, in the 9th step, overlapping region (5) are evenly distributed on the both sides of seam (3), and overlapping region (5) are 10 pixels at the width on each limit of seam (3).
6, the method of magnetic resonance image seamless splicing according to claim 1, it is characterized in that, in the 2nd step, described label is the test tube that fills with copper-bath, the internal diameter of copper-bath test tube is 4.25 millimeters, in the 14th step, the pixel signal intensities of all pixels of the topography of first image (1) in overlapping region (5) is carried out vertically the processing of fading in, the pixel signal intensities of all pixels of the topography of second image (2) in overlapping region (5) is carried out vertically the processing of fading out, its pixel signal intensities became SImg_1* (1-factor) after pixel Img_1 handled through fading in, its pixel signal intensities became factor*SImg_2 after pixel Img_2 handled through fading out, and pixel Img_1 is SImg_st=SImg_1* (1-factor)+factor*SImg_2=SImg_1-factor*SImg_1+factor*SImg_2 with pixel signal intensities after pixel Img_2 overlaps.
7, the method of magnetic resonance image seamless splicing according to claim 1, it is characterized in that, in the 9th step, overlapping region (5) is evenly distributed on the both sides of seam (3), overlapping region (5) is 10 pixels at the width on each limit of seam (3), in the 14th step, the pixel signal intensities of all pixels of the topography of first image (1) in overlapping region (5) is carried out vertically the processing of fading in, the pixel signal intensities of all pixels of the topography of second image (2) in overlapping region (5) is carried out vertically the processing of fading out, its pixel signal intensities became SImg_1* (1-factor) after pixel Img_1 handled through fading in, its pixel signal intensities became factor*SImg_2 after pixel Img_2 handled through fading out, and pixel Img_1 is SImg_st=SImg_1* (1-factor)+factor*SImg_2=SImg_1-factor*SImg_1+factor*SImg_2 with pixel signal intensities after pixel Img_2 overlaps.
8, the method of magnetic resonance image seamless splicing according to claim 1, it is characterized in that, in the 2nd step, described label is the test tube that fills with copper-bath, the internal diameter of copper-bath test tube is 4.25 millimeters, in the 9th step, overlapping region (5) is evenly distributed on the both sides of seam (3), overlapping region (5) is 10 pixels at the width on each limit of seam (3), in the 14th step, the pixel signal intensities of all pixels of the topography of first image (1) in overlapping region (5) is carried out vertically the processing of fading in, the pixel signal intensities of all pixels of the topography of second image (2) in overlapping region (5) is carried out vertically the processing of fading out, its pixel signal intensities became SImg_1* (1-factor) after pixel Img_1 handled through fading in, its pixel signal intensities became factor*SImg_2 after pixel Img_2 handled through fading out, and pixel Img_1 is SImg_st=SImg_1* (1-factor)+factor*SImg_2=SImg_1-factor*SImg_1+factor*SImg_2 with pixel signal intensities after pixel Img_2 overlaps.
CN200710047752A 2007-11-02 2007-11-02 Magnetic resonance image seamless splicing method Expired - Fee Related CN101422364B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200710047752A CN101422364B (en) 2007-11-02 2007-11-02 Magnetic resonance image seamless splicing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200710047752A CN101422364B (en) 2007-11-02 2007-11-02 Magnetic resonance image seamless splicing method

Publications (2)

Publication Number Publication Date
CN101422364A true CN101422364A (en) 2009-05-06
CN101422364B CN101422364B (en) 2010-05-19

Family

ID=40613406

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200710047752A Expired - Fee Related CN101422364B (en) 2007-11-02 2007-11-02 Magnetic resonance image seamless splicing method

Country Status (1)

Country Link
CN (1) CN101422364B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105551042A (en) * 2015-12-15 2016-05-04 沈阳东软医疗系统有限公司 Determination method and apparatus for mark point positions of scanning bed
CN105989620A (en) * 2015-03-16 2016-10-05 西门子公司 Compensating for image artifacts in medical imaging
CN111239658A (en) * 2020-02-19 2020-06-05 上海康达卡勒幅医疗科技有限公司 Method for reducing PROPELLER sequence artifacts

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1063625C (en) * 1994-08-04 2001-03-28 深圳安科高技术有限公司 Technology of reforming magnetic resonance presentation
CN100367047C (en) * 2004-11-16 2008-02-06 联想(北京)有限公司 Method and device for eliminating seam between spliced display screens

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105989620A (en) * 2015-03-16 2016-10-05 西门子公司 Compensating for image artifacts in medical imaging
CN105989620B (en) * 2015-03-16 2019-03-22 西门子公司 Compensate the image artifacts in medical imaging
CN105551042A (en) * 2015-12-15 2016-05-04 沈阳东软医疗系统有限公司 Determination method and apparatus for mark point positions of scanning bed
CN105551042B (en) * 2015-12-15 2018-09-18 沈阳东软医疗系统有限公司 A kind of scanning bed mark point method for determining position and device
CN111239658A (en) * 2020-02-19 2020-06-05 上海康达卡勒幅医疗科技有限公司 Method for reducing PROPELLER sequence artifacts
CN111239658B (en) * 2020-02-19 2022-07-19 上海康达卡勒幅医疗科技有限公司 Method for reducing PROPELLER sequence artifacts

Also Published As

Publication number Publication date
CN101422364B (en) 2010-05-19

Similar Documents

Publication Publication Date Title
US8488860B2 (en) Magnetic resonance imaging apparatus
CN110554334B (en) Method and system for coil selection in magnetic resonance imaging
CN101406396B (en) BO field drift correction in a temperature map generated by magnetic resonance tomography
CN106056537B (en) A kind of medical image joining method and device
EP2020203B1 (en) Medical guiding system
CN104484659B (en) A method of to Color medical and gray scale image automatic identification and calibration
CN101408976B (en) Enhanced system and method for volume based registration
CN101744618B (en) One-equipment room transmission PET/CT/MR image collection, registering and imaging system and method
CN105455849A (en) Mammary gland volume ultrasonic imaging device and method
CN101511273A (en) Ultrasonic breast diagnostic system
CN100559397C (en) The affine registration method and system of two dimensional image and operation forward three-dimensional viewing in the operation
CN102239495A (en) A method, apparatus, and computer program product for acquiring medical image data
KR101629058B1 (en) Ultrasonic diagnosis apparatus and program for controlling the same
CN1572246A (en) Parallel imaging based on expanded unfolding technique
US20200138382A1 (en) Methods and systems for whole body imaging
CN101422364B (en) Magnetic resonance image seamless splicing method
CN102499701B (en) Geometrical calibrating method for X-ray and fluorescent double-mode living body imaging system
CN101194845A (en) Synchronous real time dynamic tracing solid locating method and device based on dual-X-ray machine
CN105361955A (en) Medical imaging unit, medical imaging device with a medical imaging unit, and method for detecting a patient movement
CN105741303A (en) Medical image acquisition method
CN112946547A (en) System and method for generating diagnostic scan parameters from calibration images
CN107527339A (en) A kind of magnetic resonance scanning method, apparatus and system
CN110084805B (en) FOV parameter setting method and device and image processing equipment
CN108027413A (en) One kind is used for the method for calibrating magnetic resonance imaging (MRI) body mould
CN111493932A (en) Ultrasonic imaging method and system

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20100519

Termination date: 20191102

CF01 Termination of patent right due to non-payment of annual fee