<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>Hi Jan,</div><div><br></div><div>Thanks for sharing.</div><div>This does seem useful to me, but I'm not sure if I understand your method correctly.</div><div><br></div><div>For your lower and upper integration limit, do you mean the limit values for the angular range that you are "expecting"?<br>For example if you are expecting a 0-180 deg scan (although the first and last angles might not be 0 and 180 due to sampling), lower/upper integration limit would be 0 and 180 deg?</div><div><br></div><div>And why is the division 2 needed there?</div><div>I thought in <span style="font-family:arial,sans-serif;font-size:14px">rtkFDKWeightProjectionFilter.</span><span style="font-family:arial,sans-serif;font-size:14px">txx, the gap value used for the weighting is "nextAngle - previousAngle" for a certain projection.</span></div><div><span style="font-family:arial,sans-serif;font-size:14px">In this case I would expect Gap_first to be</span></div><div><span style="font-family:arial,sans-serif;font-size:14px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:14px">Gap_first = second_angle - lower_integration_limit</span></div><div><span style="font-family:arial,sans-serif;font-size:14px">(As the lower integration limit is kind of like the "virtual angle" </span><font face="arial, sans-serif"><span style="font-size:14px">preceding the first angle?)</span></font></div><div><font face="arial, sans-serif"><span style="font-size:14px"><br></span></font></div><div><font face="arial, sans-serif"><span style="font-size:14px">Thanks for your help :)</span></font></div><div><br></div><div>Cheers,</div><div>Andy</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
Date: Fri, 24 Oct 2014 17:21:27 +0200<br>
From: Jan Hoskovec <<a href="mailto:jean.hoskovec@gmail.com">jean.hoskovec@gmail.com</a>><br>
To: Andy Shieh <<a href="mailto:hsieandy@gmail.com">hsieandy@gmail.com</a>><br>
Cc: "<a href="mailto:rtk-users@public.kitware.com">rtk-users@public.kitware.com</a>" <<a href="mailto:rtk-users@public.kitware.com">rtk-users@public.kitware.com</a>><br>
Subject: Re: [Rtk-users] itFirstAngle and itLastAngle in<br>
        rtkParkerShortScanImageFilter.txt<br>
Message-ID:<br>
        <CANtP0QSnh70uETrdyTjg=u3HaUth4kRwDVfhMmKL=<a href="mailto:DhwrwzNLg@mail.gmail.com">DhwrwzNLg@mail.gmail.com</a>><br>
Content-Type: text/plain; charset=UTF-8<br>
<br>
Hi Andy,<br>
<br>
I was recently dealing with a similar problem in a different context<br>
(180? backprojection with irregular sampling and how to handle the<br>
first and last gaps) and what worked for me was<br>
<br>
Gap_first = (second angle - first angle) / 2 - lower integration limit<br>
<br>
and, analogically,<br>
<br>
Gap_last = upper integration limit - (last angle - second last angle) / 2<br>
<br>
with the integration limits being arbitrary set where I wanted them.<br>
The idea behind this was that a continuous projection value we are<br>
miming in the discrete integral should always be represented by the<br>
closest projection we have, with a known angular segment to cover.<br>
<br>
However, that was a DBP-type algorithm, for which the exact<br>
integration limits are extremely important, it may be different in the<br>
context of a short scan. But just in case you might find this<br>
useful...<br>
<br>
Cheers,<br>
<br>
Jan<br></blockquote></div></div></div>