<div dir="ltr">Hi Martin,<br><div><br></div><div>(back to cc of list)<br></div><div><br></div><div>I just got a chance to test and also got the "invalid zone" message issue:</div><br><a href="http://v.info">v.info</a> LT_AOI_3<br> <span style="font-family:monospace">+----------------------------------------------------------------------------+<br> | Name:            LT_AOI_3                                                  |<br> | Mapset:          PERMANENT                                                 |<br> | Location:        sirgas_2000                                               |<br> | Database:        /home/mneteler/grassdata                                  |<br> | Title:                                                                     |<br> | Map scale:       1:1                                                       |<br> | Name of creator: mneteler                                                  |<br> | Organization:                                                              |<br> | Source date:     Tue Aug 11 11:44:55 2020                                  |<br> | Timestamp (first layer): none                                              |<br> |----------------------------------------------------------------------------|<br> | Map format:      native                                                    |<br> |----------------------------------------------------------------------------|<br> |   Type of map: vector (level: 2)                                           |<br> |                                                                            |<br> |   Number of points:       5               Number of centroids:  0          |<br>...<br> |   Projection: UTM (zone invalid)                                           |<br></span><div><span style="font-family:monospace">...</span></div><div><br></div><div>This looks ok:</div><div><br></div><div><span style="font-family:monospace">g.region vector=LT_AOI_3 -p<br>projection: 1 (UTM)<br>zone:       -23<br>datum:      sirgas2000<br>ellipsoid:  grs80<br>north:      7397570.18085686<br>south:      7397525.93235069<br>west:       303308.13155207<br>east:       303364.49017572<br>nsres:      1.00564787<br>ewres:      1.00640399<br>rows:       44<br>cols:       56<br>cells:      2464</span></div><div><br></div><div>I checked the source code and it "simply" a bug in the printing logic in <a href="http://v.info">v.info</a>, as far as I see.</div><div>Means: all fine, just <a href="http://v.info">v.info</a> gets confused with negative UTM zone values:</div><div><a href="https://github.com/OSGeo/grass/blob/719725c22c9b2ceadbd009ae2f8acdeeef5a0fc9/vector/v.info/print.c#L497">https://github.com/OSGeo/grass/blob/719725c22c9b2ceadbd009ae2f8acdeeef5a0fc9/vector/v.info/print.c#L497</a></div><div><br></div><div>I have opened an issue for this:</div><div><a href="https://github.com/OSGeo/grass/issues/892">https://github.com/OSGeo/grass/issues/892</a></div><div><br></div><div>(with a GitHub account you can "subscribe" to it be remain informed.<br></div><div><br></div><div>You may simply ignore this "invalid zone" message for now.<br></div><div></div><div><br></div><div>Best,</div><div>Markus</div><div><br></div><div><br></div><div><br></div><div><br></div></div><br><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jul 7, 2020 at 2:18 PM Martin Bittens <<a href="mailto:martin.bittens@gmx.de">martin.bittens@gmx.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello Markus,<br>
<br>
what is the edit button?<br>
Is that in GRASS?<br>
<br>
Sorry, that I have confused you a bit. Attached I am sending 2<br>
screenshots for a better explanation (GRASS_Edit_Button.jpg,<br>
GRASS_Edit_Metadata.jpg)<br>
<br>
Just to summarize our discussion:<br>
<br>
1.) Creating a location/mapset by:<br>
<br>
# create new location from SHAPE file:<br>
grass78 -c LT_AOI_3.shp ~/grassdata/sirgas_2000<br>
Starting GRASS GIS...<br>
Creating new GRASS GIS location <sirgas_2000>...<br>
<br>
2.) g.proj -w is showing beside other information the correct<br>
projection: "SIRGAS 2000 / UTM zone 23S"<br>
<br>
3.) Import of LT_AOI_3.shp into the GRASS location/mapset which was<br>
created in 1.)<br>
<br>
4.) Checking the metadata of the imported shape file by using (in GRASS):<br>
<br>
- the edit button (explanation see attached screenshots) results in:<br>
Projection UTM (zone invalid);<br>
<br>
- the <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> command results in: Projection UTM (zone invalid);<br>
<br>
- the <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> -w command results in: Projection UTM, zone -23.<br>
<br>
<br>
That's where we are now.<br>
<br>
<br>
Regards<br>
<br>
<br>
Martin<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Regards<br>
<br>
Martin<br>
<br>
<br>
<br>
<br>
Am 7/2/2020 um 5:53 PM schrieb Markus Neteler:<br>
> Hi Martin,<br>
><br>
> On Thu, Jul 2, 2020 at 10:03 PM Martin Bittens <<a href="mailto:martin.bittens@gmx.de" target="_blank">martin.bittens@gmx.de</a>> wrote:<br>
>> Hello Markus,<br>
>><br>
>><br>
>> that is exactly what I got (see one of my previous emails).<br>
> I read them several times but I am still lost:<br>
><br>
>> But...:-):<br>
>><br>
>> If you import the shape-file you used for the setup of the<br>
>> location/mapset into this location/mapset and then check the metadata of<br>
>> the imported shape-file by<br>
>><br>
>> -using the edit button;<br>
> what is the edit button?<br>
> Is that in GRASS?<br>
><br>
> Sorry, I am doing so many things in parallel, I need a bit of more detail.<br>
><br>
> Cheers<br>
> Markus<br>
><br>
>> -the <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> command;<br>
>><br>
>> -the <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> -w command;<br>
>><br>
>> what do you get?<br>
>><br>
>><br>
>> I get:<br>
>><br>
>> -using the edit button: Projection UTM (zone invalid);<br>
>><br>
>> -the <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> command: Projection UTM (zone invalid);<br>
>><br>
>> -the <a href="http://v.info" rel="noreferrer" target="_blank">v.info</a> -w command: Projection UTM, zone -23<br>
>><br>
>><br>
>> That is what I did not figure out, why GRASS provides different meta<br>
>> data for the same shape-file although the location/mapset meta did not<br>
>> change and remains "SIRGAS 2000 / UTM zone 23S".<br>
>><br>
>><br>
>> Regards<br>
>><br>
>><br>
>> Martin<br>
>><br>
>><br>
>><br>
>> Am 7/2/2020 um 3:26 PM schrieb Markus Neteler:<br>
>>> Hi Martin,<br>
>>><br>
>>> On Thu, Jul 2, 2020 at 1:57 PM Martin Bittens <<a href="mailto:martin.bittens@gmx.de" target="_blank">martin.bittens@gmx.de</a>> wrote:<br>
>>>> Hello Markus,<br>
>>>><br>
>>>> ok, let us try.<br>
>>>><br>
>>>> Attached I send (as zip files):<br>
>>>><br>
>>>> 1.) QGIS_LT_AOI_3.zip with a shape layer created in QGIS 2.14.0. It<br>
>>>> contains contains 5 sample points extracted from an existing site (EPSG<br>
>>>> 31983);<br>
>>>><br>
>>>> 2.) Test_SG.zip with the files from the grassdata directory created by<br>
>>>> GRASS 7.8.3 after establishing the location "Test_SG" (EPSG 31983) and<br>
>>>> importing the shape layer into "Test_SG".<br>
>>> I tried and got this:<br>
>>><br>
>>> # create new location from SHAPE file:<br>
>>> grass78 -c LT_AOI_3.shp ~/grassdata/sirgas_2000<br>
>>> Starting GRASS GIS...<br>
>>> Creating new GRASS GIS location <sirgas_2000>...<br>
>>> ...<br>
>>> Welcome to GRASS GIS 7.8.4dev (0798df3d0)<br>
>>> ...<br>
>>><br>
>>> GRASS 7.8.4dev (sirgas_2000):/tmp > g.proj -w<br>
>>> PROJCS["SIRGAS 2000 / UTM zone 23S",<br>
>>>       GEOGCS["SIRGAS 2000",<br>
>>>           DATUM["Sistema_de_Referencia_Geocentrico_para_las_AmericaS_2000",<br>
>>>               SPHEROID["GRS 1980",6378137,298.257222101,<br>
>>>                   AUTHORITY["EPSG","7019"]],<br>
>>>               TOWGS84[0,0,0,0,0,0,0],<br>
>>>               AUTHORITY["EPSG","6674"]],<br>
>>>           PRIMEM["Greenwich",0,<br>
>>>               AUTHORITY["EPSG","8901"]],<br>
>>>           UNIT["degree",0.0174532925199433,<br>
>>>               AUTHORITY["EPSG","9122"]],<br>
>>>           AUTHORITY["EPSG","4674"]],<br>
>>>       PROJECTION["Transverse_Mercator"],<br>
>>>       PARAMETER["latitude_of_origin",0],<br>
>>>       PARAMETER["central_meridian",-45],<br>
>>>       PARAMETER["scale_factor",0.9996],<br>
>>>       PARAMETER["false_easting",500000],<br>
>>>       PARAMETER["false_northing",10000000],<br>
>>>       UNIT["metre",1,<br>
>>>           AUTHORITY["EPSG","9001"]],<br>
>>>       AXIS["Easting",EAST],<br>
>>>       AXIS["Northing",NORTH],<br>
>>>       AUTHORITY["EPSG","31983"]]<br>
>>><br>
>>> To me it looks as expected.<br>
>>><br>
>>> My software versions:<br>
>>><br>
>>> g.version -rge<br>
>>> version=7.8.4dev<br>
>>> date=2020<br>
>>> revision=0798df3d0<br>
>>> build_date=2020-07-02<br>
>>> build_platform=x86_64-pc-linux-gnu<br>
>>> build_off_t_size=8<br>
>>> libgis_revision=8e4dab571<br>
>>> libgis_date=2020-05-05T09:05:27+00:00<br>
>>> proj=6.3.2<br>
>>> gdal=3.0.4<br>
>>> geos=3.8.1<br>
>>> sqlite=3.32.2<br>
>>><br>
>>> Might it be that your PROJ installation is somehow broken?<br>
>>><br>
>>> best,<br>
>>> Markus<br>
>> --<br>
>> This email has been checked for viruses by Avast antivirus software.<br>
>> <a href="https://www.avast.com/antivirus" rel="noreferrer" target="_blank">https://www.avast.com/antivirus</a><br>
>><br>
><br>
<br>
<br>
--<br>
This email has been checked for viruses by Avast antivirus software.<br>
<a href="https://www.avast.com/antivirus" rel="noreferrer" target="_blank">https://www.avast.com/antivirus</a><br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature">Markus Neteler, PhD<br><a href="https://www.mundialis.de" target="_blank">https://www.mundialis.de</a> - free data with free software<br><a href="https://grass.osgeo.org" target="_blank">https://grass.osgeo.org</a><br><a href="https://courses.neteler.org/blog" target="_blank">https://courses.neteler.org/blog</a></div>