[GRASS-user] Import Shape Files

Markus Neteler neteler at osgeo.org
Tue Aug 11 03:06:16 PDT 2020


Hi Martin,

(back to cc of list)

I just got a chance to test and also got the "invalid zone" message issue:

v.info LT_AOI_3

+----------------------------------------------------------------------------+
 | Name:            LT_AOI_3
   |
 | Mapset:          PERMANENT
  |
 | Location:        sirgas_2000
  |
 | Database:        /home/mneteler/grassdata
   |
 | Title:
  |
 | Map scale:       1:1
  |
 | Name of creator: mneteler
   |
 | Organization:
   |
 | Source date:     Tue Aug 11 11:44:55 2020
   |
 | Timestamp (first layer): none
   |
 |----------------------------------------------------------------------------|
 | Map format:      native
   |
 |----------------------------------------------------------------------------|
 |   Type of map: vector (level: 2)
  |
 |
   |
 |   Number of points:       5               Number of centroids:  0
   |
...
 |   Projection: UTM (zone invalid)
  |
...

This looks ok:

g.region vector=LT_AOI_3 -p
projection: 1 (UTM)
zone:       -23
datum:      sirgas2000
ellipsoid:  grs80
north:      7397570.18085686
south:      7397525.93235069
west:       303308.13155207
east:       303364.49017572
nsres:      1.00564787
ewres:      1.00640399
rows:       44
cols:       56
cells:      2464

I checked the source code and it "simply" a bug in the printing logic in
v.info, as far as I see.
Means: all fine, just v.info gets confused with negative UTM zone values:
https://github.com/OSGeo/grass/blob/719725c22c9b2ceadbd009ae2f8acdeeef5a0fc9/vector/v.info/print.c#L497

I have opened an issue for this:
https://github.com/OSGeo/grass/issues/892

(with a GitHub account you can "subscribe" to it be remain informed.

You may simply ignore this "invalid zone" message for now.

Best,
Markus






On Tue, Jul 7, 2020 at 2:18 PM Martin Bittens <martin.bittens at gmx.de> wrote:

> Hello Markus,
>
> what is the edit button?
> Is that in GRASS?
>
> Sorry, that I have confused you a bit. Attached I am sending 2
> screenshots for a better explanation (GRASS_Edit_Button.jpg,
> GRASS_Edit_Metadata.jpg)
>
> Just to summarize our discussion:
>
> 1.) Creating a location/mapset by:
>
> # create new location from SHAPE file:
> grass78 -c LT_AOI_3.shp ~/grassdata/sirgas_2000
> Starting GRASS GIS...
> Creating new GRASS GIS location <sirgas_2000>...
>
> 2.) g.proj -w is showing beside other information the correct
> projection: "SIRGAS 2000 / UTM zone 23S"
>
> 3.) Import of LT_AOI_3.shp into the GRASS location/mapset which was
> created in 1.)
>
> 4.) Checking the metadata of the imported shape file by using (in GRASS):
>
> - the edit button (explanation see attached screenshots) results in:
> Projection UTM (zone invalid);
>
> - the v.info command results in: Projection UTM (zone invalid);
>
> - the v.info -w command results in: Projection UTM, zone -23.
>
>
> That's where we are now.
>
>
> Regards
>
>
> Martin
>
>
>
>
>
>
>
> Regards
>
> Martin
>
>
>
>
> Am 7/2/2020 um 5:53 PM schrieb Markus Neteler:
> > Hi Martin,
> >
> > On Thu, Jul 2, 2020 at 10:03 PM Martin Bittens <martin.bittens at gmx.de>
> wrote:
> >> Hello Markus,
> >>
> >>
> >> that is exactly what I got (see one of my previous emails).
> > I read them several times but I am still lost:
> >
> >> But...:-):
> >>
> >> If you import the shape-file you used for the setup of the
> >> location/mapset into this location/mapset and then check the metadata of
> >> the imported shape-file by
> >>
> >> -using the edit button;
> > what is the edit button?
> > Is that in GRASS?
> >
> > Sorry, I am doing so many things in parallel, I need a bit of more
> detail.
> >
> > Cheers
> > Markus
> >
> >> -the v.info command;
> >>
> >> -the v.info -w command;
> >>
> >> what do you get?
> >>
> >>
> >> I get:
> >>
> >> -using the edit button: Projection UTM (zone invalid);
> >>
> >> -the v.info command: Projection UTM (zone invalid);
> >>
> >> -the v.info -w command: Projection UTM, zone -23
> >>
> >>
> >> That is what I did not figure out, why GRASS provides different meta
> >> data for the same shape-file although the location/mapset meta did not
> >> change and remains "SIRGAS 2000 / UTM zone 23S".
> >>
> >>
> >> Regards
> >>
> >>
> >> Martin
> >>
> >>
> >>
> >> Am 7/2/2020 um 3:26 PM schrieb Markus Neteler:
> >>> Hi Martin,
> >>>
> >>> On Thu, Jul 2, 2020 at 1:57 PM Martin Bittens <martin.bittens at gmx.de>
> wrote:
> >>>> Hello Markus,
> >>>>
> >>>> ok, let us try.
> >>>>
> >>>> Attached I send (as zip files):
> >>>>
> >>>> 1.) QGIS_LT_AOI_3.zip with a shape layer created in QGIS 2.14.0. It
> >>>> contains contains 5 sample points extracted from an existing site
> (EPSG
> >>>> 31983);
> >>>>
> >>>> 2.) Test_SG.zip with the files from the grassdata directory created by
> >>>> GRASS 7.8.3 after establishing the location "Test_SG" (EPSG 31983) and
> >>>> importing the shape layer into "Test_SG".
> >>> I tried and got this:
> >>>
> >>> # create new location from SHAPE file:
> >>> grass78 -c LT_AOI_3.shp ~/grassdata/sirgas_2000
> >>> Starting GRASS GIS...
> >>> Creating new GRASS GIS location <sirgas_2000>...
> >>> ...
> >>> Welcome to GRASS GIS 7.8.4dev (0798df3d0)
> >>> ...
> >>>
> >>> GRASS 7.8.4dev (sirgas_2000):/tmp > 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"]],
> >>>               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"]]
> >>>
> >>> To me it looks as expected.
> >>>
> >>> My software versions:
> >>>
> >>> g.version -rge
> >>> version=7.8.4dev
> >>> date=2020
> >>> revision=0798df3d0
> >>> build_date=2020-07-02
> >>> build_platform=x86_64-pc-linux-gnu
> >>> build_off_t_size=8
> >>> libgis_revision=8e4dab571
> >>> libgis_date=2020-05-05T09:05:27+00:00
> >>> proj=6.3.2
> >>> gdal=3.0.4
> >>> geos=3.8.1
> >>> sqlite=3.32.2
> >>>
> >>> Might it be that your PROJ installation is somehow broken?
> >>>
> >>> 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
>


-- 
Markus Neteler, PhD
https://www.mundialis.de - free data with free software
https://grass.osgeo.org
https://courses.neteler.org/blog
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20200811/09268787/attachment-0001.html>


More information about the grass-user mailing list