\u00c9 muito claro que a informa\u00e7\u00e3o que estamos procurando \u00e9 esta “NAD_1983_StatePlane_Oregon_South_FIPS_3602_Feet_Intl”, mas qual n\u00famero “SRID” devemos usar?<\/p>\n
Para descobrir basta ir para ao site Prj2EPSG (http:\/\/prj2epsg.org<\/a>), e procurar pela informa\u00e7\u00e3o acima, e voc\u00ea vai obter a resposta que \u00e9 2270.<\/p>\n Agora, usando o shp2pgsql podemos carregar os dados em uma tabela chamada buildings<\/em>:<\/p>\n Os nossos dados LIDAR est\u00e3o todos em coordenadas geogr\u00e1ficas (EPSG: 4326) e n\u00f3s vamos ter de integr\u00e1-los com a camada buildings (EPSG:2270) para envitar problemas, sendo assim vamos utilizar a fun\u00e7\u00e3o ST_Transform do PostGIS pra isso:<\/p>\n Para agilizar a nossa an\u00e1lise, vamos eliminar todos os edif\u00edcios que n\u00e3o est\u00e3o contidos na \u00e1rea do nosso dado LiDAR.<\/p>\n Agora, \u00e9 publicar esta camada no GeoServer. Fa\u00e7a a publica\u00e7\u00e3o da forma tradicional, apenas na aba Cache desmarque a op\u00e7\u00e3o “Criar camada de cache para esta camada” e na aba Publishing altere as informa\u00e7\u00f5es das propriedade do KML confirme a figura abaixo:<\/p>\n Salve, e ent\u00e3o temos agora uma camada vis\u00edvel! Voc\u00ea pode visualiz\u00e1-la no Google Earth usando o KML reflector do GeoServer<\/a>.<\/p>\n http:\/\/localhost:8080\/geoserver\/wms\/kml?layers=opengeo:buildings&format_options=lookatbbox:bbox=-122.8808,42.3311,-122.8806,42.3313&mode=refresh&kmscore=50<\/a><\/p>\n Quando voc\u00ea ampliar a imagem no Google Earth, vai notar algo estranho sobre nossos edif\u00edcios: eles s\u00e3o “lisos”! Queremos edif\u00edcios em 3D, como podemos obt\u00ea-los? \u00c9 isso que iremos ver nos pr\u00f3ximos posts.<\/p>\n","protected":false},"excerpt":{"rendered":" Agora vamos adicionar ao mapa uma camada de edif\u00edcios. Para isso descompacte o arquivo BuildingFootprints.zip que baixado anteriormente. N\u00f3s vamos carregar os edif\u00edcios no banco de dados como uma tabela espacial, mas antes disso, n\u00f3s temos de descobrir o qual… Continue Reading →<\/a><\/p>\n","protected":false},"author":275,"featured_media":6072,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[24],"tags":[268],"class_list":["post-6070","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-gis","tag-lidar"],"_links":{"self":[{"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/posts\/6070","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=6070"}],"version-history":[{"count":9,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/posts\/6070\/revisions"}],"predecessor-version":[{"id":6205,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/posts\/6070\/revisions\/6205"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/media\/6072"}],"wp:attachment":[{"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/media?parent=6070"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/categories?post=6070"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.fernandoquadro.com.br\/html\/wp-json\/wp\/v2\/tags?post=6070"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}shp2pgsql -s 2270 -D BuildingFootprints.shp buildings | psql -d lidar<\/pre>\n
-- Update SRID and transform all geoms\r\nALTER TABLE buildings\r\nALTER COLUMN geom\r\nTYPE geometry(MultiPolygon,4326)\r\nUSING ST_Transform(geom, 4326);\r\n-- Rename to area column to area\r\nALTER TABLE buildings\r\nRENAME COLUMN shape_st_1\r\nTO shape_area;\r\n-- Index the table\r\nCREATE INDEX buildings_gix ON buildings USING GIST (geom);<\/pre>\n
-- Find the LIDAR extent\r\nSELECT st_extent(pa::geometry) FROM medford;\r\n-- BOX(-122.8874999 42.3125,-122.8749998 42.325)\r\n-- Delete unneeded building polygons\r\nDELETE FROM buildings\r\nWHERE NOT ST_Contains(\r\n ST_MakeEnvelope(-122.8874999, 42.3125, -122.8749998, 42.325, 4326),\r\n geom);<\/pre>\n
![\"gs_kmlconfig\"](\"http:\/\/www.fernandoquadro.com.br\/html\/wp-content\/uploads\/2016\/07\/gs_kmlconfig.jpg\")
<\/p>\n