如题，我用ITK自带的互信息配准，梯度下降优化器，速度太慢了，
用demos配准基本没什么效果，简单的位移图像都没法配过来，
可设置的参数过少，都不知道配不出来的原因（哭），
求一个C++的3D配准算法，最好基于ITK，或者自己写的库，谢谢
（如有大神指导ITK编写，感激不尽）
需求：
1 C++编写的3D图像配准程序，用于CT和PET图像
2 itk中的registration组件中的优化算法是否可以用powell等，如果是，该怎么用
3 matlab中的imregister函数互信息测度，移植到C++程序的方法也行
