画图出现问题

zhu

这是绘制程序:
x=0:2.5:357.5;  %经度范围
y=-90:2.5:90;  %维度范围
x=x';
for k=1:73
    kk=(k-1)*144+1;
    R3(:,k)=R2(kk:(kk+143));    %把数据排列成矩阵形式，144行(经度），73列（纬度）
    P3(:,k)=P2(kk:(kk+143));    %把数据排列成矩阵形式，144行(经度），73列（纬度）
end
R4=R3';   %矩阵转置:73行，144列
P4=P3';   %矩阵转置:73行，144列
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(1)
v1=-0.1:-0.1:min(min(R4));
v1=fliplr(v1);
v2=0:0.1:max(max(R4));
[c,h]=contour(x,y,R4,v2); %绘制等高线
set(h,'ShowText','on','Linestyle','-','Color','b')
hold on
[c,h]=contour(x,y,R4,v1); %绘制等高线
set(h,'ShowText','on','Linestyle',':','Color','r')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
hold on
[c,h]=contourf(x,y,-P4,[-0.05 0]); %绘制等高线
title('June(same year)')
%%%%%%%%%%%%%%%%%%%%%%%
%mapping world coast
hold on
load coast
a=length(lat);
b=length(long);
for i=600:b/2.25
    if long(i)<0   
        long(i)=long(i)+360;
    end
end
geoshow(lat,long);
grid on
xx=0:30:330;
yy=-90:30:90; 3.tif (2.55 MB, 下载次数: 2)

2014-7-2 09:27 上传

点击文件名下载附件

set(gca,'Xlim',[0, 357.5]);
set(gca,'Ylim',[-90, 90]);
set (gca,'XTick',xx,'XTickLabel',{'0','30E','60E','90E','120E','150E','180E',...
    '150W','120W','90W','60W','30W'});
set (gca,'YTick',yy,'YTickLabel',{'90S','60S','30S','EQ','30N','60N','90N'});  

晓风残月

好像不止南北美洲轮廓，从东京180往东的海岸轮廓都没有，感觉像是循环哪里出问题了，话说绘制这种等值线图为什么用matlab，用grads很快很方便

无径之林

估计你coast里long范围是【-180 180】
b=length(long);是long个数
为什么下边是for i=600:b/2.25，这样你循环肯定考虑到long截尾
通过读你的程序，我的理解是：
long是【-180 180】，但是前600之后才是西经，所以你才会从600开始
你想把西经变为正数，没那么麻烦，不用循环
直接index=find(long<0) ;  long(index)=long(index)+360;

无径之林

估计你coast里long范围是【-180 180】
b=length(long);是long个数
为什么下边是for i=600:b/2.25，这样你循环肯定考虑到long截尾
通过读你的程序，我的理解是：
long是【-180 180】，但是前600之后才是西经，所以你才会从600开始
你想把西经变为正数，没那么麻烦，不用循环
直接index=find(long<0) ;  long(index)=long(index)+360;

striver_zhu

晓风残月 发表于 2014-7-2 10:14
好像不止南北美洲轮廓，从东京180往东的海岸轮廓都没有，感觉像是循环哪里出问题了，话说绘制这种等值线图 ...

windows下可以在C++中调用grads画图吗

晓风残月

striver_zhu 发表于 2014-7-2 14:34
windows下可以在C++中调用grads画图吗

这个没有听说过，没必要混编啊，直接编写grads就好。或者matlab程序再改改也行，我觉得问题就在coast设置的long的循环里。

striver_zhu

晓风残月 发表于 2014-7-3 09:55
这个没有听说过，没必要混编啊，直接编写grads就好。或者matlab程序再改改也行，我觉得问题就在coast设置 ...

我在windows下做编程开发，想把这个做第三方插件使用。。。

沙颖凯

你把地形数据由-180 180转换成0 360的过程出错了，海岸线数据因为不是规则的网格所以事实上没有完美的转换方案。单纯把小于0的值加上360是用不了的。楼主可以选择把要画图的数据x, y转换成-180 180的，或者启用专门的地理绘图环境来画。

二爷的屁股

帮定{:5_213:}{:5_213:}