PSO和Powell混合算法在医学图像配准中的应用研究

基于互信息的图像配准方法具有自动化程度高、配准精度高等优点,已被广泛应用于医学图像的配准.但是,基于互信息的目标函数经常是不光滑的,存在许多局部极值,给问题的求解带来了很大的困难.本文讨论了互信息函数的多极值特性,并提出了一种粒子群优化算法(particle swarm optimization,PSO)和Powell混合优化方法.经检验,这种方法能有效地克服互信息函数的局部极值,大大地提高了配准精度,达到亚像素级.     

A hybrid framework for registration of carotid ultrasound images combining iconic and geometric features

Stroke is the third major cause of death worldwide behind heart disease and cancer. Carotid atherosclerosis is the most frequent cause of ischemic stroke. Early diagnosis of carotid plaque and serial monitoring of its size with the help of imaging modalities can help to prevent the atherosclerotic complications. The main difficulty is inevitable variability of patient’s head positions during image acquisitions. The time series registration of carotid images helps to improve the monitoring, characterization, and quantification of the disease by suppressing the global movements of the patient. In this work, a novel hybrid registration technique has been proposed and evaluated for registration of carotid ultrasound images taken at different times. The proposed hybrid method bridges the gap between the feature-based and intensity-based registration methods. The feature-based iterative closest point algorithm is used to provide a coarse registration which is subsequently refined by the intensity-based algorithm. The proposed framework uses rigid transformation model coupled with mutual information (MI) similarity measure and Powell optimizer. For quantitative evaluation, different registration approaches have been compared using four error metrics: visual information fidelity, structural similarity index, correlation coefficient, and MI. Qualitative evaluation has also been done using the visual examination of the registered image pairs.     

混合遗传算法在医学图像配准中的应用

针对遗传算法和梯度算法各自的优点和缺陷,提出了遗传算法和梯度算法相结合的松散耦合和紧密耦合两种混合遗传算法,并利用基于排挤的小生境技术对其性能进行了改进,使其在搜索能力和收敛速度上都有了较大提高。利用此算法成功实现了以互信息为相似性测度的医学图像配准。计算实例结果表明,本算法能够有效克服局部极值的影响,具备较强的抗噪能力,是一个精度高、稳定性好的医学图像配准算法。     

基于Contourlet变换的CT和锥形束CT图像配准算法

目的提出一种基于Contourlet变换,用于放射治疗定位的CT与锥形束CT（cone beam CT,CBCT）图像配准的方法。方法利用Contourlet变换多尺度多方向的分辨特性,将待配准图像进行Contourlet变换分解,分解后的高频方向子带合成梯度图像,采用归一化互信息作为相似性测度,把梯度图像与低频方向子带以加权函数结合,进行临床医学图像的刚性配准,有效弥补了互信息配准中缺少空间信息的不足。结果通过已知空间变换参数图像的配准结果验证了算法的准确性。配准后lO幅图像变换参数的误差极小,且均方根误差接近于0。结论该图像配准算法精确度高,并具有很好的鲁棒性,有助于提高图像引导放射治疗（image guid edradiation therapy,IGRT）中解剖组织结构和靶区的定位精度。     

基于混合优化算法的医学图像配准方法

基于互信息的图像配准算法计算复杂度高,配准速度慢。针对这一问题,本文提出一种基于改进遗传算法和Powell算法相结合的医学图像配准方法。首先针对传统遗传算法收敛速度慢、易早熟的缺陷,本文对遗传操作中的交叉运算过程提出了改进策略,并将改进的遗传算法与Powell算法相结合,充分利用遗传算法的全局搜索能力与Powell算法的局部搜索能力。与Powell算法和未改进的遗传算法相比,本文提出的算法极大地缩短了图像配准所用的时间,同时提高了算法的抗噪性。     

基于调幅-调频与熵图多模态医学图像的配准

