[GRASS-user] Import Shape Files

[Markus Neteler]

(please keep the list in CC)

On Thu, Jun 25, 2020 at 1:57 PM Martin Bittens <martin.bittens at gmx.de> wrote:
>
> Hello Markus,
>
> thank you for your reply.
>
> Usually I am creating locations and map sets with the GRASS start-up
> screen. After opening my 31983 location/mapset the command  "g.proj -w"
> results in:
> PROJCS["SIRGAS 2000 / UTM zone 23S",
>      GEOGCS["SIRGAS 2000",
> DATUM["Sistema_de_Referencia_Geocentrico_para_las_AmericaS_2000",
>              SPHEROID["GRS 1980",6378137,298.257222101,
>                  AUTHORITY["EPSG","7019"]],
>              TOWGS84[0,0,0,0,0,0,0],
>              AUTHORITY["EPSG","6674"]],
>          PRIMEM["Greenwich",0,
>              AUTHORITY["EPSG","8901"]],
>          UNIT["degree",0.0174532925199433,
>              AUTHORITY["EPSG","9122"]],
>          AUTHORITY["EPSG","4674"]],
>      PROJECTION["Transverse_Mercator"],
>      PARAMETER["latitude_of_origin",0],
>      PARAMETER["central_meridian",-45],
>      PARAMETER["scale_factor",0.9996],
>      PARAMETER["false_easting",500000],
>      PARAMETER["false_northing",10000000],
>      UNIT["metre",1,
>          AUTHORITY["EPSG","9001"]],
>      AXIS["Easting",EAST],
>      AXIS["Northing",NORTH],
>      AUTHORITY["EPSG","31983"]]
>
> This is almost the same result you sent. Just line appears in addition:
> TOWGS84[0,0,0,0,0,0,0].
>
>
> Using the command "grass78 -c epsg:31983 ~/grassdata/utm31983" like you
> did "g.proj -w" results in:
>
> PROJCS["SIRGAS 2000 / UTM zone 23S",
>      GEOGCS["SIRGAS 2000",
> DATUM["Sistema_de_Referencia_Geocentrico_para_las_AmericaS_2000",
>              SPHEROID["GRS 1980",6378137,298.257222101,
>                  AUTHORITY["EPSG","7019"]],
>              AUTHORITY["EPSG","6674"]],
>          PRIMEM["Greenwich",0,
>              AUTHORITY["EPSG","8901"]],
>          UNIT["degree",0.0174532925199433,
>              AUTHORITY["EPSG","9122"]],
>          AUTHORITY["EPSG","4674"]],
>      PROJECTION["Transverse_Mercator"],
>      PARAMETER["latitude_of_origin",0],
>      PARAMETER["central_meridian",-45],
>      PARAMETER["scale_factor",0.9996],
>      PARAMETER["false_easting",500000],
>      PARAMETER["false_northing",10000000],
>      UNIT["metre",1,
>          AUTHORITY["EPSG","9001"]],
>      AXIS["Easting",EAST],
>      AXIS["Northing",NORTH],
>      AUTHORITY["EPSG","31983"]]
>
> This looks exactly like the result you sent.

Ok, that looks fine to me.

> But anyhow, it does not
> matter how I am creating locations/mapsets imported shape files are
> showing in the meta data results "Projection UTM (zone invalid).

Then maybe the .proj file of te shape file(s) is incomplete?

Can you make one file available (even off-list)?

Best
Markus