#3060: Create location from command line using "r.in.gdal" and similar
-------------------------+-------------------------
Reporter: mankoff | Owner: grass-dev@…
Type: enhancement | Status: new
Priority: normal | Milestone: 7.0.5
Component: Startup | Version: unspecified
Keywords: | CPU: Unspecified
Platform: Unspecified |
-------------------------+-------------------------
I find that I am continually creating some new dummy location just so I
can then create the actual location I want using {{{r.in.gdal}}}. I could
maintain a global dummy location and use that, but it makes the code less
portable when I share it with others. Therefore, I often do this:
{{{
rm -fR grass
mkdir grass
grass70 -c EPSG:4326 ./grass/latlon
r.in.gdal -c input=/path/to/somefile.TIF location=foo output=tmp
exit
grass70 ./grass/foo/PERMANENT
}}}
I find it difficult to know what EPSG code I should use for somefile.TIF,
because it reports many (see below), therefore it would be nice to create
the location from outside of GRASS, but using GRASS. For example,
{{{grass70 -c input=/path/to/somefile importer=r.in.gdal}}}
{{{
Coordinate System is:
PROJCS["WGS 84 / UTM zone 33N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","32633"]]
Origin = (444892.500000000000000,8690407.500000000000000)
Pixel Size = (15.000000000000000,-15.000000000000000)
}}}
--
Ticket URL: <https://trac.osgeo.org/grass/ticket/3060>
GRASS GIS <https://grass.osgeo.org>