Commit 12a2e47a authored by Franco Mele's avatar Franco Mele
Browse files

per fare una mappa con qualità delle localizzazioni

parent 4ee89b5b
*.swp
*.history
*.DS_Store
#
# GMT-SYSTEM 4.1.3 Defaults file
#
#-------- Plot Media Parameters -------------
PAGE_COLOR = 255/255/255
PAGE_ORIENTATION = landscape
PAPER_MEDIA = a0
#-------- Basemap Annotation Parameters ------
ANNOT_MIN_ANGLE = 20
ANNOT_MIN_SPACING = 0
ANNOT_FONT_PRIMARY = Courier-Bold
# ANNOT_FONT_SIZE = 10p
ANNOT_FONT_SIZE = 8p
ANNOT_OFFSET_PRIMARY = 0.075i
ANNOT_FONT_SECONDARY = Helvetica
# ANNOT_FONT_SIZE_SECONDARY = 16p
ANNOT_FONT_SIZE_SECONDARY = 14p
ANNOT_OFFSET_SECONDARY = 0.075i
DEGREE_SYMBOL = ring
HEADER_FONT = Helvetica
# HEADER_FONT_SIZE = 18p
HEADER_FONT_SIZE = 16p
HEADER_OFFSET = 0.1875i
LABEL_FONT = Helvetica
# LABEL_FONT_SIZE = 24p
LABEL_FONT_SIZE = 20p
LABEL_OFFSET = 0.1125i
OBLIQUE_ANNOTATION = 1
PLOT_CLOCK_FORMAT = hh:mm:ss
PLOT_DATE_FORMAT = yyyy-mm-dd
PLOT_DEGREE_FORMAT = +ddd:mm:ss
Y_AXIS_TYPE = hor_text
#-------- Basemap Layout Parameters ---------
BASEMAP_AXES = WESN
BASEMAP_FRAME_RGB = 0/0/0
BASEMAP_TYPE = plain
FRAME_PEN = 1.25p
FRAME_WIDTH = 0.075i
GRID_CROSS_SIZE_PRIMARY = 0i
GRID_CROSS_SIZE_SECONDARY = 0i
GRID_PEN_PRIMARY = 0.25p
GRID_PEN_SECONDARY = 0.5p
MAP_SCALE_HEIGHT = 0.075i
TICK_LENGTH = 0.075i
POLAR_CAP = 85/90
TICK_PEN = 0.5p
X_AXIS_LENGTH = 9i
Y_AXIS_LENGTH = 6i
X_ORIGIN = 1i
Y_ORIGIN = 1i
UNIX_TIME = FALSE
UNIX_TIME_POS = -0.75i/-0.75i
#-------- Color System Parameters -----------
COLOR_BACKGROUND = 0/0/0
COLOR_FOREGROUND = 255/255/255
COLOR_NAN = 128/128/128
COLOR_IMAGE = adobe
COLOR_MODEL = rgb
HSV_MIN_SATURATION = 1
HSV_MAX_SATURATION = 0.1
HSV_MIN_VALUE = 0.3
HSV_MAX_VALUE = 1
#-------- PostScript Parameters -------------
CHAR_ENCODING = Standard
DOTS_PR_INCH = 300
N_COPIES = 1
PS_COLOR = rgb
PS_IMAGE_COMPRESS = none
PS_IMAGE_FORMAT = ascii
PS_LINE_CAP = butt
PS_LINE_JOIN = miter
PS_MITER_LIMIT = 0
PS_VERBOSE = FALSE
GLOBAL_X_SCALE = 1
GLOBAL_Y_SCALE = 1
#-------- I/O Format Parameters -------------
D_FORMAT = %lg
FIELD_DELIMITER = tab
GRIDFILE_SHORTHAND = FALSE
GRID_FORMAT = nf
INPUT_CLOCK_FORMAT = hh:mm:ss
INPUT_DATE_FORMAT = yyyy-mm-dd
IO_HEADER = FALSE
N_HEADER_RECS = 1
OUTPUT_CLOCK_FORMAT = hh:mm:ss
OUTPUT_DATE_FORMAT = yyyy-mm-dd
OUTPUT_DEGREE_FORMAT = +D
XY_TOGGLE = FALSE
#-------- Projection Parameters -------------
ELLIPSOID = WGS-84
MAP_SCALE_FACTOR = default
MEASURE_UNIT = inch
#-------- Calendar/Time Parameters ----------
TIME_FORMAT_PRIMARY = full
TIME_FORMAT_SECONDARY = full
TIME_EPOCH = 2000-01-01T00:00:00
TIME_IS_INTERVAL = OFF
TIME_INTERVAL_FRACTION = 0.5
TIME_LANGUAGE = us
TIME_SYSTEM = j2000
TIME_UNIT = d
TIME_WEEK_START = Sunday
Y2K_OFFSET_YEAR = 1950
#-------- Miscellaneous Parameters ----------
HISTORY = TRUE
INTERPOLANT = akima
LINE_STEP = 0.01i
VECTOR_SHAPE = 0
VERBOSE = FALSE
This diff is collapsed.
IN DATA 3-nov-2020 per la REVISIONE:
chiedono di spostare le 3 legende nell'angolo in alto a destra.
Modifico make_qual2 in make_qual3
Per il bollettino del 2015 uso:
$ ./make_qual3 20150101_20151231_quality.txt file_2015.ps RdYlBu_10.cpt
# GMT palette RdBu_7.cpt
#
# This product includes color specifications and designs
# developed by Cynthia Brewer (http://colorbrewer.org/).
#
# Converted to the cpt format by J.J.Green
# Diverging palette with 7 colours
#
# COLOR_MODEL = RGB
-3.50 178 24 43 -2.50 178 24 43
-2.50 239 138 98 -1.50 239 138 98
-1.50 253 219 199 -0.50 253 219 199
-0.50 247 247 247 0.50 247 247 247
0.50 209 229 240 1.50 209 229 240
1.50 103 169 207 2.50 103 169 207
2.50 33 102 172 3.50 33 102 172
# GMT palette RdGy_10.cpt
#
# This product includes color specifications and designs
# developed by Cynthia Brewer (http://colorbrewer.org/).
#
# Converted to the cpt format by J.J.Green
# Diverging palette with 10 colours
#
# COLOR_MODEL = RGB
-5.00 103 000 031 -4.00 103 000 031
-4.00 178 024 043 -3.00 178 024 043
-3.00 214 096 077 -2.00 214 096 077
-2.00 244 165 130 -1.00 244 165 130
-1.00 253 219 199 0.00 253 219 199
0.00 224 224 224 1.00 224 224 224
1.00 186 186 186 2.00 186 186 186
2.00 135 135 135 3.00 135 135 135
3.00 077 077 077 4.00 077 077 077
4.00 026 026 026 5.00 026 026 026
# GMT palette RdGy_11.cpt
#
# This product includes color specifications and designs
# developed by Cynthia Brewer (http://colorbrewer.org/).
#
# Converted to the cpt format by J.J.Green
# Diverging palette with 11 colours
#
# COLOR_MODEL = RGB
-5.50 103 000 031 -4.50 103 000 031
-4.50 178 024 043 -3.50 178 024 043
-3.50 214 096 077 -2.50 214 096 077
-2.50 244 165 130 -1.50 244 165 130
-1.50 253 219 199 -0.50 253 219 199
-0.50 255 255 255 0.50 255 255 255
0.50 224 224 224 1.50 224 224 224
1.50 186 186 186 2.50 186 186 186
2.50 135 135 135 3.50 135 135 135
3.50 077 077 077 4.50 077 077 077
4.50 026 026 026 5.50 026 026 026
# GMT palette RdYlBu_10.cpt
#
# This product includes color specifications and designs
# developed by Cynthia Brewer (http://colorbrewer.org/).
#
# Converted to the cpt format by J.J.Green
# Diverging palette with 10 colours
#
# COLOR_MODEL = RGB
-5.00 165 000 038 -4.00 165 000 038
-4.00 215 048 039 -3.00 215 048 039
-3.00 244 109 067 -2.00 244 109 067
-2.00 253 174 097 -1.00 253 174 097
-1.00 254 224 144 0.00 254 224 144
0.00 224 243 248 1.00 224 243 248
1.00 171 217 233 2.00 171 217 233
2.00 116 173 209 3.00 116 173 209
3.00 069 117 180 4.00 069 117 180
4.00 049 054 149 5.00 049 054 149
# GMT palette RdYlBu_11.cpt
#
# This product includes color specifications and designs
# developed by Cynthia Brewer (http://colorbrewer.org/).
#
# Converted to the cpt format by J.J.Green
# Diverging palette with 11 colours
#
# COLOR_MODEL = RGB
-5.50 165 000 038 -4.50 165 000 038
-4.50 215 048 039 -3.50 215 048 039
-3.50 244 109 067 -2.50 244 109 067
-2.50 253 174 097 -1.50 253 174 097
-1.50 254 224 144 -0.50 254 224 144
-0.50 255 255 191 0.50 255 255 191
0.50 224 243 248 1.50 224 243 248
1.50 171 217 233 2.50 171 217 233
2.50 116 173 209 3.50 116 173 209
3.50 069 117 180 4.50 069 117 180
4.50 049 054 149 5.50 049 054 149
# GMT palette PiYG_10.cpt
#
# This product includes color specifications and designs
# developed by Cynthia Brewer (http://colorbrewer.org/).
#
# Converted to the cpt format by J.J.Green
# Diverging palette with 10 colours
#
# COLOR_MODEL = RGB
#
# 2008-01-22 Created By Franco Mele from PiYG_10.cpt
#
-5.00 142 001 082 -4.00 142 001 082
-4.00 197 027 125 -3.00 197 027 125
-3.00 222 119 174 -2.00 222 119 174
-2.00 241 182 218 -1.00 241 182 218
-1.00 253 224 239 0.00 253 224 239
0.00 230 245 208 1.00 230 245 208
1.00 184 225 134 2.00 184 225 134
2.00 127 188 065 3.00 127 188 065
3.00 077 146 033 4.00 077 146 033
4.00 039 100 025 5.00 039 100 025
19.8 48.0 0.9 c
19.8 47.3 0.625 c
19.8 46.8 0.4 c
19.8 46.4 0.223 c
19.8 46.0 0.1 c
16.9 48.1 0.4 t
16.9 47.65 0.4 s
16.9 47.20 0.4 n
16.9 46.75 0.4 h
16.9 46.30 0.4 g
16.9 45.85 0.4 a
#x y texti
19.8 48.7 Mag
20.5 48.0 6
20.5 47.3 5
20.5 46.8 4
20.5 46.4 3
20.5 46.0 2
17.7 48.7 Depth (km)
17.8 48.1 ( 0; 5]
17.8 47.65 ( 5; 10]
17.9 47.20 (10; 30]
17.9 46.75 (30; 60]
17.9 46.30 (60;100]
17.9 45.85 >100
21.8 48.3 best
21.8 45.1 worst
21.8 48.7 Quality
#!/bin/sh
MIN_ARGS=3
if (( $# != $MIN_ARGS ));
then
echo "Usage:"
echo "$0 GMT_DATA_FILE GMT_PS_FILE GMT_CPT_FILE"
exit
fi
# Data about psxy output
GMT_CPT_FILE=seis.cpt
GMT_CPT_FILE_MAP1=seis.cpt
GMT_DATA_FILE=$1
GMT_MAP=$2
GMT_CPT_FILE=$3
# GMT_CPT_FILE=GMT_sies_10.cpt
GMT_CPT_FILE_MAP1=$3
# GMT_CPT_FILE_MAP2=seis_map2.cpt
# Legend file
GMT_LEGEND_CIRCLE=legend_circle_quality
GMT_LEGEND_TEXT=legend_quality
GMT_TITLE=' '
# Data about original (Italy) map
ITALY_MAP=italy.ps
#ITALY_GRD=italy.grd
#ZONE='6.4565/18.948167/36.670833/47.329167'
ZONE='4/23/34/49'
BAR_POSITION_1='15.75c/15.0c/3.5c/.35c'
#BAR_POSITION_2='.6c/4.9c/4c/.5c'
# posizione X e Y del foglio in inch.
POSX=3.0
POSY=1.5
# SCALA dell'immagine (dimensione in inch per ogni 100 km)
SCAL=0.9
# italy.ps (italy map) create it
PALETTE=palette.cpt
#RESAMPLED=resampled.grd
## crea PALETTE topografica
makecpt -Crelief -I -Z -T-100/2000/100 > $PALETTE
psbasemap -R$ZONE -Jm$SCAL -X$POSX -Y$POSY -Ba2/a2:."$GMT_TITLE": -K > $ITALY_MAP
#/chroot/gmt/bin//grdsample $ITALY_GRD -R$ZONE -I0.01/0.01 -G$RESAMPLED -V
#/chroot/gmt/bin//grdimage $RESAMPLED -V -Y$POSY -X$POSX -R$ZONE -C$PALETTE -P -Jm$SCAL -K > $ITALY_MAP
pscoast -Jm -R$ZONE -Df -S100/200/255 -N3 -I/0.25p/0/200/255 -L7/34.8/34.8/100 -V -O -K >> $ITALY_MAP
psscale -D$BAR_POSITION_1 -O -C$GMT_CPT_FILE_MAP1 -B1::/:" ": -K >> $ITALY_MAP
# psscale -D$BAR_POSITION_2 -O -C$GMT_CPT_FILE_MAP2 -Ba200f100::/:: -K >> $ITALY_MAP
psxy -R -Jm$SCAL -O $GMT_LEGEND_CIRCLE -S -W0.25p -K >> $ITALY_MAP
pstext -F+f10p $GMT_LEGEND_TEXT -O -R -Jm$SCAL -K >> $ITALY_MAP
# Deleting tmp files
#rm $RESAMPLED $PALETTE
# Making a work copy of italian map...
cp $ITALY_MAP $GMT_MAP
# ...adding earthquakes
psxy -R$ZONE -Jm$SCAL -O -C$GMT_CPT_FILE $GMT_DATA_FILE -S -W0.25p >> $GMT_MAP
# and convert the map in a web accepted format
# /usr/X11R6/bin//convert -trim -rotate 90 $GMT_MAP `basename $GMT_MAP ps`jpg
-100 white 0 253.48/250.43/225.29
0 253.48/250.43/225.29 100 252.48/245.86/195.57
100 252.48/245.86/195.57 200 252/241.57/169.86
200 252/241.57/169.86 300 251.9/237.81/150.76
300 251.9/237.81/150.76 400 250.38/234.76/130.95
400 250.38/234.76/130.95 500 250/231.71/111.14
500 250/231.71/111.14 600 232.67/212.67/93.333
600 232.67/212.67/93.333 700 193.05/173.05/78.095
700 193.05/173.05/78.095 800 153.43/133.43/62.857
800 153.43/133.43/62.857 900 108.1/104.76/50
900 108.1/104.76/50 1000 235.14/253.1/238.24
1000 235.14/253.1/238.24 1100 205.43/249.29/204.86
1100 205.43/249.29/204.86 1200 175.71/245.48/172.1
1200 175.71/245.48/172.1 1300 114.67/213/178.67
1300 114.67/213/178.67 1400 49.143/176.86/190.86
1400 49.143/176.86/190.86 1500 0/136.67/185.71
1500 0/136.67/185.71 1600 0/83.333/128.57
1600 0/83.333/128.57 1700 0/30/71.429
1700 0/30/71.429 1800 0/7.619/38.095
1800 0/7.619/38.095 1900 0/3.8095/19.048
1900 0/3.8095/19.048 2000 black
B white
F black
N white
The first quality rating is based on errors and goodness-of-fit:
Score
1.5 A. RMS < 0.45 sec and ERH <2.0 km and ERZ <4.0 km
0.5 B. RMS < 0.90 sec and ERH <5.0 km and ERZ <10.0 km
-0.5 C. RMS < 1.50 sec and ERH <10.0 km
-1.5 D. Worse than above
The second quality rating is based on station geometry:
Score
3.0 A. NWR > 6 and MAXGAP < 90 and either DMIN < DEPTH or DMIN <10.0
1.0 B. NWR > 6 and MAXGAP < 135 and either DMIN < 2*DEPTH or DMIN <20
-1.0 C. NWR > 6 and MAXGAP < 180 and DMIN <100
-3.0 D. Worse than above
\ No newline at end of file
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment