The gap_heatmap.py program allows quantitative monitoring of the evolution of the coverage of the National Seismological Network of Colombia in the areas of greatest interest to the Group of Evaluation and Monitoring of Seismic Activity of the Colombian Geological Survey, by calculating the maximum azimuthal GAP average in these areas and the generation of heat and contour maps.
The main script (gap_heatmap.py) has 3 commands that allow to perform the following tasks:
- [grids] Compute azimuthal GAP over a grid (in Colombia) taking into account the distribution of the seismological network (in this particular case, the Colombian National Seismological Network)
- [heatmaps] With those grids created in the previous step, generates a heatmap and a contour plot over Colombian geographical map.
- [g-h] Perform the 2 above.
- python 3.11
For non anaconda python installations:
$ apt-get install libproj-dev proj-data proj-bin
$ apt-get install libgeos-dev
$ pip install cython
Clone or download this repository and then:
$ pip install -r requirements.txt
If you have an anaconda installation:
$ conda install -c scitools cartopy
$ pip install -r requirements.txt
The program has 3 commands that allow you to carry out the following specific tasks:
- grids: Calculate the azimuthal gap in a grid over Colombia taking into account the distribution of the stations that are entered through a .csv file with their coordinates or by default, with the stations of the National Seismological Network of Colombia for the years 2015-2020.
- heatmaps: With the grids calculated in the previous step, generate a heat map and a contour map on the geographic map of Colombia.
- g-h: Perform 1 and 2.
To get the command names, type:
$ python gap_heatmap.py --help
To get help about each command, type:
$ python gap_heatmap.py [command] --help
To get distance heatmap and gap isolines, open gap-distance-plot.py edit the parameters at the end of the file and run it.
$ python gap-distance-plot.py