<div dir="ltr">Hi Jan and Simon,<div><br></div><div>Your suggestions and comments are very inspiring and helpful. Thank you for your help :)</div><div><br></div><div>Cheers,</div><div>Andy</div></div><div class="gmail_extra"><br><div class="gmail_quote">2014-10-28 0:34 GMT+11:00 Jan Hoskovec <span dir="ltr"><<a href="mailto:jean.hoskovec@gmail.com" target="_blank">jean.hoskovec@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Andy,<br>
<br>
I'm sorry not to have replied earlier. For the precisions you wanted,<br>
<br>
1) your intuition is right, the upper and lower integration limits are<br>
the values you are "expecting", the values delimiting the angular<br>
range you want to cover (whatever is the actual sampling).<br>
<br>
2) The division by two is there because in my approach, the "zone of<br>
influence" of each sampled projections begins and ends halfway between<br>
the sampled value of the gantry angle and the next / preceding sample.<br>
The weights of other than first and last samples simplify to " (next<br>
angle - previous angle) / 2" for me.<br>
<br>
Hope this would help (unless, of course, you've found a better option<br>
during the weekend :-) ).<br>
<br>
Cheers,<br>
<br>
Jan<br>
<div class="HOEnZb"><div class="h5"><br>
2014-10-25 15:41 GMT+02:00 Andy Shieh <<a href="mailto:hsieandy@gmail.com">hsieandy@gmail.com</a>>:<br>
> Hi Jan,<br>
><br>
> Thanks for sharing.<br>
> This does seem useful to me, but I'm not sure if I understand your method<br>
> correctly.<br>
><br>
> For your lower and upper integration limit, do you mean the limit values for<br>
> the angular range that you are "expecting"?<br>
> For example if you are expecting a 0-180 deg scan (although the first and<br>
> last angles might not be 0 and 180 due to sampling), lower/upper integration<br>
> limit would be 0 and 180 deg?<br>
><br>
> And why is the division 2 needed there?<br>
> I thought in rtkFDKWeightProjectionFilter.txx, the gap value used for the<br>
> weighting is "nextAngle - previousAngle" for a certain projection.<br>
> In this case I would expect Gap_first to be<br>
><br>
> Gap_first = second_angle - lower_integration_limit<br>
> (As the lower integration limit is kind of like the "virtual angle"<br>
> preceding the first angle?)<br>
><br>
> Thanks for your help :)<br>
><br>
> Cheers,<br>
> Andy<br>
><br>
><br>
>><br>
>> 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>
>><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>
><br>
><br>
</div></div><div class="HOEnZb"><div class="h5">> _______________________________________________<br>
> Rtk-users mailing list<br>
> <a href="mailto:Rtk-users@public.kitware.com">Rtk-users@public.kitware.com</a><br>
> <a href="http://public.kitware.com/mailman/listinfo/rtk-users" target="_blank">http://public.kitware.com/mailman/listinfo/rtk-users</a><br>
><br>
</div></div></blockquote></div><br></div>