<br>Hi Osc,<br><br>I'll go in a minute to your second email <br>(where you posted the debugging information)<br><br><br>but, just to answer your question:<br><br><div class="gmail_quote">The reason why the ITK image inverts that matrix<br>
is to make possible a rapid computation of <br><br> PhysicalPointToIndex and<br> IndexToPhysicalPoint<br><br>These are computations that are usually applied to many<br>pixels, and therefore it is convenient to cache the two <br>
Matrices that perform the transformation<br>(one matrix being the inverse of the other).<br><br>A change in the image spacing requires a recomputation<br>of these two matrices.<br><br><br> Regards,<br><br><br> Luis<br>
<br><br>----------------------------------------------------------------------------<br>On Fri, Apr 9, 2010 at 8:09 PM, Osc@r Math <span dir="ltr"><<a href="mailto:olguna@cimat.mx">olguna@cimat.mx</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;">
<div link="blue" vlink="purple" lang="ES-MX">
<div>
<p class="MsoNormal"><span lang="ES-TRAD">Hi All…!!!</span></p>
<p class="MsoNormal"><span lang="ES-TRAD"> </span></p>
<p class="MsoNormal"><span lang="EN-US">Recently I took back an ITK-VTK-BCB5
project, but as I know ITK & VTK improves, I decided to update those
toolkits (by the way, congrats!!! For this big effort).</span></p>
<p class="MsoNormal"><span lang="EN-US">Actually, compiled ITK 3.16 & VTK 5.4.2
Releases</span></p>
<p class="MsoNormal"><span lang="EN-US">BCB 5 with Borland C++ compiler 5.5.</span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">Today I get a runtime error that it’s
pretty mysterious (at least for me)…</span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">The application reads a 3D image by VTK
reader, show it by VTK renderization, send it to ITK pipeline, calculates
partial derivatives & magnitude gradient, then save Dx, Dy, Dz & |V|.</span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">The image has a spacing, an origin, among
others…</span></p>
<p class="MsoNormal"><span lang="EN-US">When spacing is 1.0, 1.0, 1.0 there’s
no problem, everything works fine,</span></p>
<p class="MsoNormal"><span lang="EN-US">but if spacing changes anyone of those
dimensions, application crashes.</span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">The application is capable to read image,
show it, but when tries to connect VTKReader to ITKCaster (is first action with
ITK) send an “Abnormal Program Termination”</span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">In Debug mode I realize that crashes trying
to execute </span></p>
<p class="MsoNormal"><span lang="EN-US">~vnl_svd() {}</span></p>
<p class="MsoNormal"><span lang="EN-US">In vnl_svd.h file and mark an error like:</span></p>
<p class="MsoNormal"><span lang="EN-US">… Access violation at address XXXX …</span></p>
<p class="MsoNormal"><span lang="EN-US">If continue debugging, crashes too in:</span></p>
<p class="MsoNormal"><span lang="EN-US">~vnl_diag_matrix() {}</span></p>
<p class="MsoNormal"><span lang="EN-US">In vnl_diag_matrix.h file with similar
error</span></p>
<p class="MsoNormal"><span lang="EN-US">And finally throws an “Abnormal Program
Termination”…</span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">I didn´t know why tries to invert a matrix
that contains</span></p>
<p class="MsoNormal"><span lang="EN-US">xSpacing 0 0</span></p>
<p class="MsoNormal"><span lang="EN-US"> 0 ySpacing 0</span></p>
<p class="MsoNormal"><span lang="EN-US"> 0 0 zSpacing</span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">¡¡¡ Someone could help me !!!</span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">TNX in advance</span></p>
</div>
<br><br>__________ Información de ESET Smart Security, versión de la base de firmas de virus 5014 (20100409) __________<br><br>ESET Smart Security ha comprobado este mensaje.<br><br><a href="http://www.eset.com" target="_blank">http://www.eset.com</a><br>
</div>
<br>_____________________________________<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></blockquote></div><br>