数据为课程资料:cjk_ts.dat
包含:经纬度、高度信息(应该是吧)
绘图效果:
title 与 xlabel、ylabel暂且不写上
详细代码如下
%%等高线
data = load('cjk_ts.dat');
lon = data(:, 1);
lat = data(:, 2);
t= data(:, 3);
hold on
%%插值 为绘制等高线做准备
xlon = min(lon):0.2:max(lon);
xlat = min(lat):0.2:max(lat);
xt = griddata(lon, lat, t, xlon, xlat')%%别忘记转置一列,选择最合适的插值方式
%%contour(xlon, xlat, xt)
%%调整图像
%%更改颜色
%%contour(xlon, xlat, xt)
%%等值线的填充
%%contourf(xlon, xlat, xt, 'color', 'r')
%%更改等值线的间隔
%%contourf(xlon, xlat, xt, [10:0.5:15], 'color', 'r')
%%pcolor与contour同级
pcolor(xlon, xlat, xt)
shading interp
colorbar;
%%控制颜色
%%caxis([9, 17])
colormap%%设置绘图的着色方案
%%倒置颜色
%%colormap(flipud(jet))
%%默认几个颜色,数量max = 256
%%colormap(jet(2))
%%标记
%%clabel(c, h, 'color', 'g', 'fontisize', '9')%%manual代表手动加入标记#后期精调
bound;%%插入底图,来自课程资料
%%设置等距坐标轴
axis equal;
%%限制范围
axis([120 125 28 34])
%%四周都有坐标轴
box on
%%发现遮挡
set(gca, 'layer', 'top')
结束