SeisMan

合集

Fri, 21 Jun 2013 00:00:00 +0000

本页面收集了本博客中常更新的博文或系列性博文。地学学术期刊 SAC参考手册全球地形起伏数据总结个人偏好列表 Hi-net系列博文 CentOS7 配置指南 GMT

捐赠

Wed, 05 Feb 2014 00:00:00 +0000

如果您觉得本博客的原创博文或我维护的开源项目对您有用，可以考虑通过捐赠来支持本博客。免责声明本博客致力于知识和技术分享，不以盈利为目的，是个

关于

Sat, 22 Jun 2013 00:00:00 +0000

关于博主男，90年生，据说是典型摩羯座，USTC地球物理学博士关于SeisMan 博客名: SeisMan = Seismology Manual 博主名: SeisMan = A Man Learning Seismology 博文原则博客内容以原创为

参考文献管理最佳实践

Sat, 12 May 2018 00:00:00 +0000

科研过程中文章有两大用途：日常阅读和写文章时引用。本文讲讲我在写文章引用文献过程中遇到的问题以及我的解决办法。纯属经验之谈，望其他受此困扰的

制作 Linux USB 安装镜像

Sat, 07 Apr 2018 00:00:00 +0000

在下载 Linux ISO 镜像文件后，还需要制作 USB 启动盘方能安装。在不同的操作系统，制作所用的工具、过程不同。 macOS 下制作 USB 安装镜像将安装文件转换成 dmg 格式： $ cd

在Google Earth上绘制震源球

Sat, 03 Mar 2018 00:00:00 +0000

Google Earth可以方便地放大缩小区域范围。将震源球画在Google Earth上可以方便地查看不同地区地震的震源机制。本文介绍如何把震源球显示在Google Earth上。

原理其实很简单，

利用 obspy.imaging.beachball 绘制震源球并保存为PNG格式的图片
利用 simplekml 生成KML文件，其中包含一堆点，每个点对应一个地震，点的图标设置为地震对应的震源球图片即可

用 Mac 打造合适的科研环境

Mon, 29 Jan 2018 00:00:00 +0000

本文已迁移至 https://seismo-learn.org/seismology101/computer/macos-setup/

HinetPy: Hi-net 数据申请与处理模块

Sat, 01 Apr 2017 00:00:00 +0000

HinetPy 是一个用于从日本 NIED Hi-net 网站下载地震波形数据并进行数据格式转换的Python模块。

项目主页： https://github.com/seisman/HinetPy
官方文档： https://seisman.github.io/HinetPy

计算任意深度处下表面反射波走时

Fri, 31 Mar 2017 00:00:00 +0000

地震学中的某些研究会用到 SS 前驱波的理论走时，比如常见的 S410S、S660S。进一步推广到任意深度的下表面反射震相，比如 S1000S 代表在1000 km深度处的界面的下表面反射的S波，比如 P1200P 代表在1200 km深度处的界面的下表面反射的P波。

地震学中的坐标系

Thu, 12 Jan 2017 00:00:00 +0000

地震学中不同的作者定义了不同的坐标系，搞清楚各个坐标系的定义可以避免使用过程中出现错误。

SOD 入门教程

Sat, 24 Dec 2016 00:00:00 +0000

SOD，全称是 Standing Order for Data，是一个可以自动筛选、下载以及预处理地震数据的工具，其具有高度可定制的特点，可以满足日常科研中数据申请的几乎全部需求。

官方网站： http://www.seis.sc.edu/sod/
源码： https://github.com/crotwell/sod
源码语言： Java
最新版本： 3.2.8

本文是 SOD 的入门教程，目的是使读者对 SOD 有初步的认识，详细教程及用法见官方网站。

文件管理与备份

Wed, 24 Aug 2016 00:00:00 +0000

本文已迁移至 https://seismo-learn.org/seismology101/best-practices/file-organization/

使用 Web Services 获取地震数据

Sun, 14 Aug 2016 00:00:00 +0000

地震数据通常包含三类：地震目录、台站信息、波形数据。当然，还有一些其他的地震数据，比如震源机制等等。世界上最大的地震数据中心是 IRIS DMC。除了

SAC 文件读写模块 sacio_Fortran

Tue, 19 Jul 2016 00:00:00 +0000

sacio_Fortran 是一个用 Fortran 90 写的 SAC 读写模块。在自己的程序中调用该模块，即可实现 SAC 二进制文件的读写。项目主页： http://wangliang1989.github.io/sacio_Fortran/ 文件列表下载后，有如下文件: sacio.f9

订阅学术期刊

Mon, 04 Jul 2016 00:00:00 +0000

科研工作者的一大日常就是读文章。阅读最新的文章，可以了解一个学科最新的前沿进展。常见地学期刊列表见 https://link.seisman.info/journals/ 本文介绍我知道的可以及时获取最新文章的几

