<div dir="ltr">Hi everyone,<div><br></div><div style>I'm attempting some multi-resolution demons, as input I use an initial affine transform. When monitoring convergence however, I get peaks at every pyramid level (see attached diagram). Is this normal behaviour</div>
<div style><br></div><div style> my pseudocode looks something like:<br></div><div style><br></div><div style>//Define Demons filter</div><div style><br></div><div style> RegistrationFilterType::Pointer filter = RegistrationFilterType::New();<br>
</div><div style><br></div><div style>//Create initial Affine transform and convert to Displacement field, called AffineField</div><div style><br></div><div style><div> AffineField=converter->GetOutput();</div><div> </div>
</div><div style><div><div> filter->SetInput(AffineField);</div><div> filter->SetStandardDeviations( 0.8);</div></div><div><br></div><div style>//use multiResolution</div><div style><br></div><div style><div> MultiResRegistrationFilterType::Pointer multires = MultiResRegistrationFilterType::New();</div>
<div><br></div><div> multires->SetRegistrationFilter( filter );</div><div> multires->SetNumberOfLevels( 4 );</div><div> multires->SetFixedImage( FixedImage );</div><div> multires->SetMovingImage( MovedImage );</div>
<div> unsigned int nIterations[4] = {100, 100, 50, 50 };</div><div> multires->SetNumberOfIterations( nIterations );</div><div style> multires->Update();</div><div><br></div></div><div style>Am I correct in setting the initial affine field in this manner, does this set the affine field for the first stage of the pyramid? or is it a case of the affine field being re-applied at every pyramid value and hence peaks in the metric</div>
<div style><br></div><div style>Thoughts appreciated</div><div style><br></div><div style>Thanks</div><div style><br></div><div style>Emma</div><div><br></div><div><br></div><div><br></div></div></div>