VisitMatlab: read2DScalarVTK.m

function data=read2DScalarVTK(filename)
fid=fopen(filename);
% search for POINTS Section
line=fsearch(fid,'POINTS');
npoints=sscanf(line,'POINTS %d');

% Read in points
data.Vertices=zeros(3,npoints);
data.Vertices(:)=fscanf(fid,'%f',npoints*3);

% Search for POLYGONS section
line=fsearch(fid,'POLYGONS');
npolygons=sscanf(line,'POLYGONS %d');

% read in POLYGON data
data.Faces=zeros(4,npolygons,'int64');
data.Faces(:)=fscanf(fid,'%d',npolygons*4);
data.Faces=data.Faces(2:4,:); %remove the leading number 3 - assuming mesh is trianglular
data.Faces=data.Faces+1; %matlab complains otherwise

% Search for Scalar data
line=fsearch(fid,'CELL_DATA');
ncells=sscanf(line,'CELL_DATA %d');
class=fscanf(fid,'%s',1);
var=fscanf(fid,'%s',1);
type=fscanf(fid,'%s',1);
line=fsearch(fid,'LOOKUP_TABLE');

% Read in Scalar date
data.FaceVertexCData=zeros(1,ncells);
data.FaceVertexCData(:)=fscanf(fid,'%f',ncells);

% Transpose things to order expected by patch command
data.Faces=data.Faces';
data.Vertices=data.Vertices';
data.FaceVertexCData=data.FaceVertexCData';
end

function line=fsearch(fid,str)
while true
    line=fgetl(fid);
    if numel(line) >= numel(str)
        if (strcmp(line(1:numel(str)),str))
           break
        end
    end
end
end