<div dir="ltr">Hi Pierre,<div>Thanks for the link to the add-on and the sql statements.</div><div><br></div><div>Karl<br><br>Den onsdag den 27. november 2013 14.21.46 UTC-3 skrev Pierre Racine:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Hi Karl,
<br>
<br>I do this with the ST_ExtractToRaster() function from the PostGIS Add-ons:
<br>
<br><a href="https://github.com/pedrogit/postgisaddons/releases/tag/1.21" target="_blank" onmousedown="this.href='https://www.google.com/url?q\75https%3A%2F%2Fgithub.com%2Fpedrogit%2Fpostgisaddons%2Freleases%2Ftag%2F1.21\46sa\75D\46sntz\0751\46usg\75AFQjCNER-uAEBczQbEUCrVBoeK5luCnYrg';return true;" onclick="this.href='https://www.google.com/url?q\75https%3A%2F%2Fgithub.com%2Fpedrogit%2Fpostgisaddons%2Freleases%2Ftag%2F1.21\46sa\75D\46sntz\0751\46usg\75AFQjCNER-uAEBczQbEUCrVBoeK5luCnYrg';return true;">https://github.com/pedrogit/<wbr>postgisaddons/releases/tag/1.<wbr>21</a>
<br>
<br>I just added a SUM_OF_LENGTHS method to ST_ExtractToRaster().
<br>
<br>Look at the ST_ExtractToRaster() function in the postgis_addons.sql file for an example. In your case it should look like this:
<br>
<br>SELECT ST_ExtractToRaster(ST_AddBand(<wbr>ST_MakeEmptyRaster(rast), '32BF'), 'public', 'geomtable', 'geom', null, 'SUM_OF_LENGTHS') rast
<br>FROM refrastertable;
<br>
<br>You should have a reference raster already existing in a table and your lines should be indexed.
<br>
<br>Pierre
<br>
<br>> -----Original Message-----
<br>> From: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="8TLgZGj90-MJ" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">postgis-us...@lists.<wbr>osgeo.org</a> [mailto:<a href="javascript:" target="_blank" gdf-obfuscated-mailto="8TLgZGj90-MJ" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">postgis-users-</a>
<br>> <a href="javascript:" target="_blank" gdf-obfuscated-mailto="8TLgZGj90-MJ" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">bou...@lists.osgeo.org</a>] On Behalf Of Karl Zinglersen
<br>> Sent: Wednesday, November 27, 2013 11:30 AM
<br>> To: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="8TLgZGj90-MJ" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">postgi...@googlegroups.com</a>
<br>> Cc: PostGIS Users Discussion
<br>> Subject: Re: [postgis-users] Simple Line Density
<br>> 
<br>> Hi,
<br>> I am dealing with the same types of data and issues as Jeff.
<br>> My procedure for data clean up is:
<br>> 1) CSV files into PostGIS as points
<br>> 1.1) ST_TRANSFORM to projected spatial reference (here EPSG:32621 /
<br>> UTM zone 21 N WGS84)
<br>> 2) ST_MAKELINE to "sub"-lines via PARTITION BY
<br>> 3) Create a gridded polygon layer in same srid
<br>> 4) Run equivalent to select cell_id, line_id from cells, lines where
<br>> intersects(lines.geom, cells.geom
<br>> - but my postgis runs out of memory after a while. And I am sorry to say
<br>> Spatial Analyst doesn't.
<br>> Is there more memory cheap way to do it in PostGIS (e.g. via the raster
<br>> functions - although I haven't found i yet).
<br>> 
<br>> Karl
<br>> 
<br>> Den onsdag den 26. december 2012 09.58.43 UTC-3 skrev Jeff Adams -
<br>> NOAA Affiliate:
<br>> 
<br>>         Hi Brent,
<br>> 
<br>>         Thanks for the response. I don't think your original response made
<br>> it into my inbox via the list, but now I see it on the website thread. That
<br>> seems like it might be a viable alternative, I just worry about those vessel
<br>> tracks that would be clearly outside of a particular grid cell, but whose
<br>> buffer would bleed over. How did you handle this type of situation?
<br>> 
<br>>         Jeff
<br>> 
<br>> 
<br>>         On Mon, Dec 24, 2012 at 7:22 PM, <<a>pcr...@pcreso.com</a>
<br>> <javascript:> > wrote:
<br>> 
<br>> 
<br>> Jeff,
<br>> 
<br>> Did you see my reply using vessel tracklines as vectors & a grid in Postgis to
<br>> do exactly what you describe?
<br>> 
<br>> Given we were looking at benthic impact, we buffered the tracklines to
<br>> create polygons representing the swept area of the deployed fishing gear.
<br>> These were clipped by the cells, & we could generate statistics suca as the
<br>> cumulative swept area of all tracks with each cell, number of times each cell
<br>> was crossed, & given the tracklines have a timestamp associated with them,
<br>> we could also look at the temporal pattern of tracks crossing cells, for
<br>> things like seasonal impacts & variation between seasons.
<br>> 
<br>> Cheers,
<br>> 
<br>>   Brent Wood
<br>> 
<br>
<br>______________________________<wbr>_________________
<br>postgis-users mailing list
<br><a href="javascript:" target="_blank" gdf-obfuscated-mailto="8TLgZGj90-MJ" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">postgi...@lists.osgeo.org</a>
<br><a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users" target="_blank" onmousedown="this.href='http://www.google.com/url?q\75http%3A%2F%2Flists.osgeo.org%2Fcgi-bin%2Fmailman%2Flistinfo%2Fpostgis-users\46sa\75D\46sntz\0751\46usg\75AFQjCNFq7nnxngR7sOb6icF73DL5HQVrbA';return true;" onclick="this.href='http://www.google.com/url?q\75http%3A%2F%2Flists.osgeo.org%2Fcgi-bin%2Fmailman%2Flistinfo%2Fpostgis-users\46sa\75D\46sntz\0751\46usg\75AFQjCNFq7nnxngR7sOb6icF73DL5HQVrbA';return true;">http://lists.osgeo.org/cgi-<wbr>bin/mailman/listinfo/postgis-<wbr>users</a>
<br></blockquote></div></div>