IRIS 的理论地震图生成引擎

Thu, 26 May 2016 12:41:18 +0000

通过看理论地震图有助于我们理解实际数据中的很多现象，但自己算理论地震图却不是一个简单的事情，既需要有相关代码还需要理解理论地震图计算的原理，

识别 TauP 输出中的多个 P 震相

Sun, 22 May 2016 11:46:41 +0000

在使用 TauP 计算 浅源近震 的 P 波走时时，输出中经常会出现很多个 P 震相。这个事情让人很困扰，因为搞不清这些 P 震相到底谁是谁。如果需要的只是 P 波初至，那就无所谓了，直接用最早的到时即可。有些时候需要计算某个特定 P 震相的走时，这就需要清楚地了解每个 P 震相所代表的真正路径了。

TauP 自定义模型

Sat, 21 May 2016 00:00:00 +0000

TauP 内置了常见的一维速度模型，比如 PREM、ak135、iasp91 等，这些标准速度模型位于 TauP 软件包的 StdModels 目录下。日常科研中，这些标准速度模型是不能满足需求的，因而需要自定义速度模型。

征集令

Wed, 13 Apr 2016 00:00:00 +0000

维护者征集下面列出了我发起的几个开源项目。挖了坑但是没有时间去填，所以需要征集志愿者参与维护。志愿者越多越好，项目开源的目的就是让所有人都可

用 psxy 绘制铁轨线

Wed, 30 Mar 2016 00:00:00 +0000

想要绘制黑白相间的铁轨线。既然是画线，肯定要用 psxy 了，但是 psxy 命令自身是不支持画这种黑白相间的线的，因而就需要使用一些技巧。具体技巧如下例：

判断 SAC 数据是否已去除仪器响应

Sun, 13 Mar 2016 00:00:00 +0000

有些时候，波形数据拿到手了，但是却不知道波形数据是否有被处理过，尤其是有没有去除仪器响应。

先说结论，判断一个 SAC 数据是否已经去除仪器响应的准则是：

若一个 SAC 数据中每个数据点的值都是整数，则这个数据的仪器响应没有被去除的概率为 99.9%

在极坐标下绘制深度剖面

Sat, 12 Mar 2016 00:00:00 +0000

在做大区域层析成像时，通常需要绘制深度剖面。一个简单的示例如下：

这样的深度剖面图实际上是在极坐标下绘制的。先直接给出绘制上图所使用的命令:

gmt psbasemap -R60/80/5571/6371 -JPa15c/70z -Bx5 -By200 -BWNse > test.ps

Hinet/Fnet 台站选择工具

Tue, 08 Mar 2016 00:00:00 +0000

Hinet 大概有不到 800 个台站，Fnet 大概有 150 个左右的台站。在从 Hinet 网站上下载数据时，默认是会下载所有台站的波形数据的。若只需要申请其中部分台站的数据，则需要使用 Hinet 网站自带的台站选择工具。

Hinet 网站自带的台站工具有如下优缺点：

图形界面，直观得显示了每个台站的位置
无法通过某些条件进行筛选，一次只能选中或删除一个台站
似乎需要一个 Java 的浏览器工具才能使用

因此，写了一个可以用于选择台站的 Python 工具，即 StationSelector.py，其是 HinetScripts 项目下中的一个脚本。

Hinet 配置检查脚本

Mon, 07 Mar 2016 00:00:00 +0000

HinetDoctor.py 是 HinetScripts 项目中的一个 Python 脚本。该脚本主要有如下功能：检查 Python 版本是否大于 3.3 检查是否已安装所有依赖检查配置文件是否可正常读取检查账号密码是否正确，或

用 SAC 快速拾取震相

Fri, 19 Feb 2016 00:00:00 +0000

震相拾取，是地震学的一个基本工作。SAC 提供了 plotpk 命令，可以用于手动拾取震相到时，用起来还算是比较简单方便。

关于博客无法访问的说明

Wed, 20 Jan 2016 00:00:00 +0000

最近，博客出现了无法访问的情况。就我自己遇到的情况来说，同一个实验室里，台式机连接了交换机，无法访问博客，笔记本连接了无线路由，则可以正常访

绘制颜色渐变的线条

Mon, 18 Jan 2016 00:00:00 +0000

想要绘制一条颜色变化的线段，下面是演示代码：

#!/bin/bash
# GMT v5.2.1
gmt makecpt -Crainbow -T-2/2/1 > lines.cpt
gmt psxy -JX15c/4c -R0/6/0/4 -B1 -Clines.cpt -W2p > test.ps << EOF
> -Z-1.5
1 2
2 2
> -Z-0.5
2 2
3 2
> -Z0.5
3 2
4 2
> -Z1.5
4 2
5 2
EOF

