Hello,<br><br>I ran DeformableRegistration8.cxx using the following command:<br>(I altered the reading of files suitable for DICOM)<br><br>./DeformableRegistration8 subjectout.dcm atlasout.dcm outputImagefile differenceOutputfile differenceBeforeRegistration deformationField<br>
<br>Starting Registration<br>ExceptionObject caught !<br><br>itk::ExceptionObject (0x4f51fd0)<br>Location: "void itk::MattesMutualInformationImageToImageMetric<TFixedImage, TMovingImage>::GetValueAndDerivative(const typename itk::ImageToImageMetric<TFixedImage, TMovingImage>::ParametersType&, typename itk::ImageToImageMetric<TFixedImage, TMovingImage>::MeasureType&, typename itk::ImageToImageMetric<TFixedImage, TMovingImage>::DerivativeType&) const [with TFixedImage = itk::Image<short int, 3u>, TMovingImage = itk::Image<short int, 3u>]"<br>
File: /trumpet/downloads/3DSlicer/Slicer3-lib/Insight/Code/Review/itkOptMattesMutualInformationImageToImageMetric.txx<br>Line: 1034<br>Description: itk::ERROR: MattesMutualInformationImageToImageMetric(0x1db9020): Joint PDF summed to zero<br>
<br>Does anyone know what this error means?<br><br>john<br><br><br><div class="gmail_quote">On Thu, Oct 8, 2009 at 5:19 PM, John Drozd <span dir="ltr"><<a href="mailto:john.drozd@gmail.com">john.drozd@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Sorry, my oversight. The parameters file was for DeformableRegistration1.cxx.<br>The parameters file just happened to be listed in the ITK Software Guide directly after the code listing for DeformableRegistration8.cxx, and that's why I thought a parameter file was needed for DeformableRegistration8.cxx. Looking at the code for DeformableRegistration8.cxx, I see that only a fixedImageFile, movingImageFile and outputImageFile are required.<br>
<br>Thanks,<br><font color="#888888">john</font><div><div></div><div class="h5"><br><br><div class="gmail_quote">On Thu, Oct 8, 2009 at 3:14 PM, John Drozd <span dir="ltr"><<a href="mailto:john.drozd@gmail.com" target="_blank">john.drozd@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hello,<br><br>I have two brain volumes, that I successfully loaded into ITK from two different sets of dicom series, namely,<br><br>image 1 data: (header info that I obtained when I saved my first volume as a nrrd file in Slicer)<br>
type: short<br>dimension: 3<br>space: left-posterior-superior<br>sizes: 256 256 124<br>space directions: (0,0.9375,0) (0,0,-0.9375) (-1.3,0,0)<br>kinds: domain domain domain<br>endian: little<br>encoding: gzip<br>space origin: (-79.9,-125.9,121.8)<br clear="all">
<br>image 2 data: (header info that I obtained when I saved my second volume as a nrrd file in Slicer)<br>type: unsigned short<br>dimension: 3<br>space: left-posterior-superior<br>sizes: 256 256 166<br>space directions: (0,0,-0.93905) (0,0.93905,0) (1.2024,0,0)<br>
kinds: domain domain domain<br>endian: little<br>encoding: gzip<br>space origin: (198.396,-0,-0)<br><br>I want to use a deformable registration such that I will deform the first image into the the second image. Because the images will have different modalities, I am planning to use a modified version of either DeformableRegistration8.cxx (which uses itk::BSplineDeformableTransform, itk::MattesMutualInformationImageToImageMetric with the LBFGSB Optimizer). I may also try using DeformableRegistration14.cxx which uses a gradientsteepestdescent optimizer. As a final refinement, maybe a multiresolution scheme of these routines to speed up the computation.<br>
<br>I noticed that in the parameter file for DeformableRegistration8.cxx, only one size Nx, Ny, Nz is specified for both images, but my images have different sizes and different space directions, and different space origins as listed in the nrrd headers above for these images.<br>
<br>I will try on my own, but I was wondering if someone has some helpful insight or hints as to whether there are many modifications that I have to do, to refine DeformableRegistration8.cxx or DeformableRegistration14.cxx and their associated classes to make this workable for my situation.<br>
<br>Any suggestions are welcome.<br><br>Thank you.<br><br>john<br><br>-- <br>John Drozd<br>Postdoctoral Fellow<br>Imaging Research Laboratories<br>Robarts Research Institute<br>Room 1256<br>100 Perth Drive<br>London, Ontario, Canada<br>
N6A 5K8<br>(519) 661-2111 ext. 25306<br>
</blockquote></div><br><br clear="all"><br>-- <br>John Drozd<br>Postdoctoral Fellow<br>Imaging Research Laboratories<br>Robarts Research Institute<br>Room 1256<br>100 Perth Drive<br>London, Ontario, Canada<br>N6A 5K8<br>
(519) 661-2111 ext. 25306<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>John Drozd<br>Postdoctoral Fellow<br>Imaging Research Laboratories<br>Robarts Research Institute<br>Room 1256<br>100 Perth Drive<br>London, Ontario, Canada<br>
N6A 5K8<br>(519) 661-2111 ext. 25306<br>