[ITK] (beginner question) How to get transform matrix then apply it to another volume

[Dženan Zukić]

Hi Neo,

it sounds like you did not change the OutputPixelType from unsigned char to
short, which is required to represent the entire range of CT numbers (-1000
to +3000).

Regards,
Dženan

On Sat, Aug 20, 2016 at 12:37 AM, Xiaodong Zhao <xzhao21 at lsu.edu> wrote:

> Hi Dzenan,
>
> Thanks for your reply!
>
> I tried to edit the example ImageRegistration20.cxx a bit. Basically I
> just input the the fixed MR, the moving MR, and the CT that’s
> pre-registered with the moving MR, and after the affine transformation is
> done b/w the 2 MR volumes, I just applied the same transformation matrix to
> the CT, or at least that’s what i was trying to do(i just copied the
> resampling part and changed some variables)... The MR looks good, but I got
> all these streaking artifacts from the CT,  and the CT numbers seems to be
> rescaled to from 0 to around 250, which happens to be about the same range
> of the difference image(the min and max of the IntensityRescaler in line
> 451). I attached a screenshot of one CT slice, and the code(excuse my
> POOOOOR programming skills!). Did I use a wrong resampling filter or…?
>
> Thanks again,
> Neo
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/community/attachments/20160820/a3af0bc8/attachment.html>