Yes, I noticed the huge difference as well when projecting to another coordinate system. So you reckon I don't need to do that and just follow your last mail on <br>how to make an empty raster and use map algebra?<br><br>
<div class="gmail_quote">On Fri, Mar 23, 2012 at 4:11 PM, Pierre Racine <span dir="ltr"><<a href="mailto:Pierre.Racine@sbf.ulaval.ca">Pierre.Racine@sbf.ulaval.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">> Is this ok? Once I have made the raster, I can then transform it back to srid 4326<br>
> before doing PixelAsPolygons right?<br>
<br>
</div>You definitely have to be able to view your grid and your data. So you can quickly make sure both extents are ok.<br>
<br>
This seems to be very big comparing with what I got the other day (1020x24798). You will end up with way too many polygons (361932 x 423152 = 153152249664)... Even if you cut by two to remove the ones that fall in the water. You will have to do everything in raster space (the ST_MapAlgebraFct technique) and that would be too slow.<br>

<br>
I think the best approach in this case would go like this:<br>
<br>
1) Tile your grid (smaller the better)<br>
<br>
2) Write a function burning a rasterization of each road in the right tile using ST_MapAlgebra(rast, ST_AsRaster(roadgeom)) taking the actual value of the tile into account<br>
<br>
That might work fast enough.<br>
<span class="HOEnZb"><font color="#888888"><br>
Pierre<br>
</font></span></blockquote></div><br>