domingo, 23 de junho de 2019

proj4js (geomatica.no.sapo.pt) revisitado

Transformação pontual de coordenadas entre os principais sistemas de coordenadas utilizados em Portugal Continental. 

Nas transformações que impliquem mudança de datum, são utilizados os 7 parâmetros de Helmert (Bursa-Wolf) calculados para todo o território nacional. 

EPSG:4326 - WGS84

EPSG:20790 - Hayford-Gauss Datum Lisboa (Militar)

EPSG:20791 - Hayford-Gauss Datum Lisboa (IGP)  

EPSG:27493 - Hayford-Gauss Datum 73 (IGP)

EPSG:3763 - ETRS89 / Portugal TM06

EPSG:4207 - Geográficas Datum Lisboa

EPSG:4274: Geográficas Datum 73


Nota: não utilizo o código EPSG: 5018 (Datum Lisboa/ Hayford-Gauss). Prefiro manter o EPSG:20791 o qual aparece logo a seguir ao EPSG:20790. Como só se diferenciam por uma translação devem estar juntos. E,... também prefiro continuar a utilizar IGP em vez de DGTerritório.



For more information on Proj4js and to report bugs, please visit the Proj4js Trac and Wiki at OSGeo.





quarta-feira, 17 de abril de 2019

Mata Atlântica: diferenças entre SIRGAS2000 e SAD69 (parâmetros de transformação)

Experimente clicar no mapa e irão aparecer 3 janelas: SIRGAS2000, SAD69 com transformação de WGS84 (Doppler), e SAD69  com transformação de  WGS 84 (G1150).




Pode voltar a clicar!

São as diferenças encontradas transformando de SIRGAS (considerando igual ao WGS84) para SAD69 com os seguintes parâmetros de transformação, mas com os sinais de +towgs84 opostos:
  
'EPSG:4618': "+title=SAD69 RP2005 +proj=longlat +ellps=aust_SA +towgs84=-67.35,3.88,-38.22 +units=degrees +no_defs"

 'EPSG:4291': "+title=SAD69 RP1989 +proj=longlat +ellps=aust_SA +towgs84=-66.87,4.37,-38.52 +units=degrees +no_defs"



Com os parâmetros de transformação sem modificar os sinais de +towgs84:



quinta-feira, 11 de abril de 2019

quarta-feira, 10 de abril de 2019

KML ao limite

Como as Fusion Tables vão deixar de funcionar em Dez. 2019, é uma boa oportunidade para começar a utilizar mais vezes Vector Tiles no Google Maps.

Mas primeiro vamos fazer umas experiências com kml para ver até onde é que isto aguenta. Normalmente utilizo geoxml3, e o kml-split da Mapbox. Mas neste caso não vou utilizar o kml-split.

Então temos a nossa Fusion Table que estamos a consumir numa FusionTablesLayer . Reparem que estão uns polígonos realçados a Sul, e que se fizerem zoom in aparecem outros a Norte. É neste tipo de queries dinâmicas que Fusion Tables dava jeito, não só, mas também, para definição de estilos.





Na Fusion Table é possível descarregar um ficheiro e vamos obter um kml de 16 MB. No exemplo anterior não era possível clicar num polígono e obter todos os atributos associados a cada condado. No exemplo seguinte, com kml, é, mas vai demorar mais a carregar. E as queries dinâmicas estão lá na mesma:

      








Podem dar uma vista de olhos no código-fonte através da Page Source do vosso navegador. O código não está minificado.


 

domingo, 31 de março de 2019

geomatica.no.sapo.pt revisitado

Clicar no mapa. Sistemas de Coordenadas Altitude SRTM.



Sistemas de coordenadas, via proj4js:
- ETRS89
- ETRS89 / Portugal TM06
- Hayford-Gauss Datum 73 (IGP)
- Hayford-Gauss Datum Lisboa (IGP)
- Hayford-Gauss Datum Lisboa (Militar)

Altitudes obtidas a partir de SRTM 3 segundos, reamostrado para resolução espacial de 80 metros, através de Web Service publicado pelo próprio. SRTM 1 segundo, reamostrado para resolução espacial de 25 metros, para breve. Fonte dos dados: https://www.fc.up.pt/pessoas/jagoncal/

terça-feira, 26 de setembro de 2017

Para os interessados, a página http://geomatica.no.sapo.pt foi descontinuada.

http://geomatica.no.sapo.pt/