abel数值反演的matlab实现,Abel逆变换的MATLAB实现（未完成版）

绪论

Abel逆变换常用于对于三维物体的重建中。Abel变换(Abel Transform)的意义为一个轴对称函数到向一个平面的投影，而Abel逆变换则是通过投影反演出原轴对称函数。其实际意义在于对某些轴对称的光源进行记录时可以通过平面记录器件还原出原有光源光强的三维分布，对于医学重建、火焰重建等有重要意义。

研究现状

Abel逆变换是一个典型的不适定问题，也就是说微小的数据误差在经过算法处理后，最终还原结果往往会出现巨大误差。这可以从Abel逆变换的定义式上理解一二。Abel逆变换定义式如下(摘自Wikipedia)

这个式子表示对于函数F进行Abel逆变换求得函数f(r)。同时在实际实验中，由于记录条件的限制，F往往是离散的，上式需要写成离散的格式，如下

其中上标从无穷变为n，是因为现实中的研究对象往往可认为在一定范围内，超出该范围则F值为0，这里将n看作这个范围，可以理解成CCD相机上该列值不为0的像素数目。可以看到其中存在对于F的差分运算，直接反演将产生可观的误差，且误差会被求和运算累积到无法接收的程度。因此Abel逆变换存在许多数值解法，包括：积分算子法、离散正则化方法、三次样条函数法、傅里叶-汉克尔变换方法和基函数展开(BASEX)法。在接下来的内容中，博主会对这些方法的基本思路进行介绍，尽量总结各自特点，并挑选其中较为常用的方法，给出其算法实现。

基本思路与特点

在阅读几篇文献后，对于以上数值解法有了一些了解，对其基本思路进行总结，并提出对于其特点(优缺点)的一些个人见解。对于算法的介绍按照我认为比较容易接收的逻辑进行。

三次样条函数法(插值\拟合方法)

由Abel变换的定义式可知，如果能得到像函数F的导数，则只需要进行简单的积分(求和)计算即可得到最终的结果。但是实际测量中像函数往往是离散值，因此导出了定义式的离散格式，其中微分运算用差分运算代替了，无疑带来了更大的误差。那么能否利用已知的离散像函数点还原出连续的像函数再进行求导呢？通过已知离散点还原连续函数，