去噪实验研究利器——教你画局部相似图（附程序）
研究地震去噪的同学们都知道，检验地震去噪效果最直观的方法就是画图来看。
有一种专业评估去噪性能的方法，就是使用局部相似图。
01 什么是局部相似图
局部相似图（Local similarity map），是通过显示去噪结果与去除噪声之间的局部相似性来评估信号泄漏的程度，那么一些高相似性异常的区域就表明去噪后的结果在相应的位置有信号泄漏。
就比如上图，c图相比a和b两图色块较少，其局部相似性的低值表示良好的正交化，具有更小的信号泄漏，说明该方法能够充分地保留信号并有效地衰减噪声，因此具有良好的去噪性能。
大家很关心到底怎样画出这样的图呢？我们一般都用Madagascar来画的。Madagascar即马达加斯加，是一个用于多维数据分析和可重复计算实验的开源软件包。它的任务是为地球物理学和相关领域的数字图像和数据处理研究人员提供一个方便、强大的环境和一个方便的技术工具。
Madagascar具体介绍大家可以见网站：
https://reproducibility.org/wiki/Main_Page
02 实战案例教你画局部相似图
下面，带领大家从0开始来实现局部相似图的绘制。
1.安装Madagascar运行环境。
作为开源的软件库，Madagascar的安装教程在网上很多。
Madagascar需要在linux操作系统安装，要求大家要有一些的动手能力。安装步骤主要包括：
一是下载 Madagascar 软件包。
二是解压文件。
三是进入解压文件目录并且进行安装环境设置。
四是编译并且安装软件。
五是用户环境变量设置。
六是测试Madagascar能否正常运行。
此处省略具体安装细节5000字，大家可以参考附件《Madagascar3.0稳定版安装指南》一步步安装。
建议大家安装xwindow，这样最后出图才能正常显示。
2.绘制局部相似图
首先，我们准备了segy地震数据。一个是含噪数据noisy.segy，一个是去噪后的数据denoise.segy，第三个是去掉的噪声noise.segy。显示效果是这样的：
然后，将noisy.segy和noise.segy数据上传到运行环境，对segy数据制作rsf配置文件。
在剖面数据对应目录运行程序：
这样就生成了rsf文件。
第三，编辑SConstruct文件，编写程序。
以上关键程序主要功能就是显示了局部相似图。
第四，显示局部相似图的效果。
先运行scons进行程序编译。
这里编译结果，res越大表示信号泄露越小，效果就比较好。
再运行scons view显示局部相似图。
效果就是这样的：
从这个效果来看，信号泄露非常小。
如果想在Madagascar3中显示地震信号，也可以用下面的命令：
sfgrey < noise.rsf color=e> noise.vpl
sfpen noise.vpl
以上就给大家讲解了画局部相似图的全流程方法，相信能够帮助地震去噪研究的同学画出比较专业的图表。如果对安装指南和完整程序感兴趣，可以联系我交流。再见。
扩展阅读：
