[GRASS-user] Import Shape Files

[Markus Neteler]

Hi Martin

On Wed, Jun 24, 2020 at 3:04 PM Martin Bittens <martin.bittens at gmx.de> wrote:
>
> Hello,
>
> I am trying to import shape-files from QGIS 3.12.3 to GRASS GIS 7.8.3).
> CRS  QGIS: EPSG 31983; CRS GRASS: EPSG 31983. EPSG 31983 is: SIRGAS
> 2000/UTM zone 23S.
>
> The import is working fine, the projection fits but when I check the
> meta data of the imported files I receive the message "Projection: UTM
> (zone invalid)". In order to narrow down that error I did the following
> in GRASS:
>
> 1.) I created 2 fresh locations/map sets: Loc1/Wells and Loc2/Wells
> using the same CRS (EPSG 31983);

I just tried:

grass78 -c epsg:31983 ~/grassdata/utm31983
Starting GRASS GIS...
Creating new GRASS GIS location <utm31983>...
...
Welcome to GRASS GIS 7.8.4dev (0798df3d0)
...

GRASS 7.8.4dev (utm31983):~/ec_cloud_expert > g.proj -w
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"]]


Does the `g.proj -w` output look like this for you?

All the other steps basically depend on this.

Best
Markus


> 2.) I did the import of a shape file from QGIS (EPSG 31983) to GRASS
> Loc1/Wells;
>
> 3.) Checking the meta data of the imported file results in "Projection:
> UTM (zone invalid)";
>
> 4.) I changed the location/map set from Loc1/Wells to Loc2/Wells;
>
> 5.) In Loc2/Wells I "transferred" the imported shape file from
> Loc1/Wells to Loc2/Wells using the command "v.proj";
>
> 6.) Checking the meta data now the message shows the correct projection
> "Projection: UTM (zone 23)".
>
>
> Are there any ideas why the shape file import in step 3 results in the
> error message "Projection: UTM (zone invalid)"?
>
>
> Thank you for your help.
>
>
> Regards
>
>
> Martin
>