关于seed文件的读取方法

最新推荐文章于 2024-10-15 19:19:44 发布

原创

最新推荐文章于 2024-10-15 19:19:44 发布 · 2.4k 阅读

3 ·

CC 4.0 BY-SA版权

文章标签：

#算法 #go语言 #seed文件读取

本文介绍了如何使用Golang读取地震Seed文件的头信息，包括关键数据块如010、011、012、050、052、070和074等，这些信息对于正确解析和后续的数据内容读取至关重要。通过顺序读取和关键字定位，可以确定数据块的起始位置和长度。种子文件的头信息包含台站标识码、时间信息和数据记录范围等，为地震测震和强震数据的处理提供基础。

本文是我捣鼓了一个多月的总结，目的是自己的备忘录，也可供系统内有需求的同仁以参考。

这里的seed文件的读取方法只限于地震测震和强震数据，其他应用领域的数据没有测试过所以也不知道是否适用。

参考资料《地震波形数据交换格式》和 jopens软件程序

程序实例语言：golang

这里的seed文件是指以seed为后缀名的文件，该格式广泛应用与数字地震波形数据存储和数据交换。其编码结构在《地震波形数据交换格式》中有详细说明，在这里不做赘述，有兴趣者可以下载《地震波形数据交换格式》自行研究。本文只介绍在seed文件读取解码过程中的关键点和方法步骤（只少我是认为这些是比较关键的）

seed文件一般是由两部分组成，一部分是头文件内容（使用ASCII格式存储），一部分是数据内容（使用二进制格式存储）；头文件信息的正确读取直接决定着后面的数据内容是否可以正常读取。

头文件内容是由多个数据块组成，数据块长度不固定，其数据格式大体是“数据块类型”-》“数据块长度”-》“数据块存储的各类信息” 如下所示：

卷台站头段索引子块 [11]
0110021001JXH  000005
解析：
011   子块类型
0021  子块长度
001   台站数
JXH   台站标识码（这里是5位有两个空格）
000005  台站头端顺序号   （S标志）

卷时间片索引子块 [12]
012006300012023,101,16:00:00.0000~2023,102,16:00:00.0000~000012
解析：
012   子块类型
0063  子块长度
0001  卷中时间片数量
2023,101,16:00:00.0000   时间片的开始时间
~     分隔符
2023,102,16:00:00.0000   时间片的结束时间
~     分隔符
000012                   时间片头段顺序号（T标志）

数据块和数据块之间一般没有间隔，数据块的识别我是通过文件顺序读取检索关键字的方式处理的。代码如下：

// 按照关键字查找返回信息字段起始位置和信息长度
// file 文件指针；p_start 游标开始位置; keyword 关键字(关键字后三位必须是子块类型码)
func (s *HeadInfo) GetInfoSL(file *os.File, p_start int, keyword string) (int, int, error) {
	scanner := bufio.NewScanner(file)

	file.Seek(int64(p_start), 0)
	str_length := 0
	p := 0
	buf := make([]byte, 0, 64*1024)
	scanner.Buffer(buf, 1024*1024)
	for scanner.Scan() {

		str := scanner.Text()
		ide_s := strings.Index(str, keyword)
		if ide_s != -1 {

			p += ide_s
			buffer := make([]byte, 4)
			n, err := file.ReadAt(buffer, int64(p_start+p+len(keyword)))
			if err != nil {

				//panic(err)
				return 0, 0, err

			}
			str_length, err = strconv.Atoi(string(buffer[:n]))
			if err != nil {
               log.Println("keyword:",keyword)
				//panic(err)
				return 0, 0, err

			}