CN1627095A - Method for registering non-rigid brain image based on non-homogeneous rational base spline - Google Patents
Method for registering non-rigid brain image based on non-homogeneous rational base spline Download PDFInfo
- Publication number
- CN1627095A CN1627095A CN 200310120540 CN200310120540A CN1627095A CN 1627095 A CN1627095 A CN 1627095A CN 200310120540 CN200310120540 CN 200310120540 CN 200310120540 A CN200310120540 A CN 200310120540A CN 1627095 A CN1627095 A CN 1627095A
- Authority
- CN
- China
- Prior art keywords
- image
- nurbs
- rigid body
- control vertex
- brain image
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 210000004556 brain Anatomy 0.000 title claims abstract description 21
- 238000005457 optimization Methods 0.000 claims abstract description 14
- 238000013507 mapping Methods 0.000 claims description 12
- 230000005489 elastic deformation Effects 0.000 claims description 2
- 230000001915 proofreading effect Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 3
- 230000009466 transformation Effects 0.000 description 6
- 238000003745 diagnosis Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 210000001259 mesencephalon Anatomy 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 230000014616 translation Effects 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
This invention puts forward a very new non-rigid brain image registration method based on NURBS, which can either realize registrations among small deformation images or images of large differences. This invention combines FFD with NURBS to seak the non-rigid deformation parameters of the image by multiplayer optimization characterizing that, first of all, NURBS replaces the uniform B sample function to increase its flexibility, secondly, we utilize the reverse image technology of the image process to increase the accuracy of the registration result.
Description
Technical field
The present invention relates to medical MRI (magnetic resonance image (MRI)) technical field of image processing, particularly a kind of non-rigid body brain image method for registering based on NURBS (non-uniform rational basic spline).Be mainly used in medical diagnosis midbrain images match problem, belong to intelligent information processing technology.
Background technology
Registration problems is a focus of studying in the Medical Image Processing.In general, registration can be divided into two classes: rigid body registration and non-rigid body registration.The rigid body registration is sought six unknown quantitys exactly, that is, and and three translations and three rotations.Its feature is to have affine unchangeability.Be that two points above the piece image are mapped on the other piece image, distance therebetween remains unchanged; Parallel segment is through keeping parallelism still after the conversion.Obviously, if conversion is non-rigid body between the image, the rigid body registration just can not correctly be described the relation between the image.At this time just need non-rigid body registration, this also is the emphasis of current people's research.
The non-rigid body Study of Registration of medical image a lot of method (J.V.Hajnal occurred through the development in nearly ten years, D.L.G.Hill, and D.J.Hawkes, editors.Medical ImageRegistration.CRC Press, 2001.), in these methods, Free Transform image registration (D.Rueckert of Rueckert proposition more effectively based on even basic spline, L.I.Sonoda, D.L.G.Hill, M.O.Leach, and D.J.Hawkes, " NonrigidRegistration Using Free-Form Deformations:Application to BreastMR Images; " IEEE Trans.Med.Imag., vol.18, no.8, pp.712-721, Aug.1999.).Generally speaking, their method can achieve satisfactory results.But when we found that this method diversity ratio between image is bigger, for example: the registration results between the different people was poor.By analyzing, we think that the reason of registration failure is because the FFD (Free Transform) that they use only is confined to the Uniform B-spline basis function, in order to address this problem, we have proposed a kind of more effective non-rigid body brain image registration based on NURBS (non-uniform rational basic spline).
Summary of the invention
The objective of the invention is to deficiency at existing method for registering, a kind of brand-new non-rigid body brain image method for registering based on NURBS is proposed, this method not only can realize registration between the image of small deformation, and also can be quick for the bigger situation of image difference, accurately registration.Method of the present invention can be assisted diagnosis, works out the enforcement of operation plan and operation guiding system etc., simultaneously, has reduced the human intervention in the case diagnosis process, has bigger medical value.In order to make the result more accurate, the sub-separating method of NURBS is introduced into wherein, has promptly increased degree of freedom, is used in combination with the method for multi-level optimization simultaneously.At last, we when the initialization control vertex, have considered the characteristics of brain image according to the observation to brain image, for registration provides a reasonable initial value.
For realizing such purpose, the present invention seeks the non-rigid body deformation parameter of image by multi-level optimization with FFD and NURBS based on this, has obtained reasonable experimental result.We are main innovate point.At first, NURBS has replaced even basic spline basis function.Because NURBS can adopt reference mark heterogeneous and knot vector to distribute.Increased dirigibility.Secondly, in the computational transformation image, we have utilized the reverse mapping techniques of Flame Image Process, rather than resemble the mapping of previous method employing forward.Improved the accuracy of registration results.The 3rd, more accurate in order to make the result, the sub-separating method of NURBS is introduced into wherein, has promptly increased degree of freedom, is used in combination with the method for multi-level optimization simultaneously.At last, we when the initialization control vertex, have considered the design feature of human brain according to the observation to brain image, for registration provides a reasonable initial value.
Non-rigid body brain image method for registering based on NURBS proposed by the invention comprises Free Transform, oppositely mapping, steps such as NURBS branch and multi-level optimization:
1) Free Transform
Behind the brain image process rigid body registration, the image that alignd on the whole, at this time just remaining local elastic deformation, non-rigid body registration is used for proofreading and correct these resilient bias.
If adopt the Free Transform based on Uniform B-spline, basis function and control vertex can only be confined to even distribution so.Some dirigibilities have just been lost.In order to overcome this problem, the present invention has adopted the method based on NURBS.Because NURBS not only has all character of B batten, and it also has the unevenness of basis function and control vertex.This just makes us when the search deformation parameter, can the more control summit be set in the big place of distortion, and in the zone that has matched fewer control vertex be set.After promptly adopting NURBS to replace the Uniform B-spline basis function, distortion has more dirigibility.
Here V represents control vertex, and W is a weight, and B represents the B spline base function, and p, q, r are the number of control vertex in three directions.Being arranged on the image of control mesh each point like this can be represented by this deformation formula.As mentioned above, according to the NURBS related properties, we can be provided with control vertex and knot vector according to the situation of real image.
2) oppositely mapping
Method for registering in the past adopts the forward mapping more when the computational transformation image.But we know, adopt the forward mapping, several pixels that a pixel of source images may corresponding target image, and the pixel above the target image might not find corresponding point on source images.That is to say, in target image, may produce the cavity.Method in the past normally adopts linear interpolation to come these cavities of interpolation.But when image transformation spacing relatively more violent or control vertex was excessive, the image that interpolation is come out can produce discontinuous.The present invention then uses reverse mapping to avoid this problem.At first, FFD is applied in the changing image on each voxel, obtains the position of this voxel in source images, adopts the B spline interpolation to obtain the gray-scale value of this voxel then in source images.Because the method that has adopted reverse mapping and B spline interpolation to combine, the uncontinuity of image in result images, avoided, and then improved result's accuracy.
3) NURBS branch combines with multi-level optimization
When the gap ratio of control vertex is bigger, that is to say that the degree of freedom of describing distortion is not enough, is necessary to add some new control vertexs more deformation parameter is provided.And in Computer-aided Geometric Design, various parametric lines, curved surface can reach the purpose of accurate description distortion by the way that son divides.Equally, nurbs curve can increase deformability by insert new summit in the middle of control vertex.Under the 1D situation, P is original control vertex, and P ' is new control vertex, P
iCorresponding to new control vertex P '
2i
This seed separating method of the present invention expands among the 3D entity.Obtain eight sons and divide expression.So whenever, do a second son branch, the control vertex number increases nearly octuple.
Simultaneously, we sample to image.At first, on minimum image in different resolution, the control vertex corresponding to given is optimized calculating, has calculated after one deck, increases image resolution ratio, simultaneously control vertex is carried out son and divides, until reach satisfied effect.That is to say the method that we have adopted image to combine with the control vertex multi-level optimization.Not only avoid optimization to be absorbed in local minimum, and obtained better result.
The present invention has utilized FFD and NURBS way of combining to seek the image change parameter, can be quick, accurately brain image is carried out registration.The present invention has clinical medical value, can be used for the lesion detection of medical image, the comparison of focus between the patient, normal person's growth comparison and the inspection of surgery recovery situation.The orientation problem etc. that also can be used for simultaneously, operation guiding system.
Description of drawings
Fig. 1 is the non-rigid body brain image registration Algorithm schematic diagram based on NURBS of the present invention.
Fig. 2 is an experimental result picture of the present invention.
Embodiment
For understanding technical scheme of the present invention better, be further described below in conjunction with accompanying drawing and specific embodiment.
The non-rigid body brain image method for registering principle that the present invention is based on NURBS as shown in Figure 1.Among the figure, given source images S and target image T, our target is to find deformation parameter, makes image S to change to image T by this continuous transformation.At first, according to concrete conditions such as image sizes, initial control vertex is set.Then image S is carried out Free Transform based on NURBS, in the middle of this process, we sample to image, and the number of times that specifically adopts the number of plies and NURBS to divide is corresponding.From lowermost layer one by one to top, on each layer, adopt oppositely that mapping comes the computational transformation image, the mutual information between computational transformation image and the target image again is when this layer calculates end, judge whether to reach best, if do not reach, just carry out NURBS and divide, increase image resolution ratio, continuing to calculate on one deck down, until reaching optimum.
Embodiment
In order to verify our algorithm performance, we experimentize on real people's brain image.Our image is T
1Weight term [9], size 256 * 256 * 64, bed thickness 3mm, size is 1 * 1mm in the layer, and the gray scale nonuniform field is 20%, and the initial control vertex that we get is 9 * 9 * 7.3 layers of optimization have been carried out altogether.We have chosen 7 different people and have experimentized, will be wherein; One personal accomplishment target image, all the other six respectively as source images.Fig. 2 has provided two width of cloth images wherein, a is a source images, b is a target image, c is process rigid body images after registration, d is the difference image between image b and the c, the result that e obtains for the method with us, and f is the difference image between image b and the e, the image of g for adopting the Uniform B-spline method to obtain, h is the difference image between b and the g.As can be seen from the figure our method improves a lot on accuracy, mismatches accurate a large amount of the minimizing.
For quantized result, we have compared result images, and table 1 has been listed the related coefficient (CC) between them.Have 6 class values altogether, and their and method of using even FFD are compared.Our method all has raising to each patient as can be seen from the table.
Table 1 related coefficient (CC)
Patien ???t | ????1# | ????2# | ????3# | ????4# | ????5# | ????6# |
Even FFD | ?0.9487 ????6 | ?0.8857 ????9 | ?0.9316 ????5 | ?0.9495 ????2 | ?0.8843 ????6 | ?0.9474 ????9 |
?NURBS | ?0.9494 ????1 | ?0.9074 ????2 | ?0.9321 ????8 | ?0.9501 ????4 | ?0.9082 ????1 | ?0.9480 ????9 |
As can be seen from the test results, compare with the Uniform B-spline method for registering, our method is effectively for the registration of different people, has also improved result's precision simultaneously.
Claims (4)
1, a kind of brand-new non-rigid body brain image method for registering based on NURBS comprises Free Transform, oppositely mapping, and NURBS divides and the multi-level optimization step:
1) Free Transform: behind the brain image process rigid body registration, the image that alignd on the whole, at this time just remaining local elastic deformation, non-rigid body registration is used for proofreading and correct these resilient bias;
2) oppositely mapping: FFD is applied in the changing image on each voxel, obtain the position of this voxel in source images, adopt the B spline interpolation to obtain the gray-scale value of this voxel then in source images, the reverse method that combines of mapping and B spline interpolation has guaranteed the continuity of result images;
3) NURBS branch combines with multi-level optimization: curved surface reaches accurate description distortion by the child branch, and nurbs curve inserts new point in the middle of control vertex.
2, according to the non-rigid body brain image method for registering based on NURBS of claim 1, it is characterized in that, in the Free Transform step, employing is based on the method for NURBS Free Transform, NURBS has all character of B batten, also has the unevenness of basis function and control vertex, during the search deformation parameter, can the more control summit be set in the big place of distortion, few control vertex be set matching or be out of shape few zone.
3, according to the non-rigid body brain image method for registering based on NURBS of claim 1, it is characterized in that divide with during multi-level optimization combines at NURBS, sub-separating method expands among the 3D entity, simultaneously, image is sampled, on minimum image in different resolution, corresponding to given control vertex, be optimized calculating, calculated after one deck, increased image resolution ratio, simultaneously control vertex has been carried out son and divide.
4, according to the non-rigid body brain image method for registering based on NURBS of claim 1 or 3, it is characterized in that, adopt image and the method that the control vertex multi-level optimization combines, avoided optimization to be absorbed in local minimum.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101205405A CN100346170C (en) | 2003-12-12 | 2003-12-12 | Method for registering non-rigid brain image based on non-homogeneous rational base spline |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101205405A CN100346170C (en) | 2003-12-12 | 2003-12-12 | Method for registering non-rigid brain image based on non-homogeneous rational base spline |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1627095A true CN1627095A (en) | 2005-06-15 |
CN100346170C CN100346170C (en) | 2007-10-31 |
Family
ID=34761536
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2003101205405A Expired - Fee Related CN100346170C (en) | 2003-12-12 | 2003-12-12 | Method for registering non-rigid brain image based on non-homogeneous rational base spline |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100346170C (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101666865B (en) * | 2009-09-09 | 2012-12-19 | 西北工业大学 | Method for registering diffusion tensor nuclear magnetic resonance image in local quick traveling mode |
CN102959584A (en) * | 2011-12-21 | 2013-03-06 | 中国科学院自动化研究所 | Function magnetic resonance image registration method |
CN103854276A (en) * | 2012-12-04 | 2014-06-11 | 株式会社东芝 | Image registration device and method, image segmentation device and method and medical image device |
CN105488804A (en) * | 2015-12-14 | 2016-04-13 | 上海交通大学 | Brain ASL (Arterial Spin Labeling), SPECT (Single-Photon Emission Computed Tomography) and MRI (Magnetic Resonance Imaging) image registration and fusion conjoint analysis method and system |
CN112102384A (en) * | 2020-10-14 | 2020-12-18 | 山东大学 | Non-rigid medical image registration method and system |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5365996A (en) * | 1992-06-10 | 1994-11-22 | Amei Technologies Inc. | Method and apparatus for making customized fixation devices |
US5377320A (en) * | 1992-09-30 | 1994-12-27 | Sun Microsystems, Inc. | Method and apparatus for the rendering of trimmed nurb surfaces |
US5800353A (en) * | 1996-02-12 | 1998-09-01 | Mclaurin, Jr.; Robert L. | Automatic image registration of magnetic resonance imaging scans for localization, 3-dimensional treatment planning, and radiation treatment of abnormal lesions |
JP2000300538A (en) * | 1999-04-23 | 2000-10-31 | Ge Yokogawa Medical Systems Ltd | Image displaying method, image displaying apparatus and mri apparatus |
US6486667B1 (en) * | 2000-03-31 | 2002-11-26 | Koninklijke Philips Electronics N.V. | Combination of fluid-attenuated inversion-recovery complex images acquired using magnetic resonance imaging |
JP3577512B2 (en) * | 2001-07-17 | 2004-10-13 | 独立行政法人情報通信研究機構 | Program and apparatus for extracting cerebral blood vessels from MRA images |
DE10203237B4 (en) * | 2002-01-28 | 2006-08-17 | Siemens Ag | Method for magnetic resonance imaging with automatic adjustment of the measuring field |
-
2003
- 2003-12-12 CN CNB2003101205405A patent/CN100346170C/en not_active Expired - Fee Related
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101666865B (en) * | 2009-09-09 | 2012-12-19 | 西北工业大学 | Method for registering diffusion tensor nuclear magnetic resonance image in local quick traveling mode |
CN102959584A (en) * | 2011-12-21 | 2013-03-06 | 中国科学院自动化研究所 | Function magnetic resonance image registration method |
WO2013091181A1 (en) * | 2011-12-21 | 2013-06-27 | 中国科学院自动化研究所 | Functional magnetic resonance image registration method |
US8965093B2 (en) | 2011-12-21 | 2015-02-24 | Institute Of Automation, Chinese Academy Of Sciences | Method for registering functional MRI data |
CN103854276A (en) * | 2012-12-04 | 2014-06-11 | 株式会社东芝 | Image registration device and method, image segmentation device and method and medical image device |
CN103854276B (en) * | 2012-12-04 | 2018-02-09 | 东芝医疗***株式会社 | Image registration and segmenting device and method, and medical image equipment |
CN105488804A (en) * | 2015-12-14 | 2016-04-13 | 上海交通大学 | Brain ASL (Arterial Spin Labeling), SPECT (Single-Photon Emission Computed Tomography) and MRI (Magnetic Resonance Imaging) image registration and fusion conjoint analysis method and system |
CN105488804B (en) * | 2015-12-14 | 2018-06-12 | 上海交通大学 | Brain ASL, SPECT and the method and system of MRI image registration fusion Conjoint Analysis |
CN112102384A (en) * | 2020-10-14 | 2020-12-18 | 山东大学 | Non-rigid medical image registration method and system |
CN112102384B (en) * | 2020-10-14 | 2024-03-15 | 山东大学 | Non-rigid medical image registration method and system |
Also Published As
Publication number | Publication date |
---|---|
CN100346170C (en) | 2007-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109389587B (en) | Medical image analysis system, device and storage medium | |
CN109166133A (en) | Soft tissue organs image partition method based on critical point detection and deep learning | |
CN113570627B (en) | Training method of deep learning segmentation network and medical image segmentation method | |
Xiong et al. | Mapping histological slice sequences to the allen mouse brain atlas without 3D reconstruction | |
CN112862805B (en) | Automatic auditory neuroma image segmentation method and system | |
Lui et al. | Shape-based diffeomorphic registration on hippocampal surfaces using beltrami holomorphic flow | |
Tian et al. | Improving in vivo human cerebral cortical surface reconstruction using data-driven super-resolution | |
Dai et al. | Locating anatomical landmarks on 2D lateral cephalograms through adversarial encoder-decoder networks | |
CN110223331B (en) | Brain MR medical image registration method | |
CN115830016A (en) | Medical image registration model training method and equipment | |
CN116650115A (en) | Orthopedic surgery navigation registration method based on UWB mark points | |
CN100346170C (en) | Method for registering non-rigid brain image based on non-homogeneous rational base spline | |
CN110517300A (en) | Elastic image registration algorithm based on partial structurtes operator | |
Lapeer et al. | An optimised radial basis function algorithm for fast non-rigid registration of medical images | |
CN116309754A (en) | Brain medical image registration method and system based on local-global information collaboration | |
CN113850710B (en) | Cross-modal medical image accurate conversion method | |
CN115937581A (en) | Classification and identification method of multi-site fMRI data | |
CN112598669B (en) | Lung lobe segmentation method based on digital human technology | |
Monasterio-Exposito et al. | Label augmentation to improve generalization of deep learning semantic segmentation of laparoscopic images | |
CN110276762A (en) | A kind of full-automatic bearing calibration of respiratory movement of the diffusion-weighted Abdominal MRI imaging of more b values | |
CN114419309A (en) | High-dimensional feature automatic extraction method based on brain T1-w magnetic resonance image | |
CN114463339A (en) | Medical image segmentation method based on self-attention | |
Marias et al. | Image analysis for assessing molecular activity changes in time-dependent geometries | |
CN111798500A (en) | Differential homoembryo non-rigid registration algorithm based on hierarchical neighborhood spectrum features | |
Mayer et al. | Direct registration of white matter tractographies and application to atlas construction |
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: 20071031 Termination date: 20171212 |
|
CF01 | Termination of patent right due to non-payment of annual fee |