[postgis-users] st_union(rast) returning null because of non alignment,

Rhys A.D. Stewart rhys.stewart at gmail.com
Wed Nov 7 15:12:06 PST 2012


ok, well that may be the issue, I was trying to union a rasterised table of
adjacent polygons. So will the corners may touch I doubt they overlap, but
here is the output of st_metadata for two of the polygons that are adjacent:

Clarendon 698486.043582936 617384.893096908 153 229 250 250 0 0 3448 4
Manchester 683177.872077352 632613.85425893 117 182 250 250 0 0 3448 4



On Wed, Nov 7, 2012 at 6:01 PM, Pierre Racine
<Pierre.Racine at sbf.ulaval.ca>wrote:

> By "on the same grid" I mean " at least one of any of the four corners of
> any pixel of one raster falls on any corner of the grid of the other
> raster" as stated in the ST_SameAlignment() doc.
>
> What is (are) the pixel size?
>
> > -----Original Message-----
> > From: postgis-users-bounces at lists.osgeo.org [mailto:postgis-users-
> > bounces at lists.osgeo.org] On Behalf Of Rhys A.D. Stewart
> > Sent: Wednesday, November 07, 2012 5:51 PM
> > To: PostGIS Users Discussion
> > Subject: Re: [postgis-users] st_union(rast) returning null because of non
> > alignment,
> >
> > Pierre,
> >
> > If by on the same grid you mean equal, then no hey are not. Here they
> are.
> >
> >
> > 1: POLYGON((766242.610093044 641994.463110607,779742.610093044
> > 641994.463110607,779742.610093044 648494.463110607,766242.610093044
> > 648494.463110607,766242.610093044 641994.463110607))
> > 2:POLYGON((698486.043582936 617384.893096908,736736.043582936
> > 617384.893096908,736736.043582936 674634.893096908,698486.043582936
> > 674634.893096908,698486.043582936 617384.893096908))
> >
> > Thanks,
> >
> > Rhys
> >
> >
> > On Wed, Nov 7, 2012 at 5:40 PM, Pierre Racine <
> Pierre.Racine at sbf.ulaval.ca>
> > wrote:
> >
> >
> >       Hi Rhys,
> >
> >       Both rasters must also have their upperlefx and upperlefty on the
> same
> > grid.
> >
> >       Could you do:
> >
> >       SELECT ST_AsText(ST_ConvexHull(rast))
> >
> >       on both rasters and post the result here?
> >
> >       Thanks,
> >
> >       Pierre
> >
> >
> >       > -----Original Message-----
> >       > From: postgis-users-bounces at lists.osgeo.org [mailto:
> postgis-users-
> >       > bounces at lists.osgeo.org] On Behalf Of Rhys A.D. Stewart
> >       > Sent: Wednesday, November 07, 2012 5:33 PM
> >       > To: PostGIS Users Discussion
> >       > Subject: [postgis-users] st_union(rast) returning null because
> of non
> > alignment,
> >       >
> >       > Greetings,
> >       >
> >       > according to the trusty docs online two rasters are aligned when
> they
> > share the
> >       > same skew, scale and srid. Created a raster using st_asraster
>  and
> > according to
> >       > st_metadata scale, skew and srid are all equal. Howvere st_unon
> > complains
> >       > saying that 'The two rasters provided do not have the same
> alignment.
> >       > Returning NULL'
> >       >
> >       > I noticed that  ticket #1336 as possibly relating to my issue,
> but I'm not
> > entirely
> >       > sure. My posgis version is:
> >       > POSTGIS="2.0.0 r9605" GEOS="3.3.3-CAPI-1.7.4" PROJ="Rel. 4.8.0, 6
> > March
> >       > 2012" GDAL="GDAL 1.9.0, released 2011/12/29" LIBXML="2.7.8"
> >       > LIBJSON="UNKNOWN" TOPOLOGY RASTER
> >       >
> >       >
> >       > What am I missing or doing incorrectly?
> >       >
> >       > Regards,
> >       >
> >       > Rhys
> >
> >       _______________________________________________
> >       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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20121107/af707a35/attachment.html>


More information about the postgis-users mailing list