From GeoBox

Importação de dados

Este é um exercício típico de importação de dados.

Cenário: Preciso de saber o que se posso fazer com um terreno que herdei de uma tia avó que não tenho, que viveu em Mangualde.

A minha tia avó deu-me as coordenadas do terreno :-)

st(st_geomfromtext('POLYGON((36750 104000,37990 103900,38000 103500,37000 103000,36750 104000))', 3763), wkb_geometry)

"POLYGON((36750 104000,37990 103900,38000 103500,37000 103000,36750 104000))"

select *, intersection(st_geomfromtext('POLYGON((36750 104000,37990 103900,38000 103500,37000 103000,36750 104000))',3763), wkb_geometry) from reserva where intersects(st_geomfromtext('POLYGON((36750 104000,37990 103900,38000 103500,37000 103000,36750 104000))',3763), wkb_geometry);

create table resultado as select ogc_fid as id, name as nome, description as descricao, intersection(st_geomfromtext('POLYGON((36750 104000,37990 103900,38000 103500,37000 103000,36750 104000))',3763), wkb_geometry) as intercepcao from reserva where intersects(st_geomfromtext('POLYGON((36750 104000,37990 103900,38000 103500,37000 103000,36750 104000))',3763), wkb_geometry);

select id, nome, descricao, st_area(intercepcao) from resultado; http://www.cmmangualde.pt/index.php?option=com_content&view=article&id=271&Itemid=278

wget http://www.cmmangualde.pt/images/stories/pdm_online/carta%20da%20reserva%20ecolgica%20nacional%20-%20vectorizao.rar --2010-11-04 22:54:36-- http://www.cmmangualde.pt/images/stories/pdm_online/carta%20da%20reserva%20ecolgica%20nacional%20-%20vectorizao.rar

unrar x carta\ da\ reserva\ ecolgica\ nacional\ -\ vectorizao.rar

mv Carta\ da\ Reserva\ Ecológica\ Nacional\ -\ Vectorização.kml reserva.kml

xmlstarlet val reserva.kml

recode latin1..utf8 reserva.kml

xmlstarlet val reserva.kml

ogr2ogr -s_srs EPSG:4326 -t_srs EPSG:3763 -f "PostgreSQL" PG:"host=localhost user=geobox dbname=geotuga password=geobox" reserva.kml -nln reserva

Importar um KML

sudo apt-get install xmlstarlet recode xmlcopyeditor

unrar x carta\ de\ ordenamento\ -\ vectorizao.rar mv Carta\ de\ Ordenamento\ -\ Vectorização.kml ordenamento.kml xmlstarlet val ordenamento.kml recode latin1..utf8 ordenamento.kml xmlstarlet val ordenamento.kml

Os dados estão em 4326

ogr2ogr -f "PostgreSQL" PG:"host=localhost user=geobox dbname=geotuga password=geobox" ordenamento.kml ogr2ogr -s_srs EPSG:4326 -t_srs EPSG:3763 -f "PostgreSQL" PG:"host=localhost user=geobox dbname=geotuga password=geobox" ordenamento.kml

Categoria:Portugal

ogr2ogr -f "PostgreSQL" PG:"host=localhost user=geobox dbname=geotuga password=geobox" ordenamento.kml

Último exercício

Tem que haver um plano