[postgis-users] Rast index : Could not get raster's convex hull

[Bborie Park]

Hi Hugues,

In looking at the convex hull code under the hood, I think the only
way for you to be getting that error is that at least one record's
raster has a width or height of zero.  So, try something like...

WITH foo AS (
SELECT rid, ST_Metadata(rast) As meta FROM rasters.slopes.alpes
)
SELECT rid, (rast).width, (rast).height FROM foo WHERE (rast).width <
1 OR (rast).height < 1

-bborie

On Thu, Nov 29, 2012 at 12:20 AM, Francois Hugues
<hugues.francois at irstea.fr> wrote:
> Hello,
>
> Here it is ! Sorry I have forgotten it the first time.
>
> "POSTGIS="2.0.1 r9979" GEOS="3.3.5-CAPI-1.7.5" PROJ="Rel. 4.8.0, 6 March 2012" GDAL="GDAL 1.9.1, released 2012/05/15" LIBXML="2.7.8" TOPOLOGY RASTER"
>
> Hugues.
>
>
> --
> Hugues FRANÇOIS
> Ingénieur recherche
>
> Irstea Grenoble
> 2 rue de la papeterie
> 38400 Saint-Martin d'Hères
>
> tel : +33 (0)4.76.76.27.44
> port : +33 (0)6.77.66.21.31
> fax : +33 (0)4.76.51.38.03
>
>
> -----Message d'origine-----
> De : postgis-users-bounces at lists.osgeo.org [mailto:postgis-users-bounces at lists.osgeo.org] De la part de Bborie Park
> Envoyé : jeudi 29 novembre 2012 02:20
> À : postgis-users at lists.osgeo.org
> Objet : Re: [postgis-users] Rast index : Could not get raster's convex hull
>
> Hey Hugues,
>
> What is your output for
>
> SELECT postgis_full_version()
>
> -bborie
>
> On 11/26/2012 08:01 AM, Francois Hugues wrote:
>> Hello,
>>
>> Here is an error I don't understand when I want to index som tiles
>> from a raster table :
>> ERROR:  RASTER_convex_hull: Could not get raster's convex hull
>>
>> Here is the query which does not work:
>> create index slopes_alpes_rast on rasters.slopes_alpes using
>> gist(st_convexhull(rast));
>>
>> Here are queries that does work without any error :
>>
>> 1. simple select give me a geometry for all of my tiles and no one is
>> null:
>>
>> select rid, st_convexhull(rast) from rasters.slopes_alpes
>>
>> 2. simple select + index on newly created table :
>>
>> create table rasters.tmp_index_gist as select rid, st_convexhull(rast)
>> geom from rasters.slopes_alpes create index tmp_idx on
>> rasters.tmp_index_gist using gist(geom)
>>
>> It's weird, isn't it ? Any idea of what's happening ? Where do I make
>> a mistake or what is wrong with my data ?
>>
>> Hugues.
>> _______________________________________________
>> postgis-users mailing list
>> postgis-users at lists.osgeo.org
>> http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users
>>
>
> --
> Bborie Park
> Programmer
> Center for Vectorborne Diseases
> UC Davis
> 530-752-8380
> bkpark at ucdavis.edu
> _______________________________________________
> postgis-users mailing list
> postgis-users at lists.osgeo.org
> http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users
> _______________________________________________
> postgis-users mailing list
> postgis-users at lists.osgeo.org
> http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users



-- 
Bborie Park
Programmer
Center for Vectorborne Diseases
UC Davis
530-752-8380
bkpark at ucdavis.edu