[GRASS-user] Import Shape Files

[Martin Bittens]

Hello Markus,

first of all thank very much for your efforts.

Unfortunately, the error message "Projection UTM (zone invalid) did
appear anyhow although I replaced the *.prj file by that one you sent. 
But now I figured out the following (after importing the shape file and
making active):

- using in the layers display the edit button (left click on it and then
selecting "Metadata"), the result is "Projection UTM (zone invalid)" as
I reported;

- selecting in the layers display from the vector menu "Reports and
statistics/Basic vector metadata" and then clicking on "Run", the result
is again: "Projection UTM (zone invalid)";

- selecting in the layers display from the vector menu "Reports and
statistics/Basic vector metadata" and additionally activating in the
"Print" tab the option (e) - "Print extended metadata info in shell
script style", the result is now "Projection=UTM, zone=-23".


I don't know why this happened but obviously GRASS provides different
metadata for the same shape-file depending on the command you are using.


Regards


Martin





Am 6/26/2020 um 11:56 AM schrieb Markus Neteler:
> Hello Martin,
>
> (again the list was not in copy)
>
>
> On Fri, Jun 26, 2020 at 1:19 PM Martin Bittens <martin.bittens at gmx.de> wrote:
>> Hello Markus,
>>
>> attached you can find a *.prj file created by QGIS.
> You may try to replace it with the one attached, generated by GRASS GIS.
>
> Best,
> Markus
>
>> Thank you very much for your help.
>>
>> Martin
>>
>>
>>
>> Am 6/25/2020 um 4:44 PM schrieb 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
>>
>> --
>> This email has been checked for viruses by Avast antivirus software.
>> https://www.avast.com/antivirus
>
>

--
This email has been checked for viruses by Avast antivirus software.
https://www.avast.com/antivirus