绘图效果如下图：

SAC 中将位移记录转换成速度记录

Tue, 12 Jan 2016 00:00:00 +0000

将位移记录转换成速度记录，本质上就是微分，在 SAC 中有三种方法可以实现。先生成测试用的数据: SAC> fg seis SAC> rmean; rtr; taper SAC> w seis.sac 直接对 diff 命令对波形做时间域微分:

在地图上绘制特定的经纬线

Fri, 11 Dec 2015 00:00:00 +0000

问题：想要在地图上只绘制某条经线或 / 和某条纬线，该如何实现？

GMT 中的 -B 选项可以用于控制是否绘制网格线以及网格线的间隔。如果精心设计 -B 选项中的网格间隔以及第一条网格线的偏移量，是可以使得整张图中只有一条经线或 / 和纬线并保证该经线 / 纬线落在合适的位置，但是这种方法相对来说比较麻烦且很难用于批量绘图。

SOFI2D 笔记

Tue, 08 Dec 2015 00:00:00 +0000

SOFI2D 是一个二维 PSV 交错网格有限差分代码。要使用该程序，至少需要阅读如下几篇参考文献：

Virieux, J. 1986. P-SV wave propagation in heterogeneous media: velocity-stress finite- difference method. Geophysics, 51(4), 889–901.
Levander, A.R. 1988. Fourth-order finite-difference P-SV seismograms. Geophysics, 53(11), 1425–1436.
SOFI2D User Guide

具体的用法及每个参数的意义在官方文档已经介绍的很详细，本文是我使用 SOFI2D 过程中的经验总结。

注：本文使用的是 2015 年 10 月 28 日发布的 SOFI2D。

绘制事件与台站间的连线

Wed, 18 Nov 2015 00:00:00 +0000

地震学中，经常需要绘制事件到台站的连线或台站与台站之间的连线。连线沿着大圆路径，理论上用 psxy 命令即可实现，实际绘图时要更复杂一些。

GMT 5.2.1 在 Linux 下的安装

Fri, 13 Nov 2015 00:00:00 +0000

**本文不再维护更新，最新版本的GMT的安装方法见 http://docs.gmt-china.org/install/linux.html **

注意： 为 Ubuntu 和 CentOS 用户写了一个快速安装脚本，见 gmt-easy-installer。有兴趣的可以试着安装，也欢迎测试各种不同的平台。

本文介绍如何在 Linux 下安装 GMT 5.2.x。

GMT 5.2.x 系列的安装与 GMT 5.1.x 系列的安装稍有不同，故而在前文基础上单独写一篇。

地震数据申请

Tue, 27 Oct 2015 00:00:00 +0000

申请地震波形数据，按照具体的需求，大致可以分为两大类，即事件波形数据和连续波形数据。

这两者的主要区别在于如何确定要申请的数据时间窗，前者需要震相到时信息，后者则不需要。

rdseed 用法笔记

Tue, 13 Oct 2015 00:00:00 +0000

rdseed 是 IRIS 提供的工具，用于读取 SEED 格式，从中提取出波形信息，并将波形数据保存为 SAC、AH、CSS、SEGY 或 ASCII 等多种数据格式。

判断点在多边形内

Sat, 10 Oct 2015 00:00:00 +0000

在处理数据或绘图时，有时会需要判断点是否在一个封闭多边形内。比如想要筛选出某个多边形区域内的所有地震事件。GMT 的 gmtselect 命令可以实现这一需求。

用 Wilber 3 申请地震波形数据

Mon, 28 Sep 2015 00:00:00 +0000

Wilber3 是 IRIS 提供的一个用于从 IRIS 申请地震事件波形数据的网页工具。

官方页面： http://ds.iris.edu/wilber3/find_event
优点：网页界面，简单直观
适用范围：仅适用于申请单个或少量地震事件的波形数据

本文将以 2015 年 5 月 12 日的尼泊尔 7.3 级地震为例，展示如何利用 Wilber3 申请该事件的波形数据。

一些说明

Sun, 27 Sep 2015 00:00:00 +0000

最初写博客只是想写着玩，没想到一写就写了四年，估计还会继续写下去。写的这些东西，对于很多地震学方向的人来说还是很有用的，所以这个博客以及我本人也受到了很多人的关注。一直以来，我都尽量避免在网络上公开自己的真实身份，这次因为参加培训的缘故，很多读者第一次接触到了现实中的我。我想，有些事情我应该说明或解释一下。

博客开通微信公众号

Fri, 25 Sep 2015 00:00:00 +0000

SeisMan 博客刚刚开通了微信公众号，欢迎大家关注。

微信公众号的名字是 “SeisMan 博客”，微信号是 “seismaninfo”。

