Not super obvious as I'm not seeing a reproduction here on a more recent proj. Anyone with proj6 able to reproduce? P
> On Apr 4, 2022, at 9:18 AM, Jeff Hoffmann <j...@propertykey.com> wrote: > > I am running PostgreSQL via Amazon RDS. Here is the version > information for the instance that exhibits the problem: > > version > ------------------------------------------------------------------------------------------------------------- > PostgreSQL 12.8 on aarch64-unknown-linux-gnu, compiled by > aarch64-unknown-linux-gnu-gcc (GCC) 7.4.0, 64-bit > (1 row) > > postgis_full_version > --------------------------------------------------------------------------------------------------------------------------------------------------------------- > POSTGIS="3.1.5 0" [EXTENSION] PGSQL="120" GEOS="3.8.2-CAPI-1.13.4" > PROJ="6.2.1" LIBXML="2.9.9" LIBJSON="0.12.99" LIBPROTOBUF="1.3.0" > WAGYU="0.5.0 (Internal)" > (1 row) > > Here is an example of a query that produces weird results: > > > select st_asewkt(st_transform('SRID=4269;POLYGON((-81.11425169432798 > 25.139365690162226,-79.88375846050846 > 25.135204935817377,-79.8759485474324 > 25.97511821217247,-81.11505127181294 > 25.97943803781699,-81.11425169432798 25.139365690162226))'::geometry, > 3087)); > > st_asewkt > ----------------------------------------------------------------------------------------------------------------------------------------------------------------------- > SRID=3087;POLYGON((690648.3378382007 129676.3544158864,Infinity > Infinity,Infinity Infinity,688388.1492655604 > 222840.53520809134,690648.3378382007 129676.3544158864)) > (1 row) > > Two of the points are coming up as (Infinty,Infinity). When I run the > same query using the proj4text found in the spatial_ref_sys table > instead of the srid I get the following output: > > select st_asewkt(st_transform('SRID=4269;POLYGON((-81.11425169432798 > 25.139365690162226,-79.88375846050846 > 25.135204935817377,-79.8759485474324 > 25.97511821217247,-81.11505127181294 > 25.97943803781699,-81.11425169432798 25.139365690162226))'::geometry, > '+proj=aea +lat_1=24 +lat_2=31.5 +lat_0=24 +lon_0=-84 +x_0=400000 > +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs')); > > st_asewkt > ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- > POLYGON((690648.3378382007 129676.3544158864,814558.219082916 > 132733.37665360732,812229.019377012 > 225881.33726899035,688388.1492655604 > 222840.53520809134,690648.3378382007 129676.3544158864)) > (1 row) > > This looks pretty close to right. What is the disconnect between > using the srid to call the transform vs the proj4text version of > st_transform? > > -- > Jeff Hoffmann > PropertyKey > j...@propertykey.com > _______________________________________________ > postgis-users mailing list > postgis-users@lists.osgeo.org > https://lists.osgeo.org/mailman/listinfo/postgis-users _______________________________________________ postgis-users mailing list postgis-users@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/postgis-users