<DIV>Hi Vtk users,</DIV>
<DIV>I wanted to convert an image into polydata and i wanted store it in .vtk format and i even wanted to open the corresponding .vtk file in Mayavi data visualizer.But iam getting the access voilation error by writing the following code.</DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV>#include <vtkImageImport.h><BR>#include <vtkImageDataGeometryFilter.h><BR>#include <fstream.h><BR>#include <vtkPolyDataWriter.h><BR>#include <iostream.h><BR>class dcmvtk<BR>{<BR>private:<BR> vtkImageImport *imageim;<BR> vtkImageDataGeometryFilter *gfil;<BR> vtkPolyDataWriter *vw;<BR> const char* s;<BR>public:<BR> dcmvtk();<BR> ~dcmvtk();<BR>};<BR>dcmvtk::dcmvtk()<BR>{<BR> unsigned char* p = NULL;<BR> if(p != NULL)<BR> {<BR> delete []p;<BR> p = NULL;<BR> }<BR> p = new unsigned char[256*256*2*9];<BR> s="D:/vtk/PET_after_preFilter.raw";<BR> ifstream is(s,ios::in | ios::out | ios::binary);<BR> is.read(p,256*256*2*9);<BR> is.close();<BR>
imageim=vtkImageImport::New();<BR> imageim->SetImportVoidPointer((void*)p,1);<BR> imageim->SetDataExtent(0,255,0,255,0,9);<BR> imageim->SetWholeExtent(0,255,0,255,0,9);<BR> imageim->SetDataSpacing(1.96,1.96,3.92);<BR> imageim->SetDataOrigin(0.0,0.0,0.0);<BR> gfil=vtkImageDataGeometryFilter::New();<BR> gfil->SetInput(imageim->GetOutput());<BR> ofstream ofile;<BR> <BR> const char* f=new char[256*256*9*2];<BR> f="d:/dcm2.vtk";<BR> ofile.open(f,ios::in | ios::out | ios::binary);<BR> vw=vtkPolyDataWriter::New();</DIV>
<DIV> vw->SetInput(gfil->GetOutput());<BR> vw->SetFileName(f);<BR> vw->Write();<BR>}<BR>dcmvtk::~dcmvtk()<BR>{<BR> imageim->Delete();<BR> gfil->Delete();<BR> vw->Delete();<BR>}<BR>int main(int argc, char *argv[])<BR>{<BR> dcmvtk dc;<BR> return 0;<BR>}</DIV><p>
<hr size=1>Do you Yahoo!?<br>
Yahoo! Mail - Find what you need with new enhanced search. <a href="http://us.rd.yahoo.com/evt=29917/*http://info.mail.yahoo.com/mail_250">Learn more.</a>