[发明专利]一种基于虚拟数据的轴对称参数分布图像重建方法

权利要求书

1.一种基于虚拟数据的轴对称参数分布图像重建方法，该方法包括以下步骤：

步骤一，获取单一视角的投影测量数据：利用单一视角投影覆盖被测参数分布场，可以是扇形束投影或平行束投影，得到m个投影测量数据；

步骤二，构建虚拟测量数据：共有n个视角投影覆盖被测参数分布场，n为奇数，其中一个视角的投影为真实投影，其他n-1个视角的投影为虚拟投影，每个视角的投影所覆盖的投影角度相同，使n-1个视角的虚拟投影所得到的虚拟投影测量值与m个真实投影测量值相同，总测量数据的个数为J，J＝m×n；

步骤三，根据系数矩阵和投影测量值构建求解参数值的方程组：将被测区域划分为N个网格，对于第j束投影，投影测量值b_j可以表示为：

bj=Σi=1NAjixi(j=1,2,···J)---(1)]]>

其中x_i为第i个像素中的参数值，A_ji表示第j束投影在第i个网格中的路径长度，则：

b＝Ax    (2)

其中A为A_ji构成的维数为J×N的系数矩阵，x为构成x_i的向量，b为构成b_j的向量；

步骤四，图像重建：根据代数重建方法求解(2)，得到x的表达式：

x(k+1)=x(k)+λbjk-AjkTx(k)||Ajk||2Ajk]]>(3)

jk=[k-Int(kJ)J+1]]]>

其中k为迭代次数，λ为松弛因子，Int代表取整算子，迭代的收敛性由相邻两次迭代计算得到x^(k)和x^(k-1)差的模相对于x^(k)模的变化Δ来监控，Δ可以表示为：

Δ=|x(k)-x(k-1)||x(k)|---(4)]]>

当Δ≤0.1%时停止迭代，根据求得的x^(k)可以重建轴对称参数分布场的图像。