Skip to main content

ORIGINAL RESEARCH article

Front. Public Health, 01 June 2022
Sec. Digital Public Health

Non-rigid Multi-Modal Medical Image Registration Based on Improved Maximum Mutual Information PV Image Interpolation Method

  • School of Computer and Information Science, Southwest University, Chongqing, China

With the continuous improvement of medical imaging equipment, CT, MRI and PET images can obtain accurate anatomical information of the same patient site. However, due to the fuzziness of medical image physiological evaluation and the unhealthy understanding of objects, the registration effect of many methods is not ideal. Therefore, based on the medical image registration model of Partial Volume (PV) image interpolation method and rigid medical image registration method, this paper established the non-rigid registration model of maximum mutual information Novel Partial Volume (NPV) image interpolation method. The proposed NPV interpolation method uses the Davidon-Fletcher-Powell algorithm (DFP) algorithm optimization method to solve the transformation parameter matrix and realize the accurate transformation of the floating image. In addition, the cubic B-spline is used as the kernel function to improve the image interpolation, which effectively improves the accuracy of the registration image. Finally, the proposed NPV method is compared with the PV interpolation method through the human brain CT-MRI-PET image to obtain a clear CT-MRI-PET image. The results show that the proposed NPV method has higher accuracy, better robustness, and easier realization. The model should also have guiding significance in face recognition and fingerprint recognition.

Introduction

Using artificial intelligence technology to collate and analyze medical images (14) can improve medical and health services (5) and alleviate the imbalance of medical and health resource (69) allocation, which is the focus of the development of medical artificial intelligence. Many previous studies have mentioned that the application of artificial intelligence in the field of smart medicine should be deepened, new models and means of artificial intelligence diagnosis and treatment (10) should be promoted, and a fast and accurate intelligent medical system (11) should be established. As the main component of smart medicine, medical image plays an increasingly important role in assisting medical diagnosis and promoting the rationalization of medical decision-making. Common medical images include X-ray, computed tomography (CT) (12), magnetic resonance imaging (MRI) (13, 14), positron emission tomography (PET) (15) and single-photon emission computed tomography (SPECT) (15). Medical image registration (1619) is a basic task in the process of medical image processing. In clinical practice, it is usually necessary to perform multiple imaging in multiple modes or in the same mode for the same patient, that is, to obtain information from several images and conduct comprehensive analysis. The image registration technology (2023) is used to fuse the above images, and the multi-faceted information of the human body is expressed on the same image at the same time. The internal structure and functional state of the human body are reflected from the medical image (2427), and the anatomical and physiological pathological information of the human body is more directly provided, thus playing the role of multi-information visualization at the same time.

After years of development, computer in-depth learning technology (28) and artificial intelligence technology (29) have made significant breakthroughs in theory and practice. They have made revolutionary progress in realizing the open sharing of medical information and using artificial intelligence to organize and analyze fragmented medical information. The combination of image fusion analysis and radiology is of great significance for understanding complex diseases and realizing accurate judgment, and is of great help for proposing the best diagnosis results and treatment schemes. In the process of imaging diagnosis (30, 31), single image data often cannot fully display the pathological structure of patients, and pathological features need to be expressed by multi-directional image features. Compared with a single sensor, multi-sensor (32, 33) can provide more extensive data for better identification of tissue or lesion details.

Some scholars focus on the improvement of fusion methods, such as Xiao and Pedrycz (34) and Xiao (35) proposed fuzzy value decision theory and extended it to the complex domain. Other scholars focus on improving interpolation algorithms, such as in 2005, Xiao (36) proposed a qualitative perturbation PV interpolation algorithm to avoid local extremum on grid points and non-grid points. In 2010, Xiao (37) proposed an improved Blackman-Harris PV (BHPV) interpolation algorithm, and introduced the Blackman-Harris window sinc function to replace the kernel function of the traditional PV interpolation. He can improve the extreme value problem of the traditional method. In 2021, Zang et al. (38) proposed an improved PV interpolation method to calculate the chord function. However, these methods do not pay much attention to the characteristics of the model itself. Therefore, how to fuse multiple image information with minimum complexity cost to obtain high reliability registration results is the primary problem to be solved in medical image registration.

