作者ExtraCrystal (政)
看板MATLAB
标题[绘图问题] toolbox:m_map 之「m_contour」失效
时间Sat May 7 11:24:42 2016
请问有人知道怎麽用m_contour画等高线吗？
可以帮我看一下我到底哪里出错了吗？
拜托拜托～～～
我已经卡好几天了，不知为何m_contour输出的都是空矩阵
＃原始资料：topo.xyz
(link:
https://www.dropbox.com/s/zlg3qpg3qt1mldj/topo.xyz?dl=0)
＃以下是指令：-------------------------------------
%load data
map_data = dlmread([ inputpath1 'topo.xyz']); %读取资料topo.xyz
map_lon = map_data(:,1); %经度资料
map_lat = map_data(:,2); %纬度资料
map_depth = map_data(:,3); %水深资料
%设定矩阵的界线(设定矩阵的网格）
xmin = 121+27/60+30/3600;
xmax = 121+28/60;
dx = (map_lon(2)-map_lon(1));%奠定LON的距离
ymin = 25+6/60+15/3600;
ymax = 25+6/60+58/3600;
dy = (map_lat(671)-map_lat(672));%奠定LAT的距离
F = scatteredInterpolant(map_lon, map_lat, map_depth, 'linear' ,'none'); %制作
depth与lat&lon的关系式，none 就会将外差的值变成NAN
[X, Y] = meshgrid(xmin:dx:xmax,ymin:dy:ymax);%制作网格点，预备填入资料
V = F(X ,Y);%计算出2D的水深植，这里为负值
%%绘图
figure(2)
m_proj('miller','long', [-121.467 -121.458],'lat',[ 25.104 25.116]);
m_grid('box','fancy','tickdir','in');
hold on
h=m_contour(X, Y, V);
-------------------------------------------
--
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.112.218.85
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1462591495.A.706.html