[GRASS-user] Import Shape Files

Markus Neteler neteler at osgeo.org
Fri Jun 26 07:56:29 PDT 2020


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



-- 
Markus Neteler, PhD
https://www.mundialis.de - free data with free software
https://grass.osgeo.org
https://courses.neteler.org/blog
-------------- next part --------------
A non-text attachment was scrubbed...
Name: LT_AOI_1.prj
Type: application/octet-stream
Size: 836 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20200626/4600ec3e/attachment.obj>


More information about the grass-user mailing list