Faults

Q&A
1

how can i plot the faults of sicily in a very simple way?
is there a way to enter a url and gmt refers to the site? I have plotted the stations taking lat and long and I would like to do such a thing if there are no other faster ways with faults …
thank you in advance

2

It is possible to use an URL.
Do you know any site where the data of the faults are available?

1 Like
3

Oh, thanks for the answer. yes. the site is this https://diss.ingv.it/ithdiss/#
now i can show you my code GMT where i want to put the sicily’s faults.
gmt begin italy pdf,png

gmt set IO_COL_SEPARATOR=tab

gmt info C:\programs\Lat_long_INGV.txt

gmt coast -R11.00000/16.62244/35.91997/39.36489 -Wthin -B -JM5i -Gbeige -Slightblue -I1

gmt grdimage @earth_relief_01s -R11.00000/16.62244/35.91997/39.36489 -I+d

gmt plot C:\programs\Lat_long_INGV.txt -St0.25c -Gblue

gmt end show

Im now on GMT. So for me it is usefull use the community.

4

I check the link. It suspect that you could download the data but I don’t find where.

5

here we have the possibility to dowload…i think.
6

Hi Anna, you don’t want to do this because it would use a
(16.62244 - 11.0) * 3600 x (39.36489 - 35.91997) *3600
=
20240 x 12404
grid size, which is HUGE. Better to pick a lower grid resolution here (5 or 6 minutes should be OK).

7

You could try the shp file, unzip it and then use gmt plot with the shp that you want to plot.
https://diss.ingv.it/images/downloads/DISS/DISS_3.3.0_shp.zip

8

Hi gioacchino, I don’t understand your disappointment about this piece of code.

what should I do ?

9

about this GMT give these errors

grdblend [ERROR]: Failed to remove C:\Users\Anna\AppData\Local\Temp/grdblend_resampled_2296_16.nc! [remove error: Permission denied]
grdblend [ERROR]: Failed to delete file C:\Users\Anna\AppData\Local\Temp/grdblend_resampled_2296_16.nc

Why?

10

thank you for solution.

I plotted all file from this site.

Thanks @Esteban82

1 Like
11

You are using a grid with a resolution of 1 second of arc (@earth_relief_01s). You could try with a 6-minutes grid (@earth_relief_06m) given the region of your map and size of the map.

12

That is an issue that was fixed after 6.3.0. And regarding the the grid resolution 5 or 6 minutes is probably too coarse. 2 minutes is likely better.
Latest GMT version lets you not specify the resolution and it will guess a good one (which turns out to be 3s)

gmt grdimage @earth_relief -R11.0/16.62244/35.91997/39.36489 -I+d -png lixo
grdblend [NOTICE]: Remote data courtesy of GMT data server oceania [http://oceania.generic-mapping-tools.org]

grdblend [NOTICE]: Earth Relief at 3x3 arc seconds tiles provided by SRTMGL3 (land only) [NASA/USGS].
grdblend [NOTICE]:   -> Download 1x1 degree grid tile (earth_relief_03s_g): N35E011
grdblend [NOTICE]:   -> Download 1x1 degree grid tile (earth_relief_03s_g): N35E012
...
13

Thanks you.
But if i use 6m or 2m the resolution is very bad.

Screenshot (739)
Screenshot (739)1920×1080 314 KB

So it is better with 1s or 2s? Why?

14

can you help me with fault shp file for India and neighbouring countries

15

Hello! What is your problem?

16

I want to plot active fault in India OVER A BASEMAP using GMT

17

Ok, do you have a file with the faults?

18

Figured it out.
Faults can be downloaded from (GitHub - GEMScienceTools/gem-global-active-faults: A global homogenised database of active faults maintained by the GEM Foundation)

faults = gpd.read_file(“gem_active_faults.shp”)
fault_lines = faults.geometry
fault_names = faults[“name”]
for geom, name in zip(fault_lines, fault_names):
x, y = geom.xy # Extract x and y coordinates from geometry
fig.plot(x=x, y=y, pen=“2p,red”)