|
发表于 2023-12-7 23:15:25
|
显示全部楼层
首先,你需要导入mat文件并获取GridInfo和EcomSalt的数据。以下是一个示例代码:
- load('Ex13_PRE_Salinity.mat');
- % 获取盐度数据
- salinity_data = EcomSalt.data;
- % 获取日期值
- time = EcomSalt.time;
- % 获取每天12点的数据索引
- index = find(mod(time,1) == 0.5);
- % 创建存放图片的文件夹
- mkdir('Ex13FigDir');
- % 循环生成每天12点的盐度分布图
- for i = 1:length(index)
- % 获取当天的盐度数据
- salinity = salinity_data(:, index(i));
-
- % 画出盐度分布图
- figure;
- contourf(GridInfo.lon, GridInfo.lat, salinity);
- colorbar;
-
- % 设置标题和标签
- title(['Salt Distribution at 12:00 on Day ' num2str(i)]);
- xlabel('Longitude');
- ylabel('Latitude');
-
- % 保存图片
- filename = sprintf('Ex13_DrawSalt%03d.png', i);
- saveas(gcf, fullfile('Ex13FigDir', filename));
-
- % 关闭当前图像窗口
- close;
- end
复制代码
这段代码首先加载mat文件并获取盐度数据和日期值。然后找到每天12点的数据索引,并循环生成每天12点的盐度分布图。最后将生成的图片保存到Ex13FigDir文件夹中,图片命名为Ex13_DrawSalt001.png、Ex13_DrawSalt002.png等等。
以上回复来自 FishC AI 小助理,如未能正确解答您的问题,请继续追问。 |
|