<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:DengXian;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@DengXian";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.25in 1.0in 1.25in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US">
<div class="WordSection1">
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Hi all,</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I have used RTK for a while, first of all, thanks for providing this great software.
</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am trying to perform some iterative CBCT reconstruction with a specifically defined regularization term recently. Since this regularization is not TV nor wavelet, I may not be able to directly use those RTK provided filters, like ‘ConjugateGradientConeBeamReconstructionFilter’.
What I want is simply use gradient descent method and calculate the gradient of (1/2) ||(Rf - P)||^2 and my regularization term. I got some problems on how to get the gradient of (1/2) ||(Rf - P)||^2 , which should be R*(Rf - p), f is the required results,
p is measured projections and R is forward projections. Hope to get some hints for you.</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Fow now, I tried to directly <span style="color:black">use ‘JosephForwardProjectionImageFilter’ and ‘JosephBackProjectionImageFilter’ to achieve this goal, however, the acquired values seems to be really large after several iterations.
I think I must made something wrong. Could you please give me some suggestions about this problem if I want to utilize RTF to get gradient? Or is there any implementation detail of CG/SIRT</span> filters that we can access (I am not sure if it is public)?
I put my current code here (https://drive.google.com/file/d/1GPuoYHk4ACe5VDnVwZlsIzK1OHrM2AdK/view?usp=sharing)</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thank you very much and best regards.</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Zhehao</p>
</div>
</body>
</html>