[postgis-devel] [PostGIS] #1213: [raster] RT_Intersects test fails

[Sandro Santilli]

On Sun, Oct 02, 2011 at 09:11:47PM -0000, PostGIS wrote:
> #1213: [raster] RT_Intersects test fails
> ----------------------------+-----------------------------------------------
>  Reporter:  strk            |       Owner:  dustymugs    
>      Type:  defect          |      Status:  assigned     
>  Priority:  medium          |   Milestone:  PostGIS 2.0.0
> Component:  postgis raster  |     Version:  trunk        
>  Keywords:                  |  
> ----------------------------+-----------------------------------------------
> 
> Comment(by dustymugs):
> 
>  Replying to [comment:15 strk]:
>  > I've updated gdal to !r23167 and the new failure (on the 32bit, with
>  postgis r7931) is:
>  > {{{
>  > --- rt_intersects_expected      2011-10-01 13:46:55.000000000 +0200
>  > +++ /tmp/pgis_reg_23053/test_58_out     2011-10-02 20:54:56.000000000
>  +0200
>  > @@ -218,7 +218,7 @@
>  >  2.6|2|21|ST_LineString|t
>  >  2.6|2|22|ST_LineString|t
>  >  2.6|2|23|ST_LineString|t
>  > -2.6|2|24|ST_LineString|t
>  > +2.6|2|24|ST_LineString|f
>  >  2.6|2|25|ST_LineString|t
>  >  2.6|2|26|ST_LineString|t
>  >  2.6|2|27|ST_LineString|t
>  > }}}
>  >
>  > We seem to be closer...
> 
>  Can you try the following on the 32-bit machine?  The query is equivalent
>  to the raster and geometry failing in the regression.
> 
>  {{{
>  SELECT ST_Intersects(
>          ST_AddBand(
>                  ST_MakeEmptyRaster(3, 3, 0, 0, 1, 1, 0, 0, 0),
>                  1, '8BUI', 1, 0
>          ),
>          ST_MakeLine(ARRAY[
>                  ST_MakePoint(-1.1, 1.1),
>                  ST_MakePoint(1.1, 1.1),
>                  ST_MakePoint(1.1, -1.1),
>                  ST_MakePoint(-1.1, -1.1),
>                  ST_MakePoint(-1.1, 1.1)
>          ])
>  , 1);
>  }}}
> 
>  The query should result in TRUE.

strk=#  SELECT ST_Intersects(
strk(#          ST_AddBand(
strk(#                  ST_MakeEmptyRaster(3, 3, 0, 0, 1, 1, 0, 0, 0),
strk(#                  1, '8BUI', 1, 0
strk(#          ),
strk(#          ST_MakeLine(ARRAY[
strk(#                  ST_MakePoint(-1.1, 1.1),
strk(#                  ST_MakePoint(1.1, 1.1),
strk(#                  ST_MakePoint(1.1, -1.1),
strk(#                  ST_MakePoint(-1.1, -1.1),
strk(#                  ST_MakePoint(-1.1, 1.1)
strk(#          ])
strk(#  , 1);
 st_intersects
---------------
 f
(1 row)

Sorry for not using trac, but firing up a web browser is really time
consuming for me. I really hope the newer trac will have mail-to-ticket
gateway (github has one).

--strk;

  ()   Free GIS & Flash consultant/developer
  /\   http://strk.keybit.net/services.html