[发明专利]基于二阶广义全变差的光声显微镜高分辨率图像重构方法

权利要求书

1.基于二阶广义全变差的光声显微镜高分辨率图像重构方法，其特征在于，包括：

S1、获取光声显微镜低分辨率图像；

S2、对所述光声显微镜低分辨率图像进行插值处理，得到插值图像；

S3、基于二阶广义全变差，确定优化所述插值图像的优化模型；

S4、基于所述优化模型，采用凸优化算法，重构出光声显微镜高分辨率图像。

2.根据权利要求1所述的方法，其特征在于，所述S2，包括：

利用双三次插值法对所述光声显微镜低分辨率图像进行插值，得到插值图像。

3.根据权利要求1所述的方法，其特征在于，在所述S3中，优化模型为：

argminIH12λ||IH-IH′||22+α1Σi=12||▿iIH-pi||+α0Σj=13||qi||]]>

s.t.q1=▿1p1,q2=▿2p2,q3=12(▿1p2+▿2p1)]]>

其中，I′_H为所述插值图像，α₀,α₁,λ为预设常数，p₁，p₂均为M×N维矩阵且初始化为零矩阵，其中，M，N根据I′_H的维数确定；

其中，I_H的初始值为所述I′_H，且

IH(t+1)=(λ(IH(t)+τ(▿1p1(t+1)+▿2p2(t+1)))+τIH′)/(λ+τ)]]>

其中，t为当前迭代次数，τ为预设常数；

(▿1p1(t+1))ij=-(p1(t+1))i-1j,i=M(p1(t+1))ij-(p1(t+1))i-1j,i=2,...,M-1(p1(t+1))ij,i=1;]]>

(▿2p2(t+1))ij=-(p2(t+1))ij-1,j=N(p2(t+1))ij-(p2(t+1))ij-1,j=2,...,N-1(p2(t+1))ij,j=1.]]>