<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN">
<html><body>
<p>Thanks, Pierre,</p>
<p>It seems ok to model the query the way you showed, however, that implies the generation of an entire intermediate binary raster, just to identify the region inside the whole mapped area. My issue is on avoiding such intermediate data by incorporating map algebra expression into the same query expression, Another example, ""... d.geom touch (ndvi2018 < ndvi2019)" where ndvixxxx are raster data.</p>
<p>Att.</p>
<p>Em 15-07-2019 14:43, Pierre Racine escreveu:</p>
<blockquote type="cite" style="padding-left:5px; border-left:#1010ff 2px solid; margin-left:5px; width:100%"><!-- html ignored --><!-- head ignored --><!-- meta ignored --><!-- meta ignored -->
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;">Hi,<!-- o ignored --></span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;"><!-- o ignored --> </span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;">That could translate into something like this in PostGIS:<!-- o ignored --></span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;"><!-- o ignored --> </span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;">WITH floodareas AS (<!-- o ignored --></span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;"> SELECT ST_MapAlgebra(rast, ‘1BB’, ‘CASE WHEN [rast] < 10 THEN 1 ELSE NULL END’) rast<!-- o ignored --></span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;"> FROM elevation<!-- o ignored --></span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;">)<!-- o ignored --></span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;">SELECT d.*<!-- o ignored --></span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;">FROM districts d, floodareas f<!-- o ignored --></span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;">WHERE ST_Intersects(g.geom, f.rast)<!-- o ignored --></span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;"><!-- o ignored --> </span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;">NULL becomes nodata values. You could also use ST_Reclass() instead of ST_Mapalgebra().<!-- o ignored --></span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;"><!-- o ignored --> </span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;">Pierre<!-- o ignored --></span></p>
<p class="MsoNormal"><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: #1f497d; mso-fareast-language: EN-US;"><!-- o ignored --> </span></p>
<div>
<div style="border: none; border-top: solid #E1E1E1 1.0pt; padding: 3.0pt 0cm 0cm 0cm;">
<p class="MsoNormal"><strong><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: windowtext;">De :</span></strong><span style="font-size: 11.0pt; font-family: 'Calibri',sans-serif; color: windowtext;"> postgis-users <postgis-users-bounces@lists.osgeo.org> <strong>De la part de</strong> Joao Cordeiro<br /><strong>Envoyé :</strong> 15 juillet 2019 10:19<br /><strong>À :</strong> 'PostGIS Users Discussion' <postgis-users@lists.osgeo.org><br /><strong>Objet :</strong> [postgis-users] raster regions<!-- o ignored --></span></p>
</div>
</div>
<p class="MsoNormal"><!-- o ignored --> </p>
<p>Hi there,<!-- o ignored --></p>
<p>I am not actually a POSTGIS user, or developer, although I am doing some research of integration of vector -raster in a dbms context. My main question s about identifying regions represented in a domain like the "hybrid raster" so that they can be referred to in a query. For example, I would like to write something like "<strong><em>select d from districts where d.geometry touches flooding</em></strong>", where "flooding" indicates an expression intended to decide on the pertinence of locations to a flooding risk area described, for instance through map algebra.<!-- o ignored --></p>
<p>Thanks for any hints or suggestions.<!-- o ignored --></p>
<p>Att.<!-- o ignored --></p>
</div>
<!-- html ignored --><br />
<pre>_______________________________________________
postgis-users mailing list
<a href="mailto:postgis-users@lists.osgeo.org">postgis-users@lists.osgeo.org</a>
<a href="https://lists.osgeo.org/mailman/listinfo/postgis-users">https://lists.osgeo.org/mailman/listinfo/postgis-users</a></pre>
</blockquote>
<p> </p>
<div> </div>
</body></html>