<div dir="ltr">Are you setting importer->SetImportVoidPointer(data);<br><br>"data" is the output from itk.<br><br>I hope you had a look at <a href="http://vtk.org/gitweb?p=VTK.git;a=blob;f=Imaging/Testing/Cxx/ImportExport.cxx">http://vtk.org/gitweb?p=VTK.git;a=blob;f=Imaging/Testing/Cxx/ImportExport.cxx</a><br>
<br>Jothy<br><br><br><div class="gmail_quote">On Wed, Nov 3, 2010 at 12:40 PM, Xiaopeng Yang <span dir="ltr"><<a href="mailto:yxp233@postech.ac.kr">yxp233@postech.ac.kr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div link="blue" vlink="purple" lang="EN-US">
<div>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);">Thank you so much for your warmhearted help! My code is kind of
huge. Let me explain briefly: the program is coded for segmentation using ITK, visualization
by VTK, and interface by QT. It consists of image loading, diffusion filtering,
then fast marching segmentation using diffusion filtering output as input. I
attach the code here and hope you guys can help me to find inappropriate part.
Thank you very much in advance!</span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: green;">//connect vtk to itk</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">template</span><span style="font-size: 10pt; font-family: "Courier New";"> <<span style="color: blue;">typename</span>
ITK_Exporter, <span style="color: blue;">typename</span> VTK_Importer></span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">void</span><span style="font-size: 10pt; font-family: "Courier New";"> ConnectPipelines(ITK_Exporter exporter,
VTK_Importer* importer)</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">{</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
importer->SetUpdateInformationCallback(exporter->GetUpdateInformationCallback());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
importer->SetPipelineModifiedCallback(exporter->GetPipelineModifiedCallback());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
importer->SetWholeExtentCallback(exporter->GetWholeExtentCallback());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> importer->SetSpacingCallback(exporter->GetSpacingCallback());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
importer->SetOriginCallback(exporter->GetOriginCallback());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
importer->SetScalarTypeCallback(exporter->GetScalarTypeCallback());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
importer->SetNumberOfComponentsCallback(exporter->GetNumberOfComponentsCallback());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
importer->SetPropagateUpdateExtentCallback(exporter->GetPropagateUpdateExtentCallback());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
importer->SetUpdateDataCallback(exporter->GetUpdateDataCallback());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
importer->SetDataExtentCallback(exporter->GetDataExtentCallback());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> importer->SetBufferPointerCallback(exporter->GetBufferPointerCallback());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
importer->SetCallbackUserData(exporter->GetCallbackUserData());</span></p>
<p class="MsoNormal"><span style="font-size: 10pt; font-family: "Courier New";">}</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal"><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: green;">//Paramter definiation</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">#define</span><span style="font-size: 10pt; font-family: "Courier New";"> N_SEEDS 11</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">typedef</span><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">float</span>
InputPixelType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">typedef</span><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">float</span>
OutputPixelType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">typedef</span><span style="font-size: 10pt; font-family: "Courier New";"> itk::Image< InputPixelType, 3 >
InputImageType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">typedef</span><span style="font-size: 10pt; font-family: "Courier New";"> itk::Image< OutputPixelType, 3 >
OutputImageType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">typedef</span><span style="font-size: 10pt; font-family: "Courier New";"> itk::ImageSeriesReader< InputImageType
> ReaderType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">ReaderType::Pointer reader = ReaderType::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">typedef</span><span style="font-size: 10pt; font-family: "Courier New";"> itk::GDCMImageIO
ImageIOType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">typedef</span><span style="font-size: 10pt; font-family: "Courier New";"> itk::GDCMSeriesFileNames
NamesGeneratorType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">ImageIOType::Pointer gdcmIO = ImageIOType::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">NamesGeneratorType::Pointer namesGenerator =
NamesGeneratorType::New();</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal" style="text-indent: 12pt;"><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">typedef</span><span style="font-size: 10pt; font-family: "Courier New";"> itk::VTKImageExport< OutputImageType >
ExportFilterType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">typedef</span><span style="font-size: 10pt; font-family: "Courier New";"> itk::RescaleIntensityImageFilter< </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> OutputImageType, InputImageType >
RescaleFilterType;</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">RescaleFilterType::Pointer rescaler =
RescaleFilterType::New();</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">typedef</span><span style="font-size: 10pt; font-family: "Courier New";"> itk::BinaryThresholdImageFilter<
InputImageType, </span></p>
<p class="MsoNormal"><span style="font-size: 10pt; font-family: "Courier New";">
OutputImageType > ThresholdingFilterType;</span><span style="font-size: 11pt; color: rgb(31, 73, 125);"></span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">SimpleView::SimpleView() </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">{</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">this</span>->ui = <span style="color: blue;">new</span> Ui_SimpleView;</span></p>
<p class="MsoNormal" style="text-indent: 12pt;"><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">this</span><span style="font-size: 10pt; font-family: "Courier New";">->ui->setupUi(<span style="color: blue;">this</span>);</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: green;"> </span><span style="font-size: 10pt; font-family: "Courier New"; color: green;">// Set up action signals and slots</span></p>
<p class="MsoNormal" style="text-indent: 12pt;"><span style="font-size: 10pt; font-family: "Courier New";">connect(<span style="color: blue;">this</span>->ui->actionOpen,
SIGNAL(triggered()), <span style="color: blue;">this</span>, SLOT(slotOpen()));</span><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style="text-indent: 12pt;"><span style="font-size: 10pt; font-family: "Courier New";">connect(<span style="color: blue;">this</span>->ui->pushButton_14,
SIGNAL(clicked()), <span style="color: blue;">this</span>, SLOT(slotContour()));</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal"><span style="font-size: 10pt; font-family: "Courier New";">};</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal"><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">void</span><span style="font-size: 10pt; font-family: "Courier New";"> SimpleView::slotOpen()</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">{</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> QString directory =
QFileDialog::getExistingDirectory(this,</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> tr("Choose
a directory"), QDir::currentPath());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: green;"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> namesGenerator->SetInputDirectory(directory); <span style="color: green;"></span></span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">const</span>
ReaderType::FileNamesContainer & filenames = </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
namesGenerator->GetInputFileNames();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">unsigned</span> <span style="color: blue;">int</span> numberOfFilenames = filenames.size();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> std::cout << numberOfFilenames << <span style="color: rgb(163, 21, 21);">" images loaded"</span> << std::endl; </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">for</span>(<span style="color: blue;">unsigned</span> <span style="color: blue;">int</span> fni = 0;
fni<numberOfFilenames; fni++)</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> {</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> std::cout << <span style="color: rgb(163, 21, 21);">"filename
# "</span> << fni << <span style="color: rgb(163, 21, 21);">" =
"</span>;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> std::cout << filenames[fni] << std::endl;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> }</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> reader->SetImageIO( gdcmIO );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> reader->SetFileNames( filenames );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: green;">//reader->Update();</span></span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: green;"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: green;"> //Diffusion filter</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">typedef</span>
itk::GradientAnisotropicDiffusionImageFilter<</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> InputImageType, OutputImageType >
FilterType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> FilterType::Pointer filter = FilterType::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> filter->SetInput( reader->GetOutput() );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: green;">// Set parameters for
diffusion filter</span></span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">const</span> <span style="color: blue;">unsigned</span> <span style="color: blue;">int</span>
numberOfIterations = 5;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">const</span> <span style="color: blue;">double</span> timeStep = 0.0625;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">const</span> <span style="color: blue;">double</span> conductance = 3;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> filter->SetNumberOfIterations(
numberOfIterations );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> filter->SetTimeStep( timeStep );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> filter->SetConductanceParameter( conductance );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> filter->Update();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> RescaleFilterType::Pointer rescaler =
RescaleFilterType::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> rescaler->SetOutputMinimum( 0 );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> rescaler->SetOutputMaximum( 255 );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style="text-indent: 12pt;"><span style="font-size: 10pt; font-family: "Courier New";">rescaler->SetInput( filter->GetOutput() );</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal"><span style="font-size: 10pt; font-family: "Courier New";">}</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal"><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: blue;">void</span><span style="font-size: 10pt; font-family: "Courier New";"> SimpleView::slotContour() {</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">long</span> <span style="color: blue;">int</span> seed_points[N_SEEDS][3] = {</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> { 186,
206, 138 },</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> { 138,
218, 124 },</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> { 226,
201, 124 },</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> { 309,
220, 124 },</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> { 147,
305, 83 },</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span><span style="font-size: 10pt; font-family: "Courier New";"> </span><span style="font-size: 10pt; font-family: "Courier New";">{ 111, 228, 83 },</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> { 171,
176, 83 },</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> { 245,
178, 83 },</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> { 202,
243, 83 },</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> { 112,
242, 45 },</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> { 98,
253, 7 } </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
};</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">const</span>
InputPixelType timeThreshold = 85;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> ThresholdingFilterType::Pointer thresholder1 =
ThresholdingFilterType::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> thresholder1->SetLowerThreshold( 0.0
);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> thresholder1->SetUpperThreshold( timeThreshold
);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> thresholder1->SetOutsideValue( 0 );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> thresholder1->SetInsideValue( 255 );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">typedef</span>
itk::GradientMagnitudeRecursiveGaussianImageFilter< </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> InputImageType, </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> InputImageType >
GradientFilterType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">typedef</span>
itk::SigmoidImageFilter<</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> InputImageType, </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> InputImageType >
SigmoidFilterType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> GradientFilterType::Pointer gradientMagnitude =
GradientFilterType::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SigmoidFilterType::Pointer sigmoid =
SigmoidFilterType::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> sigmoid->SetOutputMinimum( 0.0 ); </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> sigmoid->SetOutputMaximum( 1.0 ); </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">typedef</span>
itk::FastMarchingImageFilter< InputImageType, </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> InputImageType >
FastMarchingFilterType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> FastMarchingFilterType::Pointer fastMarching =
FastMarchingFilterType::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> gradientMagnitude->SetInput( </span><span style="font-size: 10pt; font-family: "Courier New";">rescaler</span><span style="font-size: 10pt; font-family: "Courier New";">->GetOutput() );</span><span style="font-size: 10pt; font-family: "Courier New";"> </span><span style="font-size: 10pt; font-family: "Courier New"; color: green;">// using
rescaler from diffusion filter as input </span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> sigmoid->SetInput( gradientMagnitude->GetOutput()
);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> fastMarching->SetInput( sigmoid->GetOutput()
);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> thresholder1->SetInput(
fastMarching->GetOutput() );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: green;">//thresholder1->Update();</span></span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New"; color: green;"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">const</span> <span style="color: blue;">double</span> sigma = 1.3; </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> gradientMagnitude->SetSigma( sigma );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">const</span> <span style="color: blue;">double</span> alpha = -0.05; </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">const</span> <span style="color: blue;">double</span> beta = 0.3; </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> sigmoid->SetAlpha( alpha );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> sigmoid->SetBeta( beta );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">typedef</span>
FastMarchingFilterType::NodeContainer NodeContainer;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">typedef</span>
FastMarchingFilterType::NodeType NodeType;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> NodeContainer::Pointer seeds =
NodeContainer::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> InputImageType::IndexType seedPosition,
seedPosition2; </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">const</span> <span style="color: blue;">double</span> initialDistance = 2.0;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> NodeType node, node1, node2;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">const</span> <span style="color: blue;">double</span> seedValue = -1*initialDistance;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">for</span> (<span style="color: blue;">int</span> i = 0; i<N_SEEDS; i++) {</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> seedPosition[0] = seed_points[i][0];</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> seedPosition[1] = seed_points[i][1];</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> seedPosition[2] = seed_points[i][2];</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> node.SetValue( seedValue );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> node.SetIndex( seedPosition );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> seeds->InsertElement( i, node );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> std::cout << <span style="color: rgb(163, 21, 21);">"Seed
Position "</span> << seedPosition << std::endl;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> }</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> fastMarching->SetTrialPoints( seeds );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> fastMarching->SetOutputSize( </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
reader->GetOutput()->GetBufferedRegion().GetSize() );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">const</span> <span style="color: blue;">double</span> stoppingTime = 100;</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> fastMarching->SetStoppingValue( stoppingTime
);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> ExportFilterType::Pointer itkExporter2 =
ExportFilterType::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> itkExporter2->SetInput(
thresholder1->GetOutput() );</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> vtkImageImport* vtkImporter2 =
vtkImageImport::New(); </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> ConnectPipelines(itkExporter2, vtkImporter2);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> vtkImageViewer2 *imageViewer =
vtkImageViewer2::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">this</span>->ui->qvtkWidget->SetRenderWindow(imageViewer->GetRenderWindow());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> vtkRenderWindowInteractor *iren =
vtkRenderWindowInteractor::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> imageViewer->SetupInteractor(iren);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
imageViewer->SetInput(vtkImporter2->GetOutput());</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> imageViewer->SetColorLevel(127);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> imageViewer->SetColorWindow(255);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> imageViewer->Render();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> imageViewer->GetRenderer()->ResetCamera();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> vtkSliderRepresentation2D *SliderRepres =
vtkSliderRepresentation2D::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">int</span> min =
imageViewer->GetSliceMin();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> <span style="color: blue;">int</span> max =
imageViewer->GetSliceMax();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->SetMinimumValue(min);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->SetMaximumValue(max);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->SetValue(<span style="color: blue;">static_cast</span><<span style="color: blue;">int</span>>((min + max) / 2));</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
SliderRepres->GetPoint1Coordinate()->SetCoordinateSystemToNormalizedDisplay();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
SliderRepres->GetPoint1Coordinate()->SetValue(0.92, 0.9);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";">
SliderRepres->GetPoint2Coordinate()->SetCoordinateSystemToNormalizedDisplay();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->GetPoint2Coordinate()->SetValue(0.92,
0.1);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->SetSliderLength(0.03);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->SetSliderWidth(0.03);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->SetEndCapLength(0.03);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->SetEndCapWidth(0.05);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->SetTubeWidth(0.005);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->SetLabelFormat(<span style="color: rgb(163, 21, 21);">"%3.0lf"</span>);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->SetTitleHeight(0.05);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderRepres->SetLabelHeight(0.05);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> vtkSliderWidget *SliderWidget =
vtkSliderWidget::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderWidget->SetInteractor(iren);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderWidget->SetRepresentation(SliderRepres);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderWidget->KeyPressActivationOff();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderWidget->SetAnimationModeToAnimate();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderWidget->SetEnabled(<span style="color: blue;">true</span>);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> vtkSliderCallback2 *SliderCb =
vtkSliderCallback2::New();</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderCb->SetImageViewer(imageViewer);</span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> SliderWidget->AddObserver(vtkCommand::InteractionEvent,
SliderCb); </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style=""><span style="font-size: 10pt; font-family: "Courier New";"> imageViewer->SetSlice(<span style="color: blue;">static_cast</span><<span style="color: blue;">int</span>>(SliderRepres->GetValue()));</span></p>
<p class="MsoNormal" style="text-indent: 12pt;"><span style="font-size: 10pt; font-family: "Courier New";">imageViewer->SetSliceOrientationToXY();</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal" style="text-indent: 12pt;"><span style="font-size: 10pt; font-family: "Courier New";"> </span></p>
<p class="MsoNormal" style="text-indent: 12pt;"><span style="font-size: 10pt; font-family: "Courier New";">SliderRepres->Delete();</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal" style="text-indent: 12pt;"><span style="font-size: 10pt; font-family: "Courier New";">SliderCb->Delete();</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal" style="text-indent: 12pt;"><span style="font-size: 10pt; font-family: "Courier New";">iren->Start();</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal" style="text-indent: 12pt;"><span style="font-size: 10pt; font-family: "Courier New";">iren->Delete();</span><span style="font-size: 10pt; font-family: "Courier New";"></span></p>
<p class="MsoNormal"><span style="font-size: 10pt; font-family: "Courier New";">}</span><span style="font-size: 11pt; color: rgb(31, 73, 125);"></span></p>
<p class="MsoNormal"><span style="font-size: 11pt; color: rgb(31, 73, 125);"> </span></p>
<div style="border-width: 1pt medium medium; border-style: solid none none; border-color: rgb(181, 196, 223) -moz-use-text-color -moz-use-text-color; padding: 3pt 0in 0in;">
<p class="MsoNormal"><b><span style="font-size: 10pt;" lang="ZH-CN">发件人</span></b><b><span style="font-size: 10pt;">:</span></b><span style="font-size: 10pt;"> Jothy [mailto:<a href="mailto:jothybasu@gmail.com" target="_blank">jothybasu@gmail.com</a>] <br>
<b><span lang="ZH-CN">发送时间</span>:</b> 2010</span><span style="font-size: 10pt; font-family: "Batang","serif";" lang="ZH-CN">년</span><span style="font-size: 10pt;"> 11</span><span style="font-size: 10pt; font-family: "Batang","serif";" lang="ZH-CN">월</span><span style="font-size: 10pt;"> 3</span><span style="font-size: 10pt; font-family: "Batang","serif";" lang="ZH-CN">일</span><span style="font-size: 10pt;" lang="ZH-CN"> </span><span style="font-size: 10pt; font-family: "Batang","serif";" lang="ZH-CN">수요일</span><span style="font-size: 10pt;" lang="ZH-CN"> </span><span style="font-size: 10pt; font-family: "Batang","serif";" lang="ZH-CN">오후</span><span style="font-size: 10pt;"> 9:15<br>
<b><span lang="ZH-CN">收件人</span>:</b> Xiaopeng Yang<br>
<b><span lang="ZH-CN">抄送</span>:</b> vtk<br>
<b><span lang="ZH-CN">主题</span>:</b> Re: [vtkusers] Error message about vtk</span></p>
</div><div><div></div><div class="h5">
<p class="MsoNormal"> </p>
<div>
<p class="MsoNormal" style="margin-bottom: 12pt;">You could post the code
snippet, so that ppl can look into it to help you better!<br>
<br>
Jothy</p>
<div>
<p class="MsoNormal">On Wed, Nov 3, 2010 at 11:54 AM, Xiaopeng Yang <<a href="mailto:yxp233@postech.ac.kr" target="_blank">yxp233@postech.ac.kr</a>> wrote:</p>
<div>
<div>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125);">Dear All,</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125);">About the error message I mentioned in the last email.
Did anybody know what causes the probem?</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125);"> </span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125);">Best regards,</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125);">Yang</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125);"> </span></p>
<div>
<div style="border-width: 1pt medium medium; border-style: solid none none; padding: 3pt 0in 0in; border-color: -moz-use-text-color;">
<p class="MsoNormal"><b><span style="font-size: 10pt;" lang="ZH-CN">发件人</span></b><b><span style="font-size: 10pt;">:</span></b><span style="font-size: 10pt;"> <a href="mailto:vtkusers-bounces@vtk.org" target="_blank">vtkusers-bounces@vtk.org</a> [mailto:<a href="mailto:vtkusers-bounces@vtk.org" target="_blank">vtkusers-bounces@vtk.org</a>]
<b><span lang="ZH-CN">代表 </span></b>Xiaopeng Yang<br>
<b><span lang="ZH-CN">发送时间</span>:</b> 2010</span><span style="font-size: 10pt; font-family: "Batang","serif";" lang="ZH-CN">년</span><span style="font-size: 10pt;"> 11</span><span style="font-size: 10pt; font-family: "Batang","serif";" lang="ZH-CN">월</span><span style="font-size: 10pt;"> 3</span><span style="font-size: 10pt; font-family: "Batang","serif";" lang="ZH-CN">일</span><span style="font-size: 10pt;" lang="ZH-CN"> </span><span style="font-size: 10pt; font-family: "Batang","serif";" lang="ZH-CN">수요일</span><span style="font-size: 10pt;" lang="ZH-CN"> </span><span style="font-size: 10pt; font-family: "Batang","serif";" lang="ZH-CN">오전</span><span style="font-size: 10pt;"> 10:22<br>
<b><span lang="ZH-CN">收件人</span>:</b> 'vtk'<br>
<b><span lang="ZH-CN">主题</span>:</b> [vtkusers] Error message about vtk</span></p>
</div>
</div>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Hi
</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I
got the following error when I run a segmentation algorithm. Could anyone
explain the cause of this error?</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Thanks,</p>
<p class="MsoNormal">Xiaopeng</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal"><img src="cid:image001.png@01CB7B9D.A12B9120" alt="Untitled.png" border="0" height="188" width="501"></p>
</div>
</div>
<p class="MsoNormal" style="margin-bottom: 12pt;"><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 <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the VTK FAQ at: <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a></p>
</div>
<p class="MsoNormal"> </p>
</div>
</div></div></div>
</div>
</blockquote></div><br></div>