当前位置：首页 > news >正文

【Delft3D FM数据后处理系列】1. 从Map.nc到精美网格图：Matlab与Surfer的进阶可视化实践

news 2026/7/14 9:06:28

1. 从Map.nc到科研级网格图：为什么需要后处理？

刚接触Delft3D FM的朋友们可能都有这样的困惑：明明模型跑出来了结果，为什么画出来的图总感觉差点意思？我在处理珠江口数值模拟项目时就深有体会——直接用quickplot导出的图片放在论文里，审稿人一眼就看出是"新手作品"。其实这就是原始数据与科研可视化之间的鸿沟。

Map.nc文件虽然包含了完整的网格信息，但直接绘制往往面临三个典型问题：

美学缺陷：默认线条粗细不均、颜色单调，缺乏学术图表应有的专业感
信息冗余：全区域网格展示会掩盖重点研究区域
背景单一：缺少岸线、地形等参照物，降低图件的解释性

这就是为什么我们需要Matlab+Surfer的组合拳。实测下来，这套方案有三大不可替代的优势：

灵活控制：从线条颜色到节点标记都能像素级调整
高效迭代：修改一个参数就能批量更新所有图表
多源融合：轻松叠加卫星底图、测量数据等辅助信息

2. Matlab基础处理：从数据读取到网格绘制

2.1 解剖Map.nc文件结构

先来看看我们的"原材料"Map.nc里藏着哪些宝贝。用ncdisp命令查看文件结构时，这几个关键变量需要特别关注：

ncdisp('Boluo2map.nc'); % 输出节选： % mesh2d_node_x Size: 1256x1 % mesh2d_node_y Size: 1256x1 % mesh2d_face_nodes Size: 4x2145

这里透露了两个重要信息：

节点坐标存储在一维数组里（1256个节点）
网格面由4节点定义（部分三角形网格会用NaN填充）

2.2 智能处理混合网格

实际项目中经常遇到四边形/三角形混合网格的情况。这段代码展示了如何智能分离两种网格类型：

FaceConnect = ncread(mapfilename,'mesh2d_face_nodes'); % 找出四边形网格(第4行非NaN) quad_idx = find(~isnan(FaceConnect(4,:))); % 找出三角形网格(第4行为NaN) tri_idx = find(isnan(FaceConnect(4,:)));

处理技巧在于patch函数的灵活应用——对四边形和三角形分别指定顶点数：

% 绘制四边形网格 patch('Faces',FaceConnect(1:4,quad_idx)','Vertices',[lon lat],... 'EdgeColor',[0.2 0.4 0.8],'LineWidth',0.5); % 绘制三角形网格 patch('Faces',FaceConnect(1:3,tri_idx)','Vertices',[lon lat],... 'EdgeColor',[0.2 0.4 0.8],'LineWidth',0.5);

建议将线条颜色设为RGB值而非简单的'b'，这样后期调整色系时更精准。

3. Surfer美学增强：BLN底图的魔法

3.1 制作专业级BLN文件

Surfer的BLN底图相当于图表的"画布"。以珠江口项目为例，优质BLN文件应该包含：

岸线边界（主轮廓）
重要岛屿（次级闭合环）
航道标记（特殊线型）
控制点标注（文字标签）

用Golden Surfer创建BLN时，注意设置好这些参数：

边界线类型：多边形/折线
高程值：Z列留空则为二维底图
属性块：记录数据来源和比例尺

3.2 Matlab-Surfer联合作战

这个改进版的BLN读取函数增加了智能配色和图层控制：

function plotBLN(blnfile) fid = fopen(blnfile); while(~feof(fid)) header = fgetl(fid); [num, ~] = sscanf(header,'%d'); data = fscanf(fid,'%f,%f',[2 num]); % 自动识别陆地/水域 if mean(data(2,:)) > 24.5 % 珠江口纬度阈值 fill(data(1,:),data(2,:),[0.9 0.9 0.7],'EdgeColor','none'); else fill(data(1,:),data(2,:),[0.7 0.8 0.9],'EdgeColor','none'); end hold on fgetl(fid); % 跳过行尾 end fclose(fid); end

使用时注意图层顺序——先画BLN底图，再叠加网格线：

plotBLN('PearlRiver.bln'); % 先绘制底图 plotGrid('Boluo2map.nc'); % 再叠加网格

4. 科研级出图技巧：从合格到卓越

4.1 焦点区域强化策略

审稿人最反感"大而全"的网格图。这是我总结的焦点区域处理流程：

用inpolygon函数标记兴趣区域节点
对核心区网格加密显示（线宽0.8pt）
外围网格淡化处理（线宽0.2pt，浅灰色）
添加比例尺和方位标

% 标记兴趣区域（示例为矩形区域） x_range = [113.5 114.2]; y_range = [22.1 22.7]; in_area = lon>=x_range(1) & lon<=x_range(2) & ... lat>=y_range(1) & lat<=y_range(2); % 差异化绘制 patch('Faces',FaceConnect(:,in_area),'Vertices',[lon lat],... 'EdgeColor','k','LineWidth',0.8); patch('Faces',FaceConnect(:,~in_area),'Vertices',[lon lat],... 'EdgeColor',[0.7 0.7 0.7],'LineWidth',0.2);

4.2 期刊适配性调整

不同期刊对图表有特定要求。以Elsevier系列期刊为例，需要特别注意：

字体：推荐Arial或Helvetica
线宽：最小0.5pt（印刷要求）
颜色模式：CMYK（非RGB）
DPI：至少600dpi（矢量图更佳）

保存设置可以封装成函数：

function saveJournalFigure(figHandle, filename) set(figHandle,'PaperUnits','inches'); set(figHandle,'PaperPosition',[0 0 7 5]); % 7x5英寸 print(figHandle,'-depsc2','-tiff',['-r600'],filename); end

5. 避坑指南：我踩过的那些雷

第一次投稿被拒的经历让我深刻认识到网格绘制的细节重要性。这里分享几个典型问题：

坐标轴陷阱：

问题：直接使用经纬度导致X/Y轴比例失调
解决：添加axis equal保持1:1比例

内存泄漏：

问题：频繁读取nc文件导致Matlab崩溃
解决：用persistent变量缓存已读数据

function [lon,lat] = getGridData(filename) persistent cacheData; if isempty(cacheData) || ~strcmp(cacheData.filename, filename) cacheData.lon = ncread(filename,'mesh2d_node_x'); cacheData.lat = ncread(filename,'mesh2d_node_y'); cacheData.filename = filename; end lon = cacheData.lon; lat = cacheData.lat; end

Surfer兼容性问题：