You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 6 Next »

This procedure outlines the process for downloading waveform and instrument data from GeoNet and using it to compute ground motion classifications. This is followed by computing ground motion intensity measures for the highest quality data.

Step-by-step guide

  1. Download a list of earthquakes from GeoNet in CSV format: https://quakesearch.geonet.org.nz
    1. Several lists can be appended to create a master file.
  2. Use mseed_write.py to load station response information and download mseed data.
    1. The data can be subset within the program to match specific years and/or magnitude ranges.
  3. Run GMC (the ground motion classifier) on the downloaded data.

    1. Create Konno matrices from gm_classifier directory. (example below):

      1. python gm_classifier/scripts/gen_konno_matrices.py dir_path/konno_matrices

    2. Extract features from gm_classifier directory. Because this can take a while, it is recommended to handle data by a yearly, or even montly, basis. Records that have already been processed will not be reprocessed, so interrupted operations will not result in a massive loss of progress. (example below):

      1. python gm_classifier/gm_classifier/scripts/extract_features.py --ko_matrices_dir base_dir_path/konno_matrices base_dir_path/gmc_record/features/YEAR base_dir_path/mseed/YEAR

    3. Run model prediction (example below):

      1. python gm_classifier/scripts/predict.py base_dir_path/gmc_record/features/YEAR base_dir_path/gmc_record/YEAR/results.csv

    4. Concatenate results from all YEAR directories under base_dir_path/gmc_record/ for a meta-results file.
  4. Filter out the best results from the meta file with gmc_filter.py. This will copy files from the mseed directory to an mseed_preferred directory, while will be used for GM IM.
  5. Generate the ground motion intensity measure data. This requires the IM_calculation Python package. It is recommended to install this package in a new Python environment via Conda or otherwise. Run the program using a shell wrapper script. See the example, im_runner.sh
  6. Finally, build the complete GM IM catalogue with by joining the ground motion classifier results with the ground motion intensity measures. An example for this is provided with the Python script gm_meta_builder.py.



  • No labels