✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
⛄ 内容介绍
地形特征点散点,规则网格,三角网格及等高线等多种地形数据的存储形式及数据结构特殊性很适合matlab的数据处理特点.本文研究如何在matlab中通过简单的编程,从而对地形数据进行必要的处理,以满足地形测量要求,并给出了数据处理实例.
⛄ 完整代码
% 地形海拔高度云图绘制
clear all;close all;
%================================================================
alldata=load('data.asc');
xmin=min(alldata(:,1));xmax=max(alldata(:,1));
ymin=min(alldata(:,2));ymax=max(alldata(:,2));
%================================================================
% Method1:griddata插值 linear,cubic,natural,nearest,v4
xl=linspace(xmin,xmax,101)';
yl=linspace(ymin,ymax,101)';
[X,Y,Fxy]=griddata(alldata(:,1),alldata(:,2),alldata(:,3),xl',yl,'linear'); % v4会超内存报错
%================================================================
% 绘图
[C,h]=contourf(X,Y,Fxy,20,'w-'); % 等高线云图
shading flat; % 平滑过渡
colorbar; % 调出颜色对应的bar
caxis([1050 1550]); % 设置bar的数值范围
colorbar('ytick',[1050:50:1550]); % 设置bar的刻度值
set(gca,'FontName','Times New Roman','FontSize',12);
set(gca, 'PlotBoxAspectRatio',[(xmax-xmin)/(ymax-ymin) 1 1],'XLim',[xmin xmax],'YLim',[ymin ymax],'ZLim',[0 1]);
ab=15;xlabel('{\itx}/m','FontName','Times New Roman','fontsize',ab);
ylabel('{\ity}/m','FontName','Times New Roman','fontsize',ab);
title('地形海拔高度云图','fontsize',ab);
%================================================================
⛄ 运行结果
⛄ 参考文献
[1]杨韫澜, 胡海彦, 陈虹,等. 基于Matlab的地形数据处理[J]. 测绘科学与工程, 2010(1):4.