Prof. Marcello Benigno: Análise Espacial com o PostGIS e o GvSIG

Uma das inúmeras vantagens de se trabalhar com o PostGIS é a possibilidade de criar views das análises espaciais, tornando os resultados das análises dinâmicas em função dos dados de entrada utilizados. Desta forma também evita-se a criação de novo layers em função de uma atualização da base.

Para exemplificar, vamos imaginar a seguinte situação: em uma determinada cidade foi definido que a área de influência de uma parada de ônibus é de 40 metros e que a prefeitura deve identificar quais são os lotes dentro de cada área de influência das paradas. Apesar da situação ser muito comum para quem já trabalha com SIG, surge a necessidade de se repetir a análise a cada nova parada inserida ou excluída no sistema, o que pode acabar gerando inconsistências na base em função do mal gerenciamento destes arquivos.

Com o uso do PostGIS e do GvSIG eu vou mostrar como resolver este problema sem a necessidade de ficar criando novas camadas a partir das análises desse problema, no caso, a criação do buffer de 40 metros e a consulta espacial (intersects).

No PgAdmin 3, crie um novo Banco de Dados Espacial com o nome cadastro (Figura 1), em seguida, baixe, descompacte e carregue no banco criado o arquivo SQL do link abaixo (Figura 2):

http://acessolivreconsultoria.com/curso_mapserver_gvsig_postgis/lote.zip

Figura 1 - Criação do Banco de Dados Espacial

Figura 2 - Carregando o Script SQL da tabela Lote

Agora vamos criar a tabela ponto_onibus a partir do script abaixo (Figura 3):

CREATE TABLE ponto_onibus(
 id serial PRIMARY KEY,
 observacao varchar(50)
 );

SELECT AddGeometryColumn (
 'public', -- esquema
 'ponto_onibus', -- nome da tabela
 'the_geom', -- nome da coluna
 '29185', -- SRID
 'POINT', -- tipo de geometria
 2 -- plano (2d)
 );

Figura 3 - Criação da tabela ponto_onibus no PgAdmin

Vamos criar duas views, uma para os buffers de 40 metros em torno das paradas de ônibus e a outra dos os lotes que estão dentro desta distância (função intersects).

-- Criando a view buffer_ponto_onibus
CREATE OR REPLACE VIEW buffer_ponto_onibus AS
SELECT id, buffer(the_geom, 40) AS the_geom 
FROM ponto_onibus;

-- Registrando a view na tabela de metadados geometry_columns
INSERT INTO geometry_columns(f_table_catalog, f_table_schema, f_table_name, f_geometry_column, coord_dimension, srid, "type")
VALUES ('', 'public', 'buffer_ponto_onibus', 'the_geom', 2, 29185, 'POLYGON');

-- Criando a view lotes_ponto_onibus
CREATE OR REPLACE VIEW lotes_ponto_onibus AS
SELECT l.*
FROM lote l,  buffer_ponto_onibus b
WHERE intersects(l.the_geom, b.the_geom);

-- Registrando a view na tabela de metadados geometry_columns
INSERT INTO geometry_columns(f_table_catalog, f_table_schema, f_table_name, f_geometry_column, coord_dimension, srid, "type")
VALUES ('', 'public', 'lotes_ponto_onibus', 'the_geom', 2, 29185, 'POLYGON');

Abra o GvSIG e mude a projeção do bloco para 29185 (Figura 4)

Figura 4 - Troca do sistema de projeção no GvSIG

Com o bloco aberto, clique no botão adicionar plano de informação e escolha a segunda aba (GeoDB), em seguida configure a conexão com o banco de dados "cadastro" (Figura 5).

Figura 5 - Configuração da conexão com o Banco de Dados

Configurada a conexão, adicione as tabelas "lote" e "parada_onibus" e as views "buffer_ponto_onibus" e "lotes_ponto_onibus" (Figura 6).

Figura 6 - Adicionando as camadas do BD

Inicie a edição da camada ponto_onibus clicando com o botão direito em cima desta camada e escolhendo a opção "Iniciar Edição" em seguida, adicione algumas paradas de forma aleatória. Termine a edição Salvando as alterações (Figura 7).