[postgis-tickets] [PostGIS] #4462: st_asmvtgeom fails with IllegalArgumentException: RobustDeterminant encountered non-finite numbers

[PostGIS]

#4462: st_asmvtgeom fails with IllegalArgumentException: RobustDeterminant
encountered non-finite numbers
----------------------+---------------------------
  Reporter:  autra    |      Owner:  pramsey
      Type:  defect   |     Status:  new
  Priority:  medium   |  Milestone:  PostGIS 2.5.3
 Component:  postgis  |    Version:  2.5.x
Resolution:           |   Keywords:
----------------------+---------------------------

Comment (by Algunenano):

 I've tested the first query in 3.0 and it works fine, both with and
 without wagyu (the new clipping library).

 Without wagyu:
 {{{
 POLYGON((536 4352,466 4352,486 4287,551 4170,551 4103,515 4077,479
 4068,405 4084,368 4126,366 4352,300 4352,312 4310,313 4159,305 4025,334
 3891,353 3832,408 3740,648 3524,667 3499,824 3358,842 3324,861 3249,899
 3165,991 3107,1119 3074,1210 3075,1229 3092,1237 3126,1228 3159,1136
 3201,1090 3251,1071 3309,979 3409,905 3493,867 3652,848 3710,764 3819,718
 3902,708 4003,643 4128,632 4229,614 4254,540 4346,536 4352),(497 4043,552
 4010,617 3943,654 3877,682 3810,729 3735,775 3676,794 3634,822 3517,841
 3467,832 3450,759 3466,722 3500,657 3600,528 3725,472 3816,407 3900,388
 3958,415 4017,451 4043,497 4043))
 }}}

 With wagyu:
 {{{
 POLYGON((466 4352,486 4287,551 4170,551 4103,515 4077,479 4068,405
 4084,368 4126,366 4352,300 4352,312 4310,313 4159,305 4025,334 3891,353
 3832,408 3740,648 3524,667 3499,824 3358,842 3324,861 3249,899 3165,991
 3107,1119 3074,1210 3075,1229 3092,1237 3126,1228 3159,1136 3201,1090
 3251,1071 3309,979 3409,905 3493,867 3652,848 3710,764 3819,718 3902,708
 4003,643 4128,632 4229,614 4254,540 4346,536 4352,466 4352),(8
 32 3450,759 3466,722 3500,657 3600,528 3725,472 3816,407 3900,388 3958,415
 4017,451 4043,497 4043,552 4010,617 3943,654 3877,682 3810,729 3735,775
 3676,794 3634,822 3517,841 3467,832 3450))
 }}}

 In both cases I'm using GEOS trunk:
 {{{
 # Select postgis_full_version();
 postgis_full_version

 -------------------------------------------------------------------------------------------------------------------------------------------
 -------------------------------------
  POSTGIS="3.0.0alpha4dev r17604" [EXTENSION] PGSQL="110" GEOS="3.8.0dev-
 CAPI-1.11.0 " PROJ="6.1.0" LIBXML="2.9.9" LIBJSON="0.13.1"
 LIBPROTOBUF="1.3.1" WAGYU="0.4.3 (Internal)"
 }}}


 I've also downgraded GEOS to 3.7.2 and it also works fine (Postgis 3.0
 still):
 {{{
 mvt_test=# Select postgis_full_version();
 postgis_full_version

 -------------------------------------------------------------------------------------------------------------------------------------------
 -----------------
  POSTGIS="3.0.0alpha4dev r17604" [EXTENSION] PGSQL="110"
 GEOS="3.7.2-CAPI-1.11.2 b55d2125" PROJ="6.1.0" LIBXML="2.9.9"
 LIBJSON="0.13.1" LIB
 PROTOBUF="1.3.1"
 (1 row)

 mvt_test=# select ST_AsText(st_asmvtgeom(geom,
 st_setsrid(st_makebox2d(st_point(1645442.8949999844, 8184994.56999849),
 st_point
 (1647681.2699999844,8187432.94499849)), 3949))) as mvtgeom from copy3;


                       mvtgeom


 -------------------------------------------------------------------------------------------------------------------------------------------
 -------------------------------------------------------------------------------------------------------------------------------------------
 -------------------------------------------------------------------------------------------------------------------------------------------
 -------------------------------------------------------------------------------------------------------------------------------------------
 ---------------------------------------------------
  POLYGON((536 4352,466 4352,486 4287,551 4170,551 4103,515 4077,479
 4068,405 4084,368 4126,366 4352,300 4352,312 4310,313 4159,305 4025,334
  3891,353 3832,408 3740,648 3524,667 3499,824 3358,842 3324,861 3249,899
 3165,991 3107,1119 3074,1210 3075,1229 3092,1237 3126,1228 3159,11
 36 3201,1090 3251,1071 3309,979 3409,905 3493,867 3652,848 3710,764
 3819,718 3902,708 4003,643 4128,632 4229,614 4254,540 4346,536 4352),(4
 97 4043,552 4010,617 3943,654 3877,682 3810,729 3735,775 3676,794 3634,822
 3517,841 3467,832 3450,759 3466,722 3500,657 3600,528 3725,472 3
 816,407 3900,388 3958,415 4017,451 4043,497 4043))
 (1 row)

 }}}



 Can you please confirm the exact version of postgis and GEOS that you are
 running?

 The latest release of Postgis is 2.5.2 which should have the same MVTs as
 3.0.0 without wagyu, so I'm inclined to agree with Paul and blame this on
 an outdated GEOS.

-- 
Ticket URL: <https://trac.osgeo.org/postgis/ticket/4462#comment:3>
PostGIS <http://trac.osgeo.org/postgis/>
The PostGIS Trac is used for bug, enhancement & task tracking, a user and developer wiki, and a view into the subversion code repository of PostGIS project.