Gwyddion教程：从AFM原始图到论文级图像处理全流程

第一次用原子力显微镜AFM拍出原始图的时候很多人反应都一样这哪是表面形貌明明是鬼画符。倾斜背景、扫描线伪影、尖峰噪声全叠在一起肉眼根本看不出样品本来的样子。这时候就该轮到 Gwyddion 上场了。Gwyddion 是一款免费开源、跨平台的 AFM/SPM 数据处理软件我这些年处理扫描探针数据基本每天都开着它。它支持大量厂商原始格式能做平面校正、行对齐、滤波、颗粒分析、粗糙度统计、三维显示还能把处理好的图和数据直接导出成论文需要的格式。最难得的是这个软件的学习曲线不算陡跟着完整流程走一遍当天就能上手。这篇教程我尽量按“从打开软件到出图发论文”这个顺序来写不绕弯子重点讲实操、讲参数怎么选、讲我踩过的坑让第一次摸 Gwyddion 的人少走弯路。不管你是刚进实验室的研究生还是被一堆扫描图困住的工程师这篇文章的思路都适用。1. 为什么 AFM 数据要单独用软件处理1.1 AFM 原始图为什么不能直接用AFM 靠探针在样品表面逐行扫描记录探针和样品之间的相互作用最终得到每个像素点的高度信息。听起来很简单但实际操作中原始数据里掺进去的伪影远比想象中多。首先样品表面往往不是刚好垂直于探针的扫描平面。装样的时候差一度扫描出来的原始图就是一张大斜面。这种情况在 AFM 数据里实在太常见了直接拿原始图给合作者看大概率会被问“你这个样品怎么一边高一边低”。这其实不是样品自己的形貌而是测量系统带来的倾斜背景必须靠软件扣掉。其次扫描器本身的压电陶瓷存在非线性、滞回、蠕变等因素会让图像在 X、Y、Z 方向出现漂移。加上环境振动、探针磨损、反馈参数没调好等问题原始数据里还会混进去很多伪影。比如相邻扫描行之间亮度跳变明显的 line artifact又比如图像上孤立的尖峰往往来自外界灰尘或探针瞬间抖动。还有一个特别容易被忽略的问题AFM 数据的 Z 方向单位是纳米甚至埃X、Y 方向扫描范围却是几微米。三个方向量纲差距很大所以原始显示出来的图天然就会“变形”立体感很突兀。Gwyddion 这类软件的重要作用之一就是把数据放到正确的数值坐标系里再逐项做校正。所以AFM 数据处理的本质是把设备采集到的原始信号转成能正确反映样品形貌的图像并给出一组可以写进论文的定量参数。这一步做不好后面的分析全是白搭。1.2 Gwyddion 为什么更合适我上手 SPM 软件时第一反应是优先用设备自带的处理软件。但用过几个不同厂商的官方软件后我发现它们有个共同点功能都全但是比较封闭格式绑定强用法也比较绑定自家设备。再来看 Gwyddion它完全是另一种思路。Gwyddion 是一个开源项目由一个长期活跃的团队持续维护几乎涵盖了 AFM 数据处理的所有常规需求支持大量厂商格式把自家仪器的数据直接拖进去就能识别还能处理 Nanoscope、SPM、NyMD、WSxM、MI/ASCII 等通用格式。Bruker、Park、NT-MDT 等设备导出的数据多数都能直接打开。处理功能完整倾斜校正、行匹配、滤波、阈值掩膜、颗粒分析、粗糙度统计、FFT 滤波、三维显示、数据导出该有的都有。开源免费更新及时不用为授权发愁实验室多台机器都能装批量处理也不另外收费。有批处理功能一批数据可以统一跑同一个流程处理大量扫描图时非常省事。这里我用一个表格直观对比一下 Gwyddion 和设备自带软件在实际使用中的体感差异对比项设备自带软件Gwyddion格式兼容性主要支持自家设备多厂商格式兼容好处理流程自定义一般按向导走灵活步骤可逆可调批处理能力部分有限制多功能完整成本通常随设备授权免费开源社区支持售后渠道为主文档、论坛、论文范例都有如果你只是偶尔处理一两张图用设备自带软件完全没问题。但如果你经常要处理几十张图、需要统一流程比较不同样品的粗糙度Gwyddion 会顺手很多。下面直接进入安装和实操环节。2. 安装与界面速览2.1 安装方法Gwyddion 官方提供 Windows、macOS、Linux 三种常见平台的安装包安装本身没什么难度。Windows 用户直接到官网下载区找 Windows 版安装包下载后一路 Next 就能装完。有一点小提醒新版安装包默认是 64 位如果你的系统还是老版本 32 位记得找对应的 32 位版本。装好之后建议把 Gwyddion 快捷方式固定在桌面或任务栏后面拖文件方便。Linux 下安装更简单不少发行版已经把它收进官方软件源。比如 Debian/Ubuntu 系列直接用 aptsudo apt install gwyddionFedora 用 dnf install gwyddionArch 系直接用 pacman -S gwyddion。用软件包管理器安装的好处是依赖自动解决升级也方便。macOS 用户到官网下载 dmg 文件即可。不过因为系统权限策略新版本 macOS 第一次打开可能需要右键选择“打开”这是正常现象不用慌。如果你喜欢尝鲜也可以去官网下载源码自己编译。但说实话日常使用真不建议自己从源码编译。Gwyddion 依赖 GTK、FFTW、libpng、libtiff 等一系列库环境配置比较繁琐遇到报错非常劝退。普通用户用官方编译好的稳定版就够了。2.2 界面一句话导航第一次打开 Gwyddion界面看上去有点“程序员风”但熟悉之后会发现布局非常清晰。顶部是菜单栏分类很直白文件、数据、工具、图像、图表、三维视图等。左侧是数据列表打开的一张或多张图都会列在这里。中间大区域是图像显示区默认是 2D 伪彩图亮色代表高、暗色代表低。右侧或底部还有数据信息栏能看扫描范围、像素数、Z 范围、当前鼠标位置坐标等参数。重点记住这几个默认面板数据列表Data List显示当前打开的所有图像双击可以切换当前图。工具箱Toolbox里面有选择、平移、放大、剖面线、距离测量等工具。很多人一打开找不到剖面线工具其实就在工具箱里图标是一个带斜线的尺子。数据信息Data Info显示所选数据的基本参数包括扫描尺寸、像素数、高度单位等。数据处理菜单绝大多数处理功能都集中在“数据”菜单下。倾斜校正在 Data → Level surface行对齐在 Data → Align rows尖峰去除在 Data → Remove Spikes。第一次接触不要急着把每个按钮都点一遍先记住两个关键入口就够用一个是数据列表窗口另一个是数据菜单下的处理功能。后面的流程全部围绕它们展开。3. 第一次完整跑通“清洗流程”3.1 打开数据文件与初始检查Gwyddion 打开文件有两种方式。一是点击菜单 File → Open二是把文件直接拖进主窗口。大多数 AFM 厂商的原始文件它都能识别后缀不一定标准可能是.spm、.001也可能是一堆附带文件里的某个主文件。拖进去之后软件会自动判断格式。如果打开之后发现图像内容不对先确认“数据列表”里选中的是不是高度通道Height 或 Topography而不是幅度图Amplitude或相位图Phase。很多仪器一次扫描会同时输出好几个通道Gwyddion 打开时默认显示的通道不一定是你想分析的那个。切换通道之后整个图像可能完全变成另一个样子。一个非常值得养成的操作习惯拿到原始数据之后先复制一份再处理。Gwyddion 里可以用 Data → Duplicate data也可以在数据列表右键选择 Duplicate。这样万一处理坏了原始图还在。而且 Gwyddion 本身采用非破坏式处理每个步骤都记录在数据历史里理论上可以逐步撤销。这个撤销功能非常强大但也容易让人产生“随便乱调也能恢复”的错觉后文我会专门提醒为什么还是要谨慎。打开数据后的第一步我建议先做肉眼检查。用鼠标滚轮放大图像仔细看有没有大面积坏区、明显扫不到的死点、或整条整条的亮暗线。如果坏区太大很多情况下可以直接裁掉避免影响后续拟合。裁切用工具箱里的矩形选择工具框选区域然后 Data → Crop 即可。3.2 剔除坏行和尖峰扫描线伪影几乎每个人都会遇到表现就是图上出现一条水平方向异常亮或暗的线宽度可能只有一行也可能连续好几行。这种情况通常由扫描过程中探针状态突变、样品台阶侧壁回扫或者外界振动造成。Gwyddion 处理这个问题的工具是 Data → Fix Lines。它能自动检测坏行并用周围行插值替代。操作上先用矩形框选包含坏行的区域尽量把坏行完整包围进去。打开 Fix Lines 对话框选择修复坏行的选项。预览没问题之后点执行。如果坏行太强一次修不干净可以重复执行或者考虑把图像切块分区域处理。但这里有个严重警告算法默认把突变的行当作坏行处理如果样品本身的形貌恰好是周期性条状结构比如光栅这时候直接套用 Fix Lines 会把真实形貌也抹掉。我的做法是先用剖面线工具拉一下确认是突然出现的异常跳变再修。尖峰是图像上孤立的亮点本质是 Z 方向瞬时跳变。轻微的尖峰用中值滤波就能压下去Gwyddion 有专门命令Data → Remove Spikes。执行时可以设定阈值比如超过该列中值多少倍算作尖峰。如果你在表征纳米颗粒建议把阈值设得保守一些因为颗粒本身也像“尖峰”阈值太低了容易把真实颗粒也当作噪声去掉。3.3 平面校正最关键的一步平面校正也就是常说的 Flatten 或 Level是整个 AFM 数据处理里最重要的一步。AFM 扫描出来的原始图Z 值通常带有一个大斜面甚至是一个二次曲面。这个斜面主要来自样品安装角度、扫描器运动方向差异不是样品真实形貌。Gwyddion 里做平面校正最常用的是两个命令Data → Level surface。Data → Remove background。区别在于Level surface 是用最小二乘法拟合一个平面并减掉它适合纠正整体倾斜Remove background 更强大可以扣除多项式、高斯甚至自定义背景函数。实操里最常见的流程是复制数据。先执行 Data → Level surface选“Subtract plane”最小二乘平面把主要大斜面试掉。如果表面还有明显的二次曲面比如边缘翘起、整体呈盆形再用 Data → Remove background选择二次多项式背景扣除。这里给一个特别重要的提醒做平面校正前一定先把图像里你真正关心的表面结构用 Mask 遮罩起来。比如你要分析基底上的纳米颗粒如果直接把带颗粒的图像用最小二乘平面拟合颗粒的存在会把拟合平面往上抬导致扣除后基底整体偏低颗粒高度被低估。正确做法是先用 Mask 把颗粒区域遮住让拟合计算忽略这些区域处理完再取消遮罩。Gwyddion 的 Mask 工具可以手动画也可以按高度阈值自动生成。菜单里的 Mask 子菜单包含 Add、Remove、Threshold 等多种操作。颗粒分析前后都要记得检查 Mask 状态否则统计结果会莫名其妙。平面校正做得好不好的判据不是“看起来平”而是你关心的区域形貌是否被保留。处理完薄膜样品时我经常会用剖面线拉一下处理前后的基底高度确认台阶和颗粒的高度差没有发生系统性偏移。3.4 行对齐处理扫描线伪影如果图像已经做了平面校正但扫描行之间还是存在明显的亮度不均下一步就是行对齐。行对齐的目的是找到各扫描行的 Z 基线让相邻行之间的高度尽量一致从而消除行方向的多项式伪影。操作路径是 Data → Align rows。这个工具很强大它允许你选择每一行要拟合多少阶的多项式通常是 1 阶或 2 阶然后对行内非掩膜区域做拟合再做行间匹配。常见参数设置参考参数推荐值说明拟合多项式阶次1 或 2阶次越高越容易把真实形貌也拟合掉是否使用 Mask是粒子和缺陷区域没遮罩会影响行基线校正范围整图或局部局部坏行可切割处理后合并我第一次用行对齐处理楼梯状台阶样品时误以为它一定会让图像变得更干净结果发现它把台阶高度都给“拉平”了。后来才明白行对齐本质上是假设样品大部分区域是平坦的然后用多数行的高度信息去校准少数偏离的行。如果样品本身就有大台阶或大颗粒必须先用 Mask 把那些区域保护起来否则行对齐会帮倒忙。3.5 滤波什么时候用什么时候别用很多新手一看到图有噪声第一反应就是点“高斯滤波”。但滤波是一把双刃剑它会改变数据本身导致后续粗糙度、粒径统计全部偏移。所以我的原则是能靠扫描参数解决的伪影尽量在采集端解决必须在软件里滤波的话滤波放在所有几何校正之后。Gwyddion 的滤波功能在“数据 → 滤波”菜单下常见的有高斯滤波适合去除高频噪声但会让边缘变模糊。中值滤波对脉冲噪声更有效同时保留边缘的效果比高斯好。均值滤波窗口内噪声会变小但边缘保留能力一般。FFT 滤波支持低通、高通、带通还能直接生成频谱图对周期性伪影或扫描仪噪声很有效。经验上探针尖触碰导致的“振铃”伪影用简单的空间域滤波往往不够。这时候可以用 FFT 滤波里的低通选项先看频谱图找出周期性噪声的峰值再设置截止频率。操作路径是 Data → Filters → FFT选择低通或带通类型调整频率范围并预览效果。但我不建议一上来就上 FFT。FFT 滤波的边界效应非常明显处理不好会让图像边缘出现很大的失真。常规做法是先用中值滤波压掉尖峰再用一个很小的 3×3 或 5×5 高斯滤波做轻平滑就能解决绝大多数视觉问题。如果要做定量分析最好在滤波前后各导出一组数据对比一下关键参数的变化幅度这样你才知道滤波到底改变了什么。4. 定量分析粗糙度、粒径与剖面4.1 粗糙度统计AFM 最重要的数据输出之一就是粗糙度参数。最常见的有Sa/Ra算术平均粗糙度指整个区域上表面高度偏差绝对值的平均。Ra 通常用于一维轮廓线Sa 用于二维区域不过很多文章里混着用。Sq/Rq均方根粗糙度对偏离平均面的大小更敏感。Sz/Rz最大高度差即最高峰和最低谷之间的距离。Ssk 偏度、Sku 峰度描述表面高度分布的对称性和尖锐程度。Gwyddion 里查看这些参数非常方便。选中图像数据点击菜单 Data → Statistics或者直接在左侧信息面板看统计值。统计窗口里有 Height、Roughness 等多个标签页还能看到高度直方图导出 CSV 也没问题。如果你只想统计图像中某一小块区域先用工具箱里的矩形选择工具框选再打开统计对话框Gwyddion 默认只统计当前选区否则统计的是整张图。这里有个容易踩的大坑统计粗糙度之前一定要先做扎实的平面校正和行对齐。粗糙度本质上就是高度分布的标准差如果图像带一个大的倾斜背景Ra 会被这个假坡度严重拉高。所以我的报告里都会写清楚图像经过 Level surface Align rows 之后Ra 是多少。这样别人复现时才知道数据流程是什么。另外一个常被忽略的问题是扫描分辨率。同样一片表面256×256 像素图和 1024×1024 像素图算出来的 Ra 很可能不同因为采样点越密细小特征越容易被抓到。如果你想比较不同样品的粗糙度尽量用相同的扫描范围、相同的像素数、相同的扫描速度否则对比结果说服力不足。4.2 颗粒分析与粒径统计做纳米颗粒、薄膜岛状生长或生物样品时颗粒分析是重头戏。Gwyddion 的颗粒分析模块在 Data → Grain analysis功能出乎意料地强。它的基本逻辑是通过高度阈值把图像二值化把高于阈值的区域认作颗粒然后逐颗粒统计面积、等效直径、高度、体积等参数。操作流程大致是这样先完成平面校正和必要的滤波。如果图像里有大块不需要分析的区域比如杂质、划痕先用 Mask 工具遮罩掉。打开 Data → Grain analysis。在 Threshold 选项卡里设置阈值。可以是绝对高度比如颗粒高度超过基底 5 nm也可以是相对值比如超过图像所有像素高度的某个百分位。我比较常用绝对阈值直观且好解释。调整分割参数。Gwyddion 支持分水岭分割Watershed如果颗粒挨得很近单纯阈值会把多个颗粒连成一个此时建议启用分水岭。在 Statistics 里勾选你需要的参数颗粒数、平均面积、等效直径、平均高度、最大高度、体积等。点击执行软件自动生成结果表。我第一次用颗粒分析时发现它输出的颗粒数量多到离谱。原因很常见没有屏蔽背景噪声阈值设太低软件把表面台阶和噪声起伏全当成了颗粒。正确做法是先用 Mask 把明显不是目标颗粒的特征挡掉再根据目标颗粒的高度范围调整阈值。判断阈值合不合适可以靠一个很好用的辅助功能让软件在原始图上描出识别出的颗粒轮廓。看着轮廓叠加图去调阈值比只看数字直观得多。另外如果颗粒密度很高、互相粘连分水岭分割几乎必开。它会把连在一起的峰按“山谷线”切开但分水岭参数选大了会把单个大颗粒切碎选小了又分不开。我的习惯是对一个小区域做参数扫描分别用几个相邻参数跑一遍对比颗粒数量分布曲线选一个不发生剧烈变化的中间值。4.3 剖面线分析剖面图是验证形貌最直观的手段。工具箱里的剖面线工具图标像个带箭头的斜线。点选之后在图像上拖动画一条线底部会弹出一个剖面图窗口显示线经过位置的高度变化。利用剖面线可以测量台阶高度。画一条垂直于台阶的线剖面图会清晰展示两级平台之间的高度差。测量颗粒高度和宽度。在颗粒上拉一条线看剖面峰的半高宽。检查平面校正效果。如果校正正确基底区域拉出来的线应该是平的而不是倾斜的。Gwyddion 的剖面图可以直接导出为图片或数据。在剖面图窗口点右键选择 Export或者用 CtrlS 保存图形。如果想把剖面数据导到 Origin 或 Python 里进一步处理可以在剖面窗口用 Copy Data 复制数值然后粘贴到文本编辑器或表格软件。剖面分析有一点务必提醒画线时尽可能沿样品形貌变化方向画不要斜穿一个颗粒就完事。颗粒直径和高度在不同方向可能有差异单条线会带来很大偏差。最少每条分析数据画三到五条平行剖面线取平均值后再写进结论。5. 三维显示与论文图导出5.1 三维图怎么调才好看Gwyddion 的三维视图是我非常喜欢的功能它可以把平面伪彩图变成一个可旋转的 3D 表面。操作很简单选中图像点击菜单 三维视图 → 新建三维视图3D View → New 3D View。默认生成的 3D 图往往不太好看主要问题是光照不对、Z 比例尺太大图形看起来像一个带刺的山脊。调 3D 图时主要关注这几个面板视角鼠标拖动旋转视角滚轮缩放。找个能看到表面起伏和边缘结构的视角不要正上方向下看那样和 2D 图区别不大。光源默认光源可能产生过大的阴影把光源调成从顶部斜照下来表面细节会更真实。Z 缩放Gwyddion 默认的 Z 轴比例容易把形貌放大太多。论文插图里比较常见的做法是不按真实比例显示但比例尺必须写清楚说明“Z 轴已缩放”避免读者误解真实起伏量级。颜色表选一个对比度合适的颜色映射比如冷色到暖色渐变看着更有层次感。调好角度和光照后在三维视图窗口里直接右键 → Save graphics或者 File → Export 导出图片。建议导出时选白色背景因为很多期刊系统不接受透明背景或渐变底色。为了清楚表示高度我通常会在 3D 视图里开启颜色比例条。Gwyddion 可以在三维视图设置里勾选“显示颜色比例条”。比例条这种元素后期在 AI 里自己加很难对整齐所以生成图片时带好最省事。5.2 论文级 2D 图导出2D 图导出也有讲究。直接在 File → Export 里选 PNG 格式可以得到位图适合快速预览。但投稿时位图至少需要 300 dpi否则印刷出来会模糊。Gwyddion 的导出对话框里有分辨率设置不要用默认的 100 dpi 直接导出。我一般设置成 600 dpi或者直接导出 SVG/PDF 矢量图。矢量图非常推荐怎么缩放都不会糊编辑起来也方便。如果期刊要求位图再从矢量图软件里转成 TIFF 也不迟。导出之前记得关掉不必要的网格线和标注。Gwyddion 默认显示网格线这在分析时有用但作为最终投稿图通常不需要。可以在“查看”菜单里关掉“显示网格”再导出图片。如果你的后续排版需要把多张图拼在一张 figure 里记得导出时保持相同的视角、相同的缩放比例、相同的颜色表。经验和教训告诉我子图之间比例不统一是编辑部最喜欢挑的毛病之一。5.3 数据矩阵与统计结果导出写文章时不仅需要图片往往还得把原始数据数值放上去。Gwyddion 可以把整张图像以 ASCII 矩阵格式导出也就是一行对应图像的一行像素用空格或制表符分隔。这个数据可以直接导入 Origin、MATLAB 或 Python 做进一步绘图。方式是 File → Export data选择 Data in spreadsheet/text 格式。统计表也可以导出。颗粒分析结果、粗糙度统计结果都能在对应窗口里找到 Export 或 Save As保存成 CSV 后用 Excel 直接打开。导出 ASCII 矩阵时有个小坑图像像素数太大时导出的文本文件可能非常大能到几百万行Excel 打开会直接卡死。如果只需要单一剖面数据建议用剖面图窗口复制如果需要整张图的面数据可以先把图像降采样到 512×512 或 256×256 再导出数据量小很多数值趋势不会变。6. 常见问题与实战经验6.1 高频问题速查表这里把我在实际使用中被问过很多次的问题整理成一张速查表希望帮你少走弯路。现象可能原因解决办法打开文件后图像全黑或全白选中的不是高度通道检查数据列表切换到 Height 通道图像像“屋面”高低不平倾斜背景未扣除Data → Level surface 先扣平面图像有明显水平条纹行伪影Data → Align rows 或 Fix lines颗粒分析数量异常多阈值太低或没有遮罩提高阈值Mask 掉杂质必要时开启分水岭3D 图看起来像刺猬Z 轴比例过大降低 Z 缩放系数导出图片模糊分辨率不足导出时设置 300 或 600 dpi数据矩阵在 Excel 里卡死文件太大降采样后导出或用 Python/Origin 打开还有一个很容易被忽略的问题重复处理。Gwyddion 的撤销功能强大但如果你保存项目并重新打开有些过程会合并。我在处理大量数据时流程顺序很固定原始数据 → 复制 → 裁剪 → 尖峰去除 → Level surface → Align rows → 滤波 → 分析。每一步确认无误后再进入下一步避免中途反复跳步骤导致混乱。6.2 几个实用窍门这里分享几个不太起眼但非常实用的操作技巧。技巧一善用平坦区域。做平面校正和行对齐时如果图像里保留一块“什么都没有”的平坦基底拟合精度会大幅提高。我在扫描时会特意选择包含干净基底区域的视野这样软件拟合出来的平面更可靠。反过来如果你的视野全是密集颗粒平面校正本身就很难做准尽量重新选区域扫描。技巧二颗粒分析阈值不要拍脑袋定。先用剖面线工具大致量出颗粒高度范围比如 10-15 nm然后从 8 nm 阈值开始跑逐渐上调到 12 nm观察粒径分布曲线什么时候最稳定。把这个范围记录下来写入方法部分。实验记录里保留这一步能让你在回复审稿意见时更有底气。技巧三批量处理用起来。Gwyddion 有专门的 Batch Processing 模块可以记录当前图像的处理流程并在其他图像上批量重放。如果你手上有几十张不同位置的扫描图都用同一套处理流程强烈建议用批量工具。路径大概是 Tools → Batch → Create/Edit界面不复杂。我不建议新手一上来就学 Python 脚本批量工具已经能解决大部分重复劳动。6.3 一套可复制的标准工作流最后总结一套我自己一直在用的标准流程你可以直接抄作业复制原始数据保留一份备份。裁剪掉明显坏区。用 Remove spikes 去除尖峰。Data → Level surface 扣掉平面倾斜。若存在行伪影用 Align rows 做行对齐。对需要分析的表面区域加 Mask再做一次 Remove background把基底修平。如果只是展示用再做一个轻微的 3×3 中值和 3×3 高斯平滑。统计粗糙度或进行颗粒分析。导出 2D 图、3D 图和数据矩阵。这套流程可以处理绝大多常规 AFM 图像。如果跑完流程图还是脏多半说明是采集端的问题还没解决这时候应该回头调扫描参数或重新制样而不是盲目靠软件硬救。用久了你会发现处理软件再厉害也救不了一块烂样品。Gwyddion 的意义不是把垃圾数据变成黄金而是把好数据里的系统性伪影剥掉让真实形貌自己站出来。希望你下次拿到 AFM 原始图时能对着软件淡定地一步一步处理而不是看得头大。