专利名称：一种非完全数据下的光学层析重建计算方法
技术领域：
本发明涉及一种光学层析重建计算方法，尤其涉及一种非完全数据下的光学层析 重建计算方法。
背景技术：
光学层析是一种重建流场三维热物理参数分布的检测技术。当待测流场中含有障 碍物、不透明物体或由于测试环境条件的限制，会引起投影数据的缺失。如果仅用已知投影 数据进行重建运算，得到的重建结果将与实际分布相差很大。医学层析检测中，曾提出过迭 代重建重投影及投影空间迭代重建重投影计算方法，但这两种方法都需要多方向下的大量 已知投影数据进行运算，因此并不适合以少投影方向数为特点的光学层析检测。本方法针对非完全数据下的光学层析重建计算的问题，提出一种基于拉格朗日插 值预估、重投影修正的光学层析计算方法，该方法首先对每个投影方向，采用拉格朗日插值 多项式进行缺失投影数据的预估计，得到一组缺失投影数据估计值。然后运用一种逐线乘 性迭代预重建算法得到一个粗重建图像，并对粗重建图像进行逐射线重投影，得到缺失投 影数据的另一组估计值。接着对二组估计值进行逐个数据对比，并对相差较大的估计值进 行权重叠加修正。重建与修正计算同时进行，最终得到所有缺失投影数据的无偏估计值及 精确的重建结果。本方法适用于由于待测流场中含有障碍物或测试条件的限制等原因引起 的数据缺失情况下的光学层析检测。

发明内容
本发明的目的在于提供一种非完全数据下的光学层析重建计算方法，该方法适用 于由于待测流场中含有障碍物或测试条件的限制等原因引起的数据缺失情况下的光学层 析检测。本发明是这样来实现的，其特征是方法步骤为1)对每个投影方向，基于已知投影数据，采用拉格朗日插值多项式进行缺失投影 数据的预估计，得到一组缺失投影数据估计值；2)运用一种逐线乘性迭代预重建算法得到一个粗重建图像，计算中只用到已知投 影数据；3)对粗重建图像进行逐射线重投影，得到缺失投影数据的另一组估计值；4)对二组估计值进行逐个数据对比，如果两者相差小于5%，则认为该数据为无 偏估计，并归入已知投影数据参与重建运算；如两者相差大于5%，则用它们的权重叠加来 修正该值并参与重建运算，随着重建迭代的进行，重投影估计值的权重逐渐增加；5)重建与修正同时进行，最终得到所有缺失投影数据的无偏估计值及精确的重建 结果。本发明的优点是该方法结合拉格朗日插值预估与重投影修正，并且将重建与修 正过程同时进行，以得到缺失数据的无偏估计值并解决非完全数据下的光学层析重建问
具体实施例方式本发明的具体计算方法及原理为当待测流场中含有障碍物、不透明物体或由于测试环境条件的限制，会引起投影 数据的缺失，此时投影数据分为已知投影数据Pk及缺失投影数据pu。1)假设竖直坐标用P表示，水平坐标用t表示。首先找出一个通过所有的已知投 影数据(tk，Pk)(注k= 1，2，3，...，Nk. Nk代表某一方向的已知投影数)的拉格朗日插值 多项式P = L(t)
Γ η ",、^ (t 一U){t-t2)...{t-D -办-tNk)P = L(t) = ^tPk----------
i=i (h -t^ih ~t2)-(tk-h-^ih-h+O-ih-tNt)这样可根据该Nk-I次拉格朗日插值多项式得到缺失投影数据Pu的一组估计值
La (Pu)。2)运用一种逐线乘性迭代预重建算法得到一个粗重建图像fK，计算中只用到已知 投影数据Pk η lf W'r * ° and Pk式中，i = 1，2，···，I ;I为投影数据总数；j = 1，2，...，丽；丽为总的重建点数；
Wi是投影矩阵W的第i行.3)做循环运算i = 1，2，. . .，I，对粗重建图像fK进行逐射线重投影，当i等于Pu 的标号时Re (pu) = Wi^式中，算符Re代表重投影估计，这样就得到了缺失数据Pu的第二组估计Re(Pu)4)对二组估计值进行逐个数据对比，如果I La (Pu) -Re (pu) | < 0. 05Re (pu)则认为该数据为无偏估计，并归入已知投影数据参与重建运算；如两者相差大于 5%，则La (Pu) = (1 _ α ) La (pu) + α Re (pu)即用它们的权重叠加来修正该值并参与重建运算，式中，α为重投影估计值的权 重因子，它随着重建迭代的进行逐渐增加.5)重建与修正同时进行，最终得到所有缺失投影数据的无偏估计值及精确的重建结果。
权利要求
一种非完全数据下的光学层析重建计算方法，其特征是方法步骤为1)对每个投影方向，基于已知投影数据，采用拉格朗日插值多项式进行缺失投影数据的预估计，得到一组缺失投影数据估计值；2)运用一种逐线乘性迭代预重建算法得到一个粗重建图像，计算中只用到已知投影数据；3)对粗重建图像进行逐射线重投影，得到缺失投影数据的另一组估计值；4)对二组估计值进行逐个数据对比，如果两者相差小于5％，则认为该数据为无偏估计，并归入已知投影数据参与重建运算；如两者相差大于5％，则用它们的权重叠加来修正该值并参与重建运算，随着重建迭代的进行，重投影估计值的权重逐渐增加；5)重建与修正同时进行，最终得到所有缺失投影数据的无偏估计值及精确的重建结果。
全文摘要
一种非完全数据下的光学层析重建计算方法，其特征是方法步骤为1)对每个投影方向，基于已知投影数据，得到一组缺失投影数据估计值；2)计算中只用到已知投影数据；3)得到缺失投影数据的另一组估计值；4)重投影估计值的权重逐渐增加；5)重建与修正同时进行，最终得到所有缺失投影数据的无偏估计值及精确的重建结果。本发明的优点是该方法结合拉格朗日插值预估与重投影修正，并且将重建与修正过程同时进行，以得到缺失数据的无偏估计值并解决非完全数据下的光学层析重建问题。
文档编号G01N21/85GK101923056SQ20101023160
公开日2010年12月22日 申请日期2010年7月20日 优先权日2010年7月20日
发明者万生鹏, 万雄, 伏燕军, 张志敏, 易江林, 高益庆 申请人:南昌航空大学