To solve this problem, through the analysis and research on rigid medical image registration based on the maximum mutual information PV (Partial Volume) (39) image interpolation method, this paper proposes a decision fusion model for multiple information image data, namely the non-rigid registration model based on the maximum mutual information NPV (Novel Partial Volume) image interpolation. The cubic B-spline interpolation method can obtain better decision results under various uncertain information conditions. In addition, the effect is particularly significant in image-assisted diagnosis and treatment.

The main contributions of this paper are listed below:

• The proposed NPV method combines cubic spline with DFP algorithm for the first time, and introduces image registration to make the registration result more reasonable and reliable.

• By comparing the synthetic image with three original images to be registered, the image registration method proposed in this paper has great advantages in comprehensively displaying pathological details, balancing the advantages of each image, and matching the corresponding spatial points.

• The improved method proposed in this paper is the traditional maximum mutual information PV image interpolation method based on gray information, so that the registration model does not need to segment the image. It is simple, accurate and robust.

The organizational structure of this paper is as follows. The second chapter focuses on the problems to be solved and the relevant implementation steps of the NPV interpolation method. In the third chapter, some examples are introduced to prove the effectiveness of the proposed method. In the fourth chapter, the full text is summarized.

Preliminaries

PV Image Interpolation

PV image interpolation (40) is a common registration method in medical research, which is convenient in updating the joint histogram of two images. Suppose that the reference image is R, and the image to be registered is F for any point p on image F. After transforming T, the point q in the reference image R is corresponding, that is, q = T(p). The adjacent pixels of q are n1, n2, n3, and n4, respectively, as shown in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. PV interpolation image.

Definition 1: Let the joint histogram be h, the remaining parameters are shown in Figure 1.

iWi=1,i=1,2,3,4,    (1)
ihrm[r(x),f(ni)]+=wi,    (2)
w1=(1-dx)×(1-dy),    (3)
w2=dx×(1-dy),    (4)
w3=dx×dy,    (5)
w4=(1-dx)×dy.    (6)

Definition 2: Let f(x) be the classical PV interpolation kernel function plotted in Figure 2, defined as:

