策划面向地理数据不是一个网格

5视图(30天)
我有一个列向量的数据(23038 x 1),对应于一个列向量的长度相同为经度纬度和另一个。我试图绘制这些数据(我用surfm surfacem, pcolorm),但它总是使MATLAB慢下来的崩溃。我假设这是因为matlab喜欢的面向数据的二维矩阵(网格)。这是正确的,如果是这样我怎么把绘图的数据?谢谢!

接受的答案

凸轮扎尔茨贝格
凸轮扎尔茨贝格 2017年11月2日
嘿,科林,
首先,真的看起来整洁的地图当你绘制的点数据是:
使用 griddata ,你可以网当前数据集分成网格:
latmin = min (lat_piomas);
latmax = max (lat_piomas);
lonmin = min (lon_piomas);
lonmax = max (lon_piomas);
gridSize = (100 - 200);
latvec = linspace (latmin、latmax gridSize (1));
lonvec = linspace (lonmax铂金生产商lonmin gridSize (2));
[LatQ LonQ] = ndgrid (latvec lonvec);
ThickQ = griddata (lon_piomas、lat_piomas、厚度、LonQ LatQ);
从那里,这是一个简单的数据以任何喜欢的格式:
axesm (“墨卡托”,“MapLonLimit”,360年[0])
geoshow (LatQ LonQ ThickQ,“DisplayType”,“表面”)
如果你想要真实的,你可以使用拓扑数据显示表面的高度,在使用颜色来显示“厚度”数据。这里有一个例子使用内置的地形数据的例子:
%获得示例地形数据
负载威尼斯平底渔船威尼斯平底渔船topolatlim topolonlim
Rtopo = georefcells (topolatlim topolonlim、大小(威尼斯平底渔船));
ElevQ = geointerp(威尼斯平底渔船,Rtopo、LatQ LonQ,“立方”);
%显示表面
axesm (“墨卡托”,“MapLonLimit”,360年[0])
geoshow (LatQ LonQ ElevQ,“DisplayType”,“表面”,CData的ThickQ)
daspectm (“米”,1000)%夸大地形的高度
希望这是有用的和有趣的。
凸轮

答案(1)

肖恩·德Wolski
肖恩·德Wolski 2017年11月1日
分散的数据列向量是什么,使用:
scatterm
如果这些数据代表一个网格,但面向单一的列向量,然后使用 geointerp 插入到一个二维网格,然后调用surfm或meshm。
4评论
科林塔特尔
科林塔特尔 2017年11月1日
是的,谢谢你!这是整个地理跨度。一个向量的厚度,一个用于纬度,经度。

登录置评。

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!