<div class="gmail_quote">Hi, <div><br></div><div>I tried with radian, </div><div>--------------------------------------------</div><div>.....</div><div>transform->SetScale(1.08438);<br>transform->SetAngle(<span style="font-family:arial, sans-serif;font-size:12px;border-collapse:collapse"><b><i>0.154531</i></b></span><b>)</b>;<div>
<br>
TransformType::TranslationType translation;<br>translation[0] = -26.6822;<br>translation[1] = -93.6822;</div></div><div>.....</div><div>--------------------------------------------</div><div><br></div><div>I got new value but still not '<span style="font-family:arial, sans-serif;font-size:12px;border-collapse:collapse">4551.86' (the value from example).</span></div>
<div><font face="arial, sans-serif"><span style="border-collapse:collapse"><br></span></font></div><div><font face="arial, sans-serif"><span style="border-collapse:collapse">I read from "ITK software guide ver 2.4" textbook, page 366. they explained (as i understand) that the output transform params from registration process are actually applied to FixeImage, in order to get the minimum metric value. </span></font></div>
<div><font face="arial, sans-serif"><span style="border-collapse:collapse"><br></span></font></div><div><font face="arial, sans-serif"><span style="border-collapse:collapse">It is still not clear that, the transform params pass directly to metric object are interpreted same way as in registration process (i mean they also applied to Fixed image instead of Moving Image). I guess the reason may come from this unobvious understanding. Do you have any idea about this? </span></font></div>
<div><font face="arial, sans-serif"><span style="border-collapse:collapse"><br></span></font></div><div><font face="arial, sans-serif"><span style="border-collapse:collapse"><br>
</span></font></div><div><font face="arial, sans-serif"><span style="border-collapse:collapse">Thanks, </span></font></div><div><font face="arial, sans-serif"><span style="border-collapse:collapse">Quyps </span></font></div>
<div><font face="arial, sans-serif"><span style="border-collapse:collapse"><br></span></font></div><div><font face="arial, sans-serif"><span style="border-collapse:collapse"><br>
</span></font></div><div><div class="gmail_quote">2009/10/3 Bill Lorensen <span dir="ltr"><<a href="mailto:bill.lorensen@gmail.com" target="_blank">bill.lorensen@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div></div><div>
/** Set/Get the angle of rotation in radians */<br>
----------<br>
void SetAngle(TScalarType angle);<br>
itkGetConstReferenceMacro( Angle, TScalarType );<br>
<br>
Try:<br>
transform->SetAngle(0.154531);<br>
</div></div><div><div></div><div><div><div></div><div><br>
On Fri, Oct 2, 2009 at 4:09 PM, Quy Pham Sy <<a href="mailto:phamsyquybk@gmail.com" target="_blank">phamsyquybk@gmail.com</a>> wrote:<br>
> hi,<br>
> I run ImageRegistration7 example in ITK.<br>
> and the output like following<br>
> -------------------------------------<br>
> Scale = 1.08438<br>
> Angle (radians) 0.154531<br>
> Angle (degrees) 8.85397<br>
> Center X = 170.327<br>
> Center Y = 148.556<br>
> Translation X = -26.6822<br>
> Translation Y = -93.4428<br>
> Iterations = 500<br>
> Metric value = 4551.86<br>
> ---------------------------------------<br>
> I made a program in which i use a same metric, transformation, interpolate<br>
> object type as above example, with same two input image file. I set above<br>
> value to calculate metric value directly from metric object, but the result<br>
> is not same as what the example made..<br>
> here is the code, did i miss something?<br>
> I'm running out of time! really appreciate any help!<br>
> -------------------------------------------------------------------------<br>
> int main( int argc, char * argv[] )<br>
> {<br>
> const unsigned int Dimension = 2;<br>
> typedef unsigned char PixelType;<br>
> typedef itk::Image< PixelType, Dimension > ImageType;<br>
> typedef itk::ImageFileReader< ImageType > ReaderType;<br>
> ReaderType::Pointer fixedReader = ReaderType::New();<br>
> ReaderType::Pointer movingReader = ReaderType::New();<br>
> fixedReader->SetFileName( "fixed.bmp" );<br>
> movingReader->SetFileName( "moving.bmp" );<br>
> try<br>
> {<br>
> fixedReader->Update();<br>
> movingReader->Update();<br>
> }<br>
> catch( itk::ExceptionObject & excep )<br>
> {<br>
> std::cerr << "Exception catched !" << std::endl;<br>
> std::cerr << excep << std::endl;<br>
> }<br>
> typedef itk::MeanSquaresImageToImageMetric< ImageType, ImageType ><br>
> MetricType;<br>
> MetricType::Pointer metric = MetricType::New();<br>
> typedef itk::CenteredSimilarity2DTransform < double > TransformType;<br>
> TransformType::Pointer transform = TransformType::New();<br>
> typedef itk::LinearInterpolateImageFunction< ImageType, double ><br>
> InterpolatorType;<br>
> InterpolatorType::Pointer interpolator = InterpolatorType::New();<br>
> ImageType::Pointer fixedImage = fixedReader->GetOutput();<br>
> ImageType::Pointer movingImage = movingReader->GetOutput();<br>
> typedef itk::CenteredTransformInitializer< TransformType, ImageType,<br>
> ImageType > TransformInitializerType;<br>
> TransformInitializerType::Pointer initializer =<br>
> TransformInitializerType::New();<br>
> initializer->SetTransform( transform );<br>
> initializer->SetFixedImage( fixedImage );<br>
> initializer->SetMovingImage( movingImage );<br>
> initializer->MomentsOn();<br>
> initializer->InitializeTransform();<br>
> metric->SetTransform( transform );<br>
> metric->SetInterpolator( interpolator );<br>
> metric->SetFixedImage( fixedImage );<br>
> metric->SetMovingImage( movingImage );<br>
> metric->SetFixedImageRegion( fixedImage->GetBufferedRegion() );<br>
> try<br>
> {<br>
> metric->Initialize();<br>
> }<br>
> catch( itk::ExceptionObject & excep )<br>
> {<br>
> std::cerr << "Exception catched !" << std::endl;<br>
> std::cerr << excep << std::endl;<br>
> return EXIT_FAILURE;<br>
> }<br>
> transform->SetScale(1.08438);<br></div></div>
> transform->SetAngle(<span style="font-family:arial, sans-serif;font-size:12px;border-collapse:collapse"><b><i>0.154531</i></b></span><b>)</b>;<div><br>
> TransformType::TranslationType translation;<br>
> translation[0] = -26.6822;<br>
> translation[1] = -93.6822;<br>
> MetricType::TransformParametersType transformParams =<br>
> transform->GetParameters();<br>
> double metricval = metric->GetValue(transformParams);<br>
> return EXIT_SUCCESS;<br>
> }<br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
</div></div></div><div>> _____________________________________<br>
> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
><br>
> Visit other Kitware open-source projects at<br>
> <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
><br>
> Please keep messages on-topic and check the ITK FAQ at:<br>
> <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
><br>
><br>
</div></blockquote></div><br><br clear="all"><br>-- <br>Pham Sy Quy<br>HCI Lab, Advanced Fusion Technology Department, <br>Room 1211, New Millennium Building<br>Konkuk University, Seoul, Korea<br>Mobile: +82-10-9800-8104<br>
</div>
</div><br><br clear="all"><br>-- <br>Pham Sy Quy<br>HCI Lab, Advanced Fusion Technology Department, <br>Room 1211, New Millennium Building<br>Konkuk University, Seoul, Korea<br>Mobile: +82-10-9800-8104<br>