matlab绘制水深数据

?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

clc;clear; %% %------------------------------1.查看文件信息------------------------------- %file为长波短波辐射6月文件 file = 'GEBCO_2021.nc'; ncdisp(file);                   %在界面查看nc文件内容，在命令行窗口阅读数据内容   %% %------------------------------2.读取数据信息------------------------------- lon = ncread(file,'lon');       %读取经度数据-180 to 180 lat = ncread(file,'lat');       %读取纬度数据 -90 to 90 crs = ncread(file,'crs');       %char，latitude_longitude elevation= ncread(file,'elevation');       %86400x43200，height_above_mean_sea_level clearvars file ;     %% %-----------------------------3.数据处理------------------------------------ ii = find((lon>117.0021)&(lon<123.0021));%渤海经度范围 jj = find((lat>37.0021)&(lat<41.0021));  %渤海维度范围 kk = find(elevation(ii,jj)>0);%查找水深大于0的数据   [LON,LAT] = meshgrid(lon(ii),lat(jj)); %将经纬度建成矩阵数据 elevation1 = elevation(ii,jj); elevation1(kk)= NaN;   clearvars lon lat ii jj kk crs elevation;%清除变量 whos;             %按字母顺序列出当前活动工作区中的所有变量的名称、大小和类型。 %--------------------------------------------------------------------------   %% figure;clf; hold on; %----选择投影方式:Robinson,Mercator,UTM，,Lambert,Stereographic,ortho----- m_proj('Robinson','long',[117,123],'lat',[37,41]);%初始化地图投影信息，将结果放入一个结构中  m_coast('patch',[.6 .6 .6],'edgecolor','k');     %在给定的地图上标出海岸线 m_grid('linestyle','none','tickdir','out','linewidth',1); m_plot(LON,LAT,'k','linewidth',1.5); m_pcolor(LON,LAT,elevation1') ;                        %画水深矩阵彩色图 m_contourf(LON,LAT,elevation1');                       %加等值线 colormap(jet(6))                %图像颜色，jet/autumn，可多颜色了自己选 c = colorbar;                   %插入色条 c.Label.FontSize = 12; c.Label.String = '水深（m)';  %插入色条名称及单位 xlabel('longitude');            %xy轴 ylabel('latitude'); title('渤海水深分布图','fontsize',15) %图像标题 % title('Water depth of Bohai (meter)');