{"id":8723,"date":"2020-04-06T07:30:05","date_gmt":"2020-04-06T10:30:05","guid":{"rendered":"http:\/\/www.fernandoquadro.com.br\/html\/?p=8723"},"modified":"2020-04-06T09:29:12","modified_gmt":"2020-04-06T12:29:12","slug":"postgis-3-extensao-raster","status":"publish","type":"post","link":"https:\/\/www.fernandoquadro.com.br\/html\/2020\/04\/06\/postgis-3-extensao-raster\/","title":{"rendered":"PostGIS 3: Extens\u00e3o Raster"},"content":{"rendered":"

A funcionalidade Raster no PostGIS faz parte da extens\u00e3o principal desde que foi introduzida. Com o lan\u00e7amento do PostGIS 3, se voc\u00ea quiser a funcionalidade Raster, voc\u00ea precisar\u00e1 instalar a extens\u00e3o principal (postgis) e tamb\u00e9m a extens\u00e3o postgis_raster. Veja:<\/p>\n

\r\nCREATE EXTENSION postgis;\r\nCREATE EXTENSION postgis_raster;\r\n<\/pre>\n
A separa\u00e7\u00e3o da funcionalidade Raster permite que se crie os pacotes do PostGIS de uma forma simplificada, sem a necessidade das depend\u00eancias da extens\u00e3o Raster, que incluem a biblioteca GDAL, um tanto quanto pesada. A funcionalidade Raster permanece intacta, no entanto, e voc\u00ea ainda pode fazer coisas bacanas com ela.<\/p>\n
Por exemplo, fa\u00e7a o download e importe um arquivo de “modelo de eleva\u00e7\u00e3o digital”. No meu caso, vou utilizar um arquivo da cidade de Vancouver.<\/p>\n
\"\"<\/center><\/p>\n
Para baixar essa imagem eu fiz o seguinte:<\/p>\n
\r\ncreatedb yvr_raster\r\npsql -c 'CREATE EXTENSION postgis' yvr_raster\r\npsql -c 'CREATE EXTENSION postgis_raster' yvr_raster\r\n<\/pre>\n
1. Criei uma base de dados e instalei nela as extens\u00f5es postgis e postgis_raster<\/p>\n
\r\nwget https:\/\/pub.data.gov.bc.ca\/datasets\/175624\/92g\/092g06_e.dem.zip\r\nunzip 092g06_e.dem.zip\r\n<\/pre>\n
2. Baixei os dados de DEM no site https:\/\/pub.data.gov.bc.ca\/datasets\/175624\/<\/p>\n
\r\nraster2pgsql -I -F -s 4269 -t 56x56 092g06_e.dem dem092g06e | psql yvr_raster\r\n<\/pre>\n
3. Carreguei os dados no banco.<\/p>\n
Ap\u00f3s o carregamento dos dados, temos uma tabela com o nome dem092g06e de eleva\u00e7\u00e3o de 56 por 56 pixels. Se voc\u00ea mapear a extens\u00e3o, eles se parecer\u00e3o com isso:<\/p>\n

\n\"\"
\n<\/center><\/p>\n
Imagine uma eleva\u00e7\u00e3o de 30 metros no n\u00edvel do solo (em um caso extremo, a Groenl\u00e2ndia e a Ant\u00e1rtica seriam 68 metros). Quanto de Vancouver ficaria debaixo d’\u00e1gua? \u00c9 principalmente um lugar montanhoso. Vamos descobrir:<\/p>\n
\r\nCREATE TABLE poly_30 AS \r\n  SELECT (\r\n   ST_DumpAsPolygons(\r\n    ST_SetBandNoDataValue(\r\n     ST_Reclass(\r\n      ST_Union(rast), \r\n      '0-30:1-1, 31-5000:0-0', '2BUI'),\r\n     0))).*\r\nFROM dem092g06e d\r\n<\/pre>\n
Existem muitas fun\u00e7\u00f5es aninhadas, portanto, lendo de dentro para fora, n\u00f3s temos:<\/p>\n
1. Uni\u00e3o de todas os chips em uma grande imagem Raster
\n2. Reclassifica\u00e7\u00e3o de todos os valores de 0 a 30 para 1 e todos os valores mais altos para 0
\n3. Defini\u00e7\u00e3o do valor “nodata” como 0, n\u00e3o nos importamos com coisas que est\u00e3o acima do nosso limite
\n4. Cria\u00e7\u00e3o de um pol\u00edgono vetorial para cada valor no raster (existe apenas um valor: “1”)<\/p>\n
O resultado fica assim:<\/p>\n

\n\"\"
\n<\/center><\/p>\n
Podemos pegar os dados de constru\u00e7\u00e3o em Vancouver e verificar quantos edif\u00edcios ficar\u00e3o embaixo da \u00e1gua:<\/p>\n
\r\nwget ftp:\/\/webftp.vancouver.ca\/OpenData\/shape\/building_footprints_2009_shp.zip\r\nunzip building_footprints_2009_shp.zip\r\nshp2pgsql -I -s 26910 -D building_footprints_2009 buildings | psql yvr_raster\r\n<\/pre>\n
Antes comparar os edif\u00edcios com a zona de inunda\u00e7\u00e3o, precisamos coloc\u00e1-los na mesma proje\u00e7\u00e3o que a zona de inunda\u00e7\u00e3o (SRID 4269).<\/p>\n
\r\nALTER TABLE buildings\r\nALTER COLUMN geom \r\nTYPE geometry(MultiPolygonZ, 4269)\r\nUSING ST_Transform(geom, 4269)\r\n<\/pre>\n
\"\"<\/center><\/p>\n
Agora podemos encontrar edif\u00edcios inundados:<\/p>\n
\r\nCREATE TABLE buildings_30_poly AS\r\n  SELECT b.* \r\n    FROM buildings b\r\n    JOIN poly_30 p\r\n      ON ST_Intersects(p.geom, b.geom)\r\n    WHERE ST_Intersects(p.geom, ST_Centroid(b.geom))\r\n<\/pre>\n
Existe outra maneira de encontrar edif\u00edcios abaixo de 30 metros, sem ter que construir um pol\u00edgono, assim:<\/p>\n
\r\nCREATE TABLE buildings_30_rast AS\r\n  SELECT b.*\r\n    FROM buildings b\r\n    JOIN dem092g06e d\r\n      ON ST_Intersects(b.geom, d.rast)\r\n    WHERE ST_Value(d.rast, ST_Centroid(b.geom)) < 30;\r\n<\/pre>\n
Como a constru\u00e7\u00e3o de pol\u00edgonos a partir da Raster pode ser cara, unir as camadas raster e a vetorial \u00e9 geralmente a maneira que queremos realizar essa an\u00e1lise.<\/p>\n

\n\"\"
\n<\/center><\/p>\n
Para trazer dados cont\u00ednuos (eleva\u00e7\u00f5es, temperaturas, resultados do modelo) para a an\u00e1lise GIS, a extens\u00e3o postgis_raster pode ser bastante \u00fatil.<\/p>\n
Fonte: Clever Elephant Blog<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"
A funcionalidade Raster no PostGIS faz parte da extens\u00e3o principal desde que foi introduzida. Com o lan\u00e7amento do PostGIS 3, se voc\u00ea quiser a funcionalidade Raster, voc\u00ea precisar\u00e1 instalar a extens\u00e3o principal (postgis) e tamb\u00e9m a extens\u00e3o postgis_raster. Veja: CREATE… Continue Reading →<\/a><\/p>\n","protected":false},"author":275,"featured_media":8731,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[24],"tags":[212,187],"class_list":["post-8723","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-gis","tag-postgis","tag-raster"],"_links":{"self":[{"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/posts\/8723","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/users\/275"}],"replies":[{"embeddable":true,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/comments?post=8723"}],"version-history":[{"count":3,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/posts\/8723\/revisions"}],"predecessor-version":[{"id":8801,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/posts\/8723\/revisions\/8801"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/media\/8731"}],"wp:attachment":[{"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/media?parent=8723"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/categories?post=8723"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/tags?post=8723"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}