Shape file into GMT

KRL · July 18, 2020, 7:03pm

Dear all,

I was wondering if you could help, I have several shape files that I would like to have in my GMT maps. Is it possible to have them read an ploted ? I have try exporting them as xy.txt files and plot them with psxy and plot commands but I am getting not a nice image.

gmt begin audemardf_fallas pdf,png 

gmt coast -JM30c -R-75/-55/0/15 -Wfaint,blue -Gwhite -Sazure -Baf -N1/thick,black 
gmt psxy audemardf_fallas_c1.txt -JM30c -R-75/-55/0/15 -Baf -Wfaint,red,dash


gmt end show

pwessel · July 18, 2020, 7:14pm

Have you tried giving plot the name of the shape file (*.shp)?

KRL · July 18, 2020, 7:32pm

I have just tried, is not plotting the Faults on the map
below the modification on the script and the result in the terminal

gmt begin audemardf_fallas pdf,png 

gmt coast -JM30c -R-75/-55/0/15 -Wfaint,blue -Gwhite -Sazure -Baf -N1/thick,black 
gmt plot Fallas_Cuaternarias_Venezuela_lineas.shp -JM30c -R-75/-55/0/15 -Baf -Wfaint,red,dash


gmt end show

./audemardf_fallas.sh: line 1: !#usr/bin/env: No such file or directory
ERROR 4: Unable to open Fallas_Cuaternarias_Venezuela_lineas.shx or Fallas_Cuaternarias_Venezuela_lineas.SHX. Set SHAPE_RESTORE_SHX config option to YES to restore or create it.
FAILURE:
Unable to open datasource `Fallas_Cuaternarias_Venezuela_lineas.shp' with the following drivers.
  -> `PCIDSK'
  -> `netCDF'
  -> `PDS4'
  -> `JP2OpenJPEG'
  -> `PDF'
  -> `MBTiles'
  -> `EEDA'
  -> `ESRI Shapefile'
  -> `MapInfo File'
  -> `UK .NTF'
  -> `OGR_SDTS'
  -> `S57'
  -> `DGN'
  -> `OGR_VRT'
  -> `REC'
  -> `Memory'
  -> `BNA'
  -> `CSV'
  -> `GML'
  -> `GPX'
  -> `KML'
  -> `GeoJSON'
  -> `GeoJSONSeq'
  -> `ESRIJSON'
  -> `TopoJSON'
  -> `OGR_GMT'
  -> `GPKG'
  -> `SQLite'
  -> `WAsP'
  -> `OpenFileGDB'
  -> `XPlane'
  -> `DXF'
  -> `CAD'
  -> `Geoconcept'
  -> `GeoRSS'
  -> `GPSTrackMaker'
  -> `VFK'
  -> `PGDUMP'
  -> `OSM'
  -> `GPSBabel'
  -> `SUA'
  -> `OpenAir'
  -> `OGR_PDS'
  -> `WFS'
  -> `WFS3'
  -> `HTF'
  -> `AeronavFAA'
  -> `EDIGEO'
  -> `GFT'
  -> `SVG'
  -> `CouchDB'
  -> `Cloudant'
  -> `Idrisi'
  -> `ARCGEN'
  -> `SEGUKOOA'
  -> `SEGY'
  -> `ODS'
  -> `XLSX'
  -> `ElasticSearch'
  -> `Carto'
  -> `AmigoCloud'
  -> `SXF'
  -> `Selafin'
  -> `JML'
  -> `PLSCENES'
  -> `CSW'
  -> `VDV'
  -> `MVT'
  -> `TIGER'
  -> `AVCBin'
  -> `AVCE00'
  -> `NGW'
  -> `HTTP'
plot [ERROR]: System call [ogr2ogr -mapFieldType Integer64=Integer -f "OGR_GMT" /var/folders/q4/f3vkf3lj7cg8fx8tp866cmt00000gn/T/gmt_ogr_3248.gmt Fallas_Cuaternarias_Venezuela_lineas.shp] FAILED with error 256.
plot [ERROR]: Cannot open file Fallas_Cuaternarias_Venezuela_lineas.shp
bash-3.2$

KRL · July 18, 2020, 7:45pm

trying with a diferent shape file

