clear clc clf tt = dlmread('ttimes.txt'); % 1 2 3 4 5 6 % ix iy oa xc yc zc ix = tt(:,3)' - 75; xc(ix) = tt(:,4); yc(ix) = tt(:,5); zc(ix) = tt(:,6); xx = reshape(xc, 75, 75); yy = reshape(yc, 75, 75); zz = reshape(zc, 75, 75); surf(xx, yy, zz) shading interp hold on raw = dlmread('atlanta.txt'); getStreets(raw, 1, 'r') getStreets(raw, 2, 'g') getStreets(raw, 3, 'b') getStreets(raw, 4, 'y') getStreets(raw, 5, 'c') getStreets(raw, 6, 'm') axis tight alpha(.7)