[Insight-users] Jacobian value

[Jesse]

Hi all,

 

In DeformableRegistration1, I want to get the Jacobian value, and I added some code at the end of the program

 

/////////////////////////////////////////////////////////////////

  registrationFilter->GetJacobianImage();
  
  itk::ImageFileWriter<ImageType>::Pointer writer3;
  writer3 = itk::ImageFileWriter<ImageType>::New();
  

    writer3->SetInput(registrationFilter->GetOutput() );

  writer3->SetFileName( "Jacobian.mhd" );
  writer3->Update();


/////////////////////////////////////////////////////////////////

 

but no new file occurs, what's wrong with it?

 

in addition, why not use 0 as the threshold of Jacobian value in method EnforceDiffeomorphism()?

 

Best

 

Jesse



			
---------------------------------
Do you Yahoo!?
 Check out the new Yahoo! Front Page. www.yahoo.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://public.kitware.com/pipermail/insight-users/attachments/20041111/5bb9047e/attachment.htm