cockroachdb/cockroach

geo/geomfn: implement ST_DistanceSpheroid({geometry,geometry,spheroid})

otan opened this issue · 2 comments

otan commented

Implement ST_DistanceSpheroid on arguments {geometry,geometry,spheroid}, which should adopt PostGIS behaviour.

Observers: Please react to this issue if you need this functionality.

For Geometry builtins, please do the following:

  • Ideally add a relevant helper function in pkg/geo/geomfn (parse and output related functions can go in pkg/geo). Add exhaustive unit tests here - you can run through example test cases and make sure that PostGIS and CRDB return the same result within a degree of accuracy (1cm for geography).
    • When using GEOS, you can reference the C API for which functions are available. Unfortunately, Windows is not currently supported when using GEOS.
  • Create a new builtin that references this function in pkg/sql/sem/builtins/geo_builtins.go. Note that we currently do not support optional arguments, so we define functions that have optional arguments once without the optional argument (using the default value in the optional position), and once with the optional argument.
  • Modify the tests in pkg/sql/logictest/testdata/logic_test/geospatial to call this functionality at least once. You can call make testbaselogic FILES='geospatial' TESTFLAGS='-rewrite' to regenerate the output. Tests here should just ensure the builtin is linked end to end (your exhaustive unit tests go the above mentioned packages!).
  • Ensure the documentation is regenerated by calling make buildshort. You can also play with it by calling ./cockroach demo --empty afterwards.
  • Submit your PR - make sure to follow the guidelines from creating your first PR.

You can follow #48552 for an example PR.

The following additional guidance has been issued on implementing this function:

use PROJ library to interpret spheroid

🤖 This issue was synced with a spreadsheet by gsheets-to-github-issues by otan on 2023-09-03T23:16:38Z. Changes to titles, body and labels may be overwritten.

Was it an oversight that the spheroid argument to ST_DistanceSpheroid isn't implemented (if I read your commit correctly)?

I'm seeing unknown signature: st_distancespheroid(geometry, geometry, string) for a query like:

SELECT ...
FROM "distapp_australiacity"
WHERE ST_DistanceSpheroid(
    "distapp_australiacity"."point",
    ST_GeomFromEWKB('\x0101000020e6100000c3f5285c8f6a6240f0164850fc7045c0'::bytea),
    'SPHEROID["WGS 84",6378137.0,298.257223563]'
) <= ("distapp_australiacity"."radius" * 70)
otan commented

Oops