I am trying to plot a geographic area, the area comprised of multiple regions. To understand, I tested this by plotting below geographic area, it is easy to plot this country map using GMT as attached in figure 1.
What I want to achieve is shown in below reference figure 2.
Figure 2.
Queries are;
Whether it is possible to plot and show individual boundaries of each region (1-8) in figure 1 (which is plotted using GMT)? So each region is recognizable individually.
Whether it is possible to plot individual regions without other regions in figure 1 using GMT?
I have downloaded data which may contain administrative regions of this country, the following image shows which formats (.cpg, .csv, .shp …) of data listed in downloaded data.
I tried with (.shp) data to plot the regions using GMT, but not successful. Can you provide an idea which format of this listed data can be used and which commands can help to do this? Thanks.
The problem is that a shape file is the set of all those files that it is complaining of missing. You need to have them all sit near by the .shp file. But I doubt that file will reproduce the figure you shown in first post. It will very likely be similar to (this is using the GMT.jl Julia wrapper)
And painting the polygons is another story.
using GMT
japan = gadm("JPN", children=true);
viz(japan, proj=:guess, title="States of Japan")
Note that the ogr2ogr command in the screen capture does ‘-f OGR_GMT’. Isn’t ‘-f GMT’ the correct one? We discusssed this before (couple of years ago). I always use ‘-f GMT’. I wonder why that works then.
