matlab中已知X，Y，Z坐标，绘制等高线图

本文采用的数据是2011年国赛的题目，使用每一点的x,y坐标，将其对应的每一元素的浓度作为z，绘制金属污染物的二维等高线图

以As为例，以下两个图形还使用了插值算法

load data%每一点坐标信息，共有319个点
load data1%浓度信息，有319个数据
As=data1(:,1);    %1为第一种重金属元素As，若想要其他重金属的，直接改为相应序号
X=data(:,1); %横坐标
Y=data(:,2); %纵坐标
[x y]=meshgrid(0:100:3e4,0:100:2e4); %根据坐标范围划分，作为插值使用
z=griddata(X,Y,As,x,y,'v4');%该命令实现插值，v4是插值方法
C=contour(x,y,z);
title('As')

绘制出来的图形为

还可以绘制整个城区的地形图，就是5个功能区在整个城区的分布情况。
其中，等高线的x,y,z坐标分别代表选取点的x,y,和海拔

function draw1 ()
clc,clear%一是生活区、二是工业区、三是山区、四是交通区、五是公园绿地区%画功能区
load data
load data1
[c1 d1]=find(data(:,4)==1); %即找到功能区1,返回的c1是每个数据对应的位置编号，即行标号；d1就是1
x1=data(c1,1); %功能区1的横坐标
y1=data(c1,2); %功能区1的纵坐标
plot(x1,y1,'r*');
[c2 d2]=find(data(:,4)==2); %同上
x2=data(c2,1);
y2=data(c2,2);
hold on
plot(x2,y2,'kd','markerfacecolor','k')
hold on
[c3 d3]=find(data(:,4)==3);
x3=data(c3,1);
y3=data(c3,2);
plot(x3,y3,'cs','markerfacecolor','c')
hold on
[c4 d4]=find(data(:,4)==4);
x4=data(c4,1);
y4=data(c4,2);
plot(x4,y4,'bp','markerfacecolor','b')
hold on
[c5 d5]=find(data(:,4)==5);
x5=data(c5,1);
y5=data(c5,2);
plot(x5,y5,'go','markerfacecolor','g')
grid on
legend('生活区','工业区','山区','交通区','公园绿地区');
hold on%画等高线
XYZ=data(:,1:end-1);
XYZ=XYZ';
[x,y]=meshgrid(0:100:3e4,0:100:2e4);
z=griddata(XYZ(1,:),XYZ(2,:),XYZ(3,:),x,y,'v4');%插值
C=contour(x,y,z)
hold on

绘制图形如下

matlab中已知X，Y，Z坐标，绘制等高线图相关推荐

matlab状态空间法算反馈阵,matlab中已知系统的状态方程怎样绘制系统阶跃响应曲线...
matlab中,怎样设置坐标系统的原点在左上角? axisijx=-pi:0.1:pi;y=sin(x);plot(x,y,'g:'); 谁能帮我做Matlab的作业:分析系统能否通过状态反馈实现镇定 ...
matlab已知圆心和半径怎么作圆,MATLAB中已知圆心和半径画圆
今天在用MATLAB编程的时候,用到了已知圆心和半径,画圆的程序,上网搜了一下,主要有下面两种,在这里总结一下:(这里我都是放在函数中做的,想画多个圆的话可以加个for循环调用一下函数,或者直接用向量 ...
MATLAB中已知圆心和半径画圆
今天在用MATLAB编程的时候,用到了已知圆心和半径,画圆的程序,上网搜了一下,主要有下面两种,在这里总结一下:(这里我都是放在函数中做的,想画多个圆的话可以加个for循环调用一下函数,或者直接用向量 ...
Java mysql获取行数_java – MySQL查询获取球体中的行(X,Y,Z坐标)？
我正在制作一个名为Minecraft with Bukkit API的游戏插件. 我有一个名为Reinforcements的数据库表,其中包含以下字段:x integer,y integer,z in ...
已知矩阵 matlab,在MATLAB中,已知矩阵A,那么A(:,2:end)表示
摘要: 已知供输工方.表示添资料准加剂.加及标应提原料有关的()的出国法等使用,品""进办理报检时口食.已知信息系统模型不包逻辑括(.... 已知信息系统构化中的结方法设计,矩阵 ...
matlab二维图形中确定一点的坐标系,在CAD中已知一个平面图中两个点的坐标如何建立坐标系以便求出图形中任一点的坐标?...
在CAD中已知一个平面图中两个点的坐标如何建立坐标系以便求出图形中任一点的坐标?以下文字资料是由(历史新知网www.lishixinzhi.com)小编为大家搜集整理后发布的内容,让我们赶快一起来看一 ...
[MATLAB]从已知矩阵中取出子阵
首先要知道冒号表示法,如下: EX: 如上有4×6矩阵,假设我们要取出第一行,第二行,第四行,则可以有B=A([1 2 4],:) 假设要取出第一行,第二行,第四行的1-2列和4-6列,可以有B=A( ...
matlab已知函数表达式画函数图像,怎么用matlab画已知函数表达式的一个函数图像？函数比较复杂的……...
点击查看怎么用matlab画已知函数表达式的一个函数图像?函数比较复杂的--具体信息答:要具体哦,如果一个自变量,函数图像就是曲线,用plot 如果是2个自变量,函数图像就是曲面了,如f=3*exp ...
已知三角形三点坐标, 求三角形的面积
摘自:http://blog.sina.com.cn/s/blog_4cddcaa3010008u7.html 已知三角形三点坐标, 求三角形的面积发表时间:2004-8-10 作者:未知[获得此文 ...

matlab中已知X，Y，Z坐标，绘制等高线图

matlab中已知X，Y，Z坐标，绘制等高线图相关推荐

最新文章

热门文章