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.bitt...@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, MarkusThank 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.bitt...@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 _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
