[Rtk-users] data loss in volume when reconstructing a anisotropic voxel

[Simon Rit]

Hi,
Can you clarify what is not working with the mailing list? Are you really
using rtk-users at openrtk.org (cc)?
I don't see anything wrong with this small piece of code. You should share
the full two versions to understand what could be wrong.
Regards,
Simon

On Fri, Apr 7, 2023 at 6:13 AM BALADASTAGIRI ROOPANAGUDI <
dastagiri051428 at gmail.com> wrote:

> Hi simon,
>
> Sorry for inconvenience, yesterday i tried to reach out to you through
> mailing list but i couldn't reach , could you please help me out to figure
> out this issue
>
> When i tried to reconstruct a anisotropic voxel there is loss of data in
> the top & bottom portion of volume and  the same doesnt happen if i try to
> reconstruct it with equal spacing out in all three direction
>
>    ConstantImageSourceType::Pointer outVol=ConstantImageSourceType::New();
>     ConstantImageSourceType::PointType origin_out;
>     ConstantImageSourceType::SizeType size_out;
>     ConstantImageSourceType::SpacingType spacing_out;
>     size_out[0] = ctDIM1;// 512
>     size_out[1] = ctDIM3;//120
>     size_out[2] = ctDIM2;//512
>     spacing_out[0] = 0.4238;
>     spacing_out[1] = 1;
>     spacing_out[2] = 0.4238;
>
>     origin_out[0] = -0.5*(size_out[0]-1)*(spacing_out[0]);
>     origin_out[1] = -0.5*(size_out[1]-1)*(spacing_out[1]);
>     origin_out[2] = -0.5*(size_out[2]-1)*(spacing_out[2]);
>     outVol->SetOrigin(origin_out);
>     outVol->SetSpacing(spacing_out);
>     outVol->SetSize(size_out);
>     outVol->SetConstant(0);
>
> I am not sure where i am going wrong,
>
> I have attached two images with same pixelspacing and variable
> pixelspacing for your references,
>
> With Regards,
> Dastagiri R
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.creatis.insa-lyon.fr/pipermail/rtk-users/attachments/20230407/8a88711d/attachment.htm>