<p>One of the main advantages of the widely applied <InlineEquation ID="IEq1"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="12190_2025_2546_Article_IEq1.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(l_{1}\)</EquationSource> </InlineEquation>-norm minimization is the great reduction in streak artifacts for computed tomography image reconstruction with limited-angle projections. However, the results are over smoothing, not sparse enough and lack of some details, causing edge blurred. To address these issues, we proposed a combined <InlineEquation ID="IEq2"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="12190_2025_2546_Article_IEq1.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(l_{1}\)</EquationSource> </InlineEquation>-norm and <InlineEquation ID="IEq3"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="12190_2025_2546_Article_IEq3.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(l_{0}\)</EquationSource> </InlineEquation>-norm regularization model and developed the NADA-HT algorithm to solve the optimization effectively using the nonmonotone alternating direction algorithm with a hard thresholding. In this paper, the convergence of the algorithm is analyzed and the weights between the <InlineEquation ID="IEq4"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="12190_2025_2546_Article_IEq1.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(l_{1}\)</EquationSource> </InlineEquation>-norm and <InlineEquation ID="IEq5"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="12190_2025_2546_Article_IEq3.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(l_{0}\)</EquationSource> </InlineEquation>-norm terms in the model are discussed. A lower bound of decrements of the objective function is derived in terms of the difference of the reconstructed images for sufficiently large iteration numbers under certain conditions. It is further shown that there exists a subsequence of the reconstructed image sequence converging to a local minimizer. Numerical experiments demonstrate the reconstruction results by using different weights of the <InlineEquation ID="IEq6"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="12190_2025_2546_Article_IEq3.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(l_{0}\)</EquationSource> </InlineEquation>-norm term and lead to discussions on parameter selection.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Convergence of NADA-HT algorithm for CT image reconstruction

  • Jiehua Zhu,
  • Xiezhang Li,
  • Guocan Feng

摘要

One of the main advantages of the widely applied \(l_{1}\) -norm minimization is the great reduction in streak artifacts for computed tomography image reconstruction with limited-angle projections. However, the results are over smoothing, not sparse enough and lack of some details, causing edge blurred. To address these issues, we proposed a combined \(l_{1}\) -norm and \(l_{0}\) -norm regularization model and developed the NADA-HT algorithm to solve the optimization effectively using the nonmonotone alternating direction algorithm with a hard thresholding. In this paper, the convergence of the algorithm is analyzed and the weights between the \(l_{1}\) -norm and \(l_{0}\) -norm terms in the model are discussed. A lower bound of decrements of the objective function is derived in terms of the difference of the reconstructed images for sufficiently large iteration numbers under certain conditions. It is further shown that there exists a subsequence of the reconstructed image sequence converging to a local minimizer. Numerical experiments demonstrate the reconstruction results by using different weights of the \(l_{0}\) -norm term and lead to discussions on parameter selection.