[ITK-users] ITK and VTK connection problem
Ashef Shahrior <[email protected]> Thu, 5 Jul 2018 19:52:02 +0600
| Newsgroups | gmane.comp.lib.itk.user |
|---|---|
| Message-ID | <CAB8nWDw+sb+mv76GaV9EtqOHuptjHFM6kSyZ=Qq92DO5yL615Q@mail.gmail.com> |
--===============0231704065==
Content-Type: multipart/alternative; boundary="000000000000ce574b057040dac0"
--000000000000ce574b057040dac0
Content-Type: text/plain; charset="UTF-8"
From ITK Software Guide Book, I am trying to do Region Growing - Connected
Threshold segmentation and show it using VTK. In the guide book, the
segmented image is written into another file. But I don't want to write it
but rather display it using VTK renderer. But it seems to be not working.
Can anyone see my code and help me out, please? Here's my code:
#include <iostream>
#include <string.h>
#include <windows.h>
#include <time.h>
#include <dos.h>
#include "itkConnectedThresholdImageFilter.h"
#include "itkImage.h"
#include "itkCastImageFilter.h"
#include "itkCurvatureFlowImageFilter.h"
#include "itkImageFileReader.h"
#include "itkImageFileWriter.h"
#include "itkVTKImageExport.h"
#include "itkVTKImageImport.h"
#include "itkRGBPixel.h"
#include "vtkImageData.h"
#include "vtkImageImport.h"
#include "vtkImageExport.h"
#include "vtkImageActor.h"
#include "vtkRenderer.h"
#include "vtkRenderWindow.h"
#include "vtkRenderWindowInteractor.h"
#include "vtkInteractorStyleImage.h"
using namespace std;
template <typename ITK_Exporter, typename VTK_Importer> void
ConnectPipelines(ITK_Exporter exporter, VTK_Importer* importer)
{
importer->SetUpdateInformationCallback(exporter->GetUpdateInformationCallback());
importer->SetPipelineModifiedCallback(exporter->GetPipelineModifiedCallback());
importer->SetWholeExtentCallback(exporter->GetWholeExtentCallback());
importer->SetSpacingCallback(exporter->GetSpacingCallback());
importer->SetOriginCallback(exporter->GetOriginCallback());
importer->SetScalarTypeCallback(exporter->GetScalarTypeCallback());
importer->SetNumberOfComponentsCallback(exporter->GetNumberOfComponentsCallback());
importer->SetPropagateUpdateExtentCallback(exporter->GetPropagateUpdateExtentCallback());
importer->SetUpdateDataCallback(exporter->GetUpdateDataCallback());
importer->SetDataExtentCallback(exporter->GetDataExtentCallback());
importer->SetBufferPointerCallback(exporter->GetBufferPointerCallback());
importer->SetCallbackUserData(exporter->GetCallbackUserData());
}
template <typename VTK_Exporter, typename ITK_Importer> void
ConnectPipelines(VTK_Exporter* exporter, ITK_Importer importer)
{
importer->SetUpdateInformationCallback(exporter->GetUpdateInformationCallback());
importer->SetPipelineModifiedCallback(exporter->GetPipelineModifiedCallback());
importer->SetWholeExtentCallback(exporter->GetWholeExtentCallback());
importer->SetSpacingCallback(exporter->GetSpacingCallback());
importer->SetOriginCallback(exporter->GetOriginCallback());
importer->SetScalarTypeCallback(exporter->GetScalarTypeCallback());
importer->SetNumberOfComponentsCallback(exporter->GetNumberOfComponentsCallback());
importer->SetPropagateUpdateExtentCallback(exporter->GetPropagateUpdateExtentCallback());
importer->SetUpdateDataCallback(exporter->GetUpdateDataCallback());
importer->SetDataExtentCallback(exporter->GetDataExtentCallback());
importer->SetBufferPointerCallback(exporter->GetBufferPointerCallback());
importer->SetCallbackUserData(exporter->GetCallbackUserData());
}
void main()
{
char inputImage[1000], seedx[10], seedy[10], lThreshold[100],
uThreshold[100];
cout<< "Input values: "<<endl;
cin >> inputImage;
cin>>seedx;
cin>>seedy;
cin>>lThreshold;
cin>>uThreshold;
const unsigned int Dimension = 2;
typedef float InternalPixelType;
typedef itk::Image< InternalPixelType, Dimension > InternalImageType;
typedef unsigned char OutputPixelType;
typedef itk::Image< OutputPixelType, Dimension > OutputImageType;
typedef itk::CurvatureFlowImageFilter< InternalImageType,
InternalImageType > CurvatureFlowImageFilterType;
CurvatureFlowImageFilterType::Pointer smoothing =
CurvatureFlowImageFilterType::New();
typedef itk::ConnectedThresholdImageFilter< InternalImageType,
InternalImageType > ConnectedFilterType;
ConnectedFilterType::Pointer connectedThreshold =
ConnectedFilterType::New();
typedef itk::CastImageFilter< InternalImageType, OutputImageType >
CastingFilterType;
CastingFilterType::Pointer caster = CastingFilterType::New();
typedef itk::ImageFileReader< InternalImageType > ReaderType;
typedef itk::ImageFileWriter< OutputImageType > WriterType;
ReaderType::Pointer reader = ReaderType::New();
WriterType::Pointer writer = WriterType::New();
reader->SetFileName( inputImage );
reader->Update();
//writer->SetFileName( outputImage );
smoothing->SetInput( reader->GetOutput() );
smoothing->SetNumberOfIterations( 5 );
smoothing->SetTimeStep( 0.125 );
//writer->SetInput( caster->GetOutput() );
const InternalPixelType lowerThreshold = atof( lThreshold );
const InternalPixelType upperThreshold = atof( uThreshold );
connectedThreshold->SetInput( smoothing->GetOutput() );
connectedThreshold->SetLower( lowerThreshold );
connectedThreshold->SetUpper( upperThreshold );
connectedThreshold->SetReplaceValue( 255 );
InternalImageType::IndexType index;
index[0] = atoi( seedx );
index[1] = atoi( seedy );
connectedThreshold->SetSeed( index );
connectedThreshold->Update();
//caster->SetInput( connectedThreshold->GetOutput() );
caster->SetInput(connectedThreshold->GetOutput());
typedef itk::VTKImageExport< InternalImageType > ExportFilterType; //
ExportFilterType::Pointer itkExporter = ExportFilterType::New(); //
//itkExporter->SetInput( caster->GetOutput() ); // this line gives
me error warning, so i did the following line
itkExporter->SetInput(connectedThreshold->GetOutput());
vtkImageImport* vtkImporter = vtkImageImport::New(); //
ConnectPipelines(itkExporter, vtkImporter);
typedef itk::VTKImageImport< InternalImageType > ImportFilterType; //
ImportFilterType::Pointer itkImporter = ImportFilterType::New(); //
vtkImageExport* vtkExporter = vtkImageExport::New();
ConnectPipelines(vtkExporter, itkImporter);
#if VTK_MAJOR_VERSION <= 5
vtkExporter->SetInput( vtkImporter->GetOutput() );
#else
vtkImporter->Update();
vtkExporter->SetInputData( vtkImporter->GetOutput() );
#endif
vtkImageActor* actor = vtkImageActor::New();
#if VTK_MAJOR_VERSION <= 5
actor->SetInput(vtkImporter->GetOutput());
#else
actor->SetInputData(vtkImporter->GetOutput());
#endif
vtkInteractorStyleImage * interactorStyle =
vtkInteractorStyleImage::New();
vtkRenderer* renderer = vtkRenderer::New();
vtkRenderWindow* renWin = vtkRenderWindow::New();
vtkRenderWindowInteractor* iren = vtkRenderWindowInteractor::New();
renWin->SetSize(500, 500);
renWin->AddRenderer(renderer);
iren->SetRenderWindow(renWin);
iren->SetInteractorStyle( interactorStyle );
renderer->AddActor(actor);
renderer->SetBackground(0.4392, 0.5020, 0.5647);
renWin->Render();
iren->Start();
actor->Delete();
interactorStyle->Delete();
vtkImporter->Delete();
vtkExporter->Delete();
renWin->Delete();
renderer->Delete();
iren->Delete();
}
Thanks in advance.
--000000000000ce574b057040dac0
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable
<div dir=3D"ltr">From ITK Software Guide Book, I am trying to do Region Gro=
wing - Connected Threshold segmentation and show it using VTK. In the guide=
book, the segmented image is written into another file. But I don't wa=
nt to write it but rather display it using VTK renderer. But it seems to be=
not working. Can anyone see my code and help me out,=C2=A0please? Here'=
;s my code:<div><br></div><div><div>#include <iostream></div><div>#in=
clude <string.h></div><div>#include <windows.h></div><div>#incl=
ude <time.h></div><div>#include <dos.h></div><div><br></div><di=
v>#include "itkConnectedThresholdImageFilter.h"</div><div><br></d=
iv><div>#include "itkImage.h"</div><div>#include "itkCastIma=
geFilter.h"</div><div><br></div><div>#include "itkCurvatureFlowIm=
ageFilter.h"</div><div><br></div><div>#include "itkImageFileReade=
r.h"</div><div>#include "itkImageFileWriter.h"</div><div><br=
></div><div>#include "itkVTKImageExport.h"</div><div>#include &qu=
ot;itkVTKImageImport.h"</div><div>#include "itkRGBPixel.h"</=
div><div><br></div><div>#include "vtkImageData.h"</div><div>#incl=
ude "vtkImageImport.h"</div><div>#include "vtkImageExport.h&=
quot;</div><div>#include "vtkImageActor.h"</div><div>#include &qu=
ot;vtkRenderer.h"</div><div>#include "vtkRenderWindow.h"</di=
v><div>#include "vtkRenderWindowInteractor.h"</div><div>#include =
"vtkInteractorStyleImage.h"</div><div><br></div><div>using namesp=
ace std;</div><div><br></div><div>template <typename ITK_Exporter, typen=
ame VTK_Importer> void ConnectPipelines(ITK_Exporter exporter, VTK_Impor=
ter* importer)</div><div>{</div><div>=C2=A0 =C2=A0 importer->SetUpdateIn=
formationCallback(exporter->GetUpdateInformationCallback());</div><div>=
=C2=A0 =C2=A0 importer->SetPipelineModifiedCallback(exporter->GetPipe=
lineModifiedCallback());</div><div>=C2=A0 =C2=A0 importer->SetWholeExten=
tCallback(exporter->GetWholeExtentCallback());</div><div>=C2=A0 =C2=A0 i=
mporter->SetSpacingCallback(exporter->GetSpacingCallback());</div><di=
v>=C2=A0 =C2=A0 importer->SetOriginCallback(exporter->GetOriginCallba=
ck());</div><div>=C2=A0 =C2=A0 importer->SetScalarTypeCallback(exporter-=
>GetScalarTypeCallback());</div><div>=C2=A0 =C2=A0 importer->SetNumbe=
rOfComponentsCallback(exporter->GetNumberOfComponentsCallback());</div><=
div>=C2=A0 =C2=A0 importer->SetPropagateUpdateExtentCallback(exporter-&g=
t;GetPropagateUpdateExtentCallback());</div><div>=C2=A0 =C2=A0 importer->=
;SetUpdateDataCallback(exporter->GetUpdateDataCallback());</div><div>=C2=
=A0 =C2=A0 importer->SetDataExtentCallback(exporter->GetDataExtentCal=
lback());</div><div>=C2=A0 =C2=A0 importer->SetBufferPointerCallback(exp=
orter->GetBufferPointerCallback());</div><div>=C2=A0 =C2=A0 importer->=
;SetCallbackUserData(exporter->GetCallbackUserData());</div><div>}</div>=
<div><br></div><div>template <typename VTK_Exporter, typename ITK_Import=
er> void ConnectPipelines(VTK_Exporter* exporter, ITK_Importer importer)=
</div><div>{</div><div>=C2=A0 =C2=A0 importer->SetUpdateInformationCallb=
ack(exporter->GetUpdateInformationCallback());</div><div>=C2=A0 =C2=A0 i=
mporter->SetPipelineModifiedCallback(exporter->GetPipelineModifiedCal=
lback());</div><div>=C2=A0 =C2=A0 importer->SetWholeExtentCallback(expor=
ter->GetWholeExtentCallback());</div><div>=C2=A0 =C2=A0 importer->Set=
SpacingCallback(exporter->GetSpacingCallback());</div><div>=C2=A0 =C2=A0=
importer->SetOriginCallback(exporter->GetOriginCallback());</div><di=
v>=C2=A0 =C2=A0 importer->SetScalarTypeCallback(exporter->GetScalarTy=
peCallback());</div><div>=C2=A0 =C2=A0 importer->SetNumberOfComponentsCa=
llback(exporter->GetNumberOfComponentsCallback());</div><div>=C2=A0 =C2=
=A0 importer->SetPropagateUpdateExtentCallback(exporter->GetPropagate=
UpdateExtentCallback());</div><div>=C2=A0 =C2=A0 importer->SetUpdateData=
Callback(exporter->GetUpdateDataCallback());</div><div>=C2=A0 =C2=A0 imp=
orter->SetDataExtentCallback(exporter->GetDataExtentCallback());</div=
><div>=C2=A0 =C2=A0 importer->SetBufferPointerCallback(exporter->GetB=
ufferPointerCallback());</div><div>=C2=A0 =C2=A0 importer->SetCallbackUs=
erData(exporter->GetCallbackUserData());</div><div>}</div><div><br></div=
><div><br></div><div>void main()</div><div>{</div><div>=C2=A0 =C2=A0 char i=
nputImage[1000], seedx[10], seedy[10], lThreshold[100], uThreshold[100];</d=
iv><div><br></div><div>=C2=A0 =C2=A0 cout<< "Input values: "=
;<<endl;</div><div>=C2=A0 =C2=A0 cin >> inputImage;</div><div>=
=C2=A0 =C2=A0 cin>>seedx;</div><div>=C2=A0 =C2=A0 cin>>seedy;</=
div><div>=C2=A0 =C2=A0 cin>>lThreshold;</div><div>=C2=A0 =C2=A0 cin&g=
t;>uThreshold;</div><div><br></div><div>=C2=A0 =C2=A0 const=C2=A0 =C2=A0=
=C2=A0unsigned int=C2=A0 =C2=A0 Dimension =3D 2;</div><div><br></div><div>=
=C2=A0 =C2=A0 typedef=C2=A0 =C2=A0float=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =
=C2=A0InternalPixelType;</div><div>=C2=A0 =C2=A0 typedef itk::Image< Int=
ernalPixelType, Dimension >=C2=A0 InternalImageType;</div><div><br></div=
><div>=C2=A0 =C2=A0 typedef unsigned char=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0=
=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 OutputPixel=
Type;</div><div>=C2=A0 =C2=A0 typedef itk::Image< OutputPixelType, Dimen=
sion > OutputImageType;</div><div><br></div><div>=C2=A0 =C2=A0 typedef i=
tk::CurvatureFlowImageFilter< InternalImageType, InternalImageType > =
CurvatureFlowImageFilterType;</div><div><br></div><div>=C2=A0 =C2=A0 Curvat=
ureFlowImageFilterType::Pointer smoothing =3D CurvatureFlowImageFilterType:=
:New();</div><div><br></div><div>=C2=A0 =C2=A0 typedef itk::ConnectedThresh=
oldImageFilter< InternalImageType, InternalImageType > ConnectedFilte=
rType;</div><div><br></div><div>=C2=A0 =C2=A0 ConnectedFilterType::Pointer =
connectedThreshold =3D ConnectedFilterType::New();</div><div><br></div><div=
>=C2=A0 =C2=A0 typedef itk::CastImageFilter< InternalImageType, OutputIm=
ageType > CastingFilterType;</div><div><br></div><div>=C2=A0 =C2=A0 Cast=
ingFilterType::Pointer caster =3D CastingFilterType::New();</div><div><br><=
/div><div>=C2=A0 =C2=A0 typedef=C2=A0 itk::ImageFileReader< InternalImag=
eType > ReaderType;</div><div>=C2=A0 =C2=A0 typedef=C2=A0 itk::ImageFile=
Writer<=C2=A0 OutputImageType=C2=A0 > WriterType;</div><div><br></div=
><div>=C2=A0 =C2=A0 ReaderType::Pointer reader =3D ReaderType::New();</div>=
<div>=C2=A0 =C2=A0 WriterType::Pointer writer =3D WriterType::New();</div><=
div><br></div><div>=C2=A0 =C2=A0 reader->SetFileName( inputImage );</div=
><div>=C2=A0 =C2=A0 reader->Update();</div><div>=C2=A0 =C2=A0 //writer-&=
gt;SetFileName( outputImage );</div><div><br></div><div>=C2=A0 =C2=A0 smoot=
hing->SetInput( reader->GetOutput() );</div><div>=C2=A0 =C2=A0 smooth=
ing->SetNumberOfIterations( 5 );</div><div>=C2=A0 =C2=A0 smoothing->S=
etTimeStep( 0.125 );</div><div><br></div><div>=C2=A0 =C2=A0 //writer->Se=
tInput( caster->GetOutput() );</div><div>=C2=A0 =C2=A0 const InternalPix=
elType lowerThreshold =3D atof( lThreshold );</div><div>=C2=A0 =C2=A0 const=
InternalPixelType upperThreshold =3D atof( uThreshold );</div><div><br></d=
iv><div>=C2=A0 =C2=A0 connectedThreshold->SetInput( smoothing->GetOut=
put() );</div><div>=C2=A0 =C2=A0 connectedThreshold->SetLower(=C2=A0 low=
erThreshold=C2=A0 );</div><div>=C2=A0 =C2=A0 connectedThreshold->SetUppe=
r(=C2=A0 upperThreshold=C2=A0 );</div><div><br></div><div>=C2=A0 =C2=A0 con=
nectedThreshold->SetReplaceValue( 255 );</div><div><br></div><div>=C2=A0=
=C2=A0 InternalImageType::IndexType=C2=A0 index;</div><div><br></div><div>=
=C2=A0 =C2=A0 index[0] =3D atoi( seedx );</div><div>=C2=A0 =C2=A0 index[1] =
=3D atoi( seedy );</div><div><br></div><div>=C2=A0 =C2=A0 connectedThreshol=
d->SetSeed( index );</div><div>=C2=A0 =C2=A0 connectedThreshold->Upda=
te();</div><div><br></div><div>=C2=A0 =C2=A0 //caster->SetInput( connect=
edThreshold->GetOutput() );</div><div><span style=3D"white-space:pre"> <=
/span>caster->SetInput(connectedThreshold->GetOutput());</div><div><b=
r></div><div>=C2=A0 =C2=A0 typedef itk::VTKImageExport< InternalImageTyp=
e > ExportFilterType; //</div><div>=C2=A0 =C2=A0 ExportFilterType::Point=
er itkExporter =3D ExportFilterType::New(); //</div><div><br></div><div>=C2=
=A0 =C2=A0 =C2=A0 =C2=A0 //itkExporter->SetInput( caster->GetOutput()=
); // this line gives me error warning, so i did the following line</div><=
div><span style=3D"white-space:pre"> </span>itkExporter->SetInput(connec=
tedThreshold->GetOutput());</div><div><br></div><div><br></div><div><br>=
</div><div>=C2=A0 =C2=A0 vtkImageImport* vtkImporter =3D vtkImageImport::Ne=
w();=C2=A0 //</div><div>=C2=A0 =C2=A0 ConnectPipelines(itkExporter, vtkImpo=
rter);</div><div><br></div><div>=C2=A0 =C2=A0 typedef itk::VTKImageImport&l=
t; InternalImageType > ImportFilterType; //</div><div>=C2=A0 =C2=A0 Impo=
rtFilterType::Pointer itkImporter =3D ImportFilterType::New(); //</div><div=
><br></div><div><br></div><div>=C2=A0 =C2=A0 vtkImageExport* vtkExporter =
=3D vtkImageExport::New();</div><div>=C2=A0 =C2=A0 ConnectPipelines(vtkExpo=
rter, itkImporter);</div><div><br></div><div>#if VTK_MAJOR_VERSION <=3D =
5</div><div>=C2=A0 =C2=A0 vtkExporter->SetInput( vtkImporter->GetOutp=
ut() );</div><div>#else</div><div>=C2=A0 =C2=A0 vtkImporter->Update();</=
div><div>=C2=A0 =C2=A0 vtkExporter->SetInputData( vtkImporter->GetOut=
put() );</div><div>#endif</div><div><br></div><div>=C2=A0 =C2=A0 vtkImageAc=
tor* actor =3D vtkImageActor::New();</div><div>#if VTK_MAJOR_VERSION <=
=3D 5</div><div>=C2=A0 =C2=A0 actor->SetInput(vtkImporter->GetOutput(=
));</div><div>#else</div><div>=C2=A0 =C2=A0 actor->SetInputData(vtkImpor=
ter->GetOutput());</div><div>#endif</div><div><br></div><div>=C2=A0 =C2=
=A0 vtkInteractorStyleImage * interactorStyle =3D vtkInteractorStyleImage::=
New();</div><div><br></div><div>=C2=A0 =C2=A0 vtkRenderer* renderer =3D vtk=
Renderer::New();</div><div>=C2=A0 =C2=A0 vtkRenderWindow* renWin =3D vtkRen=
derWindow::New();</div><div>=C2=A0 =C2=A0 vtkRenderWindowInteractor* iren =
=3D vtkRenderWindowInteractor::New();</div><div><br></div><div>=C2=A0 =C2=
=A0 renWin->SetSize(500, 500);</div><div>=C2=A0 =C2=A0 renWin->AddRen=
derer(renderer);</div><div>=C2=A0 =C2=A0 iren->SetRenderWindow(renWin);<=
/div><div>=C2=A0 =C2=A0 iren->SetInteractorStyle( interactorStyle );</di=
v><div><br></div><div>=C2=A0 =C2=A0 renderer->AddActor(actor);</div><div=
>=C2=A0 =C2=A0 renderer->SetBackground(0.4392, 0.5020, 0.5647);</div><di=
v><br></div><div>=C2=A0 =C2=A0 renWin->Render();</div><div>=C2=A0 =C2=A0=
iren->Start();</div><div><br></div><div>=C2=A0 =C2=A0 actor->Delete(=
);</div><div>=C2=A0 =C2=A0 interactorStyle->Delete();</div><div>=C2=A0 =
=C2=A0 vtkImporter->Delete();</div><div>=C2=A0 =C2=A0 vtkExporter->De=
lete();</div><div>=C2=A0 =C2=A0 renWin->Delete();</div><div>=C2=A0 =C2=
=A0 renderer->Delete();</div><div>=C2=A0 =C2=A0 iren->Delete();</div>=
<div>}</div></div><div><br></div><div>Thanks in advance.</div></div>
--000000000000ce574b057040dac0--
--===============0231704065==
Content-Type: text/plain; charset="us-ascii"
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Disposition: inline
The ITK community is transitioning from this mailing list to discourse.itk.org. Please join us there!
________________________________
Powered by www.kitware.com
Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html
Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.php
Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ
Follow this link to subscribe/unsubscribe:
https://itk.org/mailman/listinfo/insight-users
--===============0231704065==--