<br>
<br>
<br>
<table width=100%>
<tr valign=top>
<td width=40%><font size=1 face="sans-serif"><b>Pierre Racine <Pierre.Racine@sbf.ulaval.ca></b>
</font>
<br><font size=1 face="sans-serif">Sent by: postgis-users-bounces@postgis.refractions.net</font>
<p><font size=1 face="sans-serif">06/07/2011 07:30 PM</font>
<table border>
<tr valign=top>
<td bgcolor=white>
<div align=center><font size=1 face="sans-serif">Please respond to<br>
PostGIS Users Discussion <postgis-users@postgis.refractions.net></font></div></table>
<br>
<td width=59%>
<table width=100%>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">To</font></div>
<td><font size=1 face="sans-serif">PostGIS Users Discussion <postgis-users@postgis.refractions.net></font>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">cc</font></div>
<td>
<tr valign=top>
<td>
<div align=right><font size=1 face="sans-serif">Subject</font></div>
<td><font size=1 face="sans-serif">Re: [postgis-users] How do you use PostGIS
Raster?</font></table>
<br>
<table>
<tr valign=top>
<td>
<td></table>
<br></table>
<br>
<br>
<br><font size=2 color=#004080 face="Calibri">What are you missing to implement
any of those in PostGIS raster?</font>
<br>
<br><font size=2 color=#004080 face="Calibri">For the first thing: maybe
nothing, but I'd have to check. It would be convenient if there was a "rasterize"
function which would serve as the inverse of a ST_DumpAsPolygons().</font>
<br>
<br><font size=2 color=#004080 face="Calibri">For the second thing: Landsat
scenes are ~7000x8000 pixels; ST_DumpAsPolygons() gives a huge set of geometries.
The MODIS fire perimeters are also huge. The first step is a "global"
(image-wide) analysis. It was very slow when I tried it. We need an ST_Intersection
that returns a clipped raster; and we also need an ST_Union() or ST_Collect()
which returns an aggregated raster.</font>
<br><font size=2 color=#004080 face="Calibri"> </font>
<br><font size=2 face="Tahoma"><b>From:</b> postgis-users-bounces@postgis.refractions.net
[mailto:postgis-users-bounces@postgis.refractions.net] <b>On Behalf Of
</b>Bryce L Nordgren<b><br>
Sent:</b> Tuesday, June 07, 2011 3:12 PM<b><br>
To:</b> PostGIS Users Discussion<b><br>
Subject:</b> Re: [postgis-users] How do you use PostGIS Raster?</font>
<br><font size=3 face="Times New Roman"> </font>
<br><font size=2 face="Arial"><br>
I'm not using PostGIS raster yet. We've been using PostGIS for years to
keep track of active fire and burn scar detections from the MODIS instruments
on board the Aqua and Terra satellites. Part of our process is to resample
the point detections onto a grid which aggregates burned area over the
course of a fire season. Aggregating onto a grid has been performed using
an external program, but we'd like to be able to use PostGIS raster.</font><font size=3 face="Times New Roman">
<br>
</font><font size=2 face="Arial"><br>
More recently, we're comparing a vector fire perimeter from the MODIS instrument
 to a higher resolution Landsat image of the area. We select the areas
of the Landsat raster which are inside the perimeter vs. outside the perimeter
and calculate histograms on the separate image regions. Currently, we're
using ENVI and IDL to do this more or less manually. We'd like to be able
to perform this inside PostGIS, as PostGIS is the source of the vector
data anyway.</font><font size=3 face="Times New Roman"> <br>
</font><font size=2 face="Arial"><br>
Bryce</font><tt><font size=2>_______________________________________________<br>
postgis-users mailing list<br>
postgis-users@postgis.refractions.net<br>
</font></tt><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users"><tt><font size=2>http://postgis.refractions.net/mailman/listinfo/postgis-users<br>
</font></tt></a>
<br>