./audemardf_fallas.sh: line 1: !#usr/bin/env: No such file or directory
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Unable to write feature 434 from layer Fallas_Cuaternarias_Triangulo_2.
ERROR 1: Terminating translation prematurely after failed
translation of layer Fallas_Cuaternarias_Triangulo_2 (use -skipfailures to skip errors)
plot [ERROR]: System call [ogr2ogr -mapFieldType Integer64=Integer -f "OGR_GMT" /var/folders/q4/f3vkf3lj7cg8fx8tp866cmt00000gn/T/gmt_ogr_3440.gmt Fallas_Cuaternarias_Triangulo_2.shp] FAILED with error 256.
plot [ERROR]: Cannot open file Fallas_Cuaternarias_Triangulo_2.shp
bash-3.2$ 
bash-3.2$ ```

![Screenshot 2020-07-18 15.44.58|690x421](upload://cPlUkcT6cKA3fj8sN4Kc4QqUv9l.png)

Joaquim · July 18, 2020, 7:46pm

The .shx is the index file and should not be mandatory to read the shape file. If you can make that shapefile available we can look (see if my gdal is more permissive?). But only within two couple of hours.

KRL · July 18, 2020, 7:48pm

Okey, please may I have an e-mail to send the folder

pwessel · July 18, 2020, 7:48pm

Also fix your shebang line (1st line in the script). Should be #!/usr/bin/env bash

Joaquim · July 18, 2020, 7:49pm

If file is not too big you can zip it and post it here (drag-n-drop). Otherwise post a link.

KRL · July 18, 2020, 7:53pm

Fallas_Cuaternarias.zip (774.2 KB)

KRL · July 18, 2020, 7:57pm

these are the shapefiles I am trying to get into GMT, still, geological units are missing but I can begin with faults

Uploading: Fallas_Cuaternarias.zip…

Esteban82 · July 18, 2020, 9:13pm

In the link I explained how to convert shp to txt for GMT. In your txt file, each line (or polygon) should be separated with >.

KRL · July 18, 2020, 10:47pm

My QGIS version don’t have GMT format to export which QGIS is this version is it MacOsX?
I try usin ESRI Shapefile format , and the plot command… and is not working

gery · July 19, 2020, 12:50am

Hola K.,

This works:

for shp in *.shp
do
nombre=basename $shp .shp
ogr2ogr -f “GMT” ${nombre}.gmt $shp -skipfailures
done

It generates these ones but just ignore them:

ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
ERROR 1: Features without geometry not supported by GMT writer.
Warning 1: The output driver does not seem to natively support Integer64 type for field ID. Converting it to Real instead. -mapFieldType can be used to control field type conversion.
Warning 1: The output driver does not seem to natively support Integer64 type for field ID. Converting it to Real instead. -mapFieldType can be used to control field type conversion.
ERROR 1: Features without geometry not supported by GMT writer.

Check the enclosed image, plotted with GMT (6!!!) and QGIS.

Hope this helps,

Gery

KRL · July 19, 2020, 1:09am

Muchas gracias! lo voy a intentar!, la imagen en GMT tiene algunos detalles sigue conectando algunos polígonos que no se encuentran conectados entre ellos realmente.

Joaquim · July 19, 2020, 1:36am

Not sure why fancy is off, but another (so simpler) alternative. With GMT Julia (and mixed mode, sorry the other I don’t remember from hart)

D = gmtread("c:/v/Fallas_Cuaternarias/Fallas_Cuaternarias_Venezuela_lineas.shp");
coast(region=(-75,-55,0,15), proj=:Mercator, figsize=17, W="faint,blue", G="white", S="azure", N="1/thick,black")
plot!(D, W="faint,red,dash", fmt=:png, show=true)

gery · July 19, 2020, 1:37am

Sip, tienes razón, well, some features are not connected because it’s a quick plot, start playing with that piece of code to get the gmt files first, once you have them you need to plot them separately. Please keep us updated about your progress.

Gusto ayudarte!

KRL · July 19, 2020, 2:01am

I will try this way as well it looks so Beautiful! this is exacly how it should look like

KRL · July 19, 2020, 5:57pm

Dear Joaquim I am at the moment installing GMT Julia , which I am consulting the installation process from this page http://gmt.soest.hawaii.edu/doc/latest/julia_wrapper.html#installing

any particular reason why shows up this ERROR

Hope not to take too much of your time and thank you

KRL · July 19, 2020, 6:10pm

Joaquim · July 19, 2020, 8:05pm

In Julia, like in Python, before using a package you must load it first. But just do like is explain in the readme page, Type

] add GMT

and then before using it

using GMT