基于OpenCV模拟地震记录矢量化软件.docVIP

基于OpenCV模拟地震记录矢量化软件 　　摘 要：随着地震观测技术与手段的不断更新，数字化地震观测逐渐取代传统的模拟地震观测。为了保存历史?稻荩?如何有效提取模拟地震记录中的珍贵信息是地震工作者的研究重点。该项目提出了基于OpenCV的模拟地震记录矢量化算法，使用C#开发相应软件。该软件利用有震波形区域提取方式可以快速、准确地将模拟地震记录矢量化。 　　关键词：模拟地震记录；数字化；OpenCV；有震波形区域提取；波形追踪 　　0 引言 　　世界上第一张地震图要追溯到1889年，英国人米尔恩和尤因利用安置在德国波茨坦的现代地震仪记录下了发生在日本的一次地震。此后70年的时间里，人类观测的地震信息都是以模拟地震记录的方式存在。模拟地震记录蕴含着大量宝贵的地震信息，具有不可估量的价值。直到19世纪70年代，美国率先建成了由十个数字化台站组成的数字化台网。随后数字地震记录技术逐渐取代模拟地震记录技术，一直发展沿用至今。这造成了地震数据应用研究中，模拟地震信息与数字地震信息之间有一条明显的数据断裂。如今大多数模拟地震记录存放在仓库里经受着自然的侵蚀，蕴藏着宝贵信息的模拟地震记录无法被计算机处理。若能有效提取这些模拟地震信息，将会推动地震研究工作快速发展。 　　1 研究现状 　　目前为止，对于模拟地震记录的研究主要分为两个阶段：存储与矢量化。 　　存储是将模拟地震记录中的信息由易损坏的纸质信息转化为方便存储与备份的数字图片信息。2012年开始，哈佛大学投入大量精力进行相关模拟地震记录电子存档工作，该项计划被命名为哈佛大学地震波形归档计划。2015年，河北地震局蒋宏毅等人研发了一套模拟地震记录数字化存储管理系统，实现了模拟地震记录的数字化存储与查询。但数字存储只是将纸质模拟地震记录扫描为数字图片，没有进行矢量化，计算机仍无法识别和处理模拟地震记录，也无法挖掘记录中的宝贵信息‘5]。 　　许多科研工作者进行模拟地震记录矢量化的研究：Teves-Costa（ 1999）使用CADcore和AutoCAD开展模拟地震记录的矢量化工作；Pintore（2005）基于人工神经网络提出了Teseo方法，可以手动和自动矢量化模拟地震记录；徐涛等（2014）利用Matlab GUI开发了一个交互的模拟地震矢量化程序；王茂发等（2016）提出了CSF算法，并开发了相应的程序。虽然矢量化研究工作取得了显著成绩，但许多矢量化工作只是将一小部分含有复杂波形的图片截取出来单独进行矢量化研究，脱离了整张图片使得后期波形拼接与时间拼接变得难以实现。 鉴于跨平台计算机视觉库OpenCV在人脸识别、车牌识别、动态手势检测等方面的准确性与高效性，因此，本文基于OpenCV设计了新的模拟地震记录矢量化算法，并使用程序设计语言C#开发软件。其中包含分块提取波动区域、自动提取复杂波等特有的功能模块，能够快速、准确地处理整张模拟地震记录。 　　2 软件矢量化过程 　　2.1 流程 　　本软件实现模拟地震记录矢量化的过程共分为8个步骤（图1所示）：加载地震图纸、图纸二值化、搜索起点、有震波形区域提取、扫描波动区域、判断波形类型、拼接和反演。下面给出每一个步骤的具体描述。 　　2.2加载地震图纸 　　将整张图纸加载至矢量化软件中是模拟地震记录矢量化最基本的工作。早期的模拟监测是将一段时间（一般为1天）内3个方向（东西、南北、上下）的地震波动记录在一张图纸上。观察整张图（如图2）纸可知，波形密集，包含的信息量巨大，该项目使用C#常用控件PictureBox实现了软件的图片加载功能，可以加载、处理整张模拟地震记录。 　　2.3 图纸二值化 　　许多因素如光线、图纸的新旧程度等，会对模拟地震记录扫描而成的图片产生影响，所以必须对图纸进行二值化处理，它是在图片预处理阶段减少图片噪声非常有效的方式。一张灰度的光栅图片类似一个二维的数组，每一个元素在数组中的值在0-255之间。利用OpenCV的Threshold函数对图片进行二值化处理：设置一个阈值（目前实验获得经验数值为155），将大于阈值的像素点置为255，小于阈值的置为0，如图3所示： 　　2.4搜索起点 　　为了确保数据的一致性，软件在提取模拟地震记录中的信息时，始终在同一个坐标系下处理。以图片的左上角为原点，图片的上边缘为x轴正值方向，图片的左边缘为y轴正值方向。取原点附近的一个点c（x，0）为起点，沿着垂直方向向下搜索，每次步长为一个像素点。搜素算法具体描述如下： 　　（I）水平方向 　　水平方向x值保持不变（公式1），为了避免图片左侧空白边缘对搜索结果的影响，一般x取值稍大于0（实验取值为10）。 Xi=x（i=l，2…m） 　　（1） 　　（2）垂直方向 　　垂直方向y值由0不断变大