lizadams
commented
3 years ago Will work on obtaining openstreet map as *.osm file and then load it into QGIS and save as a shapefile. curl --connect-timeout 10 --retry 10 --retry-delay 10 --retry-max-time 500 -H "Host: overpass-api.de" -H "Content-Type: text/xml" -d 'relation["adminlevel"="2"];(.;>;); out body;' http://overpass-api.de/api/interpreter -o countries-border.osm That seemed to provide a specific boundary. I am not sure how to select or obtain a larger area.
Is your feature request related to a problem? Please describe. The current world basemap in VERDI doesn't have the equivalent to county or finer administrative boundaries, so when you get a small 3km(?) domain over Botaga Columbia, there isn't a boundary visible on the world background map. This issue was revealed when Donna requested and obtained an example of a mercator projection from James West.
For the mercator projection for Bogota use the LANDMASK variable from the GRIDCRO2D file to compare with the background map. The only basemap that we have that displayed any boundaries was the river basemap for the Bogota Domain. The Columbia Mercator file is 9km resolution, and it shows the administrative boundary lines from the World Map.
/work/MOD3EVAL/jeast/data/eqm_inner The current world map in VERDI is: worldMapFileName = "Department_of_State_Valid_QGIS.shp"; // OK as Shapefile
Describe the solution you'd like Obtain county or finer administrative boundary files for additional regions of the world, if they are available.
Describe alternatives you've considered If we can obtain publically available openstreet maps as shapefiles, we may be able to convert these high resolution shapefiles from a polygon to a polyline feature for them to be used in VERDI. (this is what Michael Geigert did for the high resolution state map that he provided to us.)
Additional context I will also add the GRIDCRO2D to the github repository and name it GRIDCRO2D_Mercator_Bogota.nc