Skip to content

SQL, extract municipio

Peter edited this page Feb 11, 2020 · 9 revisions

Exemplo de RR-BoaVista

Pelo datasets.ok.org.br/city-codes temos confirmação de que o código IBGE é 1400100 e a entrada Wikidata é Q181056, onde fica destacado que o OSM-relation-ID é 326286. Ainda assim, convém no SQL conferir se não existem outros polígonos de boundary:administrative com mesmo código IBGE. Se fossemos recuperar direto pelo OSM-relation-ID da Wikidata, convém lembrar que é negativo na representação osm_id. Outro problema é que queremos o "polígono puro", mas a relation inclui o ponto do admin_centre (falta conferir se (conforme Do zero) o osm2pgsql separou as duas coisas (geometria pura de ST_Polygon).

 SELECT osm_id, ST_GeometryType(way) as geom_type, tags 
 FROM planet_osm_polygon 
 WHERE tags->>'boundary'='administrative' and tags->>'IBGE:GEOCODIGO'='1400100';

A listagem confirma apenas um item e com o perfil esperado.