<div dir="ltr">Hello,<div><br></div><div>Unfortunatly, I don't find a way to print raster results in OpenJump. Can you explain me how I can do it ? </div><div><br></div><div>Thank you for your help !</div><div dir="ltr"><div><br></div><div>Pierre</div></div><br><div class="gmail_quote"><div dir="ltr">Le lun. 20 juin 2016 à 16:16, Pierre Racine <<a href="mailto:Pierre.Racine@sbf.ulaval.ca" target="_blank">Pierre.Racine@sbf.ulaval.ca</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I would try to understand, visually using OpenJump, how your tiles are clipped in the two cases... I guess sometimes 1 pixel wide tiles get clipped in a strange way.<br>
<br>
Pierre<br>
<br>
> -----Original Message-----<br>
> From: postgis-users [mailto:<a href="mailto:postgis-users-bounces@lists.osgeo.org" target="_blank">postgis-users-bounces@lists.osgeo.org</a>] On Behalf<br>
> Of Pierre Bails<br>
> Sent: Wednesday, June 15, 2016 11:41 AM<br>
> To: <a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a><br>
> Subject: [postgis-users] Raster tile size impact results<br>
><br>
> Hi all,<br>
><br>
> I'm looking for get all pixel's value in a raster insert thanks to raster2pgsql<br>
> intersect by a polygon.<br>
> Nevertheless, there are some things that I don't understand: depends on tile<br>
> size, I don't have the same result... The error between 2 raster is also correlate<br>
> to the pixel size (I assume it's due to the pixel number).<br>
><br>
> Can you explain me how the tile size can influence the final result ?<br>
><br>
> Here is the query :<br>
> SELECT DISTINCT SUM((px).val)<br>
> FROM(<br>
>   SELECT *<br>
>   FROM  mns INNER JOIN<br>
> st_setsrid(st_geomfromgeojson('{"type":"Polygon","coordinates":[[[x,y],[x,y],[x<br>
> ,y] ]]}'), 4326) as geom on ST_Intersects(mns.rast, geom)<br>
> ) AS clip, ST_PixelAsCentroids (ST_Clip(clip.rast, clip.geom),1) AS px;<br>
><br>
><br>
> DSM px size 70cm:<br>
> raster2pgsql -a -t 250x250 -s 4326 -F pathToDSM.tiff dsmTable (result: 46 347<br>
> 846)<br>
><br>
> raster2pgsql -a -t 50x50 -s 4326 -F pathToDSM.tiff dsmTable    : 46 382 873<br>
> Difference: 35027<br>
><br>
><br>
> DSM px size 20cm :<br>
> raster2pgsql -a -t 250x250 -s 4326 -F pathToDSM.tiff dsmTable :  567 505 888<br>
><br>
> raster2pgsql -a -t 50x50 -s 4326 -F pathToDSM.tiff dsmTable    : 596 897 610<br>
> Difference: 29 391 722<br>
><br>
><br>
> If we considere 250x250's result as reference, there is a difference of 0.07% in<br>
> the first case, and 5% in the second. If the polygon area increases, the error rate<br>
> increases.<br>
><br>
> Thank you,<br>
><br>
> Pierre<br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/postgis-users" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/postgis-users</a></blockquote></div></div>