function zz=regridmesh(data,xx,yy)
centroids=squeeze(mean(reshape(data.Vertices(data.Faces,:),[size(data.Faces,1),3,3]),2));
AA=scatteredInterpolant(centroids(:,1), centroids(:,2), data.FaceVertexCData);
zz=AA(xx,yy);
end