figure(1) set(gcf,'position',[100 100 1200 650] ,'color','w') for i=2013:2022 % suptitle([num2str(i) '年月平均大气可降水量 (mm)']) for j=1:12 temp=tp(:,:,years==i & months==j); monmean(i-2012,j)=mean(temp(:)); subplot(3,4,j) pcolor(LON,LAT,temp'); shading interp colormap(jet); % caxis([0 15]) colorbar ax = gca; title([ num2str(j) ' 月']) ax.XLabel.String = '经度(°)'; ax.XLabel.FontSize = 10; ax.YLabel.String = '纬度(°)'; ax.YLabel.FontSize = 10; ax.FontSize = 10; pause(0.0002) end suptitle([num2str(i) '年月平均大气可降水量 (mm)']) print(gcf,'-djpeg','-r300',['.\Pic\1_' num2str(i) '年月平均大气可降水量.jpg']); clf end
时间: 2024-04-25 18:26:34 浏览: 135
这是一段 MATLAB 代码,用于绘制 NetCDF 格式的数据文件中的大气可降水量数据,并将每年的月平均可降水量作为一幅图像保存到本地。
具体来说,`figure(1)` 表示创建一个名为 1 的新图形窗口,`set(gcf,'position',[100 100 1200 650] ,'color','w')` 表示设置图形窗口的位置和背景色。
`for i=2013:2022` 表示循环遍历从 2013 年到 2022 年的每一年。内部循环 `for j=1:12` 表示遍历每一年的每一个月。
`temp=tp(:,:,years==i & months==j)` 表示从原始数据文件中获取当前年份和月份对应的大气可降水量数据,`monmean(i-2012,j)=mean(temp(:))` 表示计算当前月份的平均可降水量,并将其保存到 `monmean` 数组中。
`subplot(3,4,j)` 表示在当前图形窗口中创建一个子图,`pcolor(LON,LAT,temp')` 表示在子图中绘制经纬度网格和对应的可降水量数据,`shading interp` 表示使用插值方法进行颜色填充,`colormap(jet)` 表示使用 Jet 颜色映射表,`colorbar` 表示添加颜色条,`ax = gca` 表示获取当前子图的坐标轴对象,`title([ num2str(j) ' 月'])` 表示设置子图标题,`ax.XLabel.String = '经度(°)'` 和 `ax.YLabel.String = '纬度(°)'` 分别表示设置坐标轴标签,`pause(0.0002)` 表示暂停 0.0002 秒,以便在绘制子图时能够逐步显示。
`suptitle([num2str(i) '年月平均大气可降水量 (mm)'])` 表示设置当前年份的总标题,`print(gcf,'-djpeg','-r300',['.\Pic\1_' num2str(i) '年月平均大气可降水量.jpg'])` 表示将当前图形窗口保存为 JPEG 文件,`clf` 表示清空当前图形窗口,以便绘制下一个年份的数据。
阅读全文