修改 SAC 所允许的最大文件数目

Sat, 12 Sep 2015 00:00:00 +0000

SAC 在处理数据时，一次性最多只能读入 1000 个 SAC 文件。当读入过多的 SAC 文件时，会出现如下警告: Max files: reading first 1000 files. 并只读入文件列表中的前 1000 个文件。想要突破最多

GMT 4 在 Mac 下的安装

Sat, 05 Sep 2015 00:00:00 +0000

GMT 官方没有为 Mac 提供 GMT4 的安装包，因而只能通过其他方式安装。Mac 下有很多第三方软件管理工具，比如 MacPorts 和 Homebrew。MacPorts 和 Homebrew 中自带最新版本的 GMT，可以直接安装，是比较推荐的使用方式。某些情况下可能会希望自己手动编译，比如自己修改了 GMT 源码。

因而，本文会介绍四种安装 GMT4 的方法，列举如下

用 Homebrew 安装 GMT4
用 MacPorts 安装 GMT4
基于 Homebrew 安装 GMT4
基于 MacPorts 安装 GMT4

GMT 添加断层名

Fri, 21 Aug 2015 00:00:00 +0000

在 GMT 中绘制断层很简单，只要有断层的坐标数据，直接用 psxy 即可绘制。

有时不仅要绘制断层，还需要在断层线的旁边加上断层名，最直接的想法是用 pstext 将断层名写上去。使用 pstext 添加断层名时会遇到两个问题：

text 的位置；
text 的角度（一般情况下，断层名与断层线平行）；

如果只是几条断层的话，可以手动指定并微调 text 的位置和角度，以达到想要的效果。如果要标记的断层名有点多，或者想要批量绘制不同区域的多张断层图，手动微调的方法就有些麻烦了。当然，如果真有需求的话，可以手动微调，建立一个断层名的数据库，以后用起来就方便多了。

GMT 绘制图中图

Sun, 16 Aug 2015 00:00:00 +0000

英文一般称为 inset map 或 insert map，中文还不确定有没有专门的名词，一般常称为图中图、小图、嵌入图。

所谓 insert map，即在绘制小区域地图时，在图上的某个角落里绘制一张更大范围的地图，并在大范围地图上绘制一个方框，以更清楚地表明小区域的地理位置。

若是了解 GMT 中图层的概念以及图层叠加的原理，就可以很容易地把 insert map 画出来。由于 GMT4 和 GMT5 语法上的差异，这里同时给出 GMT4 和 GMT5 两个版本的绘图脚本。二者原理是差不多的，只是语法上有些差异，因而下面只解释 GMT5 脚本中的细节，对 GMT4 脚本不做过多介绍。

GMT 区域填色

Wed, 05 Aug 2015 00:00:00 +0000

提出问题

GMT 中可以用 psxy 绘制多段数据，画出多个闭合的多边形。一个示例命令如下:

gmt psxy input.dat -R0/360/-90/90 -JX15c/10c -Bx60 -By30 -W1p -L > test.ps

安装 evt2sac

Fri, 31 Jul 2015 00:00:00 +0000

evt2sac 用于将 NEIC 的 EVT 格式的数据转换为 SAC 格式。

主页： https://seiscode.iris.washington.edu/projects/evt2sac

博客被 Google 认定为 “对用户非常有用” 啦

Wed, 29 Jul 2015 00:00:00 +0000

无聊的时候会在搜索引擎里以关键词 “seisman” 搜自己，看看在搜索引擎里的排名是怎样的。从一开始的完全搜不到，到后来 Google 搜索结果中排第一，再

地球物理学家列表

Thu, 23 Jul 2015 00:00:00 +0000

本页面已迁移至请访问: https://link.seisman.info/scientists/

pssac2 使用教程

Sun, 19 Jul 2015 00:00:00 +0000

pssac2 是 pssac 的升级版。相比于 pssac 而言增加了更多的功能，选项也更加符合 GMT 的风格。

pssac2 与 pssac 的选项之间有很多相似之处，所以本文只介绍二者之间的不同之处，读者应首先阅读《pssac 使用教程》。

在地图上绘制波形

Sat, 18 Jul 2015 00:00:00 +0000

将 SAC 波形数据绘制在地图上也算是常见的需求之一。最常见的情况是，在地图上有一些台站，需要将波形画在台站的附近。

此问题的难点在于，在画地图时，-R 指定的是区域范围，即横轴是经度，纵轴是纬度，而 pssac 在绘制波形时，横轴是时间，纵轴是振幅。因而想要将波形放在特定的位置，不是一件容易的事情。

pssac 使用教程

Fri, 17 Jul 2015 00:00:00 +0000

