[Rtk-users] using rtk in matlab : input images for filters from matlab data types

[Simon Rit]

Hi Julia,
I would have used
volume = py.GetImageFromArray(result)
instead of your last line. Does this work? If it does, then you have
to be sure that you set correctly the origin, spacing and direction
(the ITK meta information). One way to do this if you already have an
image with the correct meta information would be to use
py.volume.CopyInformation(volume)
Simon

On Thu, Mar 23, 2017 at 11:14 AM, Julia Semyakishkina
<prekrasnaya1985 at gmail.com> wrote:
> Hello.
> In matlab example from http://wiki.openrtk.org/index.php/SimpleRTK volume
> which is projected generates in SheppLoganPhantomFilter(). But i have my own
> volume, that is stored in matlab regular matrix for example.
>
> What i want to do:
> forward = py.JosephForwardProjectionImageFilter();
> forward.SetGeometry(g);
> projs = forward.Execute(projs, volume);
>
> The problem is to initialize the volume parameter.
> By this code
> data_size = size(matimg);
> transpose = matimg';
> result=py.numpy.reshape(transpose(:)', int32(data_size));
> volume = py.Image.fromarray(result);
>
> i store volume in type of py.Image, but in Execute() method i have to pass
> py.SimpleRTK.SimpleRTK.Image type. So, is there way to do such convertation
> ? Or is there some another way to init py.SimpleRTK.SimpleRTK.Image with my
> own matlab data ?
>
>
>
> _______________________________________________
> Rtk-users mailing list
> Rtk-users at public.kitware.com
> http://public.kitware.com/mailman/listinfo/rtk-users
>