matlab 展示SCS三维海底地形图

Hua_Zhang

matlab 展示SCS三维海底地形图

只为秀图而已咯，哈哈

1. clc ; close ; clear;
   
2. path = 'xx';
   
3. lon0=ncread([path,'ETOPO1_Bed_g_gmt4.grd'],'x');
   
4. lat0=ncread([path,'ETOPO1_Bed_g_gmt4.grd'],'y');
   
5. topo0=ncread([path,'ETOPO1_Bed_g_gmt4.grd'],'z');
   
6. topo0=topo0';
   
7. 
   
8. lonmin = 98;
   
9. lonmax = 130;
   
10. latmin = 0;
    
11. latmax = 30;
    
12. ix = find(lon0>lonmin & lon0<lonmax);
    
13. iy = find(lat0>latmin & lat0<latmax);
    
14. 
    
15. tlon = lon0(ix);
    
16. tlat = lat0(iy);
    
17. topo = topo0(iy,ix);
    
18. data = topo; data(data>=0)=nan;
    
19. [xx,yy] = meshgrid(tlon,tlat);
    
20. max = max(data(:));min = min(data(:));
    
21. %% ----------------plot topo by joker------------
    
22. close ;
    
23. figure
    
24. view([12 -10 50])
    
25. %surf(x',y',topog,[1 1],'-k','linewidth',2);
    
26. hold on
    
27. surf(xx,yy,data);
    
28. shading interp;
    
29. axis([lonmin lonmax latmin latmax min max])%lon:30~55(x);lat:-40~0(y);
    
30. xlabel('longitude(E)','FontSize',16);
    
31. ylabel('latitude(S)','FontSize',16);
    
32. zlabel('Depth(m)','FontSize',16);
    
33. %colormap(b2r(min,max));
    
34. mycolor=zeros(3,202);
    
35. mycolor(1,:)=[0:0.01:1 zeros(1,101)+1];
    
36. mycolor(2,:)=[0:0.01:1 1:-0.01:0];
    
37. mycolor(3,:)=[zeros(1,101)+1 1:-0.01:0];
    
38. mycolor = mycolor';
    
39. colormap(mycolor);
    
40. h = colorbar('v');
    
41. caxis([min max]);
    
42. 
    
43. set(get(h,'ylabel'),'string','Depth(m)','fontsize',21);
    
44. text(lonmin+5,latmax-5,10,'China','color','k','fontsize',25);
    
45. text(lonmax-5,latmax-15,10,'Pacific','color','k','fontsize',20);
    
46. text(lonmin+15,latmin+10,10,'SCS','color','k','fontsize',20);
    
47. 
    
48. set(gca,'fontsize',16);
    
49. disp(' !!OK!! ');
    

复制代码

xiaofeixiazyh

赞！学习一下

贫道敬孔

漂亮！！