# hk

**Repository Path**: codepanda64/hk

## Basic Information

- **Project Name**: hk
- **Description**: 使用远震P波接收函数，进行HK搜索，得到莫霍面深度和k(vp/vs)
- **Primary Language**: Python
- **License**: Apache-2.0
- **Default Branch**: master
- **Homepage**: None
- **GVP Project**: No

## Statistics

- **Stars**: 0
- **Forks**: 0
- **Created**: 2023-12-21
- **Last Updated**: 2024-05-05

## Categories & Tags

**Categories**: Uncategorized

**Tags**: None

## README

<!--
 * @Author: codepanda
 * @Date: 2023-12-21 16:11:21
 * @LastEditors: CodePanda
 * @Email: codepanda@qq.com
 * @LastEditTime: 2023-12-25 17:27:52
-->

# Anconda 或者 miniconda 的安装和使用

## 下载 Anconda 或者 miniconda

- [Anaconda3 for windows](https://repo.anaconda.com/archive/Anaconda3-2023.09-0-Windows-x86_64.exe)
- [Minicode3 for windows](https://mirrors.bfsu.edu.cn/anaconda/miniconda/Miniconda3-latest-Windows-x86_64.exe)

## 安装 Anconda 或者 miniconda

    运行安装程序下一步，具体步骤略

## 打开 Anaconda 命令行

- 打开命令行
  在 Windows 下
  开始菜单 -> Anaconda Prompt(anaconda) 或者 Anconda Powershell Prompt(anaconda)

- 增加科大源
  在命令行输入一下命令
  ```powershell
    conda config --add channels https://mirrors.ustc.edu.cn/anaconda/pkgs/free/
    conda config --add channels https://mirrors.ustc.edu.cn/anaconda/pkgs/main/
    conda config --add channels https://mirrors.ustc.edu.cn/anaconda/cloud/conda-forge/
    conda config --add channels conda-forge
    conda config --set show_channel_urls yes
  ```

## 下载程序包

```powershell
  # 需要先安装 git
  git clone https://gitee.com/codepanda64/hk.git
  cd hk
```

## 创建程序运行 Python 虚拟环境

- 创建 Python 虚拟环境

  > 如果安装的 Anaconda 打开命令行会有一个默认的 `base` 环境， 可以直接使用 `base` 环境，安装需要的 Python 包，如下：

  ```powershell
  # 如果安装 miniconda 跳过此步骤
    conda env update -f environment.yml -n base
  ```

  > 可以通过执行 `conda deactivate` 退出 `base` 环境，再创建新的环境

  ```powershell
    conda env create -f environment.yml
  ```

- 进入 Python 虚拟环境
  ```powershell
    conda activate hk
  ```

# 使用程序

`hk.py` 为入口程序， `hk.ini` 为配置文件

- hk.ini 文件说明

  - `[input]` 为输入配置信息
    - `rf_path` 为接收函数根路径，根据具体情况进行修改
    - `stations_name` 为需要处理台站的代码，需要和`rf_path`下的目录有对应，多个台站用逗号（英文）分隔，例如：`X1.53110,X1.53111`，或者使用 `*` 表示处理 `rf_path` 下的所有目录。
    - `best_dir_name` 为挑选好的接收函数，放到的目录名称，默认为 `best`
    - `rf_ext` 为接收函数的后缀名，默认为 `sac` 或 `eqr`
      > 输入目录的结构为： `rf_path`/`station_name`/`best`/`*.eqr`
  - `[hk]` 为 HK 搜索所需的参数
    - `model_name` 计算走时，射线参数使用的模型，默认为 `iasp91`
    - `hmin` 为搜索的最下深度，可根据地理位置的莫荷面深度通常的认知适当调整， 默认为 30
    - `hmax` 为搜索的最上深度，可根据地理位置的莫荷面深度通常的认知适当调整， 默认为 60
    - `kmin` 为搜索的最小 Vp/Vs， 可根据地理位置的 Vp/Vs 通常的认知适当调整, 默认为 1.5
    - `kmax` 为搜索的最大 Vp/Vs, 可根据地理位置的 Vp/Vs 通常的认知适当调整, 默认为 2.5
    - `num` 为搜索网格的数量
    - `vp` P 波的群速度, 可根据地理位置的 Vp 通常的认知适当调整，默认为 6.2
    - `pre_time` 接收函数截取的开始时间
    - `post_time` 接收函数截取的截止时间
    - `w1` 计算 Hk 时，Ps 的权重
    - `w2` 计算 Hk 时，PpPs 的权重
    - `w3` 计算 Hk 时，PsPs+PpSs 的权重
  - `[output]` 为输出参数
    - `out_root` 为输出路径
    - `hK_results_file_name` 为输出文件名
    - `hk_img_name` 为输出图像名

- 使用

  ```powershell
    python hk.py
  ```

- 输出文件
  `{station_name}_{model_name}_{hmin}-{hmax}_{kmin}-{kmax}_{num}_{vp}.txt`
  格式为： 台站名, best_H, best_H_sig, best_k, best_k_sig