f(x)={1-|x|,|x| <10,|x| 1.    (7)

In PV interpolation, the value of each gray pair in the joint histogram is the sum of several decimals with small changes, so it not only makes the calculated mutual information value more accurate, but also smooths the change of the mutual information curve when the registration has a small transformation.

FIGURE 2
www.frontiersin.org

Figure 2. Basic kernel curve.

The process of medical image registration mainly includes four key parts, graphic transformation, registration optimization method, deformation image interpolation method and similarity measure analysis. The flow chart is shown in Figure 3. First of all, it is necessary to find the appropriate graphic transformation model and select the appropriate optimization method to solve the equation. Then, the floating image is transformed on the results obtained by solving, and the deformed image is interpolated. Finally, the similarity measure of the two images is calculated, and whether the next iteration is needed according to the similarity degree.

FIGURE 3
www.frontiersin.org

Figure 3. The flowchart of PV image interpolation method.

Proposed Methods

The registration process based on image intensity is actually an optimization problem, that is, after selecting the appropriate similarity measure for the image to be registered, the similarity measure between images can be maximized as far as possible. The registration process is mainly composed of five parts, image preprocessing, image transformation, gray interpolation, objective function optimization, and similarity measure calculation. When achieving high-precision alignment of medical images, we select the reference image CT as the reference image, and select PET and MRI as floating image information sources to provide information.

In this chapter, we will focus on the five steps of the proposed method, as shown in Figure 4. First of all, we need to extract the parameter values of the graphics. The second step of image transformation is necessary in the coarse registration. In the third step, the DFP optimization algorithm is used to optimize the image. In the fourth step, the deformed graphics are interpolated. Finally, the similarity analysis of the graphics is completed and the results are output.

FIGURE 4
www.frontiersin.org

Figure 4. Procedure followed by the proposed method.

Preliminaries

Coordinate origin is defined in the gray center of the image, the coordinates are (cx, cy, cz), g(x, y, z) is the gray value of point (cx, cy, cz).

cx=xxg(x,y,z)xg(x,y,z),cy=yyg(x,y,z)yg(x,y,z),cz=zzg(x,y,z)zg(x,y,z).    (8)

For the images to be registered, in order to improve the registration accuracy, it is necessary to remove the interference such as mask and bed frame in the image. In order to improve the registration speed and balance the relationship between the registration speed and accuracy, some voxel points are usually selected as the sample set of the objective function, namely, the images to be registered are sampled. In this paper, regular grid spacing is used for sampling. Several individual pixels are selected from CT, MRI, and PET images to be registered as sample sets.

Image Transformation

The key to register multiple different modal images is to find a suitable transformation model so that all pixels on the floating image can take you through the transformation model to find the unique pixel coordinates corresponding to them in the reference image and maintain the same or similar anatomical information. In this paper, rigid transformation is used for global coarse registration and then non-rigid transformation is used for local elastic registration.

(1) Rigid transformation based on affine transformation

Rigid transformation mainly includes linear transformation and affine transformation. The affine transformation has the ability to map the floating image line to the line in the reference image space and ensure the parallel relationship between different lines, and it can be expressed as a combination of linear transformations such as rotation, scaling and translation. Therefore, this paper uses affine transformation method for coarse registration of CT, MRI, and PET images.

Assuming that the pixel point on the original image is the spatial coordinate p0(x0,y0,z0) after p0(x0, y0, z0) translation transformation, the affine transformation can be expressed as:

[x0y0z01]=[a11a12a130a21a22a230a31a32a3300001] [x0y0z01].    (9)

(2) Non-rigid transformation based on cubic B-spline

For a three-dimensional floating image Ω = {(x, y, z)|0 ≤ xX, 0 ≤ yY, 0 ≤ zZ}, the control points ϕi,j,k in the image constitute the size of the grid nx, ny, nz, and limit its movement within the range of four grid points in its neighborhood, as below:

nx=Xλx,nx=Yλy,nx=Zλz.    (10)

Among them, λx, λy, λz represents the distance of grid points on the x, y, z axis, and selects different sizes according to different registration requirements. Usually, the number that can divide the image size X, Y, Z is selected to reduce the computational complexity and ensure the stability of image transformation. The single transformation T is used to transform the spatial coordinates of the control point to (Δx, Δy, Δz), and the coordinate transformation of all points in the neighborhood of the control point (x, y, z) is fitted with the change of the control point, which is expressed as (x + Δx, y + Δy, z + Δz). The mathematical expression of the transformation is as follows.

T(x,y,z)=i=03m=03n=03Bi(u)Bm(w)ϕi+l,j+m,k+n.    (11)

Among them, i, j, k is the index of the control point coordinates (Δx, Δy, Δz) in the grid, u, v, w is the distance between the corresponding point and the integer coordinate point, l, m, n is the number of B-spline basis functions, Bl, Bm, Bn represents l, m, n basis functions, respectively. The calculation formulas of each parameter are as follows:

{i=xλx-1j=yλy-1k=zλz-1    (12)
{u=x-xv=y-yw=z-z.    (13)

where ⌊x⌋, ⌊y⌋, ⌊z⌋ represents an integer not greater than x, y, z. The expressions of zeroth to cubic B-spline basis functions are as follows.

{B0(u)=(1-u)36B1(u)=3u3-6u2+46B2(u)=-3u3+3u2+3u+16B3(u)=u36    (14)

DFP Optimization Algorithm

The transformation matrix is essentially an optimization problem for solving multi-parameter objective function. The use of DFP in the registration process can make the calibrated image more optimized.

In the registration process, the cubic B-spline transformation function is used as the objective function, and the transformation parameters are solved by DFP algorithm. The specific steps are as follows:

(1) Given parameters δ ∈ (0, 1), σ ∈ (0, 0.5), initial point x0ϵn, termination error 0 ≤ ε ≤ 1. Initially symmetric positive definite matrix H0 and let k = 0.

(2) Calculate gk = ∇f(xk), when ||gk|| ≤ ε stops, output xk as the approximate minimum point.

(3) Calculate the search direction dk = −Hkgk.

(4) Let mk be the minimum non-negative integer m satisfying the following inequality.

f(xk+δmdk)f(xk)+σδmgkTdk.    (15)

Let αk=δmk,xk+1=xk+αkdk, then calculate sk, yk;

(5) From Equation (16) calculate Hk+1;

Hk+1={Hk,             skTyk0Hk-HkykykTHkykTHkyk+skskTskTyk,ykTsk>0.    (16)

(6) Let k = k + 1, turn to Step 1.

Image Interpolation Algorithm After Deformation

Since the reference image and floating image are discrete image data, the pixel coordinates are integers, and the position of the pixel coordinates of the image obtained by graphic transformation is not necessarily exactly located in the integer coordinates. Therefore, it is inevitable to interpolate the pixels on the transformed image. Interpolation belongs to a resampling method of the image. It not only needs to consider the image accuracy after interpolation, but also needs to consider the computational complexity and the protection of the original image data. PV and trilinear interpolation methods are the most widely used in interpolation calculation based on mutual information registration. In the three-dimensional case, the objective function T = f(x) satisfies the following conditions for the point m(xi, xj, xk) of floating image MRI and the point n(Yi + ηj, Yj + ηj, Yk + ηk) of reference image CT.

(1) T(x) ≥ 0, xZ;

(2) m=-+f(m-η)=1,mZ,0η<1

Because cubic spline satisfies the above two conditions at the same time, using cubic spline as kernel function f(x) as below, plotted in Figure 5.

f(x)={ax3+bx2+cx+d,|x|<10,otherwise    (17)

And the gray contribution histogram can be obtained as follows:

h[F(xi,xj,xk),R(yi+p,yj+q,yk+r)]+=f1(p-ηi)f2(q-ηj)f3(r-ηk).    (18)

where p, q, r ∈ (0, 1) as an integer.

FIGURE 5
www.frontiersin.org

Figure 5. Flow chart of DFP optimization algorithm.

With the increase of kernel function order, the local extremum generated by interpolation is eliminated. However, the computational complexity will be more when the order is >3. Therefore, considering the registration time benefit and accuracy, cubic spline is selected for gray histogram calculation.

The kernel function obtained by cubic spline interpolation method is shown in Figure 6. Compared with the basic image interpolation method (Figure 2), the image interpolation curve obtained by the improved method is smoother, which effectively avoids the occurrence of burrs and sharp points. Especially, the proposed NPV method has better accuracy for changing smaller points.

FIGURE 6
www.frontiersin.org

Figure 6. Cubic B-spline interpolation kernel function interpolation curve.

Similarity Measure Analysis

Similarity measure is the consideration of calculating the similarity between reference image and floating image. This paper uses the similarity measure based on gray level, which includes the sum of squares of error, mean square error, mutual information and other methods. Since mutual information is an automatic calculation method, in order to improve the registration speed, this paper uses mutual information method to analyze the similarity measure.

The calculation of mutual information is based on Shannon entropy, which is defined as follows.

H=ipilog1pi=-ipilogpi.    (19)

where pi represents the probability of information appearing. The mutual information between reference image CT and floating image MRI is as follows.

H(RCT)=-rpR(r)logpR(r).    (20)
H(FMRI)=-fpF(f)logpF(f).    (21)

The joint entropy of two images is defined as follows.

H(RCT,FMRI)=-r,fpPF(r,f)logpRF(r,f).    (22)

NMI(R,F)=H(RCT) + H(FMRI)H(RCT,FMRI) is used to calculate mutual information, and the registration is completed when mutual information reaches the maximum.

Examples

The brain CT, MRI, and PET images with AD were extracted from ANDI and OASIS databases for analysis. And we used Python to analyze and extract the data. Figure 7A is the CT reference image, Figure 7B is MRI reference data, and Figure 7C is PET reference data images. Medical image preprocessing includes image gray normalization, removing interference background and extracting region of interest. Define a three-dimensional coordinate system for two images to be registered. Taking CT as reference image R and MRI as floating image F, the gray level of body pixels on two images from different devices is normalized to 256 gray level.

FIGURE 7
www.frontiersin.org

Figure 7. CT/MRI/PET reference image. (A) shows CT reference data images. (B) is MRI reference data image. (C) is PET reference data image.

First, we use CT images as reference images and MRI images as floating images for medical registration. In Figures 8A–D are the segmented registration images of CT and MRI images, respectively. In Figure 8A, the initial registration image is preprocessed. Figure 8B is the medical image after multimodal registration. Figure 8C is a medical image of affine transformation global coarse registration. Figure 8D shows the medical image after high-precision registration of non-rigid transformation.

FIGURE 8
www.frontiersin.org

Figure 8. Registration images of brain CT and MRI. (A) shows the pre-processed initial registration image. (B) medical images after multi-modal registration. (C) is affine transformation global rough registration of medical images. (D) represents medical image with high precision registration after non-rigid transformation.

Secondly, we use CT-MRI images as reference images and PET images as floating images for medical registration. In Figures 9A–D are partial registration images of CT-MRI and PET images, respectively. Among them, Figure 9A pre-processed initial registration image. Figure 9B is the medical image after multimodal registration. Figure 9C is a medical image of affine transformation global coarse registration. Figure 9D For the final high-precision registration results, we add false colors to facilitate physician recognition.

FIGURE 9
www.frontiersin.org

Figure 9. Segmental registration images of brain CT-MRI and PET images. (A) shows the pre-processed initial registration image. (B) medical images after multi-modal registration. (C) is affine transformation global rough registration of medical images. (D) represents the medical image after non-rigid transformation and high-precision registration, and the addition of false color is convenient for physicians to identify.

Discussion

In one fusion, that is, when CT images and MRI images are aligned, the mutual information between the two is 0.781. In the second fusion, namely, when the CT-MRI image and PET image are aligned, the mutual information between the two is 0.622, as shown in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Comparison of mutual information fusion methods.

It can be seen from Table 1 that the proposed NPV method has improved the mutual information of the first and second time compared with the proposed PV method, indicating that this method has higher accuracy in graph registration.

In the process of image diagnosis, it is difficult for any single image to fully display the pathological structure and case characteristics of patients. Different images have prominent advantages for pathological diagnosis, and different imaging equipment, imaging angle and image mode often lead to many differences in diagnosis results. Therefore, reasonable medical image registration is of great significance. In order to verify the new registration algorithm proposed in this paper, CT images, MRI images, and PET images are used as original images for registration. Firstly, image preprocessing (Figures 8A, 9A) and multimodal registration (Figures 8B, 9B) were performed on the above three images. Then, global coarse registration (Figures 8C, 9C), B-spline function solution, image deformation and NPV interpolation, similarity measure analysis and other steps were performed on the first two images to obtain the registration results of CT and MRI images Figure 8D). Thirdly, the registration results were re-registered with PET images, and the secondary operation was completed according to the above steps (Figures 9A–C). The final synthesis results are shown in Figure 9D.

Conclusion

Medical image registration is an important subject in medical research and is of great significance in medical diagnosis and medical image analysis. This paper presents a new NPV image registration model. The proposed image registration model preprocesses the image before registration, so that the subsequent registration results are more accurate and reliable. The cubic spline interpolation and DFP optimization algorithm are introduced to operate the image. In addition, the proposed method also strictly controls the similarity measure, and ensures that the similarity measure is optimal through multiple iterations. The method has strong operability. By comparing the synthetic image with three original images to be registered, the image registration method proposed in this paper has great advantages in comprehensively displaying pathological details, balancing the advantages of each image, and matching the corresponding spatial points, which plays an important role in clinical application and theoretical research. Through the specific human brain image, compared with the traditional PV image, the sensitivity test shows that the image registration model proposed in this paper has high sensitivity, high reliability and popularization to gray change. Compared with the traditional registration method, this method has higher accuracy, better robustness and is easy to implement. In the direction of face recognition and fingerprint recognition, this model should also have certain reference significance.

However, due to the different imaging principles of medical images in different modes, there may be great differences in attributes. This method is a kind of non-rigid registration, which often requires the design of feature method and the adjustment and optimization of parameters, and it is difficult to obtain a highly universal registration effect. Moreover, although the proposed NPV image registration model is easy to implement, it needs a large amount of calculation and high-performance computer software, which is of little help to the reality of remote areas.

Therefore, our subsequent work is mainly to further improve the algorithm and extend it to other image recognition fields, so that it has higher applicability in areas with normal health conditions.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

Author Contributions

The author confirms being the sole contributor of this work and has approved it for publication.

Conflict of Interest

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Shihui Y, Wan Y, Shaoyi D, Jun S. Review of medical image registration based on deep learning. Pattern Recognit Artif Intell. (2021) 34:287–99. doi: 10.16451/j.cnki.issn1003-6059.202104001

PubMed Abstract | CrossRef Full Text

2. Hao P, Chen Z, Jiang S, Wang Y. Research on non-rigid registration of multimodal medical images based on Demons algorithm. Biomed Eng J. (2014) 31:161–5. doi: 10.7507/1001-5515.20140032

CrossRef Full Text | Google Scholar

3. Zhou X, Wang M, Zhou L, Wu Q. Research on image denoising model improved by Demons algorithm. Physics. (2015) 64:145–54. doi: 10.7498/aps.64.024205

CrossRef Full Text

4. Penney GP, Weese J, Little JA, Desmedt P, Hill DLG, Hawkes DJ. A comparison of similarity measures for use in 2-D-3-D medical image registration. IEEE Trans Med Imaging. (1998) 17:586–55. doi: 10.1109/42.730403

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Panayides S, Amini A, Filipovic ND, Sharma A, Tsaftaris SA, Young A, et al. AI in medical imaging informatics: current challenges and future directions. IEEE J Biomed Health Inf. (2020) 24:1837–57. doi: 10.1109/JBHI.2020.2991043

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Zhu F, Zhu X, Huang Z, Ding M, Li Q, Zhang X. Deep learning based data-adaptive descriptor for non-rigid multi-modal medical image registration. Signal Process. (2021) 183:108023. doi: 10.1016/j.sigpro.2021.108023

CrossRef Full Text | Google Scholar

7. Zhou Q, Yu Z. Medical image registration based on biogeographic optimization algorithm. J Yunnan Minzu Univers. (Natural Science Edition). (2020) 29:614–8. doi: 10.3969/j.issn.1672-8513.2020.06.014

CrossRef Full Text

8. Jun S, Linlin W, Shanshan W, Yanxia C, Gan W, Dongming W, et al. A review of the application of deep learning in medical imaging. Chin J Image Graph. (2020) 25:1953–81. doi: 10.11834/jig.200255

CrossRef Full Text

9. Xue W. Research on Multimodal Medical Image Registration Based on Demons. Lanzhou Jiaotong University (2018).

10. Kong M, He F, Li L. Development status and strategy of artificial intelligence assisted diagnosis and treatment. China Eng Sci. (2018) 20:86–91. doi: 10.15302/J-SSCAE-2018.02.013

CrossRef Full Text

11. Wu X, Pi Y, Zhou J. CT/MRI hybrid registration method and its application in radiotherapy planning system. Chin J Med Phys. (2020) 37:1148–54. doi: 10.3969/j.issn.1005-202X.2020.09.013

CrossRef Full Text

12. Jin Y, Ma G. Investigation and evaluation of optimal registration for medical CT images. In: 2010 3rd International Congress on Image and Signal Processing. Yantai: IEEE (2010). p. 2794–7.

Google Scholar

13. Yang W, Zhong L, Chen Y, Lin L, Lu Z, Liu S, et al. Predicting CT image from MRI data through feature matching with learned nonlinear local descriptors. IEEE Trans Med Imaging. (2018) 37:977–87. doi: 10.1109/TMI.2018.2790962

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Nikolaev AV, de Jong L, Weijers G, Groenhuis V, Mann RM, Siepel FJ, et al. Quantitative evaluation of an automated cone-based breast ultrasound scanner for MRI-3D US image fusion. IEEE Trans Med Imaging. (2021) 40:1229–39. doi: 10.1109/TMI.2021.3050525

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Brusaferri L, Emond EC, Bousse A, Twyman R, Whitehead AC, Atkinson D, et al. Detection efficiency modeling and joint activity and attenuation reconstruction in non-TOF 3-D PET from multiple-energy window data. IEEE Trans Radiat Plasma Med Sci. (2022) 6:87–97. doi: 10.1109/TRPMS.2021.3064239

CrossRef Full Text | Google Scholar

16. Van Holen R, Vandenberghe S. Optimization of a stationary small animal SPECT system for simultaneous SPECT/MRI. In: 2013 IEEE Nuclear Science Symposium and Medical Imaging Conference (2013 NSS/MIC) Seoul: IEEE (2013). p. 1–4.

Google Scholar

17. Cui Z, Mahmoodi S, Conway J, Guy M, Lewis E, Havelock T, et al. A novel non-rigid registration method based on nonparametric statistical deformation model for medical image analysis. In: 2017 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC). Atlanta, GA: IEEE (2017). p. 1–3.

Google Scholar

18. Huang P, Guo C, Chen H, Zhang H. An overview of image registration methods based on deep learning. J Hangzhou Univ Electron Sci Technol. (2020) 40:37–44. doi: 10.13954/j.cnki.hdu.2020.06.007

PubMed Abstract | CrossRef Full Text

19. Tang Z. Algorithm Research Based on Mutual Information and Demons Medical Image Registration [Dissertation]. Liaoning Normal University (2015).

20. Kim S, Tai Y-W. Hierarchical nonrigid model for 3D medical image registration. In: 2014 IEEE International Conference on Image Processing (ICIP). Paris: IEEE (2014). p. 3562–6.

Google Scholar

21. Marstal K, Berendsen F, Dekker N, Staring M, Klein S. The continuous registration challenge: evaluation-as-a-service for medical image registration algorithms. In: 2019 IEEE 16th International Symposium on Biomedical Imaging (ISBI 2019). Venice: IEEE (2019). p. 1399–402.

Google Scholar

22. Schaefer G, Tait R, Zhu SY. Overlay of thermal and visual medical images using skin detection and image registration. Conf Proc IEEE Eng Med Biol Soc. (2006) 2006:965–7. doi: 10.1109/IEMBS.2006.259275

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zana F, Klein JC. A multimodal registration algorithm of eye fundus images using vessels detection and Hough transform. IEEE Trans Med Imaging. (1999) 18:419–28. doi: 10.1109/42.774169

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Klein S, Staring M, Murphy K, Viergever MA, Pluim JPW. elastix: a Toolbox for Intensity-based medical image registration. IEEE Trans Med Imaging. (2010) 29:196–205. doi: 10.1109/TMI.2009.2035616

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Ait-Aoudia S, Mahiou R. Medical image registration by simulated annealing and genetic algorithms. In: Geometric Modeling and Imaging (GMAI '07). Zurich: IEEE (2007). p. 145–8.

Google Scholar

26. Tang K, Li Z, Tian L, Wang L, Zhu Y. ADMIR–affine and deformable medical image registration for drug-addicted brain images. IEEE Access. (2020) 8:70960–8. doi: 10.1109/ACCESS.2020.2986829

CrossRef Full Text | Google Scholar

27. Xinyu Z, Ziyi W, Guogang C, Ying C, Yijie W, Mengxue L, et al. Chaotic brain storm optimization algorithm in objective space for medical image registration. In: 2020 5th International Conference on Intelligent Informatics and Biomedical Sciences (ICIIBMS). Okinawa: IEEE (2020). p. 81–4.

Google Scholar

28. Zhang L, Zhou L, Li R, Wang X, Han B, Liao H. Cascaded feature warping network for unsupervised medical image registration. In: 2021 IEEE 18th International Symposium on Biomedical Imaging (ISBI). Nice: IEEE (2021). p. 913–6.

Google Scholar

29. Geng M, Meng X, Yu J, Zhu L, Jin L, Jiang Z, et al., Content-noise complementary learning for medical image denoising. IEEE Trans Med Imaging. (2022) 41:407–19. doi: 10.1109/TMI.2021.3113365

PubMed Abstract | CrossRef Full Text

30. Segars WP, Tsui BMW, Cai J, Yin FF, Fung GSK, Samei E. Application of the 4-D XCAT phantoms in biomedical imaging and beyond. IEEE Trans Med Imaging. (2018) 37:680–92. doi: 10.1109/TMI.2017.2738448

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Pang S, Pang C, Zhao L, Chen Y, Su Z, Zhou Y, et al. SpineParseNet: spine parsing for volumetric MR image by a two-stage segmentation framework with semantic image representation. IEEE Trans Med Imaging. (2021) 40:262–73. doi: 10.1109/TMI.2020.3025087

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Xiao F. A new divergence measure for belief functions in D-S evidence theory for multisensor data fusion. Inf Sci. (2020) 514:462–83. doi: 10.1016/j.ins.2019.11.022

CrossRef Full Text | Google Scholar

33. Cheng C, Xiao F. A distance for belief functions of orderable set. Pattern Recogn Lett. (2021) 145:165–70. doi: 10.1016/j.patrec.2021.02.010

CrossRef Full Text | Google Scholar

34. Xie D, Xiao F, Pedrycz W. Information quality for intuitionistic fuzzy values with its application in decision making. Eng Appl Artif Intell. (2021) 109:104568. doi: 10.1016/j.engappai.2021.104568

CrossRef Full Text | Google Scholar

35. Xiao F. CEQD: a complex mass function to predict interference effects. IEEE Trans Cybernet. (2021) 1–13. doi: 10.1109/TCYB.2020.3040770

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Xiao F. CaFtR: a fuzzy complex event processing method. Int J Fuzzy Syst. (2021) 1–14. doi: 10.1007/s40815-021-01118-6

CrossRef Full Text | Google Scholar

37. Xiao F. Generalization of Dempster–Shafer theory: a complex mass function. Appl Intell. (2020) 50:3266–75. doi: 10.1007/s10489-019-01617-y

CrossRef Full Text | Google Scholar

38. Zang L, Yang S, Zhang K, Shen S. Coupling gradient mutual information and PV interpolation image registration algorithm. J Surv Mapp Sci. (2021) 46–48:68–74. doi: 10.16251/j.cnki.1009-2307.2021.04.011

CrossRef Full Text

39. Yang J, Lu Z. A shape representation based on polar vector fourier descriptors. In: 2016 12th International Conference on Computational Intelligence and Security (CIS). Wuxi: IEEE (2016). p. 507–10.

Google Scholar

40. Dave V, Mahapatra C, Manchanda R. A mathematical model of the calcium transient in urinary bladder smooth muscle cells. Annu Int Conf IEEE Eng Med Biol Soc. (2015) 2015:5359–62. doi: 10.1109/EMBC.2015.7319602

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: PV image interpolation method, Shannon entropy, non-rigid medical image registration, DFP algorithm, CT-MRI-PET image

Citation: He L (2022) Non-rigid Multi-Modal Medical Image Registration Based on Improved Maximum Mutual Information PV Image Interpolation Method. Front. Public Health 10:863307. doi: 10.3389/fpubh.2022.863307

Received: 27 January 2022; Accepted: 21 March 2022;
Published: 01 June 2022.

Edited by:

Yuanpeng Zhang, Nantong University, China

Reviewed by:

Hongru Zhao, Soochow University, China
Fuyuan Xiao, Chongqing University, China
Yanhui Zhou, Southwest University, China

Copyright © 2022 He. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Liting He, bGl0aW5naGUxQG91dGxvb2suY29t

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.