<div dir="ltr"><div>Hi Thomas, <br></div><br><div>I have no suggestion other than to create an issue here <a href="https://github.com/OSGeo/grass/issues">https://github.com/OSGeo/grass/issues</a> with a reproducible example so others can easily test. Also report the GDAL and PROJ versions you have.</div><div><br></div><div>Best, <br></div><div>Vero<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El mar, 19 oct 2021 a las 14:50, Thomas Adams (<<a href="mailto:tea3rd@gmail.com">tea3rd@gmail.com</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi all,</div><div><br></div><div>I have run into a problem with the development version of GRASS 7.8.6 and with the latest versions of PROJ and GDAL. I had a significant issue a couple of months ago and had to re-build GRASS, etc. and I have been trying to get back where I was since that time.</div><div><br></div><div>A major part of the work I do is working with US NOAA/NWS raster radar precipitation data in what's called xmrg format... a binary grid file that relies on a very obtuse 'projection' called HRAP, described here (and within linked files): <a href="https://www.weather.gov/owp/oh_hrl_distmodel_hrap" target="_blank">https://www.weather.gov/owp/oh_hrl_distmodel_hrap</a> -- it basically relies on a polar stereographic projection with a spherical Earth model with these parameters:</div><div><br></div>-PROJ_INFO-------------------------------------------------<br>name       : Stereographic<br>proj       : stere<br>a          : 1337.784777<br>es         : 0.0<br>f          : 0.0<br>lat_0      : 90.0000000000<br>lat_ts     : 60.0000000000<br>lon_0      : -105.0000000000<br>k_0        : 1.0000000000<br>x_0        : 401.0<br>y_0        : 1601.0<br>-PROJ_UNITS------------------------------------------------<br>unit       : meters<br>units      : meters<br>meters     : 1.0000000000<div><br></div><div>Previously, I had been able to re-project raster and vector maps to/from this location from other map locations with different projections. Now I cannot do this; I now get this error:</div><div><br></div><div>v.proj location=world_location mapset=teaiii input=st_us                        <br>proj_create_operations: Source and target ellipsoid do not belong to the same celestial body<br>proj_create_operations: Source and target ellipsoid do not belong to the same celestial body<br>WARNING: proj_create() failed for '(null)'<br>ERROR: Unable to initialize coordinate transformation<br><br></div><div>Any ideas what I might try?</div><div><br></div><div>Thank you,</div><div>Tom<br></div><div><br></div><div><br></div><div><br><br></div></div>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</blockquote></div>