I think: 1. It should be "DATASET UNSTRUCTURED_GRID" not "DATASET STRUCTURED_GRID" 2. In CELLS it should be NC*5 not NC*4, because this represents the total number of integer values in the CELL list. In your case it is (Number of elements*(1 + Number of nodes per element)) 3. In the CELL list you should have Number of nodes per element, node1, node2, node3, node4 ... ... Therefore, it should be: fprintf(fp,'%d %d %d %d %d\n',4,Cells(1,i)-1,Cells(2,i)-1,Cells(3,i)-1,Cells(4,i)-1);
4. Because you are displaying Quadrilateral, in the list of CELL_TYPE you should use: fprintf(fp,'%d\n',9); not fprintf(fp,'%d\n',4); 4 is for polyline see figure 3 of (5) 5. See http://www.vtk.org/VTK/img/file-formats.pdf for detail Good luck! HNG On Fri, Nov 19, 2010 at 7:45 AM, 乔彦超 <oliver.q...@gmail.com> wrote: > Hi: > > i want write my results into VTK format, then i can use the paraview. > but i have some problem. my compute model is 700KM*700KM, there are 71*71 > nodes. it's a Quadrilateral mesh. > > my matlab code as followed > % Model size, m > xsize=700000; > ysize=700000; > % Defining resolution > xnum=71; > ynum=71; > NP=xnum*ynum; > NC=xnum*ynum-2*xnum+1; > % Defining gridline positions for regular basic grid > % In horizontal direction > gridx=zeros(xnum,1); > for i=2:1:xnum > gridx(i)=gridx(i-1)+xstp; > end > % Defining gridline positions for 51x51 irregular basic grid > % In vertical direction > gridy=zeros(ynum,1); > for i=2:1:6 > gridy(i)=gridy(i-1)+20000; > end > for i=7:1:46 > gridy(i)=gridy(i-1)+1250; > end > for i=46:1:51 > gridy(i)=gridy(i-1)+10000; > end > for i=51:1:71 > gridy(i)=gridy(i-1)+25000; > end > gridy(ynum)=ysize; > > %write the points > a=0; > for i=1:1:ynum > for j=1:1:xnum > a=a+1; > Points(1,a)=gridx(j); > Points(2,a)=gridy(i); > end > end > % write the cell > nn=0; > for j=1:1:ynum > for i=1:1:xnum > nn=nn+1; > Cells(1,nn)=(j-1)*(xnum+1)+i; > Cells(2,nn)=(j-1)*(xnum+1)+i+1; > Cells(3,nn)=j*(xnum+1)+i+1; > Cells(4,nn)=j*(xnum+1)+i; > Cells(5,nn)=nn; > end > end > % write VTK file > fp=fopen('Resault.vtk','w'); > fprintf(fp,'# vtk DataFile Version 3.0\n'); > fprintf(fp,'vtk output\n'); > fprintf(fp,'ASCII\n'); > fprintf(fp,'DATASET STRUCTURED_GRID\n'); > fprintf(fp,'POINTS %d float\n',NP); > for i=1:NP > fprintf(fp,'%e %e %e \n',Points(1,i),Points(2,i),0); > end > fprintf(fp,'\n'); > fprintf(fp,'CELLS %d %d\n',NC,NC*4); > aa=0; > for i=1:NC > aa=aa+1; > fprintf(fp,'%d %d %d %d > %d\n',aa,Cells(1,i)-1,Cells(2,i)-1,Cells(3,i)-1,Cells(4,i)-1); > end > fprintf(fp,'\n'); > fprintf(fp,'CELL_TYPES %d\n',NC); > for i=1:NC > fprintf(fp,'%d\n',4); > end > fprintf(fp,'\n'); > fprintf(fp,'POINT_DATA %d\n',NP); > > when i tried to open it with paraview. > it sayed "ERROR: In > C:\berk\ReleaseBuilds\ParaView3\VTK\IO\vtkDataReader.cxx, line 2078 > ERROR: In C:\berk\ReleaseBuilds\ParaView3\VTK\IO\vtkDataReader.cxx, line > 2078 > ERROR: In > C:\berk\ReleaseBuilds\ParaView3\VTK\Parallel\vtkPDataSetReader.cxx, line 705 > vtkPDataSetReader (0AEA4D18): Expecting 'DIMENSIONS' insted of: POINTS > ERROR: In > C:\berk\ReleaseBuilds\ParaView3\VTK\Parallel\vtkPDataSetReader.cxx, line 705 > vtkPDataSetReader (0AEA4D18): Expecting 'DIMENSIONS' insted of: POINTS" > > What's the problem? could you help me? > > thank you very much. > > Best regards! > > -- > Yanchao Qiao > Ph.D student of Geodynamics > > +86-136 010 926 41 > > College of Earth Science > Graduate University of Chinese Academy of Sciences > > 19A, YuQuan Road, Beijing, China, 100049 > Email:oliver.q...@163.com <email%3aoliver.q...@163.com>; > oliver.q...@gmail.com;qiaoyancha...@mails.gucas.ac.cn > > _______________________________________________ > Powered by www.kitware.com > > Visit other Kitware open-source projects at > http://www.kitware.com/opensource/opensource.html > > Please keep messages on-topic and check the ParaView Wiki at: > http://paraview.org/Wiki/ParaView > > Follow this link to subscribe/unsubscribe: > http://www.paraview.org/mailman/listinfo/paraview > >
_______________________________________________ Powered by www.kitware.com Visit other Kitware open-source projects at http://www.kitware.com/opensource/opensource.html Please keep messages on-topic and check the ParaView Wiki at: http://paraview.org/Wiki/ParaView Follow this link to subscribe/unsubscribe: http://www.paraview.org/mailman/listinfo/paraview