注意：本文仅针对 Prof. Lupei Zhu 的 pssac 程序，其仅可在 GMT 4下使用。对于GMT5用户，GMT 在 5.3.1 之后自带了 pssac，可以使用 gmt pssac 直接调用。

GMT 绘制无刻度轴

Wed, 15 Jul 2015 00:00:00 +0000

提出问题

要绘制一张底图，要求

左边和下边有标注和刻度；
右边和上边无标注且无刻度；

最直接的想法是使用 -BWSen，绘图得到的效果图如左图所示，右边和上边虽然没有了标注，但是刻度线依然存在。实际想要实现的效果是右图。

SEED 格式转 SAC 格式

Fri, 19 Jun 2015 00:00:00 +0000

本文会介绍将 SEED、miniSEED 等格式转换为 SAC 格式所使用的工具和用法。需要注意，本文并不会详细介绍工具的具体语法，只介绍常用的功能，其他功能及选项，读者应自行阅读相关文档。

GMT5 自定义坐标轴

Sat, 06 Jun 2015 00:00:00 +0000

问题的提出

GMT 绘图时可以使用 -B 选项指定每个坐标轴的刻度、标注以及网格线。默认情况下的设置有时可能无法满足需求，比如想要在特定位置设置刻度，比如任意格式的日期、非 log10 的刻度等，就需要自己设置了。

底图边框与网格线

Tue, 19 May 2015 00:00:00 +0000

本文介绍 GMT 的 - B 选项的两种稍特殊的用法，以实现只绘制边框和只绘制网格线的目的。

两点说明：

示例采用 GMT5 语法，但方法对于 GMT4 也同样适用；
想不起来当初为什么会有只绘制网格线的需求了，所以这里只写如何实现，而不说为什么会有如此奇怪的需求。

球状模型下的 tau 积分

Fri, 15 May 2015 00:00:00 +0000

在计算震相走时的过程中，对 $\tau(p)$ 函数的积分是其中比较重要的一步，具体 $\tau(p)$ 函数如何使用，在地震学教材里都有介绍。

CPS330 的安装

Tue, 12 May 2015 00:00:00 +0000

本文已迁移至 https://seismo-learn.org/software/cps/

震相走时的计算

Sat, 02 May 2015 00:00:00 +0000

震相走时的计算是地震学的基本问题之一。

目前，计算一维分层模型下的震相走时，主要使用两个程序包：ttimes 和 TauP。 ttimes 相对来说有些古老，所以多数情况下更推荐使用 TauP。但 TauP 不支持 PKPab、 PKPbc 这种震相名，所以偶尔可能还是需要使用 ttimes。

ttimes 和 TauP 都是根据 Buland and Chapman（BSSA,1983）提出的理论框架来实现的，本文是对这篇文章的核心内容的简单总结。

调用 SAC 进行 Hilbert 变换

Wed, 29 Apr 2015 00:00:00 +0000

需要在 C 程序中对数据做 Hilbert 变换，自己写显然是不可能的啦，重复造轮子不说，写的还不一定对。找了一些代码，发现都写的好复杂，后来发现 SAC 在函数库中提供了 Hilbert 变换的接口，可以直接调用。

firtrn 可以用于对数据进行 Hilbert 变换，但是官方没有给一个比较直白的例子。下面就贴一个我写的例子啦。

GMT 5.1.2 在 Windows 下的安装

Sat, 18 Apr 2015 00:00:00 +0000

本文介绍如何在 Windows 7 下安装 GMT5.1.2。

一些产品的推广链接

Sat, 11 Apr 2015 00:00:00 +0000

本文列出了我在使用的部分产品/服务的推广链接，读者若点击我的推广链接，注册并使用相关的产品/服务，则我可以获取一定的好处。

从 3D 数据中截取剖面

Wed, 08 Apr 2015 00:00:00 +0000

本文所说的 3D 数据是指 XYZD 数据。常见的情况是数据中包含四列：经度、纬度、深度、值。

GMT 几乎是不可以处理 3D 数据的，比如与网格相关的命令只支持 2D netCDF 文件。另外，GMT 绘制 3D 图的效果也很一般，因而通常情况下都需要从 3D 数据中截取剖面，绘制剖面图。

本文试着介绍如何从 3D 数据中截取自己想要的剖面。

全球地形起伏数据：SRTM15_PLUS

Tue, 07 Apr 2015 00:00:00 +0000

SRTM15_PLUS 是 SRTM30_PLUS 的升级版，其空间分辨率为 15 弧秒，相当于 0.5km。其包含了陆地高程和海洋深度数据，具体信息参考 SRTM30_PLUS 的介绍。数据下载： README Google KMZ 14GB netCDF 数据 netCDF 数据

修改 GMT 生成的 PS 文件的标题

Sat, 14 Mar 2015 00:00:00 +0000

