NavPT foi à vida

Pois é. Segundo o próprio site:

"Por razões de ordem superior e completamente alheias à equipa do Projecto NavPT, somos obrigados a terminar o projecto."

Enquanto não soubermos mais pormenores acerca das "razões de ordem superior", resta-nos lembrar que se tratava de um projecto que permitia visualizar trajectórias, em tempo real, relativas à RIV (Região de Informação de Voo) de Lisboa. Permitia também, entre outras coisas, ter acesso às comunicações entre as aeronaves e a torre de controlo.

Ou seja, nada a que não possamos ter acesso através de sites internacionais, o que poderá levantar alguma suspeição sobre a legitimidade das alegadas "razões de ordem superior". Para mais considerações sobre este e outros aspectos, consultar por exemplo este blog.

Tile Map Service

Depois de termos os nossos tiles gerados, precisamos mesmo de um servidor de mapas? Acho que não (proof of concept).

Proj4js - transformação de coordenadas

Relativamente à aplicação para transformação pontual de coordenadas que disponibilizo aqui, desde o ano passado, foram feitas algumas (pequenas) alterações:

• Actualização dos parâmetros de Bursa-Wolf, sendo agora utilizados os parâmetros mais recentes do IGP;
• Formatação das coordenadas geográficas (graus decimais) de saída para 8 casas decimais;
• Formatação das coordenadas cartográficas (metros) de saída para 2 casas decimais.

Da página do IGP referida constam também pontos transformados os quais poderão servir para comparar resultados.

Link:
Transformação de coordenadas com Proj4js (implementação em javascript da biblioteca PROJ)

Serviços: parâmetro opcional

Todos os serviços disponibilizados em geodivagar.appspot.com admitem um parâmetro opcional callback para a chamada de uma função em javascript. Pode ser utilizado, por exemplo, o JSONScriptRequest.

Para isso, incluir o script na página.

Definir o callback na chamada ao serviço:

   function onMapClick(overlay,latlng) {

  if (latlng) {
   coords=latlng;
   request = 'http://geodivagar.appspot.com/geogauss?lat=' + latlng.y + '&lng=' + latlng.x + '&callback=onComplete';   
   aObj = new JSONscriptRequest(request);
   aObj.buildScriptTag();
   aObj.addScriptTag();   
       }
   }

Definir a função que recebe o objecto JSON, por exemplo:

 function onComplete(jData) {
   var myHtml = 'X: ' + jData.x.toFixed(2)+ '
 Y: ' + jData.y.toFixed(2);
           map.openInfoWindow(coords, myHtml);
 }

Demo completa disponível aqui.

Dica simples e eficaz

Directamente do blog do Darren Cope:

Merge A Directory of Shapefiles Using OGR

mkdir merged
for %f in (*.shp) do (
if not exist merged\merged.shp (
ogr2ogr -f “esri shapefile” merged\merged.shp %f) else (
ogr2ogr -f “esri shapefile” -update -append merged\merged.shp %f -nln Merged )
)

Mais 2 serviços JSON

Mais 2 serviços com resposta JSON:

/wmsdetails Informação detalhada acerca de um determinado serviço WMS

Exemplo:

http://geodivagar.appspot.com/wmsdetails?url=http://mapas.igeo.pt/wms/caop/continente

Parâmetros: url do serviço WMS

Resposta:

{
     layers: [
          "Distritos",
          "NUT3",
          "NUT2",
          "Freguesias",
          "CAOP",
          "Concelhos",
     ],
     operations: [
          "GetCapabilities",
          "GetMap",
          "GetFeatureInfo",
          "DescribeLayer",
          "GetLegendGraphic",
          "GetStyles",
     ],
     title: "Carta Administrativa Oficial de Portugal (CAOP - Continente) - Versão 6.0",
     abstract: "Serviço de mapas WMS: Continente - Limites Administrativos Oficiais (Limites de País, Limites de NUT 2, Limites de NUT3, Limites de Distrito, Limites de Concelho e Limites de Freguesia) para o Continente e Ilhas. A esta informação está associada a toponímia, bem como outra informação descritiva como seja a área oficial de cada circunscrição administrativa.",
     version: "1.1.1",
     type: "OGC:WMS"
}

/layerdetails Informação detalhada acerca de uma determinada layer de um serviço WMS

Exemplo:

