<div dir="ltr"><div>On Thu, Sep 15, 2016 at 1:02 AM, Rich Shepard <<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>> wrote:<br>>   The target location has these PROJ_INFO parameters defined in a Shapefile<br>> *.prj:<br>><br>> name: NAD_1983_HARN_StatePlane_Oregon_North_FIPS_3601_Feet_Intl<br>> datum: nad83harn<br>> ellps: grs80<br>> proj: lcc<br>> lat_1: 44.33333333333334<br>> lat_2: 46<br>> lat_0: 43.66666666666666<br>> lon_0: -120.5<br>> x_0: 2500000<br>> y_0: 0<br>> no_defs: defined<br>><br>><br>> The source location has these PROJ_INFO parameters defined in the<br>> metadata.xml with EPSG 2992:<br>><br>> name: NAD83 / Oregon Lambert (ft)<br>> datum: nad83<br>> ellps: grs80<br>> proj: lcc<br>> lat_1: 43<br>> lat_2: 45.5<br>> lat_0: 41.75<br>> lon_0: -120.5<br>> x_0: 399999.9999984<br>> y_0: 0<br>> no_defs: defined<br>> nadgrids: WO<br>><br>>   When I imported this OFGDB source and selected epsg 2992 a window with 6<br>> flavors appeared. I selected #6: Washington/Oregon HARM adjustment.<br><br>yes, it would apply<br>nadgrids=WO<br><br>according to the list in the location manager selection window (I just tried here with EPSG 2992).<br><br>> Apparently this was not used.<br><br>I checked:<br>- in GRASS GIS 7.2 is is used, should be ok:<br><br>g.proj -p                                                                      <br>-PROJ_INFO-------------------------------------------------<br>name       : NAD83 / Oregon Lambert (ft)<br>datum      : nad83<br>ellps      : grs80<br>proj       : lcc<br>lat_1      : 43<br>lat_2      : 45.5<br>lat_0      : 41.75<br>lon_0      : -120.5<br>x_0        : 399999.9999984<br>y_0        : 0<br>no_defs    : defined<br>nadgrids   : WO<br>-PROJ_EPSG-------------------------------------------------<br>epsg       : 2992<br>-PROJ_UNITS------------------------------------------------<br>unit       : foot<br>units      : feet<br>meters     : 0.3048<br><br>Question: do you use 7.2.svn?<br><br>- in trunk ("7.3") it is not used due to NAD file removal during the recent code sprint upon long discussions with the GDAL maintainer Even Rouault (<a href="https://trac.osgeo.org/grass/changeset/69211">https://trac.osgeo.org/grass/changeset/69211</a>). Here it should use the respective file through GDAL/PROJ4 since we don't want to keep our private copies any longer.<br><br><br>>   From the target location/mapset I tried v.proj<br>> location=ODOT2014mapset=data_source.gdb input=state_roads_2014 (I renamed<br>> the long original name to this with g.rename.) and was told there<br>> was no projection conversion table available.<br><br>Sounds to me that you are using the trunk version of GRASS GIS?<br><br>>   Should I re-import the data and try to ensure that the HARM adjustment is<br>> accepted, or is there a way of reprojecting the source to match the target?<br><br>In the first place I would use GRASS GIS 7.2.svn since the NAD management stuff in trunk is currently experimental.<br>It is a fairly complex beast. The idea is that GRASS no longer bothers about NAD and just receives related info from GDAL/PROJ4. Seems not to be complete yet (I'll update <a href="https://trac.osgeo.org/grass/ticket/2456">https://trac.osgeo.org/grass/ticket/2456</a> accordingly).<br><br>For now, here a cmd line way for GRASS GIS 7.2:<br><br></div><span style="font-family:monospace,monospace"># create dummy location just in order to fetch the list<br># of available datums for given EPSG code:<br></span><div><span style="font-family:monospace,monospace">grass72 -c epsg:2992 ~/grassdata/dummy --exec g.proj -t epsg=2992</span><span style="font-family:monospace,monospace"><span style="font-family:monospace,monospace"> datumtrans=-1<br></span>Cleaning up temporary files...<br>Creating new GRASS GIS location/mapset...<br>Executing <g.proj -t datumtrans=-1 epsg=2992> ...<br>---<br>1<br>Used in whole nad83 region<br>towgs84=0.000,0.000,0.000<br>Default 3-Parameter Transformation (May not be optimum for older datums; use this only if no more appropriate options are available.)<br>---<br>2<br>Used in Florida<br>nadgrids=FL<br>Transforms 'Old NAD83' to 'HPGN NAD83'<br>---<br>3<br>Used in Maryland<br>nadgrids=MD<br>Transforms 'Old NAD83' to 'HPGN NAD83'<br>---<br>4<br>Used in Tennessee<br>nadgrids=TN<br>Transforms 'Old NAD83' to 'HPGN NAD83'<br>---<br>5<br>Used in Wisconsin<br>nadgrids=WI<br>Transforms 'Old NAD83' to 'HPGN NAD83'<br>---<br>6<br>Used in Washington - Oregon<br>nadgrids=WO<br>Transforms 'Old NAD83' to 'HPGN NAD83'<br>Execution of <g.proj -t datumtrans=-1 epsg=2992> finished.<br>Cleaning up temporary files...<br><br><br></span></div><div><span style="font-family:monospace,monospace"># now generate location using EPSG code and related datumtransform #6:<br>grass72 -c epsg:2992 ~/grassdata/oregon2992_nad83_WO --exec g.proj -t epsg=2992 datumtrans=6 -p<br>Cleaning up temporary files...<br>Creating new GRASS GIS location/mapset...<br>Executing <g.proj -t epsg=2992 datumtrans=6 -p> ...<br>-PROJ_INFO-------------------------------------------------<br>name       : NAD83 / Oregon Lambert (ft)<br>datum      : nad83<br>ellps      : grs80<br>proj       : lcc<br>lat_1      : 43<br>lat_2      : 45.5<br>lat_0      : 41.75<br>lon_0      : -120.5<br>x_0        : 399999.9999984<br>y_0        : 0<br>no_defs    : defined<br>nadgrids   : WO<br>-PROJ_EPSG-------------------------------------------------<br>epsg       : 2992<br>-PROJ_UNITS------------------------------------------------<br>unit       : foot<br>units      : feet<br>meters     : 0.3048<br>Execution of <g.proj -t epsg=2992 datumtrans=6 -p> finished.<br>Cleaning up temporary files...<br><br></span></div><div><span style="font-family:monospace,monospace"># eventually start using the new location<br></span></div><div><span style="font-family:monospace,monospace">grass72 ~/grassdata/oregon2992_nad83_WO/PERMANENT/</span><br><br></div><div><br></div><div>But! When entering this location and running<br>g.proj -p<br>still nadgrids is not listed any more... confusing.<br>Anyone with more insights here?</div><div><br>Markus</div></div>