<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><span style="font-family: monospace;">Hi all,<br><br><br>I am faced with a problem
regarding the pasting or stitching of oblique images in a final volume.
For instance, suppose I have three images, each one oriented at Rx=[-5,
0, +5] degrees, how could I read these oriented images
and place all three in a common volume?<br><br>I am not certain that the OrientImageFilter can do these, even though it can allow us to use the following command:<br><br>orienter --> SetUseImageDirection( );</span><br><br><br>any insight on this would be appreciated.<br><br><br><br>Pascal<br><br><br><br><br><br>--- On <b>Mon, 2/2/09, Bill Lorensen <i><bill.lorensen@gmail.com></i></b> wrote:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;">From: Bill Lorensen <bill.lorensen@gmail.com><br>Subject: Re: [Insight-users] itk::OrientImageFilter<br>To: "Xavier Planes" <xplanes@cistib.upf.edu><br>Cc: "insight-users@itk.org" <insight-users@itk.org><br>Date: Monday, February 2, 2009, 10:27 AM<br><br><pre>I think the OrientImageFIlter's letter notation is the opposite of the<br>DICOM orientations.<br><br>Instead of LPS, try RAI.<br><br>Bill<br>On Mon, Feb 2, 2009 at 10:01 AM, Xavier
Planes <xplanes@cistib.upf.edu><br>wrote:<br>> Hi,<br>><br>><br>><br>> I've tried to use the filter itk::OrientImageFilter.<br>I've<br>> read a several documentation about DICOM orientation standards, tried some<br>> source code examples and read several emails in the mailing list without<br>> success… . My main goal is to be able to transform the DICOM volume<br>image<br>> that can be in any DICOM orientation (LPS, RAS, RAI, LPI, …) to an itk<br>> Volume image always oriented in the same way (for example LPS).<br>><br>><br>><br>> I've tried this code:<br>><br>><br>><br>> typedef itk::OrientedImage<dcmAPI::PixelType, 3><br>> OrientedVolumeType;<br>><br>> typedef itk::ImageSeriesReader< OrientedVolumeType ><br>> VolumeReader;<br>><br>><br>><br>> VolumeReader::Pointer
volumeReader = VolumeReader::New();<br>><br>> itk::GDCMImageIO::Pointer dicomIO =<br>itk::GDCMImageIO::New();<br>><br>> volumeReader->SetImageIO( dicomIO );<br>><br>> volumeReader->SetFileNames( sliceFilePaths );<br>><br>> volumeReader->Update();<br>><br>> OrientedVolumeType::Pointer volume =<br>> volumeReader->GetOutput();<br>><br>><br>><br>> itk::OrientImageFilter<VolumeType,VolumeType>::Pointer orienter =<br>> itk::OrientImageFilter<VolumeType,VolumeType>::New();<br>><br>> orienter->UseImageDirectionOn();<br>><br>> orienter->SetDesiredCoordinateOrientation(<br>> itk::SpatialOrientation::ITK_COORDINATE_ORIENTATION_LPS);<br>><br>> orienter->SetInput(volume);<br>><br>>
orienter->Update();<br>><br>><br>><br>> It works fine with RAI input image orientation. However,<br>it<br>> doesn't work for RSA input image orientation because the output image<br>is not<br>> in LPS orientation. I've used VolView to check that the input images<br>are<br>> correct and it works fine: the VTK output image is always in the same<br>> orientation.<br>><br>><br>><br>> I've also debugged a little bit with DICOM RSA input orientation. The<br>> itk::OrientImageFilter does this steps:<br>><br>> 1. RSA -> RAS. Permute Y And Z axis<br>><br>> 2. RAS -> LAS. Flip X axis<br>><br>> 3. LAS -> LPS. Flip Y axis<br>><br>> However the output is not correct. It seems that when input image is<br>> permuted (Y and Z axis), the orientation of the third axis is also<br>> indirectly flipped (X). So after the permute filter, the image
orientation<br>> is not RAS, as it was expected by the filter.<br>><br>><br>><br>> I was thinking to create a table for permutation and flipping parameters<br>> depending on each orientation. However I don't have sample DICOM data<br>for<br>> all the possible orientations.<br>><br>><br>><br>> Has anybody tried this filter? Has anybody tried RSA orientation input<br>> image? Do you know what input image orientations have been tried?<br>><br>><br>><br>> Thanks!<br>><br>> Xavi<br>><br>> _______________________________________________<br>> Insight-users mailing list<br>> Insight-users@itk.org<br>> http://www.itk.org/mailman/listinfo/insight-users<br>><br>><br>_______________________________________________<br>Insight-users mailing list<br>Insight-users@itk.org<br>http://www.itk.org/mailman/listinfo/insight-users<br></pre></blockquote></td></tr></table><br>