地形图的三维可视化
clear;clc
a=imread('c0808tin_Tin21.tif');%%读取tiff影像数据，这里只能读取到值，其他的坐标等信息需要后面人为添加
[m,n] = size(a); %%获取数据的行列号
a = a(22:112,1:n); %%根据需要删除数据中无效的行列号
[m,n] = size(a);
a = a(1:3:m,1:4:n); %%均匀降采样
b = find(a < -20);
a(b(1:end)) = NaN; %%将极小值赋值为无效值
c=smoothdata(a,'gaussian',9); %%高斯平滑
e = find(c == 0);
c(e(1:end)) = NaN;
meshc(c)
% % aminn=min(min(a))
% % amaxx=max(max(a))
% meshc(d)%%绘制带有等高线的网格图，这里只是用一个变量，所以绘制结果的x，y范围为0-300，若需要显示经纬度则需定义与数据量相同的x，y向量，使用meshc(x,y,c)绘图
% caxis([-0.025,0.015]);%%调整一下colorbar的范围