最近写了个画东海地形图的程序

shashagyx

%利用ETOPO1数据绘制东海海底地形图

A=importdata('G:\Kuroshio\data\geodas_ETOPO1-7376\ETOPO1-7376\ETOPO1-7376.xyz');

depth=A(:,3);

dep=reshape(depth,721,601);

dep(find(dep>=0))=NaN;

dep=dep.*(-1);

dep=dep';

load NOAAcoastline.dat;

data=NOAAcoastline;

lon=[120:1/60:132];

lat=[32:-1/60:22];

[x,y]=meshgrid(lon,lat);

m_proj('miller','lon',[120 132],'lat',[22 32]);

m_grid('linestyle','none');

hold on

m_plot(data(:,1),data(:,2),'r','linewidth',1);

[c h]=m_contour(x,y,dep ,[100,200,1000],'k');

clabel(c,h, 'manual','color','b','linewidth',1);

xlabel('经度');

ylabel('纬度');

print(gcf,'-dpng','G:\Kuroshio\data\东海海底地形图.png');

mofangbao

已经转移到matlab版块 建议加图

二爷名声在外

貌似不只是没图，数据也没传…（第一行读取数据）
仅仅这段程序是没有任何意义的啊！

西西

自己原创的？不错不错，值得鼓励
只是楼主的标题打错了，我还很奇怪这是啥地方。

沙颖凯

发到MATLAB版去再传个图。单看程序感觉有可以改进的地方，x, y和dep的大小没法保证一样，这么写可能会有问题。

小小轩轩琪

很赞啊，不错不错，牛人一枚

masterfan

没图呀

baishizhou

建议经纬度也从数据里取，用reshape形成同样维数的数组.

qhmemory

谢谢楼主，程序很不错，建议把图贴出来，不是更加直观~~

757084134

谢谢分享哈