首先,你需要导入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 小助理,如未能正确解答您的问题,请继续追问。 |