Análises Espaciais com o PostGIS Raster - Parte 1

O projeto PostGIS Raster possibilita o armazenamento e manipulação de dados matriciais em Banco de Dados Espaciais através da linguagem SQL.

Nesta série de tutoriais, serão demonstradas algumas funcionalidades desta extensão através da utilização dos dados de uso e ocupação do solo da mesorregião de Aracicaba (SP) e adjacências. Os dados podem ser obtidos através deste link: http://ge.tt/1Hye7zJ2/v/0?c

Caso você ainda não tenha o PostgreSQL e o PostGIS instalados no Ubuntu, este post descreve o processo: Instalação do PostGIS 2.0 no Ubuntu via Terminal

1) Criação do Banco de Dados Espaciais

Com o terminal aberto, entre com o comando:

$ createdb sp

Uma vez criado o Banco de Dados, deve-se então habilitar a extensão PostGIS no mesmo:

$ psql sp
psql (9.4.4)
Type "help" for help.

postgis=# CREATE EXTENSION postgis;

Digite \q para sair do Banco de Dados.

2) Conversão do arquivo Raster para SQL

Para realizar a conversão, abra o terminal e navegue até o diretório onde o arquivo uso_ocupacao.tif se encontra, em seguida, execute o seguinte comando:

$  raster2pgsql -s 4326 -I -M -C uso_ocupacao.tif -F -t 30x30 uso_ocupacao_solo > uso_ocupacao_solo.sql

Cada flag, que são as letras iniciadas com um "-" determina uma propriedade no arquivo sql que é gerado, neste caso em específico, temos:

• -s : sistema de referência espacial do raster (SRID 4326)
• -I: cria um índice GiST na coluna do raster
• -M: adiciona um "Vacuum analise" ao script
• -C: cria as constraints da tabela raster
• -F: adiciona uma coluna com o nome do arquivo
• -t 30×30: resolução espacial do raster

Outras opções sobre as flags e sobre o raster2pgsql  podem ser encontradas no link:

http://postgis.refractions.net/docs/using_raster.xml.html#RT_Loading_Rasters

Agora é possível carregar todos os arquivos no Banco de Dados (municipios.sql, hidrografia.sql e uso_ocupacao_solo.sql) através de um único comando:

$  for f in *.sql; do psql -f $f -d sp; done

3) Visualização das tabelas no QGIS:

Uma vez configurada a conexão entre o QGIS e o Banco PostGIS, através do plugin DB Manager podemos então visualizar as tabelas do tipo raster, como mostra a Figura 1:

Figura 1 - Com o botão direito do mouse é possível escolher a opção "Add to canvas"

Para melhor representação visual da tabela, clique com o botão direito no layer e em "Properties -> Style" realize o ajuste exibido na Figura 2:

Figura 2 - Ajuste das cores do mapa temático

A Figura 3 mostra como deve ficar o resultado no QGIS.

Figura 3 - Mapa de uso e ocupação do solo

No próximo post demonstrarei algumas consultas e análises espaciais úteis utilizando os dados desta tabela, através da linguagem SQL... até breve :)