http://geodivagar.appspot.com/layerdetails?url=http://mapas.igeo.pt/wms/caop/continente&layer=Distritos

Parâmetros: url do serviço WMS e layer

Resposta:

{
     boundingBoxWGS84: [
         -10.1913,
         36.8987,
         -5.7141,
         42.1887,
     ],
     styles: {
           default: {
                legend: http://mapas.igeo.pt/wms/caop/continente?version=1.1.1&service=WMS&request=GetLegendGraphic&layer=Distritos&format=image/png,
                title: "default"
           }
     },
     crsOptions: [
          "EPSG:27492",
          "EPSG:20790",
          "EPSG:4326",
          "EPSG:4258",
          "EPSG:102160",
          "EPSG:102161",
          "EPSG:102164",
          "EPSG:102165"
     ],
    title: "Distritos"
}

Foi utilizada a biblioteca OwsLib

Restantes serviços disponíveis em http://geodivagar.appspot.com/

Serviços no App Engine

Aproveitei o facto de ter criado conta no Google App Engine para disponibilizar alguns serviços com parâmetros passados na URL e com resposta em JSON (quem quiser chamar REST a isto, está à vontade).

Os serviços são os seguintes (poucos por enquanto):

/geogauss Projecção directa de coordenadas geográficas ETRS89 em coordenadas cartográficas PTTM06

Exemplo: http://geodivagar.appspot.com/geogauss?lat=38.01132704&lng=-7.87167688

Parâmetros: "lat"=latitude "lng"=longitude

Resposta: JSON

/gaussgeo Projecção inversa de coordenadas cartográficas PTTM06 em coordenadas geográficas ETRS89

Exemplo: http://geodivagar.appspot.com/gaussgeo?x=22958.6550821&y=-183907.735462

Parâmetros: "x" "y"

Resposta: JSON

/arcmer Comprimento de arco de meridiano (metros) sobre o elipsóide associado ao sistema WGS-84 entre duas latitudes

Exemplo: http://geodivagar.appspot.com/arcmer?lat1=38.0&lat2=39.0

Parâmetros: "lat1"=latitude inicial "lat2"=latitude final

Resposta: JSON

/elevation Gera uma URL do GOOGLE CHART com o perfil de terreno entre dois pontos

Example: http://geodivagar.appspot.com/elevation?lat1=38.5&lat2=39.0&lng1=-8.0&lng2=-8.0

Parâmetros: "lat1" e "lng1" coordenadas do ponto de partida, "lat2" e "lng2" coordenadas do ponto de chegada

Nota: a utilização deste serviço pressupôe a aceitação dos Termos de Serviço da API do Google Maps bem como dos Termos das Google Chart Tools

Resposta: JSON

/dms2decimal Converte graus, minutos e segundos em graus decimais

Exemplo: http://geodivagar.appspot.com/dms2decimal?d=39&m=40&s=5.73

Parâmetros: "d"=graus "m"=minutos "s"=segundos

Resposta: JSON

/decimal2dms Converte graus decimais em graus, minutos e segundos

Exemplo: http://geodivagar.appspot.com/decimal2dms?decimal=39.6682583333

Parâmetro: "decimal"= graus decimais

Resposta: JSON

/geoutil Faz parsing de coordenadas geográficas

Exemplo 1: http://geodivagar.appspot.com/geoutil?str=39%2040m%205.73s%20N%208%207m%2059s%20W

Exemplo 2: http://geodivagar.appspot.com/geoutil?str=39.5%20N%209.0%20W

Parâmetro: "str"= aceita strings como "39 40m 5.73s N 8 7m 59s W" ou "39.5 N 9.0 W"

Resposta: JSON

Ainda não implementei serviços para processamento em batch mas poderei vir a fazê-lo no futuro. Adicionalmente, poderei vir a disponibilizar mais serviços dependendo da disponibilidade de tempo.

Na sequência do 10º aniversário do fim do SA

Na sequência do 10º aniversário do fim do Acesso Selectivo (SA) o qual ocorreu dia 1 de Maio, encontrei uma página com as primeiras edições da coluna Innovation da revista GPS World:

• GPS: A multipurpose system by D. Wells and A. Kleusberg
  January/February 1990, Vol. 1, No. 1, pp. 60-63.


