Bug#1008612: closed by Sebastiaan Couwenberg <sebastic at xs4all.nl> (Re: Bug#1008612: libspatialite7: erroneous geometry conversion using ST_Transform)
Brian Miles
selimnairb at gmail.com
Tue Mar 29 16:35:12 BST 2022
Hi Bas,
I apologize that I mis-stated my original description of the bug. It is not that the results for 26910 transforms are the same as for 32610 transforms, but that transforms to 26910 are simply incorrect for Debian/Ubuntu.
As this is only happening on Debian/Ubuntu, it is not clear to me that bringing this upstream will result in resolution. Are there any other option to addressing this?
I am happy to help with testing/debugging if that helps.
Thanks,
Brian
> On Mar 29, 2022, at 11:18 AM, Debian Bug Tracking System <owner at bugs.debian.org> wrote:
>
> This is an automatic notification regarding your Bug report
> which was filed against the libspatialite7 package:
>
> #1008612: libspatialite7: erroneous geometry conversion using ST_Transform
>
> It has been closed by Sebastiaan Couwenberg <sebastic at xs4all.nl>.
>
> Their explanation is attached below along with your original report.
> If this explanation is unsatisfactory and you have not received a
> better one in a separate message then please contact Sebastiaan Couwenberg <sebastic at xs4all.nl> by
> replying to this email.
>
>
> --
> 1008612: https://bugs.debian.org/cgi-bin/bugreport.cgi?bug=1008612
> Debian Bug Tracking System
> Contact owner at bugs.debian.org with problems
>
> From: Sebastiaan Couwenberg <sebastic at xs4all.nl>
> Subject: Re: Bug#1008612: libspatialite7: erroneous geometry conversion using ST_Transform
> Date: March 29, 2022 at 11:05:11 AM EDT
> To: Brian Miles <selimnairb at gmail.com>, 1008612-done at bugs.debian.org
>
>
> Hi Brian,
>
> I'm not in a position to address this issue, please take it upstream. Either the mailing list:
>
> https://groups.google.com/g/spatialite-users
>
> Or the issue tracker:
>
> https://www.gaia-gis.it/fossil/libspatialite/index
>
> Be aware that upstream development is very slow, don't get your hopes up for a quick fix.
>
> Kind Regards,
>
> Bas
>
> --
> GPG Key ID: 4096R/6750F10AE88D4AF1
> Fingerprint: 8182 DE41 7056 408D 6146 50D1 6750 F10A E88D 4AF1
>
>
> From: Brian Miles <selimnairb at gmail.com>
> Subject: libspatialite7: erroneous geometry conversion using ST_Transform
> Date: March 29, 2022 at 10:47:18 AM EDT
> To: submit at bugs.debian.org
>
>
> Package: libspatialite7
> Version: 5.0.1-2
>
> Hello,
>
> I have discovered a bug in libspatialite7’s ST_Transform() function where transforms from WGS84 geographic coordinates (EPSG:4326) to NAD83 UTM 10N (26910) returns the result for WGS84 UTM 10N (32610). This appears to only be occurring on Debian 11 (and perhaps other versions, as well as Ubuntu 20.04, and perhaps other versions). The behavior on macOS and Alpine Linux appears to be correct. Here are the steps to reproduce:
>
> # DB Setup
>
> CREATE TABLE IF NOT EXISTS test(id INTEGER PRIMARY KEY);
> SELECT AddGeometryColumn ('test', 'point', 4326, 'POINT', 'XY', 1);
>
> INSERT INTO test (point) VALUES (GeomFromText('POINT(-122.338928 47.629273)', 4326));
> INSERT INTO test (point) VALUES (GeomFromText('POINT(-122.338923 47.62927)', 4326));
> INSERT INTO test (point) VALUES (GeomFromText('POINT(-122.33895 47.629253)', 4326));
> INSERT INTO test (point) VALUES (GeomFromText('POINT(-122.392051 47.62888)', 4326));
> INSERT INTO test (point) VALUES (GeomFromText('POINT(-122.392038 47.62888)', 4326));
>
>
> # Results from Debian 11 (ARM64 and x86_64) and Ubuntu 22.04 (ARM64):
>
> spatialite> SELECT AsText(point), AsText(ST_Transform(point, 26910)) FROM test;
> POINT(-122.338928 47.629273)|POINT(549665.036792 5275308.276363)
> POINT(-122.338923 47.62927)|POINT(549665.415273 5275307.94615)
> POINT(-122.33895 47.629253)|POINT(549663.402919 5275306.039517)
> POINT(-122.392051 47.62888)|POINT(545674.35684 5275231.954471)
> POINT(-122.392038 47.62888)|POINT(545675.333511 5275231.962125)
>
> spatialite> SELECT AsText(point), AsText(ST_Transform(point, 32610)) FROM test;
> POINT(-122.338928 47.629273)|POINT(549665.158878 5275308.476981)
> POINT(-122.338923 47.62927)|POINT(549665.53736 5275308.14677)
> POINT(-122.33895 47.629253)|POINT(549663.525012 5275306.240133)
> POINT(-122.392051 47.62888)|POINT(545674.479273 5275232.145082)
> POINT(-122.392038 47.62888)|POINT(545675.455944 5275232.152738)
> spatialite>
>
>
> # Results from macOS 12 (Homebrew or MacPorts; ARM64 and x86_64):
>
> spatialite> SELECT AsText(point), AsText(ST_Transform(point, 26910)) FROM test;
> AsText(point) AsText(ST_Transform(point, 26910))
> ---------------------------- ----------------------------------
> POINT(-122.338928 47.629273) POINT(549665.158879 5275308.47686)
> POINT(-122.338923 47.62927) POINT(549665.537361 5275308.146648)
> POINT(-122.33895 47.629253) POINT(549663.525012 5275306.240011)
> POINT(-122.392051 47.62888) POINT(545674.479273 5275232.14496)
> POINT(-122.392038 47.62888) POINT(545675.455944 5275232.152617)
>
> spatialite> SELECT AsText(point), AsText(ST_Transform(point, 32610)) FROM test;
> AsText(point) AsText(ST_Transform(point, 32610))
> ---------------------------- -----------------------------------
> POINT(-122.338928 47.629273) POINT(549665.158878 5275308.476981)
> POINT(-122.338923 47.62927) POINT(549665.53736 5275308.14677)
> POINT(-122.33895 47.629253) POINT(549663.525012 5275306.240133)
> POINT(-122.392051 47.62888) POINT(545674.479273 5275232.145082)
> POINT(-122.392038 47.62888) POINT(545675.455944 5275232.152738)
> spatialite>
>
>
> # Results from Alpine Linux 3.15 (ARM64 and x86_64):
>
> sqlite> SELECT AsText(point), AsText(ST_Transform(point, 26910)) FROM test;
> POINT(-122.338928 47.629273)|POINT(549665.158879 5275308.47686)
> POINT(-122.338923 47.62927)|POINT(549665.537361 5275308.146648)
> POINT(-122.33895 47.629253)|POINT(549663.525012 5275306.240011)
> POINT(-122.392051 47.62888)|POINT(545674.479273 5275232.14496)
> POINT(-122.392038 47.62888)|POINT(545675.455944 5275232.152617)
>
> sqlite> SELECT AsText(point), AsText(ST_Transform(point, 32610)) FROM test;
> POINT(-122.338928 47.629273)|POINT(549665.158878 5275308.476981)
> POINT(-122.338923 47.62927)|POINT(549665.53736 5275308.14677)
> POINT(-122.33895 47.629253)|POINT(549663.525012 5275306.240133)
> POINT(-122.392051 47.62888)|POINT(545674.479273 5275232.145082)
> POINT(-122.392038 47.62888)|POINT(545675.455944 5275232.152738)
>
>
> To make sure this wasn’t PROJ that was broken, I did the following PROJ queries:
>
> # Results from Debian 11 (ARM64 and x86_64) and Ubuntu 22.04 (ARM64):
>
> echo -122.338928 47.629273 | cct -c1,2 -d6 -z0 -t0 +proj=utm +zone=10 +ellps=GRS80 +datum=NAD83 +units=m +no_defs
> 549665.158878 5275308.476964 0.000057 0.0000
>
> echo -122.338928 47.629273 | cct -c1,2 -d6 -z0 -t0 +proj=utm +zone=10 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
> 549665.158878 5275308.476981 0.000000 0.0000
>
> # Results from macOS 12 (Homebrew or MacPorts; ARM64 and x86_64):
>
> echo -122.338928 47.629273 | cct -c1,2 -d6 -z0 -t0 +proj=utm +zone=10 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
> 549665.158878 5275308.476981 0.000000 0.0000
>
> echo -122.338928 47.629273 | cct -c1,2 -d6 -z0 -t0 +proj=utm +zone=10 +ellps=GRS80 +datum=NAD83 +units=m +no_defs
> 549665.158878 5275308.476964 0.000057 0.0000
>
> # Results from Alpine Linux 3.15 (ARM64 and x86_64):
>
> echo -122.338928 47.629273 | cct -c1,2 -d6 -z0 -t0 +proj=utm +zone=10 +ellps=GRS80 +datum=NAD83 +units=m +no_defs
> 549665.158878 5275308.476964 0.000057 0.0000
>
> echo -122.338928 47.629273 | cct -c1,2 -d6 -z0 -t0 +proj=utm +zone=10 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
> 549665.158878 5275308.476964 0.000057 0.0000
>
>
> So it seems that the PROJ output is the same for all three (or six if you include architecture), so the problem likely lies in libspatialite7 on Debian/Ubuntu.
>
>
> Please let me know if you need more information.
>
> Thanks,
>
> Brian
>
>
More information about the Pkg-grass-devel
mailing list