- 积分
- 31
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-10-21
- 最后登录
- 1970-1-1
|
发表于 2017-5-25 10:36:32
|
显示全部楼层
%%%%—————Code to draw 3D bathymetry——————————
%-------Created by 猛牛,May 14,2012--------------------
clear;clf;clc;
ncload data.nc
d_elev=elev(721:1600,361:1900); % Luzon Strait
lon=X(361:1900);
lat=Y(721:1600);
[xx,yy]=meshgrid(lon,lat);
%---------------------------------------------
%%
close all
view([10 -10 45])
hold on
fontsize=16
surfl(xx,yy,d_elev);
hold on;
surf(xx,yy,d_elev);
set(gca,'fontsize',fontsize-4)
shading interp
% axis([120 135 14 31 -7000 900])
xlabel('经度','FontSize',fontsize);
ylabel('纬度','FontSize',fontsize);
zlabel('水深地形(m)','FontSize',fontsize);
h=colorbar('v','fontsize',fontsize);
caxis([-7000 3000]);
set(get(h,'ylabel'),'string','水深地形(m)','fontsize',fontsize);
% text(118,22,1000,'Taiwan Island','color','k','fontsize',fontsize,'fontweight','bold');
% text(118,19,0,'Luzon Strait','color','k','fontsize',fontsize,'fontweight','bold');
% text(118,fontsize,100,'Luzon Island','color','k','fontsize',fontsize,'fontweight','bold');
% title('Bathymetry of Luzon Strait','FontSize',31,'fontweight','bold');
load a.mat
a=a/255;
a=flipud(a);
colormap(a)
%alpha(0.5) %设置透明度
|
|