本文整理自： http://gmt.soest.hawaii.edu/boards/1/topics/1820

需要说明的一点，本文是要修改 PS 文件的标题，而不是 PS 文件的文件名。

用编辑器打开任意一个 GMT 生成的 PS 文件，查看文件的前几行，其中一行大概如下:

%%Title: GMT v5.1.1 (r12968) [64-bit] Document from psxy

这就是所谓的 PS 文件的标题。当用 evince 或 gv 打开 PS 文件时，该标题会在软件的某个位置显示。

修改 Y 轴的坐标标注的方向

Wed, 11 Mar 2015 00:00:00 +0000

在用 GMT 绘制底图时，通常 Y 轴的标注都是水平的，比如如下绘制底图的命令:

gmt psbasemap -R0/50/-10/10 -JM10c -Bx10 -By10 > test.ps

效果如下图

有些时候想要让 Y 轴的标注逆时针旋转 90 度，即 Y 轴的标注方向与 Y 轴平行而不是与 X 轴平行。想要实现的效果如下图。

征集《SAC 参考手册》维护者

Sat, 07 Mar 2015 00:00:00 +0000

如题，希望能够征集到多名志愿者参与《SAC 参考手册》的维护。

fk 用法笔记

Sat, 28 Feb 2015 00:00:00 +0000

fk 是 Prof. Lupei Zhu 写的一个用于计算 水平分层模型 下的理论格林函数并合成理论地震图的代码包。代码是开源的，可以直接编译使用。

慢度与射线参数

Thu, 05 Feb 2015 00:00:00 +0000

本文介绍几个专业名词：

慢度
水平慢度 / 垂直慢度
射线参数
视速度 / 视慢度

这些名词其实很简单，但是在大三一直到研一的几年间，这几个名词在我脑中的概念是有些混乱的，我把混乱的原因归结于不同课本对同一个专业名词乱用了不同的词汇。故而在这里总结一下。

不同比例尺曲线的画法

Sun, 01 Feb 2015 00:00:00 +0000

本文由 Joe Wang（cjmuqiao@163.com）投稿。

目标效果图

强大的跨平台 PDF 处理工具：cpdf

Tue, 27 Jan 2015 00:00:00 +0000

日常生活中对 PDF 的最基本的操作大概就是合并和切割了。之前一直在用强大的 PDFtk，但是由于一些底层依赖的问题，PDFtk 在短时间内是不可能出现在 CentOS7 下了，只能寻找新的 PDF 替代工具了。

Coherent PDF，简称 cpdf，是个用于处理 PDF 的命令行工具。该软件支持 Windows、Linux 和 Mac，且对于非商业用途是免费的。

主页： http://community.coherentpdf.com/

用 taup_time 计算震相走时及射线信息

Sat, 24 Jan 2015 00:00:00 +0000

本文已迁移至 https://seismo-learn.org/software/taup/taup_time/

GMT 中添加注释和标注

Sun, 28 Dec 2014 00:00:00 +0000

在日常绘图时经常会需要在已经画好的图上添加一些注释或标注，比如在某个位置加文字、箭头或者连线。一般情况下，若注释或标注出现在单个底图的内部，则可以很容易通过 pstext 和 psxy 来实现；若需要在底图外添加注释或标注，则会复杂一些。本文试图解决这个问题。

Windows 下使用 GMT 的正确姿势

Wed, 10 Dec 2014 00:00:00 +0000

警告：由于某些未知的更新，本文所推荐的notepad++终端插件无法在终端中显示GMT输出的各种信息，因而请勿完全按照本文所说的进行操作！！！

我几乎是完全在 Linux 下工作的。有过几次在 Windows 下使用 GMT 的经历，个人觉得用户体验是非常糟糕的，Windows 自带的工具太简陋，用来画 GMT 图太繁琐，浪费好多时间。虽然未来我也几乎不会在 Windows 下写 bat 脚本，但还是想将自己的一些经验写下来，希望能够提高 Windows 下 GMT 的用户体验。

在正式开始本文之前，首先要认清几个基本的事实。如果你无法认同这几个事实，那么就没有再读下去的必要了：

GMT 最初是在 Linux 下开发的 ** 纯命令行 ** 工具。在 Linux 下使用 GMT 要比在 Windows 下方便很多。所以，如果能在 Linux 下使用，就不要在 Windows 下使用；
大多数图都不可能用一个 GMT 命令就完成，因而需要将一系列 GMT 命令写到 bat 脚本中；
将 GMT 命令写入 bat 脚本的目的有两个：首先是记录 GMT 绘图的过程，其次才是批量绘图；
写 bat 脚本的过程，大部分时间和精力都是花在命令的调试和图形的微调上；
花半个小时时间折腾一下工具，就可以换取更好的用户体验以及更高的调试效率，这是一笔很划算的交易；

