<div dir="ltr">Hello. I want to compute the Normalized Mutual information of two images when I give a specific transformation. Below is my code.<div><br></div><div><div>#include "itkMutualInformationImageToImageMetric.h"</div><div>#include "itkRandomImageSource.h"</div><div>#include "itkTranslationTransform.h"</div><div>#include"itkCenteredRigid2DTransform.h"</div><div>#include "itkCenteredTransformInitializer.h"</div><div>#include "itkLinearInterpolateImageFunction.h"</div><div>#include"itkImageFileReader.h"</div><div>#include"itkNormalizedMutualInformationHistogramImageToImageMetric.h"</div><div>#include "itkMattesMutualInformationImageToImageMetric.h"</div><div>#include "itkMutualInformationImageToImageMetric.h"</div><div><br></div><div># include <cstdlib></div><div># include <iostream></div><div># include <iomanip></div><div># include <fstream></div><div># include <iomanip></div><div># include <cmath></div><div># include <ctime></div><div># include <cstring></div><div># include <random></div><div><br></div><div><br></div><div><br></div><div><br></div><div>using namespace std;</div><div><br></div><div><br></div><div><br></div><div>int main(int argc, char *argv[]);</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>const    unsigned int    ImageDimension = 2;</div><div>typedef  unsigned char    PixelType;</div><div>typedef itk::Image<PixelType, ImageDimension> FixedImage;</div><div>typedef itk::Image<PixelType, ImageDimension> MovingImage;</div><div><br></div><div>typedef itk::ImageFileReader<FixedImage> FixedImageTypeReader;</div><div>typedef itk::ImageFileReader<MovingImage> MovingImageTypeReader;</div><div><br></div><div>typedef itk::CenteredRigid2DTransform<double> RigidTransformType;</div><div><br></div><div><br></div><div>typedef itk::NormalizedMutualInformationHistogramImageToImageMetric<FixedImage, MovingImage> MetricType;</div><div>//typedef itk::MattesMutualInformationImageToImageMetric<FixedImage, MovingImage> MetricType;</div><div>//typedef itk::MutualInformationImageToImageMetric<FixedImage, MovingImage> MetricType;</div><div><br></div><div><br></div><div><br></div><div>typedef itk::CenteredTransformInitializer<RigidTransformType, FixedImage, MovingImage> InitializerType;</div><div><br></div><div>RigidTransformType::Pointer rigidtransform;</div><div><br></div><div><br></div><div><br></div><div><br></div><div>MetricType::MeasureType test(double *arg , MetricType::Pointer metric){</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span>MetricType::MeasureType v=0;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span>int i;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre"> </span>double a, e, f;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">        </span>a = arg[0] ;</div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>e = arg[1];</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">    </span>f = arg[2];</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">    </span></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">       </span>rigidtransform->SetAngle(a);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">        </span>RigidTransformType::OutputVectorType t = rigidtransform->GetTranslation();</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>t[0] = e;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">      </span>t[1] = f;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">      </span>rigidtransform->SetTranslation(t);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>RigidTransformType::ParametersType p;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>p.SetSize(5);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>p[0] = a;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">      </span>p[1] = 255.5;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>p[2] = 255.5;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>p[3] = e;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">      </span>p[4] = f;</div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">     </span></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">       </span>metric->UseAllPixelsOn();</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span>metric->SetTransform(rigidtransform);</div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">      </span>cout<<p<<endl;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">     </span>return metric->GetValue(p);</div><div>}</div><div><br></div><div><br></div><div><br></div><div><br></div><div>int main( int argc, char *argv[] )</div><div>{</div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span></div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">      </span>FixedImageTypeReader::Pointer fixedreader = FixedImageTypeReader::New();</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">       </span>MovingImageTypeReader::Pointer movingreader = MovingImageTypeReader::New();</div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span>fixedreader->SetFileName(argv[1]);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>movingreader->SetFileName(argv[2]);</div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">        </span>fixedreader->Update();</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">      </span>movingreader->Update();</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">     </span>FixedImage::Pointer fixedImage = fixedreader->GetOutput();</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>MovingImage::Pointer movingImage = movingreader->GetOutput();</div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">      </span></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">       </span>MetricType::Pointer metric;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">    </span>rigidtransform = RigidTransformType::New();</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">    </span></div><div>  </div><div><span class="gmail-Apple-tab-span" style="white-space:pre">    </span>MetricType::HistogramSizeType h;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">       </span>h.SetSize(2);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>h[0] = h[1]=64;</div><div><br></div><div><br></div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">     </span>metric = MetricType::New();</div><div>  </div><div><span class="gmail-Apple-tab-span" style="white-space:pre"> </span>metric->SetTransform(rigidtransform);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">       </span>metric->SetHistogramSize(h);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">        </span>metric->SetNumberOfFixedImageSamples(10000);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">        </span>metric->SetNumberOfSpatialSamples(10000);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span>metric->SetFixedImage(fixedImage);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>metric->SetMovingImage(movingImage);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">        </span>metric->SetFixedImageRegion(fixedImage->GetLargestPossibleRegion());</div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">    </span>const unsigned int numberOfParameters = rigidtransform->GetNumberOfParameters();</div><div>  </div><div>  </div><div><span class="gmail-Apple-tab-span" style="white-space:pre">      </span></div><div>  </div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span></div><div>  </div><div><span class="gmail-Apple-tab-span" style="white-space:pre">    </span>std::cout<<"Spacing: "<<movingImage->GetSpacing()<<", Origin: "<<movingImage->GetOrigin()<<", Direction: "<<movingImage->GetDirection() <<std::endl;</div><div>  </div><div> </div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">     </span>itk::LinearInterpolateImageFunction<MovingImage, double>::Pointer interpolator = itk::LinearInterpolateImageFunction<MovingImage, double>::New();</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>interpolator->SetInputImage(movingImage);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span>metric->SetInterpolator(interpolator);</div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">     </span>InitializerType::Pointer initializer = InitializerType::New();</div><div><span class="gmail-Apple-tab-span" style="white-space:pre"> </span>initializer->SetTransform(rigidtransform);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>initializer->SetFixedImage(fixedreader->GetOutput());</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">    </span>initializer->SetMovingImage(movingreader->GetOutput());</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>initializer->GeometryOn();</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">  </span>initializer->InitializeTransform();</div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">        </span>std::cout << fixedreader->GetOutput()->GetLargestPossibleRegion().GetSize();</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">       </span>RigidTransformType::InputPointType mc = rigidtransform->GetCenter();</div><div>  </div><div><span class="gmail-Apple-tab-span" style="white-space:pre">     </span></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">       </span>MovingImage::SizeType c = movingreader->GetOutput()->GetLargestPossibleRegion().GetSize();</div><div><br></div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">     </span>rigidtransform->SetCenter(mc);</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">      </span>std::cout << "parameters: " << rigidtransform->GetParameters() << std::endl;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">    </span>std::cout << "pcenter: " << mc << std::endl;</div><div>  </div><div><span class="gmail-Apple-tab-span" style="white-space:pre">    </span>double p[3];</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span>p[0] = 0.52;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span>p[1] = -138;</div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span>p[2] = 52;</div><div> </div><div><span class="gmail-Apple-tab-span" style="white-space:pre">    </span>std::cout << "Value: " << test(p,metric) << std::endl;</div><div><br></div><div><br></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">   </span></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">       </span></div><div><span class="gmail-Apple-tab-span" style="white-space:pre">       </span>return EXIT_SUCCESS;</div><div>}</div></div><div><br></div><div>Whenever I execute the code, I get a strange message (see attached image). What could be the problem?</div></div>