背景：基于传统互信息量的多模态医学图像配准方法配准时需要利用二维直方图或者Parzen窗函数的方法估计概率密度分布,进而计算互信息量,这种方式计算速度慢,而且只考虑了图像的灰度信息,容易出现误配。 目的：针对目前主流的配准方法鲁棒性差、耗时的缺点,提出了一种新的基于调幅-调频(AM-FM)特征互信息量的快速配准方法。 方法：该方法考虑了图像的空间和结构信息;首先通过AM-FM模型对图像进行分解,得到图像的AM-FM特征,与图像的灰度特征一起组成高维特征;然后利用熵图和最小生成树加快AM-FM特征互信息量的计算,从而实现了医学图像的快速配准。 结果与结论：对20组磁共振T1-T2加权图像、CT/正电子发射计算机断层成像图像进行了实验,结果表明该方法在图像空间分辨率较低,有噪声影响等情况下均可以达到较好的结果,且配准精度优于国际上的主流方法,具有计算速度快,精度高,鲁棒性强的特点,适于临床应用。     

基于形态学梯度和互信息的医学图像配准方法

基于互信息的图像配准方法,已被广泛用于医学图像的配准.但是该方法计算量较大且无法处理图像空间信息,导致运行时间较长且易陷入局部极值.为解决此问题,本研究提出了一种基于形态学梯度和互信息相结合的医学图像配准新方法,该方法充分利用图像的灰度信息和空间几何形状,可节省运行时间且有效改善传统互信息方法中的局部极值问题.实验结果表明,该方法的配准精度和运行速度明显优于传统方法.     

引入高斯函数的互信息法多模态图像配准

目的：最大互信息作为相似度测量在医学图像配准中已被广泛应用。在计算图像互信息时,为了避免引入新的灰度值一般采用部分体积插值统计两幅图像的联合直方图。但用该方法计算中,当图像平移整数点时,统计联合直方图会出现缺陷,使目标函数出现局部极值,从而造成误配准。方法：将高斯函数引入到直方图统计中,选取适当的邻域,用高斯函数计算邻域内各点像素对联合直方图的贡献。利用高斯函数的平滑性,避免了在互信息计算过程中统计图像联合直方图时出现误差。使用Powell优化方法,寻找最佳的优化参数,实现图像的最佳配准。结果：采用CT-PET数据进行实验,该方法平滑了目标函数,有效地消除了局部极值,提高了多模态图像配准的精确性,并且,对噪音图像配准也产生很好的效果。结论：该方法适用于多模态医学图像配准,克服了传统互信息计算时的不足,提高了配准的正确率和精确度。     

基于互信息与边缘梯度相结合的多模医学图像配准方法

针对传统互信息配准方法未利用图像空间信息的缺点,本文研究了图像边缘信息的梯度相似性.首先采用小波模极大值边缘检测提取出图像边缘,提出将边缘图像的梯度相似性系数与传统的互信息相乘作为图像配准的目标函数.然后通过使用Powell优化算法对目标函数进行寻优,得出配准变换参数.最后在互信息的基础上引入图像边缘梯度信息,突出了全局最优解.实验结果表明,该方法可以得到精确、有效的配准结果.     

一种基于样本抽样性质的图像配准方法

基于互信息的医学图像配准是当前常用的方法,但互信息的计算量大,对此可采用欠采样技术来减少计算量,但欠采样会降低配准的精度。本文中,我们基于样本抽样均值分布定理,提出了利用多次采样的改进方法,它兼具较高的精度和较快的配准速度。文中的数值计算结果证明了这一点。     

基于伪超声和互信息的磁共振与超声图像配准

手术中超声图像与术前磁共振图像的配准在手术导航系统中非常重要。我们利用磁共振和超声成像原理,提出了基于伪超声和互信息,并结合多分辨率技术与Powell优化算法对两种模态图像进行配准的方法,该方法可以有效降低寻优过程中陷入局部极值收敛的概率,提高两种模态图像的配准精度。实验结果表明,我们提出的基于伪超声和互信息的配准方法比目前手术导航系统中普遍采用的标记点方法具有更高的配准精度。     

基于ITK多图谱配准的脑部皮质下核团的分割研究