GMT 边框发虚的原因及解决办法

Tue, 09 Dec 2014 00:00:00 +0000

GMT 绘制的 PS 格式的图片一向以高精度著称。在画图的时候偶尔会遇到，PS 图片中其他地方精度都很高，唯独边框看上去有些发虚，精度比较低。

上次遇到这个问题的时候，因为着急画图，所以也没细想，直接把边框发虚的图用上了。今天再次遇到这个问题，认真的想了想，找到了原因所在，也找到了解决办法。

CentOS 7 下的软件安装方法及策略

Sun, 23 Nov 2014 00:00:00 +0000

2010 年开始正式接触 Linux，入门发行版是 Ubuntu 10.10，后来过渡到 Ubunu 11.04，这其中也尝试了很多其他主流的发行版。进入实验室之后，开始用 CentOS 5，然后是 CentOS 6，现在进化到 CentOS 7。

使用了四年的 Linux，前三年都是在瞎折腾，浪费了不少时间，也得到了不少经验与教训。现在可能是真的老了，已经不愿意折腾了，只希望配置好一个系统之后，就一直使用下去。

利用 taup_setsac 为 SAC 文件标记理论到时

Mon, 10 Nov 2014 00:00:00 +0000

本文已迁移至 https://seismo-learn.org/software/taup/taup_setsac/

利用 taup_pierce 计算射线穿透点

Sat, 08 Nov 2014 00:00:00 +0000

taup_pierce 是 TauP 提供的一个命令行工具，用于计算各种震相在不同深度的穿透点信息。本文通过几个例子介绍一个 taup_pierce 的使用方法。

TauP 的安装

Wed, 08 Oct 2014 00:00:00 +0000

本文已迁移至 https://seismo-learn.org/software/taup/install/

rdseed 的安装

Tue, 07 Oct 2014 00:00:00 +0000

rdseed 用于将 SEED 格式的数据转换为 SAC 格式。

CentOS 7 下安装 WPS Office

Wed, 01 Oct 2014 00:00:00 +0000

在 Microsoft Office 盛行的今日，Linux 下一直没有一个能与之完美兼容的 Office 套件。大多数发行版默认安装的 OpenOffice 或 LibreOffice，也只是刚刚达到能用的地步，界面以及功能的差异还是很明显的。

或许 WPS Office for Linux 是个不错的选择。

中文官网： http://linux.wps.cn/
英文官网： http://www.wps.com/linux/

CentOS 7 安装 fcitx 中文输入法

Sat, 20 Sep 2014 00:00:00 +0000

本文不再更新

目前不确定是否可行，具体参考： https://github.com/1dot75cm/myrepo

前面 CentOS7.0 下安装小小输入法一文中已经介绍了如何在 CentOS7 下安装小小输入法来输入中文，实际用起来还是有不少不如意的地方。

今天恰好发现了 CentOS 7 下安装 fcitx 输入法的简单方法，如下。

Linux 下合并 PDF

Thu, 18 Sep 2014 00:00:00 +0000

经常需要将多个 PDF 合并为一个文件，这里收集了网上给出的一些方案，并总结一下。

推荐使用cpdf

将博文以 PDF 格式导出

Sun, 14 Sep 2014 00:00:00 +0000

在每篇博文的标题下方最右侧，有一个 PDF 的图标，点击即可下载当前博文的 PDF 版本。

对于每篇博文，利用 pandoc 将 markdown 格式的源码转换为 PDF 格式，并将其上传到七牛云中，供读者下载。

Hi-net 生成仪器响应 SAC PZ 文件

Sat, 13 Sep 2014 00:00:00 +0000

前文 Hinet 仪器响应中介绍了 Hi-net 的仪器响应的详细信息。据此，用 Python 写了一个脚本，将 Hi-net 的 Channel Table 文件转换为常用的 SAC PZ 文件。

脚本位于：ch2pz.py。

Hi-net WIN32 转 SAC 脚本

Fri, 12 Sep 2014 00:00:00 +0000

该脚本用于实现 Hinet WIN32 格式到 SAC 格式的转换，脚本位于 rdhinet.py。

Hi-net win32tools

Sun, 07 Sep 2014 00:00:00 +0000

从 Hi-net 获取的波形数据为不常见/用的WIN32 格式，Hi-net 提供了 win32tools 来实现波形数据的基本处理和转换。

Hi-net 的仪器响应

Sat, 06 Sep 2014 00:00:00 +0000

前文介绍了 channel table 文件的格式以及每一列的含义，其中 7-13 列记录了 channel 的仪器信息。由于 Hi-net 并没有像 IRIS 那样提供 RESP 或 SAC_PZ 文件，因而需要用户根据 7-13 列的信息手动推导仪器响应。