• The limitations of GPS by A. Kleusberg and R.B. Langley
  March/April 1990, Vol. 1, No. 2, pp. 50-52.


• Why is the GPS signal so complex? by R.B. Langley
  May/June 1990, Vol. 1, No. 3, pp. 56-59.


• Electronic charts and GPS by M.J. Casey and P. Kielland
  July/August 1990, Vol. 1, No. 4, pp. 56-59.


• The issue of selective availability by Y. Georgiadou and K.D. Doucet
  September/October 1990, Vol. 1, No. 5, pp. 53-56.


• Comparing GPS and GLONASS by A. Kleusberg
  November/December 1990, Vol. 1, No. 6, pp. 52-54.


• The GPS receiver: An introduction by R.B. Langley
  January 1991, Vol. 2, No. 1, pp. 50-53.


• Precise, real-time dredge positioning by S.R. DeLoach
  February 1991, Vol. 2, No. 2, pp. 43-45.


• The orbits of GPS satellites by R.B. Langley
  March 1991, Vol. 2, No. 3, pp. 50-53.


• Ionospheric effects on GPS by J.A. Klobuchar
  April 1991, Vol. 2, No. 4, pp. 48-51.


• GPS vehicle location and navigation by E.J. Krakiwsky
  May 1991, Vol. 2, No. 5, pp. 50-53.


• Continuous monitoring of crustal deformation by Y. Bock
  June 1991, Vol. 2, No. 6, pp. 40-47.


• The mathematics of GPS by R.B. Langley
  July/August 1991, Vol. 2, No. 7, pp. 45-50.


• GPS in civil aviation by K.D. McDonald
  September 1991, Vol. 2, No. 8, pp. 52-59.


• GPS – satellites of opportunity for ionospheric monitoring by D. Coco
  October 1991, Vol. 2, No. 9, pp. 47-50.

ECW Compressor

A ERMapper sempre disponibilizou no seu site utilitários gratuitos nomeadamente a famosa ferramenta de compressão ECW Compressor (com GUI e tudo, mas também com o limite de 500 MB por imagem). Depois da empresa ter sido adquirida pela ERDAS deixei de ver o compressor na secção de downloads no site desta última, apesar de lá constarem outros produtos herdados da ERMapper, nomeadamente o visualizador gratuito ER Viewer.

Assim, para quem não queira utilizar o binómio GDAL/ECW ou directamente o ECW SDK (o qual neste momento também não se encontra disponível) ou outra ferramenta que utilize o SDK, disponibilizo para download a versão ECW Compressor 2.6. Actualmente não tenho bem a certeza acerca do licenciamento desta ferramenta por isso a utilização da mesma fica ao critério de cada um (por outras palavras, não culpem o mensageiro).

LINK

Pesquisa por códigos postais de 7 dígitos

Normalmente, um código postal de 7 dígitos (cp7) corresponde a uma frente de quarteirão (pelo menos em zonas urbanas). Consequentemente, a uma rua podem estar associados vários cp7 tanto do mesmo lado da rua como do lado oposto. Por isso, a pesquisa por cp7, na maioria dos casos, não chega para definir de forma inequívoca uma morada sendo por isso necessário completar a informação com o respectivo número de policia. Em princípio o nome do arruamento estará subentendido no cp7.

No entanto, existe quem insista na pesquisa por códigos postais de 7 dígitos.

Por exemplo, a pesquisa por códigos postais dos Mapas do Sapo deixa muito a desejar. Quem pesquisar por meia dúzia de códigos postais, constata rapidamente que o que é feito é uma comparação entre o código postal e a descrição de cada ponto de interesse (POI) constante da base de dados. Uma abordagem ad hoc, portanto. Se não houver nenhum POI cuja morada contenha o código postal pretendido, o resultado será errático (por exemplo, pesquisando pelo código postal do Fórum Picoas, 1050-996, vamos obter um resultado próximo da Fundação Calouste Gulbenkian, o que não tem nada a ver) .

Abordagem mais interessante encontramos em geocoder.pinguimcomfrio.net. O autor construiu uma base de dados com os códigos postais e respectivos arruamentos. Quando se pesquisa por um determinado código postal, a respectiva morada (rua e localidade) constante da base de dados é enviada para o geocoder do Google. Com esta solução ao menos temos a garantia que o resultado se encontra na rua pretendida, o que já é muito bom.