大量研究表明,阿尔茨海默症(AD)的病变与大脑皮质下核团的萎缩息息相关,某些核团的萎缩(如海马)可能成为AD疾病早期诊断的标志,而皮质下核团的分割是研究核团萎缩模式的重要前提。基于AD患者和正常人各30例3DT1W-MR图像,先结合直方图分析和三维形态学分析方法对图像进行脑组织提取,后采用ITK配准算法将10个脑图谱图像经两阶段分别配准到提取脑组织后的图像空间。第一阶段实现基于均方差的仿射配准,第二阶段实现基于互信息的B样条形变配准,两阶段的配准均采用线性插值法和梯度下降的优化搜索方法。最后采用STAPLE融合算法,对配准后得到的10个目标图像进行图像融合,得到最终的分割结果。结果表明:除尾状核外,分割得到的其余6对核团的体积与常用的FSL-FIRST算法的分割结果无统计学差别(P>0.05);AD患者的右侧伏核和双侧海马发生萎缩(P<0.05)。因此,基于ITK配准框架的多图谱配准分割方法能有效分割MR图像上边界不明确的皮质下核团。     

基于互信息的医学图像配准技术研究进展

基于互信息的图像配准方法具有自动化程度高、配准精度高等优点,已被广泛应用于医学图像的配准。详细地介绍了互信息的基本概念、归一化互信息、互信息计算中常用的变换模型、插值方法、优化算法及提高配准速度策略等。同时对互信息在配准技术中存在的问题及解决方法进行了系统地介绍,最后对今后互信息的研究工作作了展望。     

基于互信息的医学图像配准中互信息的计算

基于互信息的配准方法是医学图像配准领域的重要方法.互信息是图像配准中常用的相似性度量,具有鲁棒、精度高等优点,但基于互信息的配准计算量大,制约了它的实际应用.我们采用基于多分辨率和混合优化策略的配准方法,在图像的不同灰度等级数下进行配准,分析了互信息的计算量与灰度等级数的关系,并用人头部的MRI图像和CT图像做了二维的单模模拟实验和多模实际配准实验,结果显示在灰度等级数为32和64时,与灰度等级数为256时相比,配准精度没有明显改变,而计算量下降显著.     

Development and validation of a CT-3D rotational angiography registration method for AVM radiosurgery

In this paper a novel technique is proposed and validated for radiosurgery treatment planning of arteriovenous malformations (AVMs). The technique was developed for frameless radiosurgery by means of the CyberKnife, a nonisocentric, linac-based system which allows highly conformed isodose surfaces to be obtained, while also being valid for other treatment strategies. The technique is based on registration between computed tomography (CT) and three-dimensional rotational angiography (3DRA). Tests were initially performed on the effectiveness of the correction method for distortion offered by the angiographic system. These results determined the registration technique that was ultimately chosen. For CT-3DRA registration, a twelve-parameter affine transformation was selected, based on a mutual information maximization algorithm. The robustness of the algorithm was tested by attempting to register data sets increasingly distant from each other, both in translation and rotation. Registration accuracy was estimated by means of the "full circle consistency test." A registration quality index (expressed in millimeters) based on these results was also defined. A hybrid subtraction between CT and 3DRA is proposed in order to improve 3D reconstruction. Preprocessing improved the ability of the algorithm to find an acceptable solution to the registration process. The robustness tests showed that data sets must be manually prealigned within approximately 15 mm and 20 degrees with respect to all three directions simultaneously. Results of the consistency test showed agreement between the quality index and registration accuracy stated by visual inspection in 20 good and 10 artificially worsened registration processes. The quality index showed values smaller than the maximum voxel size (mean 0.8 mm compared to 2 mm) for all successful registrations, while it resulted in much greater values (mean 20 mm) for unsuccessful registrations. Once registered, the two data sets can be used for CyberKnife treatment planning. Target delineation is performed on 3DRA while dose calculation and DRR generation are performed on CT. In conclusion, a method was developed for using 3DRA images for AVM frameless radiosurgery treatment planning. The method proved to be feasible, robust, and accurate for clinical use. 3DRA can be performed at different times or locations compared to standard, frame based stereotactic angiography. Unlike two-dimensional angiography, 3DRA allows examination of the shape of the AVM and of the surrounding target from any arbitrary point of view during treatment planning. The method can be applied to any case of intermodality registration, is operator-independent, and allows estimation of registration quality. Further research is desirable to improve time resolution in order to distinguish between feeding and draining vessels.     

Targeted 2D/3D registration using ray normalization and a hybrid optimizer

