[Insight-users] Crash in ShapeDetectionLevelSetImageFilter

Nils Hanssen hanssen at caesar . de
Tue, 8 Jul 2003 09:42:32 +0200


This is a multi-part message in MIME format.

------=_NextPart_000_0003_01C34535.41758570
Content-Type: text/plain;
	charset="iso-8859-1"
Content-Transfer-Encoding: 7bit

Hi,

I am trying to understand the behaviour of the ShapeDetectionLSIF.
Therefore, I set the propagation-weighting to zero. By setting the
curvature-weighting to a value of one I would expect that the inital surface
is shrinking to a point (I set the MaxRMSError to zero) and the number of
iterations very high.

The filter is crashing in
SegmentationLevelSetFunction<TImageType, TFeatureImageType>
::PropagationSpeed(const NeighborhoodType &neighborhood, const
FloatOffsetType &offset) const
[...]
--> else return (
static_cast<ScalarValueType>(m_SpeedImage->GetPixel(idx)) ); // crashing
here
[...]

When I set the propagation-weighting to - for instance - 0.0001 the filter
is not crashing, but the contour is not shrinking according to the mean
curvature.

Is that the correct behaviour of the filter?

I would appreciate any suggestions. Thanks!


Regards,
Nils

-------------------------
Nils Hanssen
Surgical Systems Laboratory
research center c ae sa r
Ludwig-Erhard-Allee 2
53175 Bonn
fon: +49-228-9656-197
fax: +49-228-9656-117
http://www . caesar . de/ssl





------=_NextPart_000_0003_01C34535.41758570
Content-Type: text/html;
	charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META HTTP-EQUIV=3D"Content-Type" CONTENT=3D"text/html; =
charset=3Diso-8859-1">


<META content=3D"MSHTML 6.00.2600.0" name=3DGENERATOR></HEAD>
<BODY>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2>Hi,</FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2></FONT></SPAN>&nbsp;</DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial size=3D2>I am =
trying to=20
understand the behaviour of the ShapeDetectionLSIF. </FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial =
size=3D2>Therefore, I set the=20
propagation-weighting to zero. By setting the curvature-weighting to a =
value of=20
one I would expect that the inital surface is shrinking to a point (I =
set the=20
MaxRMSError to zero) and the number of iterations very high.=20
</FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2></FONT></SPAN>&nbsp;</DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial size=3D2>The =
filter is=20
crashing in </FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2>SegmentationLevelSetFunction&lt;TImageType,=20
TFeatureImageType&gt;<BR>::PropagationSpeed(const NeighborhoodType=20
&amp;neighborhood, const FloatOffsetType &amp;offset) =
const</FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2>[...]</FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial size=3D2>--&gt; =
else return (=20
static_cast&lt;ScalarValueType&gt;(m_SpeedImage-&gt;GetPixel(idx)) ); // =

crashing here</FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2>[...]</FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2></FONT></SPAN>&nbsp;</DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial size=3D2>When I =
set the=20
propagation-weighting to&nbsp;- for instance - 0.0001 the filter is not=20
crashing, but the contour is not shrinking according to the mean=20
curvature.</FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2></FONT></SPAN>&nbsp;</DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial size=3D2>Is =
that the correct=20
behaviour of the filter?</FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2></FONT></SPAN>&nbsp;</DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial size=3D2>I =
would appreciate=20
any suggestions. Thanks!</FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2></FONT></SPAN>&nbsp;</DIV>
<DIV><SPAN class=3D765043007-08072003></SPAN>&nbsp;</DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2>Regards,</FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2>Nils</FONT></SPAN></DIV>
<DIV><SPAN class=3D765043007-08072003><FONT face=3DArial=20
size=3D2></FONT></SPAN>&nbsp;</DIV>
<P><FONT face=3D"Courier New" =
size=3D2>-------------------------</FONT><BR><FONT=20
face=3D"Courier New" size=3D2>Nils Hanssen</FONT><FONT=20
face=3D"Times New Roman"><BR></FONT><FONT face=3D"Courier New" =
size=3D2>Surgical=20
Systems Laboratory</FONT><BR><FONT face=3D"Courier New" =
size=3D2>research center c=20
ae sa r</FONT><FONT face=3D"Times New Roman"><BR></FONT><FONT =
face=3D"Courier New"=20
size=3D2>Ludwig-Erhard-Allee 2<BR></FONT><FONT face=3D"Courier New" =
size=3D2>53175=20
Bonn</FONT><FONT face=3D"Times New Roman"><BR></FONT><FONT =
face=3D"Courier New"=20
size=3D2>fon: +49-228-9656-197</FONT><FONT face=3D"Times New =
Roman"><BR></FONT><FONT=20
face=3D"Courier New" size=3D2>fax: +49-228-9656-117</FONT><FONT=20
face=3D"Times New Roman"><BR><U></U></FONT><U><FONT face=3D"Courier New" =

color=3D#0000ff size=3D2><A href=3D"http://www . caesar . de/ssl"=20
target=3D_blank>http://www . caesar . de/ssl</A></FONT></U> </P><BR>
<DIV>&nbsp;</DIV></BODY></HTML>

------=_NextPart_000_0003_01C34535.41758570--