<DIV>Hi friends,</DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV>I wanted to open the .vtk file using mayavi visualizer.But iam getting no scalar data represented erro when i was trying to watch the isosurface in the modules of mayavi.What changes i should made? </DIV>
<DIV> </DIV>
<DIV>#include <vtkImageImport.h><BR>#include <vtkImageDataGeometryFilter.h><BR>#include <fstream.h><BR>#include <vtkPolyDataWriter.h></DIV>
<DIV>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,8);<BR> imageim->SetWholeExtent(0,255,0,255,0,8);<BR> imageim->SetDataSpacing(1.96,1.96,3.92);<BR> imageim->SetDataOrigin(0.0,0.0,0.0);<BR> imageim->SetDataScalarTypeToUnsignedShort();<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->SetFileTypeToBinary();<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>Yahoo! Mail Mobile<br>
<a href="http://us.rd.yahoo.com/mail_us/taglines/mobile/*http://mobile.yahoo.com/learn/mail">Take Yahoo! Mail with you!</a> Check email on your mobile phone.