Hi-net Channel Table 文件

Fri, 05 Sep 2014 00:00:00 +0000

从 Hi-net 下载连续波形数据会得到 ZIP 压缩文件，其中包含了两个后缀为 ch 的文件，即 Channel Table 文件。两个文件内容是相同的，区别在于一个是 euc 编码，一个是 sjis 编码。只需要用其中一个即可，这里选择 euc 编码的文件。

Channel Table 文件描述了每个 channel 的仪器信息，某种程度上等效于常见的 RESP 或 SAC PZ 仪器响应文件。但是实际上 Channel Table 有一定的局限性，无法理解 Hi-net 为什么会用这种格式来指定仪器响应信息。

尤其需要注意的是，对于 Hinet 而言，channel table 可以完整的描述仪器响应，而对于 Fnet 而言，channel table 给出的信息不足。因而关于 channel table 的内容，仅适用于 Hinet，不适用于 Fnet，不确定是否适用于其他台网。

Hi-net WIN32 格式

Thu, 04 Sep 2014 00:00:00 +0000

对于连续波形数据而言，从 Hi-net 下载的文件为 zip 压缩文件，其中包含了多个后缀为 cnt 的文件。每个 cnt 文件中包含了多台站多分量的一分钟的波形数据。这些 cnt 文件

Hi-net 事件波形数据

Wed, 03 Sep 2014 00:00:00 +0000

Hi-net 基于 JMA Unified Hypocenter Catalog 提供的事件信息，提供了事件波形数据。事件波形申请地址事件波形下载地址查看网页源代码可以发现，申请的实质是调用另一个 openRequest 函数，其中

Hi-net 走时和震源机制数据

Tue, 02 Sep 2014 00:00:00 +0000

Hi-net 基于 JMA Unified Hypocenter Catalog，提供了到时数据以及震源机制解目录。

数据下载地址数据格式官方文档

本文不谈数据格式，只讨论数据的下载问题。

JMA 统一地震目录

Mon, 01 Sep 2014 00:00:00 +0000

JMA 提供了日本及其周边地区的区域地震目录。JMA首先提供了一个初级的地震目录，并不定期发布最终版的地震目录。 Hi-net 网站中可以直接获取 JMA 初级地震目录

Hinet 连续波形数据申请及下载的脚本实现

Sat, 30 Aug 2014 00:00:00 +0000

前文已经说到，Hi-net 连续波形数据申请的本质可以简化为如下几步：

指定参数，构建 query string
向 Hi-net 发送数据请求
等待数据准备
下载数据

Hi-net 连续数据申请的源码分析

Fri, 29 Aug 2014 00:00:00 +0000

上一篇博文介绍了 Hi-net 数据申请的流程以及其中的一些限制。从表面上看，数据申请的过程其实就是鼠标点点点。

Hi-net 连续波形数据

Wed, 27 Aug 2014 00:00:00 +0000

本文将介绍如何从 Hi-net 申请并下载连续波形数据，并介绍 Hi-net 数据的一些细节。

连续波形数据申请地址： https://hinetwww11.bosai.go.jp/auth/download/cont/?LANG=en

日本高密度地震台网 Hi-net

Mon, 25 Aug 2014 00:00:00 +0000

Hi-net，全称为 High Sensitivity Seismograph Network，是 NIED 下的一个覆盖日本全境的高质量、高密度的短周期地震台网。

Python 发送邮件到 BREQ_FAST

Sat, 23 Aug 2014 00:00:00 +0000

利用 Python 的标准库 smtplib 实现了通过脚本发送邮件的功能。

源码地址：sendmail.py

我所使用的软件/服务列表

Tue, 05 Aug 2014 00:00:00 +0000

本文罗列我个人正在使用或曾经使用过的一些软件或服务。

博客托管的一些调整

Thu, 17 Jul 2014 00:00:00 +0000

博客托管一直是一个比较头疼的问题。作为一个中文博客，访问者必然主要来自于国内，也有少量访问者来自于美国、加拿大等国家的中国人。如果将博客托管在国内主机，可以保证国内大部分人的访问速度，但是托管在国内服务器上的网站必须进行备案，想想就觉得很麻烦；如果将博客托管在国外主机，则大部分国内读者访问网站的速度会很慢，甚至可能由于墙的问题而完全无法访问。

彻底卸载 Intel Parallel Studio

Wed, 16 Jul 2014 00:00:00 +0000

Intel 的开发工具默认安装在 /opt/intel 目录下，由于一时手贱，用 sudo rm -rf /opt/intel 把整个 intel 文件夹给删除了。在重新安装 Intel Parallel Studio 的时候一直提示 “已安装该产品”，进而导致无法

CentOS 7 配置指南