Re: [GRASS-user] Which EPSG code for location default?
On Thu, 23 Sep 2021, Jeff McKenna wrote: - goto https://epsg.org Jeff, Thanks for excellent advice. I decided to be consistent and use the NAD83(HARN)/meters EPSG for Washington South to be consistent with what I've used in other projects and other states. That EPSG code is 2856. Now I'll create a new location with that projection and reproject all imported maps to that location and rename it to the project's location name. I think it was in the late 1980s-early 1990s that I used UTM zones before switching to State Plane Coordinates. On this project the first data set I imported used UTM Zone 10 and I stuck with it ... until now. Regards, Rich ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Which EPSG code for location default?
On Thu, 23 Sep 2021, Micha Silver wrote: IMHO you should *not* do that. Just create a new Location. Micha, I thought that would be the answer. Then I need to reproject all data at least twice, then rename the new location. Regards, Rich ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Which EPSG code for location default?
On 9/23/21 3:38 PM, Rich Shepard wrote: Yet, ... I have a question: is there a way in grass to change the PROJ_INFO in the PERMANENT mapset of a location to a new projecion? Or, do I manually edit that file? IMHO you should *not* do that. Just create a new Location. Thanks in advance, Rich ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user -- Micha Silver Ben Gurion Univ. Sde Boker, Remote Sensing Lab cell: +972-523-665918 https://orcid.org/-0002-1128-1325 ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Which EPSG code for location default?
On Thu, 23 Sep 2021, Jeff McKenna wrote: On 2021-09-23 8:55 a.m., Jeff McKenna wrote: NAD83 / UTM zone 10N + NAVD88 It could also be that no one yet has submitted this to the EPSG geodetic committee, for inclusion (UTM zone 10N with that vertical datum). *hint* *maybe you should* (PS I've submitted several, just use the mindset "well someone has to") Jeff, Interesting. Actually, since this project (like all of them) have data sources using different projections I don't care which is used as long as I can reproject data from the source to a single appropriate one. I could use Washington State Plane (South) and that would work just as well, So.all of my answers to you are very unhelpful, ha. Not at all. But that's all I've got for you now. And that's good enough. Yet, ... I have a question: is there a way in grass to change the PROJ_INFO in the PERMANENT mapset of a location to a new projecion? Or, do I manually edit that file? Thanks in advance, Rich ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Which EPSG code for location default?
On Thu, 23 Sep 2021, Jeff McKenna wrote: Taking a look at your output, my brain said "good question", ha. Jeff, There are so many value-specific codes I didn't see an overall code. Usually the key line to look for is the "AUTHORITY:" line, which yours lacks. Ah-ha. That's good information. So the usual next step is to start searching on https://epsg.io for the text "UTM zone 10N". But more recently the official EPSG homepage includes a search tool that does a live search on its own database (giving guaranteed accurate/uptodate responses): - goto https://epsg.org - click the top menu: EPSG Dataset/Text Search - enter: UTM zone 10N and press "Go" - gives a response of 12 items EPSG:26910 does look somewhat similar :) I'll check for myself; I had the impression that all EPSG codes were four digits. Hope that gives you some hints... (at least on the official EPSG code perspective) It certainly does. Thanks very much, Rich ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Which EPSG code for location default?
On 2021-09-23 8:55 a.m., Jeff McKenna wrote: NAD83 / UTM zone 10N + NAVD88 It could also be that no one yet has submitted this to the EPSG geodetic committee, for inclusion (UTM zone 10N with that vertical datum). *hint* *maybe you should* (PS I've submitted several, just use the mindset "well someone has to") So.all of my answers to you are very unhelpful, ha. But that's all I've got for you now. -jeff -- Jeff McKenna GatewayGeo: Developers of MS4W, MapServer Consulting and Training co-founder of FOSS4G http://gatewaygeo.com/ ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Which EPSG code for location default?
Hi Rich, Taking a look at your output, my brain said "good question", ha. Usually the key line to look for is the "AUTHORITY:" line, which yours lacks. So the usual next step is to start searching on https://epsg.io for the text "UTM zone 10N". But more recently the official EPSG homepage includes a search tool that does a live search on its own database (giving guaranteed accurate/uptodate responses): - goto https://epsg.org - click the top menu: EPSG Dataset/Text Search - enter: UTM zone 10N and press "Go" - gives a response of 12 items EPSG:26910 does look somewhat similar :) Hope that gives you some hints... (at least on the official EPSG code perspective) -jeff -- Jeff McKenna GatewayGeo: Developers of MS4W, MapServer Consulting and Training co-founder of FOSS4G http://gatewaygeo.com/ On 2021-09-22 7:49 p.m., Rich Shepard wrote: Applying gdalsrsinfo to a shapefile's .prj produced these results: PROJ.4 : +proj=utm +zone=10 +datum=NAD83 +units=m +vunits=m +no_defs OGC WKT2:2018 : COMPOUNDCRS["NAD83 / UTM zone 10N + NAVD88 height", PROJCRS["NAD83 / UTM zone 10N", BASEGEOGCRS["NAD83", DATUM["North American Datum 1983", ELLIPSOID["GRS 1980",6378137,298.257222101, LENGTHUNIT["metre",1]], ID["EPSG",6269]], PRIMEM["Greenwich",0, ANGLEUNIT["Degree",0.0174532925199433]]], CONVERSION["UTM zone 10N", METHOD["Transverse Mercator", ID["EPSG",9807]], PARAMETER["Latitude of natural origin",0, ANGLEUNIT["Degree",0.0174532925199433], ID["EPSG",8801]], PARAMETER["Longitude of natural origin",-123, ANGLEUNIT["Degree",0.0174532925199433], ID["EPSG",8802]], PARAMETER["Scale factor at natural origin",0.9996, SCALEUNIT["unity",1], ID["EPSG",8805]], PARAMETER["False easting",50, LENGTHUNIT["metre",1], ID["EPSG",8806]], PARAMETER["False northing",0, LENGTHUNIT["metre",1], ID["EPSG",8807]], ID["EPSG",16010]], CS[Cartesian,2], AXIS["(E)",east, ORDER[1], LENGTHUNIT["metre",1, ID["EPSG",9001]]], AXIS["(N)",north, ORDER[2], LENGTHUNIT["metre",1, ID["EPSG",9001, VERTCRS["NAVD88 height", VDATUM["North American Vertical Datum 1988"], CS[vertical,1], AXIS["gravity-related height (H)",up, LENGTHUNIT["metre",1, ID["EPSG",9001] The data are the Columbia River and southern Washington state. It's been a couple of decades since I last dug into proj and I'd appreciate being pointed to a document that will teach me which of the several EPSG codes above is the one to which I want to apply to PROJ_INFO in the PERMANENT mapset. TIA, Rich ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Which EPSG code for location default?
Applying gdalsrsinfo to a shapefile's .prj produced these results: PROJ.4 : +proj=utm +zone=10 +datum=NAD83 +units=m +vunits=m +no_defs OGC WKT2:2018 : COMPOUNDCRS["NAD83 / UTM zone 10N + NAVD88 height", PROJCRS["NAD83 / UTM zone 10N", BASEGEOGCRS["NAD83", DATUM["North American Datum 1983", ELLIPSOID["GRS 1980",6378137,298.257222101, LENGTHUNIT["metre",1]], ID["EPSG",6269]], PRIMEM["Greenwich",0, ANGLEUNIT["Degree",0.0174532925199433]]], CONVERSION["UTM zone 10N", METHOD["Transverse Mercator", ID["EPSG",9807]], PARAMETER["Latitude of natural origin",0, ANGLEUNIT["Degree",0.0174532925199433], ID["EPSG",8801]], PARAMETER["Longitude of natural origin",-123, ANGLEUNIT["Degree",0.0174532925199433], ID["EPSG",8802]], PARAMETER["Scale factor at natural origin",0.9996, SCALEUNIT["unity",1], ID["EPSG",8805]], PARAMETER["False easting",50, LENGTHUNIT["metre",1], ID["EPSG",8806]], PARAMETER["False northing",0, LENGTHUNIT["metre",1], ID["EPSG",8807]], ID["EPSG",16010]], CS[Cartesian,2], AXIS["(E)",east, ORDER[1], LENGTHUNIT["metre",1, ID["EPSG",9001]]], AXIS["(N)",north, ORDER[2], LENGTHUNIT["metre",1, ID["EPSG",9001, VERTCRS["NAVD88 height", VDATUM["North American Vertical Datum 1988"], CS[vertical,1], AXIS["gravity-related height (H)",up, LENGTHUNIT["metre",1, ID["EPSG",9001] The data are the Columbia River and southern Washington state. It's been a couple of decades since I last dug into proj and I'd appreciate being pointed to a document that will teach me which of the several EPSG codes above is the one to which I want to apply to PROJ_INFO in the PERMANENT mapset. TIA, Rich ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user