<DIV><BR>hi,users,</DIV>
<DIV> I wanted to use ITK read PNG image and segement it,then through VTK and show on the VC's dialog.my program can through complie but can't show the result.Next It is my code .<EM>Does anybody could give me some help?</EM><I> thanks .</I></DIV>
<DIV><I> </I><I> sincerely mexiao</I></DIV>
<DIV>class CVTK <BR>{<BR>public:<BR> vtkWin32RenderWindowInteractor *iren;<BR> vtkImageViewer *viewer;</DIV>
<DIV> typedef itk::Image< float, 2 > InternalImageType1;<BR> <BR> typedef itk::Image< unsigned char, 2 > OutputImageType1;<BR> typedef itk::ImageFileReader< InternalImageType1 > ReaderType1;<BR> typedef itk::ImageFileWriter< OutputImageType1 > WriterType1;<BR> typedef itk::CastImageFilter< InternalImageType1, OutputImageType1 > CastingFilterType;<BR> typedef itk::CurvatureFlowImageFilter< InternalImageType1, InternalImageType1 > CurvatureFlowImageFilterType;<BR> typedef itk::ConnectedThresholdImageFilter< InternalImageType1, InternalImageType1> ConnectedFilterType;<BR> float lowerThreshold ;<BR> float upperThreshold ;<BR> InternalImageType1::IndexType index;<BR> typedef itk::Image< unsigned char, 2 > ImageType;<BR> typedef itk::ImageToVTKImageFilter<ImageType> ConnectorType;</DIV>
<DIV>public:<BR> vtkCellPicker * picker;<BR> static void MyExecute(vtkObject *caller, unsigned long eid, void *clientdata, void *calldata);<BR> void BeginRenderOn(CStatic* aStatic);<BR> CVTK();<BR> virtual ~CVTK();<BR>};</DIV>
<DIV>//////////////////////////////////////////////////////////////////////////////////////////</DIV>
<DIV>CVTK::CVTK()<BR>{ <BR> lowerThreshold = 150.0;<BR> upperThreshold = 180.0;<BR> index[0] =60;<BR> index[1] =116;<BR> ReaderType1::Pointer reader = ReaderType1::New();<BR> reader->SetFileName("G:\\ITKVTK\\test\\test\\myProject\\Debug\\BrainProtonDensitySlice.png");<BR> WriterType1::Pointer writer = WriterType1::New();<BR> writer->SetFileName("G:\\ITKVTK\\test\\test\\myProject\\Debug\\aa.png");<BR> CurvatureFlowImageFilterType::Pointer smoothing = CurvatureFlowImageFilterType::New();<BR> ConnectedFilterType::Pointer connectedThreshold = ConnectedFilterType::New();<BR> CastingFilterType::Pointer caster = CastingFilterType::New();<BR> smoothing->SetInput( reader->GetOutput() );<BR> connectedThreshold->SetInput( smoothing->GetOutput() );<BR> caster->SetInput( connectedThreshold->GetOutput() );<BR> writer->SetInput( caster->GetOutput() );<BR><BR> smoothing->SetNumberOfIterations( 5 );<BR> smoothing->SetTimeStep( 0.125 );<BR> connectedThreshold->SetLower( lowerThreshold );<BR> connectedThreshold->SetUpper( upperThreshold );<BR> connectedThreshold->SetReplaceValue( 255 );<BR> connectedThreshold->SetSeed( index );<BR> ConnectorType::Pointer connector = ConnectorType::New();<BR> connector->SetInput( caster->GetOutput()); <BR> <BR>// writer->Update();<BR><BR> viewer = <SPAN class="t_tag" onclick="tagshow(event)" href="tag.php?name=vtk">vtk</SPAN>ImageViewer::New();<BR> viewer->SetInput(connector->GetOutput());<BR> viewer->SetColorWindow(256);<BR> viewer->SetColorLevel(128);<BR><BR>iren = vtkWin32RenderWindowInteractor::New();<BR> viewer->SetupInteractor(iren);<BR> // connector->Update();<BR>picker=vtkCellPicker::New();<BR> //picker->SetTolerance(0.00001);<BR><BR>vtkCallbackCommand * cmd = vtkCallbackCommand::New();<BR>cmd->SetCallback(CVTK::MyExecute);<BR>cmd->SetClientData( this );<BR>picker->AddObserver(vtkCommand::EndPickEvent, cmd);<BR>iren->SetPicker(picker);<BR>cmd->Delete(); <BR><BR>}</DIV>
<DIV>/////////////////////////////////////////////////////<BR>BOOL CDlgVTKITKDlg::OnInitDialog()<BR>{<BR> CDialog::OnInitDialog();<BR> SetIcon(m_hIcon, TRUE); // Set big icon<BR> SetIcon(m_hIcon, FALSE); // Set small icon<BR> // TODO: Add extra initialization here<BR> this->m_vtk.BeginRenderOn(&(this->m_RenderArea));//have been defined CVTK m_vtk;<BR> UpdateData(FALSE);<BR> <BR> return TRUE; // return TRUE unless you set the focus to a control<BR>}<BR>///////////////////////////////////////////////////////////////////////<BR>void CVTK::BeginRenderOn(CStatic * aStatic)<BR>{<BR>CRect rect;<BR>aStatic->GetClientRect(&rect);<BR>viewer->SetSize(rect.Width(),rect.Height());<BR>viewer->SetParentId(aStatic->m_hWnd); </DIV>
<DIV>viewer->Render(); </DIV>
<DIV>// iren->Initialize();<BR>}</DIV>
<DIV>/////////////////////////////////////////////////////////////////<BR>void CVTK::MyExecute(vtkObject *caller, unsigned long eid, void *clientdata, void *calldata)<BR>{<BR>CVTK *m_viewer = reinterpret_cast<CVTK*>(clientdata );<BR> if (m_viewer->picker->GetCellId()<0)<BR>{<BR> AfxMessageBox("CellID < 0");<BR>}<BR> else<BR>{<BR> AfxMessageBox("CellID > 0");<BR> double ScreenCoordinate[3],globalCoordinate[3]; </DIV>
<DIV> m_viewer->picker->GetSelectionPoint(ScreenCoordinate); <BR> m_viewer->picker->GetPickPosition(globalCoordinate); </DIV>
<DIV>};<BR>}<BR></DIV><br><!-- footer --><br><span title="neteasefooter"/><hr/>
<a href="http://www.yeah.net">网易邮箱，中国第一大电子邮件服务商</a>
</span>