Volume constrained 2-phase segmentation method utilizing a linear system solver based on the best uniform polynomial approximation of

详细信息    查看全文

• 作者：Stanislav Harizanov ; ^{sharizanov@parallel.bas.bg" class="auth_mail" title="E-mail the corresponding author} ; ; Svetozar Margenov ^{margenov@parallel.bas.bg" class="auth_mail" title="E-mail the corresponding author} ; Pencho Marinov ^{pencho@bas.bg" class="auth_mail" title="E-mail the corresponding author} ; Yavor Vutov ^{yavor@parallel.bas.bg" class="auth_mail" title="E-mail the corresponding author}
• 关键词：65F10 ; 41A10 ; 68W10 ; 68U10
• 刊名：Journal of Computational and Applied Mathematics
• 出版年：2017
• 出版时间：15 January 2017
• 年：2017
• 卷：310
• 期：Complete
• 页码：115-128
• 全文大小：1628 K

文摘

Volume constrained image segmentation aims at improving the quality of image reconstruction via incorporating physical information for the underline object of interest into the mathematical modeling of the segmentation problem. In this paper, we develop a general framework for 3D 2-phase image segmentation, based on constrained class="mathmlsrc">class="formulatext stixSupport mathImg" data-mathURL="/science?_ob=MathURL&_method=retrieve&_eid=1-s2.0-S0377042716302928&_mathId=si45.gif&_user=111111111&_pii=S0377042716302928&_rdoc=1&_issn=03770427&md5=4c1cd42b0e0dbb2387784bbb6957b15c" title="Click to view the MathML source">ℓ²class="mathContainer hidden">class="mathCode">${\ell }^2$ minimization of a non-local regularizer, the Euler–Lagrange derivative of which is the discrete graph-Laplacian of a weighted graph, associated with the image voxels. It involves a convenient change of basis in the image domain, for which the optimization function is decomposed element-wise. Using univariate polynomial approximation techniques, we show that the transformation matrix does not need to be explicitly computed and its action is well approximated by a suitable matrix polynomial. The error is independent of the domain size, thus our approach is applicable to high resolution data. The model allows for adding arbitrary linear terms into the optimization function in order to increase the control on the output, in particular to add its scalar product with another, already known, segmentation vector. Such a “hybrid” process may significantly improve on the individual quality of each of the involved segmentations.