<div dir="ltr">Dear proj list,<div>I'm having a bit of trouble with proj and I'm hoping someone can tell me where I've gone wrong. For context, my primary issue is converting a raster from NAD83 state plane to NAD27 state plane using GDAL, which uses proj under the hood. It looks like it's proj that is causing me problems. I used conda to install a fresh environment with proj8 and proj-data on a RHEL7 machine.</div><div>Here's the conversion that I think is giving the wrong answers (EPSG 6570 is NAD83(2011) South Carolina State Plane (ft) and EPSG 32033 is NAD27 South Carolina South):</div><div><br></div><div>cs2cs EPSG:6570 EPSG:32033 << EOF</div><div>2320610.000 353815.000 0</div><div>EOF</div><div><br></div><div>which returns:</div><div>2320626.52      353751.41 0.00<br><div><div dir="ltr" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Using the NGS NCAT online program, that should be:</font></div><div><span style="color:rgb(51,51,51);font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13.2px">2320568.869   </span><span style="color:rgb(51,51,51);font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13.2px">353687.842</span><font face="arial, helvetica, sans-serif"><br></font></div><div><span style="color:rgb(51,51,51);font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13.2px"><br></span></div><div><span style="color:rgb(51,51,51);font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13.2px">I did the transform in parts and found that the issue is in the reference frame transform between NAD83(2011) lat lon and NAD27 lat lon, where no change happens, but it should. That is:</span></div><div>cs2cs EPSG:6318 EPSG:4267<span style="color:rgb(51,51,51);font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13.2px"><br></span></div><div>with inputs of </div><div>32d48'4.701"N   79d57'23.546"W 0.000<br></div><div>has outputs of</div><div>32d48'4.701"N   79d57'23.546"W 0.000<br></div><div><br></div><div>If I use NAD83 (EPSG 4269) instead of NAD83(2011), things look right and I get output of:</div><div>32d48'4.066"N   79d57'24.227"W 0.000<br></div><div><br></div><div>Trying to change from NAD83(2011) to NAD83 (EPSG 6318 to 4269) gives no change, so I'm assuming PROJ doesn't know that transform. I'm guessing that because it doesn't know that transform, it can't figure out that it should do the NAD83 to NAD27 part, so instead of being a couple meters off, it's off by 50 meters or so. If I do the transforms in parts, SP83 to NAD83(2011), NAD83(2011) to NAD83, NAD83 to NAD27, NAD27 to SP27, I get the right answer within the couple meters due to the missing transform.</div><div><br></div><div>Is the transform from NAD83(2011) to NAD83 something I just missed installing? Without it, it looks like my GDAL process will be painful. Is there a better way to handle it?</div><div><br></div><div><font face="arial, helvetica, sans-serif">Kirk Waters, PhD                     | NOAA Office for Coastal Management<br></font></div><div><font face="arial, helvetica, sans-serif">Applied Sciences Program      | 2234 South Hobson Ave<br>843-740-1227 (empty office)   | Charleston, SC 29405</font></div><div><font face="arial, helvetica, sans-serif">843-324-2203 (cell during COVID)    </font><div><font face="arial, helvetica, sans-serif"><a href="http://coast.noaa.gov/digitalcoast" target="_blank">coast.noaa.gov/digitalcoast</a></font><br></div></div><div><br></div></div></div></div></div></div></div></div></div></div></div>