[Insight-users] Question concerning level set term.

[tomtom]

Hello everyone,

I am trying to implement an own level set term.
The evolution equation that I have looks like: phi_t = A(x,y) * N,
where A(x,y) is real and N is the inwards directed unit normal vector.

So, as far as I understood things, I can use the backwardGradient, which
should be N, if I normalize it with its length.

Is it right that the Value()-method should look like:

....Value(const LevelSetInputIndexType& iP)
{
LevelSetOutputRealType outputValue = itk::NumericTraits<
LevelSetOutputRealType >::Zero;

LevelSetOutputRealType AValue = ComputeA(x,y);

LevelSetGradientType backwardGradient =
this->m_CurrentLevelSetPointer->EvaluateBackwardGradient( iP );
double lengthBackwardGradient = ...;

for(unsigned int dim = 0; dim < ImageDimension; dim++)
outputValue+=backwardGradient[dim] * AValue / lengthBackwardGradient;

return outputValue;
}

Or am I doing something wrong.

Thanks for any remarks or hints.



--
View this message in context: http://itk-insight-users.2283740.n2.nabble.com/Question-concerning-level-set-term-tp7426557p7426557.html
Sent from the ITK Insight Users mailing list archive at Nabble.com.