Hi again<br><br>I have some problems regarding the use of Level Set Methods. I have read the manual and I am now trying to use one of the example files, LaplacianSegmentationLevelSetImageFilter.cxx.<br><br>With the 2D test image,
BrainProtonDensitySlice.PNG and VentricleModel.PNG they use in the manual I get the same results as they do. However when I try with a 3D volume it does not go to well. <br>Images and code can be found at <a href="http://www.idi.ntnu.no/~rieck/dipp/">
http://www.idi.ntnu.no/~rieck/dipp/</a><br><br>speed.PNG is the model I give the filter, the initial segmentation for it to enhance. in.PNG is the image I want to segment. Levelset.PNG is the result of the LaplacianSegmentation. From code:
<br> <br>diffusion->SetInput( image );<br> laplacianSegmentation->SetInput( model );<br> laplacianSegmentation->SetFeatureImage( diffusion->GetOutput() );<br> thresholder->SetInput( laplacianSegmentation->GetOutput() );
<br><br> //debug<br> //ImageType::Pointer debug = laplacianSegmentation->GetOutput();<br> //debug->Update();<br> //return debug;<br><br>debug is here levelset.PNG.<br><br>This is sent trough a thresholding filter to get the final result, seen in
result.PNG.<br><br><br>I have done the same thing with the 2D test image, and the corresponding laplacianSegmentation->GetOutput (VentricleLevelSet.png) does not look anything alike levelset.PNG.<br><br>I am doing something fundamentaly wrong? I have tried with a variety of parameters as well, without any luck. The parameters I used for these example images were: file1 file2 5 1 1 -2 0.
<br><br><br>regards,<br>Christian.<br><br>