hi Luis,<br><br>fixed image size: [215, 172, 48]<br>fixed image origin: [44.0625, 0, 0]<br>fixed image spacing: [0.3125, 0.3125, 1.9808]<br>fixed image direction: <br>1 0 0<br>0 1 0<br>0 0 1<br><br>moving image size: [512, 512, 52]<br>
moving image origin: [0, 0, 0]<br>moving image spacing: [0.3125, 0.3125, 1.92911]<br>moving image direction: <br>1 0 0<br>0 1 0<br>0 0 1<br><br>image mask size: [215, 172, 48]<br>image mask origin: [44.0625, 0, 0]<br>
image mask spacing: [0.3125, 0.3125, 1.9808]<br>image mask direction: <br>1 0 0<br>0 1 0<br>0 0 1<br><br>when GeometryOn() is used: Initial Transform Parameters: [0, 0, 0, 2.34375, 53.125, 2.64351]<br><br><br>When MomentsOn() is used: Initial Transform Parameters: [0, 0, 0, 7.70064, 49.5835, -0.0834114]<br>
<br><br><br>cheers, <br><br>Michael<br><br><br><div class="gmail_quote">On Mon, Mar 1, 2010 at 11:11 PM, Luis Ibanez <span dir="ltr"><<a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.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;">Hi Michiel,<br>
<br>
It is very unlikely that the CenteredTransformInitializer<br>
will fail to overlap the two images.<br>
<br>
<br>
DO NOT assign any value manually to the Transform.<br>
That's what the Initializer is for.<br>
<br>
<br>
Please do the following:<br>
<br>
1) Set the initializer to use GeometryOn();<br>
<br>
2) Run the code again and<br>
post to this mailing list the following information:<br>
<br>
a) The Transform just as it is initialized by the<br>
initializer ( no need to run registration iterations).<br>
<br>
b) Origin, Spacing, Orientation and number of pixels<br>
of both the Fixed Image and Moving Images.<br>
<br>
<br>
Thanks<br>
<br>
<br>
Luis<br>
<br>
<br>
<br>
----------------------------------------------------<br>
On Mon, Mar 1, 2010 at 11:40 AM, michiel mentink<br>
<div><div></div><div class="h5"><<a href="mailto:michael.mentink@st-hughs.ox.ac.uk">michael.mentink@st-hughs.ox.ac.uk</a>> wrote:<br>
> spoke too soon :(<br>
><br>
> Now the moving image dimensions are zeroed too.<br>
><br>
> help?<br>
><br>
> On Mon, Mar 1, 2010 at 4:30 PM, michiel mentink<br>
> <<a href="mailto:michael.mentink@st-hughs.ox.ac.uk">michael.mentink@st-hughs.ox.ac.uk</a>> wrote:<br>
>><br>
>> ah solved it.<br>
>><br>
>> For any that are interested:<br>
>><br>
>> TransformType::OutputVectorType translation_vector;<br>
>> translation_vector[0] = 0.0;<br>
>> translation_vector[1] = 0.0;<br>
>> translation_vector[2] = 0.0;<br>
>> transform->SetTranslation(translation_vector);<br>
>><br>
>> and then comment out<br>
>> initializer->InitializeTransform();<br>
>><br>
>> Still doesn't solve my problem though.<br>
>><br>
>> cheers, Michael<br>
>><br>
>><br>
>> On Mon, Mar 1, 2010 at 3:39 PM, michiel mentink<br>
>> <<a href="mailto:michael.mentink@st-hughs.ox.ac.uk">michael.mentink@st-hughs.ox.ac.uk</a>> wrote:<br>
>>><br>
>>> Dear all,<br>
>>><br>
>>> is there anyway of setting the initial iteration of the<br>
>>> CentredTransformInitializer manually?<br>
>>><br>
>>> The CentredTransformInitializer uses either GeometryOn() or momentsOn()<br>
>>> that calculates<br>
>>> the center of the image, but that produces an initial parameter that<br>
>>> fails the registration even<br>
>>> before it starts, because its errors:<br>
>>> Description: itk::ERROR: MeanSquaresImageToImageMetric(0x9dfa3f0): All<br>
>>> the points mapped to outside of the moving image<br>
>>><br>
>>> Alternatively, can I just overwrite the parameters that are produced by<br>
>>> the centredTransformInitializer? How?<br>
>>><br>
>>> I'm using the VersorRigid3DTransform.<br>
>>><br>
>>> cheers, Michael<br>
>>><br>
>><br>
><br>
><br>
</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>
> Kitware offers ITK Training Courses, for more information visit:<br>
> <a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.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>
</blockquote></div><br>