forked from alirezaniki/OBSTransformer
-
Notifications
You must be signed in to change notification settings - Fork 2
/
build_json.sh
39 lines (32 loc) · 1.65 KB
/
build_json.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
#!/bin/bash
# May 18 2023
# Alrieza Niksejel ([email protected])
# provides the required station_list.json file
# station info file needs to follow the below format:
# longitude(deg) latitude(deg) network station channel depth(km)
# Full path to the dataset
data_dir="/media/dalquake/T7/3k_OBS_validation_dataset/Github/dataset"
# Full path to the station info file
station_info="/media/dalquake/T7/3k_OBS_validation_dataset/Github/stations.dat"
pickhome=`pwd`
[[ ! -d $data_dir ]] || [[ ! -f $station_info ]] && echo "requirement(s) missing" && exit
[[ -d json ]] && rm -r json/* 2> /dev/null || mkdir json
cd $data_dir
n=`grep -cE '^[1-9]|^-' $station_info`
i=1
for sta in `find ./ -type f -name "*mseed" | cut -f 3 -d / | cut -f 2 -d . | sort | uniq`; do
net=`grep "$sta" $station_info | awk '{print $3}' | head -1`
cha=`grep "$sta" $station_info | awk '{print $5}' | head -1`
lat=`grep "$sta" $station_info | awk '{print $2}' | head -1`
lon=`grep "$sta" $station_info | awk '{print $1}' | head -1`
dep=`grep "$sta" $station_info | awk '{print $6}' | head -1`
if [[ $i = 1 ]]; then
echo "{\"$sta\": {\"network\": \"$net\", \"channels\": [\""$cha"1\", \""$cha"2\", \""$cha"Z\"], \"coords\": [$lat, $lon, $dep]}," >> $pickhome/json/station_list.json
elif [[ $i = $n ]]; then
echo "\"$sta\": {\"network\": \"$net\", \"channels\": [\""$cha"1\", \""$cha"2\", \""$cha"Z\"], \"coords\": [$lat, $lon, $dep]}}" >> $pickhome/json/station_list.json
else
echo "\"$sta\": {\"network\": \"$net\", \"channels\": [\""$cha"1\", \""$cha"2\", \""$cha"Z\"], \"coords\": [$lat, $lon, $dep]}," >> $pickhome/json/station_list.json
fi
echo "done with station $sta"
i=$((i+1))
done