X-ray images are often used to guide minimally invasive procedures in interventional radiology. The use of a preoperatively obtained 3D volume can enhance the visualization needed for guiding catheters and other surgical devices. However, for intraoperative usefulness, the 3D dataset needs to be registered to the 2D x-ray images of the patient. We investigated the effect of targeting subvolumes of interest in the 3D datasets and registering the projections with C-arm x-ray images. We developed an intensity-based 2D/3D rigid-body registration using a Monte Carlo-based hybrid algorithm as the optimizer, using a single view for registration. Pattern intensity (PI) and mutual information (MI) were two metrics tested. We used normalization of the rays to address the problems due to truncation in 3D necessary for targeting. We tested the algorithm on a C-arm x-ray image of a pig's head and a 3D dataset reconstructed from multiple views of the C-arm. PI and MI were comparable in performance. For two subvolumes starting with a set of initial poses from +/-15 mm in x, from +/-3 mm (random), in y and z and +/-4 deg in the three angles, the robustness was 94% for PI and 91% for MI, with accuracy of 2.4 mm (PI) and 2.6 mm (MI), using the hybrid algorithm. The hybrid optimizer, when compared with a standard Powell's direction set method, increased the robustness from 59% (Powell) to 94% (hybrid). Another set of 50 random initial conditions from [+/-20] mm in x,y,z and [+/-10] deg in the three angles, yielded robustness of 84% (hybrid) versus 38% (Powell) using PI as metric, with accuracies 2.1 mm (hybrid) versus 2.0 mm (Powell).     

CT/MRI混合配准方法及其在放疗计划系统中的应用

旨在研究放疗中图像配准方法,特别是针对放疗中常用的CT、MRI,提出基于混合框架的配准方法,该方法主要包 括两个方面:（1）采用掩膜（Mask）提取感兴趣区域、形态学运算等图像处理方法以及CPU多线程并行技术,大幅度提高配 准速度;（2）采用由全局到局部的混合配准策略,首先利用基于仿射变换的刚性配准整体配准,以防止图像间偏差过大,在 此基础上针对感兴趣区域采用B样条弹性配准,调整局部形变。通过实验表明,采用预处理及加速策略的刚性配准,在保 持其精度的情况下,提速比可达10倍,测试结果已达到临床需求;此外,采用基于GPU加速的混合配准策略,其配准速度 提至约4 min。     

基于互信息医学刚性图像配准方法的研究

目的:基于互信息的配准方法是医学图像配准领域的重要方法,具有鲁棒性,精度高等优点。本文探究医学刚性图像配准的有效算法和关键技术。方法:基于互信息配准方法,利用Powell多参数算法和改进的PV插值算法,得到两幅图像之间的最大互信息和最佳配准参数。结果:二维磨牙CT图像配准实验表明,配准速度快,精度提高,验证了插值方法的有效性。结论:方法和算法可提高配准速度,能有效抑制互信息目标函数的局部极值。     

A study on medical image registration by mutual information with pyramid data structure

The registration method based on mutual information is currently a popular technique for medical image registration, but the computation of the mutual information is complex and the registration speed is slow. In this work, a new slice accumulation pyramid (SAP) data structure was proposed to expedite the registration process. A numerical comparative study between the new data structure and the existing wavelet pyramid (WP) data structure was given, and the results confirmed that the new pyramid data structure was superior to the WP in both the calculation efficiency and the optimizing performance. Finally, SAP was applied to remove the artifacts between CT and MRI data sets, and the results showed the validation of SAP to registration of mulmodality images.     

A computationally efficient method for automatic registration of orthogonal x-ray images with volumetric CT data

The paper presents a computationally efficient 3D-2D image registration algorithm for automatic pre-treatment validation in radiotherapy. The novel aspects of the algorithm include (a) a hybrid cost function based on partial digitally reconstructed radiographs (DRRs) generated along projected anatomical contours and a level set term for similarity measurement; and (b) a fast search method based on parabola fitting and sensitivity-based search order. Using CT and orthogonal x-ray images from a skull and a pelvis phantom, the proposed algorithm is compared with the conventional ray-casting full DRR based registration method. Not only is the algorithm shown to be computationally more efficient with registration time being reduced by a factor of 8, but also the algorithm is shown to offer 50% higher capture range allowing the initial patient displacement up to 15 mm (measured by mean target registration error). For the simulated data, high registration accuracy with average errors of 0.53 mm +/- 0.12 mm for translation and 0.61 +/- 0.29 degrees for rotation within the capture range has been achieved. For the tested phantom data, the algorithm has also shown to be robust without being affected by artificial markers in the image.