Planet Georezo » Autres planets

Un Planet est un site Web dynamique qui agrège le plus souvent sur une seule page, le contenu de notes, d'articles ou de billets publiés sur des blogs ou sites Web afin d'accentuer leur visibilité et de faire ressortir des contenus pertinents aux multiples formats (texte, audio, vidéo, Podcast). C'est un agrégateur de flux RSS. Il s'apparente à un portail web.
Vous pouvez lire le billet sur le blog La Minute pour plus d'informations sur les RSS !

Canaux

6472 éléments (0 non lus) dans 116 canaux

Dans la presse

Du côté des éditeurs

Toile géomatique francophone

Géomatique anglophone

Blogs GeoRezo.net

Autres planets

Planet OSGeo

Abonnez vous à cette page (RSS) !

Autres planets

Sauter à la page: 1 2 3 ... 12

21 janvier 2019
8:49
Stefano Costa: Usare SIGECweb su macOS o GNU/Linux
sur Planet OSGeo
SIGECweb è la piattaforma del Sistema informativo generale del catalogo accessibile via web.

Tramite SIGECweb tutti gli utenti abilitati (sia interni al Ministero, sia esterni) possono accedere al patrimonio informativo delle schede di catalogo per effettuare nuove catalogazioni, digitalizzazioni delle schede cartacee esistenti o campagne di revisione delle schede già inserite in banca dati, nonché ricerche estese. Si tratta di uno strumento fondamentale per la tutela del patrimonio culturale italiano.

Da alcuni anni, per una serie di problemi tecnici, SIGECweb è accessibile solo tramite una specifica versione del browser Firefox, la 45, che l’Istituto Centrale per il Catalogo e la Documentazione (ICCD) mette a disposizione tramite il proprio sito. In pratica questa versione personalizzata non sostituisce il browser esistente ma viene installata esclusivamente per l’uso del SIGECweb.

Il file di installazione fornito è disponibile solo per sistemi Windows, e una volta installato crea un nuovo profilo utente separato da quello principale (non troveremo la cronologia dei siti visitati, le password salvate, i segnalibri, le opzioni etc).

Se abbiamo esigenza di usare SIGECweb da macOS o GNU/Linux, è possibile comunque seguire attentamente i suggerimenti riportati nell’articolo Install an older version of Firefox e scaricare la stessa versione 45 per il nostro sistema operativo dall’archivio ufficiale di Mozilla:
- Firefox ESR 45.9.0 per macOS
- Firefox ESR 45.9.0 per GNU/Linux
La versione per GNU/Linux è direttamente funzionante una volta estratti i file dall’archivio compresso. È tuttavia necessario creare un diverso profilo utente, analogamente alla versione predisposta da ICCD. In questo profilo utente sarà anche salvata l’opzione di disattivare gli aggiornamenti, che dobbiamo limitare al solo Firefox per SIGECweb.

Il modo più semplice per creare il nuovo profilo utente è eseguire Firefox dal terminale dalla directory dove è stato scaricato:
```
» ./firefox -P
```
Schermata dell’archivio con la versione 45.9.0 di Firefox La finestra di selezione del profilo utente. Da qui, solo la prima volta, dobbiamo creare il nuovo profilo SIGECweb e deselezionare la voce “Utilizza il profilo selezionato senza chiedere all’avvio”.
Una volta creato il profilo utente dedicato, è necessario disattivare gli aggiornamenti automatici, che purtroppo non possiamo installare. Visto che è importantissimo mantenere aggiornata la versione di uso generale, l’uso di due profili separati ci permette di limitare questa disattivazione alla versione “parallela” dedicata a SIGECweb.

Adesso siamo pronti per continuare a usare il normale Firefox (debitamente aggiornato alla versione più recente!) per tutte le nostre attività, riservando la versione 45 a SIGECweb.

Chiudiamo tutte le finestre di Firefox, e poi apriamo nuovamente la versione normale. A questo punto comparirà nuovamente la finestra di selezione del profilo: possiamo riattivare l’utilizzo del profilo utente default senza chiedere all’avvio della versione Firefox normale, mentre useremo il profilo “SIGECweb” in modo dedicato avviando dal terminale.
Ripristinare il profilo di default per la versione aggiornata di Firefox.
Sempre dalla directory dove è stata scaricata la versione 45, possiamo usare questo comando:
```
» ./firefox -P SIGECweb --no-remote
```
L’opzione -P indica quale profilo caricare (senza far comparire la finestra di scelta) mentre l’opzione --no-remote permette di eseguire in parallelo entrambe le versioni di Firefox senza che quella dedicata a SIGECweb prenda il sopravvento.

Ora possiamo dedicarci alla catalogazione di beni culturali anche se il nostro sistema operativo non è Windows.
19 janvier 2019
13:21
Paulo van Breugel: Update release GRASS GIS 7.4.4
sur Planet OSGeo

For those who missed it, a new update release GRASS GIS 7.4.4 is available since the 4th of January. It mainly brings bugfixes, but it also includes an important new function, the module r.mapcalc.simple. This module is especially important for a better integration with QGIS. It therefore has already been dubbed the “QGIS friendship” release …
Continue reading Update release GRASS GIS 7.4.4
1:11
GeoTools Team: GeoTools 20.2 Released
sur Planet OSGeo

The GeoTools team is happy to announce the release of GeoTools 20.2: geotools-20.2-bin.zip geotools-20.2-doc.zip geotools-20.2-project.zip geotools-20.2-userguide.zip maven repository This release is a stable release and is recommend for new development and production systems. This release is made in conjunction with GeoServer 2.14.2. Improvements Add the possibility to specify the target
1:08
GeoServer Team: GeoServer 2.14.2 Released
sur Planet OSGeo
We are happy to announce the release of GeoServer 2.14.2. Downloads are provided (zip|war|exe) along with docs (html|pdf) and extensions.

This is a stable release of the GeoServer 2.14 series and is recommended for all production systems. Users of prior releases of GeoServer are encouraged to upgrade.

This release is made in conjunction with GeoTools 20.2 and GeoWebCache 1.14.2. Thanks to all who contributed to this release.

For more information please see our release notes (2.14.2 | 2.14.1|2.14.0|2.14-RC).
Improvements and Fixes
This release includes a number of new features and improvements:
- WMTS Restful binding is now accessible to all users (used to be limited to admins) and works with workspace specific services
- gs:DownloadEstimator (almost always) returns true when estimating full raster downloads at native resolution
- Cannot create jp2k coverage through rest (IndexOutOfBounds)
- KML ignores sortBy parameter when querying records
- NullPointerException when using env() function with LIKE operator in CSS filters
- Can’t modify existing GWC blobstore via UI without renaming
- NPE if a Jiffle Rendering Transformation is used with Channel Selection
- OpenLayers2 preview does not trigger automatically on IE8
- Bad rendering with JAI-EXT and Input/Output TransparentColor options
- Complex MongoDB generated properties are not correctly handlded in SLDs
- Move the GeoServer ENV Parametrization documentation to a more general Section
- Allow expressions in ColorMapEntry labels for GetLegendGraphic
About GeoServer 2.14 Series
Additional information on the GeoServer 2.14 series:
- New MongoDB extension added
- Style editor improvements including side by side editing
- Nearest match support for WMS dimension handling
- Upgrade notes documenting REST API feature type definition change
- State of GeoServer 2.14 (SlideShare)
- GeoServer Ecosystem (SlideShare)
- GeoServer Developers Workshop (SlideShare)
17 janvier 2019
18:58
Martin Davis: Fun wit JEQL: Hilbert Curves
sur Planet OSGeo

Hilbert Curve of order 4:

Hilbert Curve of order 6:Code is in the JEQL script repo.

import jeql.std.function.HashFunction;

hilbertOrder = 6;
side = Val.toInt( Math.pow(2, hilbertOrder) );
count = side * side;

radius = 1;

t = select * from Generate.sequence( 0, count-2 );

t = select i, geom: Geom.buffer(hilbertEdge, 0.4)
let
hilbertPt1 = HashFunction.hilbertPoint(hilbertOrder, i),
hilbertPt2 = HashFunction.hilbertPoint(hilbertOrder, i+1),
hilbertEdge = Geom.createLineFromPoints( hilbertPt1, hilbertPt2 )
from t;

t1 = select *,
styleFill: clr, styleStroke: clr, styleStrokeWidth: 1
let
clr = Color.toRGBfromHSV(Val.toDouble(i) / count, 1, 1)
from t;
Mem t1;

The function hilbertPoint uses the efficient algorithm from� [threadlocalmutex.com] Code is on Github.
16 janvier 2019
13:31
Fernando Quadro: O usuário root no GeoServer
sur Planet OSGeo
Se você já usou ou está ambientado a utilizar o GeoServer, deve lembrar que ao realizar a instalação, no primeiro acesso, você utiliza o usuário admin (que é o administrador do sistema).

Porém, não é desse usuário que irei falar nesse post, e sim do usuário “root” que muitos desconhecem.

Devido a natureza altamente configurável, a segurança do GeoServer pode ocasionar que um administrador interrompa a autenticação normal, desabilitando todos os usuários, inclusive as contas administrativas.

Por esse motivo, o subsistema de segurança do GeoServer contém uma conta “root” que está sempre ativa, independentemente do estado da configuração de segurança. Assim, no estilo UNIX, essa conta fornece o status de “super usuário” e destina-se a fornecer um método de acesso alternativo para corrigir problemas de configuração. É importante ressaltar que seu nome não pode ser alterado e a senha da conta “root” é a senha mestra.

Por padrão, a senha mestra é gerada e armazenada em um arquivo denominado security/masterpw.info usando texto simples (não encriptado). Ao atualizar o diretório de dados de um GeoServer existente (versões 2.1.x e inferiores), o algoritmo tenta descobrir a senha de um usuário com a permissão ROLE_ADMINISTRATOR. Se tal senha for encontrada e o tamanho da mesma tiver no mínimo 8 caracteres, o GeoServer usa essa senha como senha mestra.

Desta forma, o arquivo security/masterpw.info é um risco de segurança. O administrador deve ler este arquivo, verificar a senha mestra registrando-se no GeoServer como usuário “root” e depois disso este arquivo deve ser removido.

Por padrão, o login na interface administrativa e nas APIs REST usando a senha mestra está desativado. Para ativá-lo, você precisará alterar manualmente o Master Password Provider no arquivo config.xml, (geralmente localizado em security/masterpw/default/config.xml), adicionando a seguinte instrução:
```
<loginenabled>true</loginenabled>
```
Fonte: GeoServer Documentation
15 janvier 2019
20:41
Stefano Costa: I libri che ho letto nel 2018
sur Planet OSGeo

Nel 2018 ho fatto una scelta piuttosto originale per guidare le mie letture: ho recuperato a casa dei miei genitori un pacco di libri rimasti impilati negli ultimi anni, generalmente arrivati in regalo, che per qualche motivo non mi ero mai portato via. Libri che non avevo letto, ecco. Libri che avevo dimenticato di leggere. Uno magari potrebbe decidere di leggere solo una certa autrice per tutto l’anno, solo libri pubblicati nel 1965, o qualche criterio del genere.

Ho finito la pila di libri non letti e poi mi sono tuffato dentro 4 3 2 1 come una lontra nel fiume.
Daniele De Silva, Non avevo capito niente
Questo non era rimasto impilato ma mi ha dato l’idea di una annata a tema. Un po’ sconnesso nell’incedere ma è stata una bella lettura, forse troppo veloce.
Annamaria Fassio, I giorni del Minotauro
Un giallo piemontese edito da Frilli. Ben architettato e ambientato.
Chinua Achebe, Non più tranquilli
Se con Le cose crollano eravamo di fronte alla tragedia di una società antica, di un suo protagonista inizialmente invincibile, il salto di due generazioni ci porta in un tempo dell’Africa quasi contemporaneo, eppure ancora legato a quella società antica, ai suoi legami indissolubili da cui è difficile sciogliersi anche per chi è apparentemente molto brillante.
Antonio Pennacchi, Canale Mussolini
A me Canale Mussolini non è piaciuto. Ho faticato per leggerlo. Ho sperato a lungo che abbandonasse la prosa dal fare omerico via via che la storia si dipanava, ripetendosi sempre uguale, e invece niente. Ci ho sperato un po’ perché varie persone erano rimaste entusiaste. Di sicuro è scritto per conquistare il lettore. Di sicuro non avevo mai letto una giustificazione così lunga del fascismo, così assolutoria per chi ne è stato protagonista al minuto e così, diciamo, controversa nella figura del narratore. È difficile sospendere il giudizio mentre si legge questa epopea familiare fatta di personaggi tragici, immobili.

Come dice, scusi? Canale Mussolini avrebbe vinto il Premio Strega e io non capisco un’acca della letteratura italiana contemporanea? Ma io sto unicamente raccontando il libro come l’ho trovato io, che guarda caso poi risulta scritto da un prete, e sempre quest’anno mi son trovato con un altro illustre premiato ancor più penoso, e in fin dei conti se per far rinascere l’epica italiana dobbiamo sorbirci un revisionismo palloso e ripetitivo, allora lasciamola nella tomba. Questa è la mia versione dei fatti, poi vedete voi e andate in malora.
Enrico Giannichedda, Quasi giallo
Non conosco molti archeologi che scrivono gialli (ad eccezione, forse, di Fred Vargas che non è proprio un’archeologa). Quindi non sono preoccupato di dire che questo quasi giallo non mi è piaciuto moltissimo. Forse perché le parti archeologiche, che abbondano e quasi debordano, mi sono note in buon dettaglio sia per studio sia per averle sentite proprio dalla voce di Enrico Giannichedda, e in questa cornice sembrano paradossalmente meno interessanti. Forse perché le parti gialle sono abbastanza solidamente nel solco del genere, senza però essere veramente nel solco giallo della copertina, magari più noir o 10YR 2/1. Altri che hanno letto il libro non sono rimasti soddisfatti dal finale, che però tutto sommato a me piace.
Michele Serra, Ognuno potrebbe
Michele Serra è un membro di quella categoria di persone che svolgono il meschino dovere di partorire ogni giorno un pensierino da scolaro delle elementari per la pubblica fruizione. La forma libresca non migliora l’esito e questa storia di un precario, scritta da uno che il precario non sa nemmeno recitarlo sotto forma di macchietta, non piace proprio.
Timur Vermes, Lui è tornato
Tanto inquietante (oggi ancora più di quando è stato scritto, solo nel 2011) quanto scopppiettante, questo è il libro che più mi è piaciuto leggere nel 2018, carico di continui rilanci che solo un buonsenso ormai intorpidito può considerare assurdi. Non ho visto il film che ne è stato tratto, in cui il protagonista non è Hitler bensì Mussolini, ma ne ho parlato con alcuni che lo hanno visto: il discorso è andato a parare sul “messaggio” che l’autore trasmette. Io penso che sia un messaggio molto implicito, che l’autore abbia saputo mostrare in modo eccellente come funzionano, come possono funzionare certi meccanismi psicologici e sociali che hanno conseguenze rapidamente irreparabili, senza bisogno di manifestare una “ovvia” valutazione negativa che avrebbe reso molto meno incisivo il suo autentico messaggio. In ogni caso, lui è tornato già parecchie volte negli ultimi anni e sembra che non siano state vendute abbastanza copie di questo libro.
Edoardo Nesi, Storia della mia gente
Chi abbia dato un premio a questo libro, dovrebbe almeno giustificare il vistoso errore grammaticale del titolo, poiché questa è una storia individuale e solipsistica, vissuta e narrata alla prima persona singolare, da parte di un autore che, del suo passato rimpianto da rampollo fallito di famiglia operosa, ricorda i nomi propri dei macchinari (femminili, ovviamente) ma non quelli degli operai. Uno che dedica varie pagine al Martini nell’ambito del tracollo economico della piccola media impresa manifatturiera italiana. Uno che è felice di scendere in piazza a manifestare perché lo fa stare bene.

Il motivo del blasone è presto detto, essendo l’autore stato precedentemente inserito nella shortlist ha deciso di titillare il premio stesso con continui rimandi ad esso, rendendolo (immaginiamo a propria insaputa) co-protagonista di questa farsa che sa rendersi lucida solo in forma onirica – rivelandosi anche genuinamente impregnata di razzismo.
Paul Auster, 4 3 2 1
Questo libro, questi libri, sono anzitutto una smisurata forma di devozione, una lunghissima dichiarazione d’amore verso la scrittura e la letteratura, verso New York e Parigi e forse anche verso l’essere ebrei negli Stati Uniti. Non per caso la prima parte mi ha ricordato fortissimamente Middlesex, così come il rimando continuo tra storia personale e storia collettiva, non un semplice sfondo ma un palcoscenico.

Il volume è imponente e ha richiesto una certa disciplina nella lettura, evitando tassativamente di leggere più di un capitolo al giorno per non andare in confusione, ma la trama è certo uno degli elementi meno portanti del capolavoro, come prevedibile. Trama che è composta anche di frammenti, mattoni autobiografici composti in modi sempre diversi.

Ci sono capoversi lunghissimi che vorresti non finissero mai. Ci sono liste, ma che liste, di libri, di film, di poesie, di musica. C’è tantissimo sesso, muoiono molte persone e sono sempre le stesse persone i personaggi che vivono 4, 3, 2 vite leggermente o completamente diverse dall’una che tiene il filo.

Si ride, si gode e si soffre moltissimo con questo libro.
Paolo Lazzarin, Patagonia
Questo che mi ha regalato Elisa il 26 dicembre è un delizioso resoconto fotografico di viaggio. La Patagonia, al di là delle frasi da guida turistica, è difficile da spiegare, e la vastità degli spazi soverchia la vista, la mente. Questo viaggio si è svolto principalmente lungo la catena andina, via terra in direzione sud e via mare in direzione nord. Il nostro viaggio, più spezzato e incoerente, era stato certamente meno lento ma comunque sfogliando le pagine ho ritrovato quelle montagne, quelle strade e quei guanachi — perché poi sono uno dei compagni di viaggio più memorabili. Grazie.

Nel 2019 ho deciso che leggerò solo libri scritti da autrici.
14:42
gvSIG Team: Working with historical maps: Georeferencing images. GIS applied to archaeology
sur Planet OSGeo

A new module of the free course about GIS for archaeologists is now available, and we are now halfway through the course. At this module we will see how we can use Geographic Information Systems to georeference historical maps, so we can use them for our projects and cross them with current cartography.

During the exercise, we will proceed to georeference a scanned image with historical information of a city using reference cartography.

If you have not done the previous exercises, we recommend you to start with the first module.

We remind you that the data to follow the course can be downloaded in a compressed .zip file from this link.

Remember that for doubts about software use, carrying out the exercises, and problems that you find when carrying out the course, you can use the gvSIG users mailing list.

Module 4. Managing historical maps: Georeferencing images

Previous modules:
Module 1: Preparation of an archaeological project

Module 2: Digitalizing prospection areas

Module 3: Analysis of prospection results: Vector geoprocessing
14:42
gvSIG Team: Trabajando con mapas antiguos: Georreferenciación de imágenes. SIG aplicados a Arqueología
sur Planet OSGeo

Os traemos un nuevo módulo con el que llegamos al ecuador del curso gratuito de SIG aplicados a Arqueología. En este tema veremos como podemos utilizar los Sistemas de Información Geográfica para georreferenciar mapas antiguos, pudiendo así utilizarlos para nuestros proyectos y cruzarlos con cartografía actual.

Durante el ejercicio, a partir de una imagen escaneada con información histórica de una ciudad, procederemos a georreferenciarla utilizando cartografía de referencia.

Si aún no has realizado los ejercicios anteriores, te recomendamos que empieces por el primer módulo.

Os recordamos que los datos para poder seguir el curso los podéis descargar comprimidos en un fichero .zip en este enlace.

Para dudas en el manejo del software, realización de los ejercicios y problemas que puedan surgir en la realización del curso, siempre podéis utilizar la lista de usuarios de gvSIG.
Tema 4. Gestionando mapas antiguos: georreferenciación de imágenes
Temas anteriores
Tema 1: Preparación de un proyecto arqueológico: vistas, capas y tablas

Tema 2: Digitalizar zonas de prospecciones

Tema 3: Análisis de resultados de prospecciones: Geoprocesos vectoriales
8:54
OSGeo-fr: Les vidéos des rencontres QGIS 2018 sont en ligne!
sur Planet OSGeo

Les vidéos des conférences sont en ligne désormais. Les liens sont disponibles sur cette page
8:52
OSGeo-fr: Les vidéos des rencontres QGIS 2018 sont en ligne!
sur Planet OSGeo

Les vidéos des conférences sont en ligne désormais. Les liens sont disponibles sur cette page
7:00
Blog 2 Engenheiros: Como desenhar uma Área de Estudo no QGIS 3.4 considerando uma Bacia Hidrográfica?
sur Planet OSGeo

Todo e qualquer trabalho deve ter um limite espacial, e não é diferente para os estudos ambientais.

É comum utilizarmos como referência os limites das bacias hidrográficas para determinar nossas áreas de estudo. Porém, como realizamos esse procedimento dentro do QGIS 3.4?

Confira nossa postagem e veja diferentes formas de desenhar sua área de estudo.
Desenho livre no QGIS
Um dos caminhos é desenhar manualmente os limites da sua bacia hidrográfica utilizando um modelo digital de terreno (MDT) e suas curvas de nível.

Para realizar esse procedimento, iremos baixar um MDT da JAXA (Agência Espacial Japonesa), a qual disponibiliza gratuitamente, mediante cadastro, esses dados para o mundo inteiro com resolução de 30 metros (1 arco segundo).

O registro é realizado fornecendo dados simples, como nome, instituição, nacionalidade, propósito e e-mail. Clique aqui para realizá-lo. Logo após o registro, você receberá um e-mail para confirmar sua conta e em seguida, outro com o ID e senha de acesso.

Ao acessar o banco de dados da JAXA, você encontrará um mapa índice, no qual você deverá selecionar o local da sua área de interesse. Em seguida, uma nova janela irá abrir para que você dê mais um zoom na sua área de interesse. Por fim, após essa ultima seleção, você conseguirá fazer o download do MDT.
Mapa índice apresentado no site da JAXA para download do modelo digital do terreno.
Após realizar o download o seu MDT, adicione ele ao QGIS utilizando o gerenciador de fonte de dados (“Open Data Source Manager”), especificamente na guia Raster. Ou ainda, clicando em Camada (Layer), adicionar camada (Add Layer) e adicionar camada raster (Add Raster Layer).

Como o MDT baixado abrange uma área muito extensa e para evitar processamento desnecessário de curvas de nível, iremos recortar ele usando a ferramenta “Clip raster by extent”.

Essa ferramenta esta disponível pelo GDAL, dentro de Processing Tools (Processing > Toolbox), ou pelo atalho Ctrl + Alt + T.

Nesta ferramenta, você irá determinar qual é o arquivo raster a ser recortado (1), onde ele será recordado (2) e onde o resultado será salvo (3). No nosso caso, iremos demarcar um retângulo em uma bacia hidrográfica qualquer, mas você pode abranger sua área de interesse e acrescentar alguns quilômetros a este recorte.
Ferramenta para recortar arquivos matriciais no QGIS.
Com um MDT menor, podemos executar a ferramenta “Contour” para extrairmos as curvas de nível do nosso raster. Esta ferramenta também encontra-se no Processing Toolbox e ao abrir ela, você terá que fornecer dados como o MDT (Input Layer), intervalo entre as curvas de nível (Interval between contour lines) e onde as curvas de nível serão salvas.
Ferramenta para criação de curvas de nível no QGIS.
Lembrando que este procedimento pode ser realizado sem as curvas de nível, porém, elas auxiliam bastante na delimitação de bacias hidrográficas, pois os divisores de águas são sempre perpendiculares às linhas das curvas de nível.

Desta forma, verifique a localização do seu empreendimento (ou área de interesse) e lance ele no seu mapa (ou abra o shapefile de localização dele) e em seguida, marque um ponto a jusante dele (de preferência, usando shapefile da rede hidrográfica local).

A partir desse ponto a jusante, iremos desenhar nossa área de estudo.

Para criar um novo shapefile (e desenhar nossa área de estudo), iremos clicar em Camadas (Layer), Criar Camada (Create Layer) e Nova Camada Shapefile (New Shapefile Layer). Na janela que será aberta, você deverá fornecer informações como o nome do arquivo, geometria (neste caso, deve ser polígono), sistema de coordenadas e os campos da tabela de atributos.
Criação de um novo shapefile no QGIS 3.4.
Agora, para realizar o desenho, você deve clicar no ícone com um lápis amarelo enquanto o shapefile criado esteja selecionado. Ao realizar esse procedimento, você irá ligar o modo de edição.

Logo ao lado deste botão (lápis amarelo), há um ícone para criação de novos polígonos, clique nele para começar a desenhar sobre o MDT e as curvas de nível.

Existem várias formas de começar a delimitar sua bacia hidrográfica (neste caso, nossa área de estudo). É possível começar marcando os topos (partes mais altas) no MDT e a partir do ponto a jusante, iremos ligar esses pontos com linhas perpendiculares às das curvas de nível. O vídeo a seguir demonstra esse procedimento.
Como delimitar uma bacia hidrográfica manualmente.
Há outros vídeos, como este aqui, porém, são mais extensos (mas em português).

Após finalizar a delimitação da sua bacia hidrográfica, clique com o botão direito e o QGIS irá criar um novo polígono. Não se esqueça de clicar no disquete para salvar as modificações realizadas no seu shapefile.

O resultado deste procedimento é apresentado na figura a seguir para o MDT que baixamos do site da JAXA.
Área de estudo delimitada usando MDT e curvas de nível.
Desta forma, você pode criar inúmeras áreas de estudos, seja contornando os divisores de águas a partir de um MDT ou contornando a área industrial por meio de uma imagem aérea.
Usando limites já definidos de BH
Você pode baixar os limites de bacias hidrográficas de fontes oficiais como IMA (SC) e EPAGRI. Após realizar o download destes arquivos, basta copiar e colar a sua bacia hidrográfica de interesse para delimitar sua área de estudo.

Esse procedimento é realizado criando-se um novo arquivo vetorial (como apresentado acima). Liga-se o modo edição para o arquivo contendo as bacias hidrográficas já delimitadas, seleciona-se o polígono da bacia hidrográfica de interesse e copia-se ela.

Agora, vá até o shapefile criado (em branco) e cole o polígono recém copiado.

Existe ainda uma terceira forma de delimitar bacias hidrográficas de forma automática, porém, deixaremos essa opção para outras postagem. Caso tenha ficado com alguma dúvida, deixe ela nos comentários que estaremos respondendo assim que possível.
14 janvier 2019
11:54
Fernando Quadro: Sorteio: 01 bolsa de estudos da Geocursos
sur Planet OSGeo

Promoção: De volta às Aulas com a GEOCURSOS

A Geocursos lançou uma nova promoção que vai sortear uma bolsa de estudos de 100% no dia 15 de fevereiro. O ganhador poderá escolher entre os cursos de:

– OpenLayers 4 (EAD)
– PostgreSQL (EAD)
– DBA PostGIS (Ao vivo)

Veja como participar, lendo o regulamento logo abaixo.

REGULAMENTO DO SORTEIO

– SIGA o perfil da @geocursos no instagram
– CURTA a foto oficial da promoção (no instagram)
– MARQUE um amigo(a) nos comentário da publicação do instagram (não vale fake, empresas ou famosos)
– Pensamento POSITIVO

Participe quantas vezes quiser, quanto mais comentar, mais chances de ganhar, desde que usando nomes de amigos diferentes.

Informações sobre os cursos é só falar com a @geocursos no direct.

O Resultado será divulgado no dia 15/02 nos stories.

PS1: O perfil do ganhador deve estar público/aberto.
PS2: O ganhador tem até dia 18/02 para se manifestar via direct (@geocursos), caso contrário, realizaremos novo sorteio.
PS3: A inscrição no curso (OpenLayers ou PostgreSQL) se dará 3 dias após divulgação do resultado, e o ganhador terá 3 meses para realizar o curso.
PS4: O Curso DBA PostGIS, é online e ao vivo, e tem início no dia 16/03 (sábado).

Boa sorte!
11 janvier 2019
3:04
Just van den Broecke: Cheers to 2018
sur Planet OSGeo
Never too late to look back. Inspired by the concise bullet-point-style of Tom Kralidis’ Cheers to 2018, my personal highlights of living and working in the Open Source Geospatial and OSGeo(.nl|org)-world in 2018. I’ll post 2019 reso’s later. Not just personal projects, but also events from the OSGeo.nl Local Chapter in which I am involved since its beginnings.
January
- One of my 2018-resolutions was to focus more on SensorWeb/Spatiotemporal inspired by working since 2015 on the Smart Emission project. Started a blog-post-series “Emits”.
- With OSGeo.nl we always start the year (since 2011) with a joint OSGeo.nl and OpenStreetMap NL new-years party at Cafe Dudok in Hilversum. Many of the 2018-plans discussed there became reality. Not at least the two QGIS-Days and FOSS4G-NL 2018 conference (see below).
- The first Dutch QGIS Day hosted by the IHE Delft Institute for Water Education. Very well attended (160?) by both regular QGIS-users and IHE-students. Big thanks to Erik Meerburg (what would OSGeo.nl be without Erik?) and Hans van der Kwast (IHE).
- Started working on a contract for the European Commission Joint Research Center (Ispra, Italy) to deploy five AirSensEUR air-quality stations within the Netherlands, connecting these to the Smart Emission Platform. More in this blog-post.
February
- As first-timer went to TheThingsNetwork Conference. What an amazing event! For some years already into Lora(WAN) and maintaining a TTN Gateway, I learned so much, also from the workshops there, that I wrote a separate blog-post about this event.
- Started contracting (eventually until sept 30, 2018) at Dutch Kadaster to migrate the entire Smart Emission (SE) Platform to Kadaster infrastructure. Initially just migrating the SE Docker Images, normally two-weeks (to be honest, a day’s) work, BUT, as a (nice) surprise learned we had to move SE “To The Cloud”, that is to Azure Kubernetes (AKS), well that kept me and one Kadaster-colleague (learned a lot from you Gerwin!) busy till end of september at least 2-3 days a week!
March
- As first-timer went to the OSGeo 2018 Code Sprint in Bonn. Mainly to work on GeoHealthCheck. (GHC). Another highlight! Not the least since I finally met the great Tom Kralidis! I had been working on GHC with Tom a couple of years but we somehow never met IRL. Spent most of week on a single PR (per-Resource scheduling) which required a re-architecture of GHC. But more importantly socializing with Tom and other folks there, well, Bonn, need I say more…Thanks organizers: Astrid Emde et al. for this great week!
- In that same codesprint week I also attended several talks at the yearly FOSSGIS Conference there. The FOSSGIS Conference is a major event in Europe as it involves all German-speaking (DACH) countries, plus their OpenStreetMap communities. Even with my high-school German language knowledge very well to follow. As an exception to the rule I was grateful to be allowed to give a GeoHealthCheck presentation in English. Thanks to CCC, FOSSGIS always videotapes all talks, so you can watch some of my stuttering here.
April
- With OSGeo.nl we and our community-members more and more try to organize small, hands-on events that are lightweight to organize. The NLExtract Hacking Day was a success! Thanks to organizer Frank Steggink (who also wrote the report), and a great space offered by WebMapper (Edward MacGillavry) in the Social Impact Factory Utrecht.
May
- What happened in this month? Luckily GitHub tracks my activity. Many commits for the Smart Emission Kubernetes challenge. I see a week without commits in May: made a great hiking trip with friends in Ireland: the Kerry Way. Rain? Only one evening when walking to the pub (always live music there)!
June
- More work on map5.nl, a subscription service for Dutch topographic, historical- and embellished hillshade and arial maps I started to host in 2015. Plus ongoing Smart Emission stuff.
- Gradually, a heat-wave was overtaking The Netherlands. Climate change?
July
- Heat-wave more and more building up.. But that did not stop the great OSGeo.nl organizing team to host the second FOSS4G-NL in Almere. July 10-11. Kudos to Erik Meerburg and team for two great days of workshops and conf!
- At the FOSS4G-NL I had the honor to open the conference and present the national Dutch GeoForum: geoforum.nl. We jointly (Kadaster, WebMapper, OSGeo.nl) extended the existing Dutch Kadaster PDOK Forum. As OSGeo.nl we did not want to launch yet another forum but to bundle forces. Since that launch hundreds of folks registered.
- Week after: off to FOSS4G-EU in Guimarães,Portugal. FOSS4G-EU is always something special. Folks gather from literally every country in Europe. From the Baltic States to Croatia. In a cosy atmosphere Jorge Gustavo Rocha and the LOC made this into a great event. Great talks during the day and every evening (and sometimes nights) in this small-town university city it was easy to find and gather with the OSGeo-folks on the main square Largo da Oliveira.
August
- Started to seriously deploy the five AirSensEUR (ASE) stations. First at an RIVM air quality station in Breukelen next to the main A2 highway. With the great help of Jan Vonk of RIVM we managed to attach and digitally connect these stations at their premises. To have an idea, read my blog post on this effort. The ASE stations will gather raw, uncalibrated, Air Quality data there for two months and use the onsite RIVM measurements as reference for later calibration.
September
- Again a small gap in my GitHub commit history: two weeks with my family (and our great Labrador Jamie) in rural France: Auvergne and Ardeche. I always love to set out GPS-hikes at home, hacking some of the French IGN topo-maps for personal use, and see how these turn out on the spot.
Oktober
- Contract with Dutch Kadaster ended on sept 30. Just in time we migrated the entire Smart Emission stack to Azure Kubernetes (AKS). Well, IMHO Kubernetes is just awesome, especially after a few years of ad-hoc Docker deployment. Yes, a learning curve and devops-paradigm shift. Though for smaller projects in 2019 I am now combining Ansible with Docker (Compose).
- After two months of gathering data for calibration, deployed the five AirSensEUR (ASE) stations at their target locations within The Netherlands. At each location challenges: an ASE needs 3 things: 1) to be firmly/rain-protected mounted physically, 2) 220V current and 3) WIFI. That combination turned out to be a challenge at each spot. But finally all five ASEs are now humming at their spots.
- And with OSGeo.nl we had the second 2018 QGIS-Day at Ordina Utrecht. Thanks to the LOC: Jonas, Coen et al. Observation: all hands-on workshops were overcrowded!
November
- More work on map5.nl and working on a next SAAS: GeoHealthCheck in the Cloud, as a subscription service, under the name geoqos (“Geo-Cause”) on geoqos.com. Hopefully that makes the GeoHealthCheck project sustainable.
- Presented Smart Emission at the Geo Users Festival by Geonovum. Bit of beamer issues, my slides on Slideshare.
- What’s that day after the GIS Day? Yes, the PostGIS Day! That’s exactly what the LOC team (Steven M. Ottens and Tom van Tilburg) organized. Deliberately keeping it low-profile to attract high-profile PostGIS hackers. And how they did! The format was an unconference. Even as a considered PostGIS veteran (starting with a PostGIS workshop from Paul Ramsey in Lausanne 2006) my brain was crushing hearing stuff like from Martijn Meijers – Indexing AIS data / point clouds with a Space Filling Curve.
December
- Through my various projects thinking: All GeoData is Spatio-Temporal. Presented in Eindhoven at a Ruimteschepper event: LoRaWAN TheThingsNetwork makers event. Slides are here. For one thing attempted to bridge the world of IoT with OGC SensorWeb via the SensorThings API in a PoC.
- That’s almost it of 2018. A last brainwave on GeoHealthCheck usability made me to introduce interactive Plotly.js based interactive history Run-graphs. Together with and thanks to Tom Kralidis reviewing that was my final GitHub commit in 2018. Ok finalized on jan 3, 2019 but what a year!
10 janvier 2019
12:25
Fernando Quadro: Curso DBA PostGIS – Turma 5
sur Planet OSGeo

Prezado leitor,

Você tem interesse em aprender a trabalhar com banco de dados espacial, e possui conhecimento prévio em algum banco de dados? Então esta é a sua oportunidade!

A GEOCURSOS acaba de lançar a Turma 5 do Curso DBA PostGIS. Este curso online oferece uma visão completa que vai desde uma revisão sobre o PostgreSQL até tópicos avançados do PostGIS, apresentando como trabalhar em sua totalidade com esta poderosa extensão espacial do banco PostgreSQL.

Este curso é formado pelos nossos cursos PostGIS Básico (16 horas online) + PostGIS Avançado (20 horas online) e acontecerá entre os dias 16 de março e 19 de maio (aos sábados).

Se você fosse comprar os cursos separadamente sairia pelo valor de R$ 798,00. Porém o curso esta com uma super promoção, e está saindo por apenas R$ 549,00.

Para maiores informações e para ver a ementa completa do curso, acesse:

[www.geocursos.com.br]
9:42
FOSS4G 2019 Bucharest: Registration for FOSS4G 2019 Bucharest is now open!
sur Planet OSGeo

We are honored to announce the international geospatial community that the registration for the geospatial event of the year FOSS4G 2019 is now open! And more…

We’ve made a promise to the community and we are delighted to stand before you today and say that we kept our word. We promised to do our best to make this global FOSS4G the most affordable of all times and it is! Joining the geospatial international community for one week packed with dynamic, valuable, timely events in a vibrant European capital city is more at hand than ever before. The EarlyBird for FOSS4G 2019 Bucharest is at the historical level of 330 EUR and it offers you everything you could ever need for such a tech gathering: access to plenary session (28-30 august), access to conference tracks (28- 30 august), access to 300+ talks, access to mobile app, ability to contact attendees, access to Ice Breaker Party (29 august), access to Gala Dinner (29 august), access to Pub Race Offers (28-30 august), access to Academic papers post event, access to presentations post event, conference kit, email support.

Furthermore, FOSS4G 2019 #historicalearlybird comes with a surprise not to be missed! For the first 100 participants, we have prepared 100 special FOSS4G 2019 T-shirts! This way, we – the Bucharest Local Organizing Committee – show our appreciation for your early, comforting support for the conference. If this you might have expected, our next surprise will definitely.. surprise you! We are deeply thankful for the opportunity to invite you in our home, Romania and we would like to offer you a small, captivating Romanian Passport. Curious? Be among the first 100 FOSS4G 2019 participants and you’ll receive a wonderful collection of interesting facts of Romania all packed in the Curious Passport.

No FOSS4G would be possible without the valuable input of sponsors. They are the conference’s backbone contributing to the values, principles and progress of our community. Support us in our work to organize FOSS4G2019 Bucharest as one of the most memorable FOSS4G events. The call for sponsors is now open. Please, take a look at our suggestions and find the one that best suits your needs and preferences.

Sponsoring is not yet possible but you want to actively contribute to the success of FOSS4G 2019 Bucharest? Become a volunteer! Besides valuable input in making FOSS4G 2019 happen, you will receive a free ticket to the conference (including the ice breaker and gala dinner). Interested?

Sponsorship is not for you, neither volunteering but you still want to help FOSS4G 2019, then Spread the Word!

Of course, FOSS4G 2019 is much much more than T-shirts and wonderful new places to visit. Program, Workshop and Academic committees are almost fully formed and ready for action! We are honored to say that the FOSS4G 2019 program will be shaped by remarkable people in the geospatial international community. Check out the Committee webpage to see who will review your FOSS4G 2019 contribution as Call for Papers and Workshop will open next week!

Follow us in the following months and find out all about the geospatial event of the year!
09 janvier 2019
10:30
gvSIG Team: Analysis of prospection results: Vector geoprocessing. GIS applied to Archeology
sur Planet OSGeo

A new module of the free GIS course applied to Archeology is now available. With this module we are already starting to explore some of the advanced tools of Geographic Information Systems.

We will see how to apply vector geoprocesses to analyze the geographic information related to an archaeological site. We will also work with the information obtained from a GPS in a GPX file.

If you have not done the previous exercises, we recommend you to start with the first module.

We remind you that the data to follow the course can be downloaded in a compressed .zip file from this link.

Remember that for doubts about software use, carrying out the exercises, and problems that you find when carrying out the course, you can use the gvSIG users mailing list.
Module 3. Analysis of prospection results
Previous modules:
Module 1: Preparation of an archaeological project

Module 2: Digitalizing prospection areas. GIS applied to Archeology
9:57
gvSIG Team: Análisis de resultados de prospecciones: Geoprocesos vectoriales. SIG aplicados a Arqueología
sur Planet OSGeo

Os traemos un nuevo módulo del curso gratuito de SIG aplicados a Arqueología. Con este tema ya empezamos a explorar algunas de las herramientas avanzadas de los Sistemas de Información Geográfica.

Veremos como aplicar geoprocesos vectoriales para analizar la información geográfica relacionada con un yacimiento. También se trabajará con la información obtenida con un GPS en un fichero GPX.

Si aún no has realizado los ejercicios anteriores, te recomendamos que empieces por el primer módulo.

Os recordamos que los datos para poder seguir el curso los podéis descargar comprimidos en un fichero .zip en este enlace.

Para dudas en el manejo del software, realización de los ejercicios y problemas que puedan surgir en la realización del curso, siempre podéis utilizar la lista de usuarios de gvSIG.
Tema 3. Análisis de resultados de prospecciones
Temas anteriores
Tema 1: Preparación de un proyecto arqueológico: vistas, capas y tablas

Tema 2: Digitalizar zonas de prospecciones
08 janvier 2019
12:58
Markus Neteler: GRASS GIS 7.4.4 released: QGIS friendship release
sur Planet OSGeo
We are pleased to announce the GRASS GIS 7.4.4 release

What’s new in a nutshell

The new update release GRASS GIS 7.4.4 is release with a few bugfixes and the addition of r.mapcalc.simple esp. for QGIS integration. An overview of the new features in the 7.4 release series is available at New Features in GRASS GIS 7.4.

As a stable release series, 7.4.x enjoys long-term support.

Binaries/Installer download:
- winGRASS 7.4.4: 32bit standalone installer | 64bit standalone installer
- winGRASS 7.4.4 OSGeo4W: 32bit OSGeo4W installer | 64bit OSGeo4W installer
- EPEL6 (CentOS6, RHEL6, …)
- Fedora and EPEL7 (CentOS7, RHEL7, …)
- Max OS X (fully bundled binaries)
- … further binary packages for other Linux distributions will follow shortly, please check at software downloads.
Source code download:
- [https:]]
- [https:]]
- To get the GRASS GIS 7.4.4 source code directly from SVN, see here.
More details:

See also our detailed announcement:
- [https:]]
- [https:]] (overview of new 7.4 stable release series)
- [https:]] (list of available addons)
  First time users may explore the first steps tutorial after installation.
About GRASS GIS

The Geographic Resources Analysis Support System ( [https:]] ), commonly referred to as GRASS GIS, is an Open Source Geographic Information System providing powerful raster, vector and geospatial processing capabilities in a single integrated software suite. GRASS GIS includes tools for spatial modeling, visualization of raster and vector data, management and analysis of geospatial data, and the processing of satellite and aerial imagery. It also provides the capability to produce sophisticated presentation graphics and hardcopy maps. GRASS GIS has been translated into about twenty languages and supports a huge array of data formats. It can be used either as a stand-alone application or as backend for other software packages such as QGIS and R geostatistics. It is distributed freely under the terms of the GNU General Public License (GPL). GRASS GIS is a founding member of the Open Source Geospatial Foundation (OSGeo).

The GRASS Development Team, January 2019

The post GRASS GIS 7.4.4 released: QGIS friendship release appeared first on GFOSS Blog | GRASS GIS and OSGeo News.
7:59
Marco Bernasocchi: New Year’s present – QField 1.0 RC1
sur Planet OSGeo

It was a long and winding road but we are very excited to announce the general availability of QField 1.0 Release Candidate 1.

Packed with loads of useful features like online and offline features digitizing, geometry and attributes editing, attribute search, powerful forms, theme switching, GPS support, camera integration and much more, QField is the powerful tool for those who need to edit on the go and would like to avoid standing in the swamp with a laptop or paper charts.

With a slick user interface, QField allows using QGIS projects on tablets and mobile devices. Thanks to the QGIS rendering engine, the map-results are identical and come with the full range of styling possibilities available on the desktop.

We ask you to help us test as much as possible this Release Candidate so that we can iron out as many bugs as possible before the final release of QField 1.0.

You can easily install QField using the playstore ( [qfield.org] ), find out more on the documentation site (http://qfield.org) and report problems to our issues tracking system ( [qfield.org] )

QField, like QGIS, is an open source project. Everyone is welcome to contribute to make the product even better – whether it is with financial support, enthusiastic programming, translation and documentation work or visionary ideas.

If you want to help us build a better QField or QGIS, or need any services related to the whole QGIS stack don’t hesitate to contact us.
6:02
Blog 2 Engenheiros: Como salvar seu projeto SIG e manter o Layout do mapa no ArcGIS e QGIS?
sur Planet OSGeo
Quando estamos trabalhando em um projeto ambiental, utilizamos vários dados de diferentes compartimentos, seja do meio físico, biótico ou socioeconômico. E quando essa informação pode ser representada espacialmente, usamos um sistema de informações geográficas (SIG) para visualizá-las.

Após finalizar a interpretação desses dados, exportamos os resultados no formato de mapas, os quais devem atender um padrão definido, caso contrário, teríamos vários mapas (do mesmo projeto) com aparências (i.e. layouts) diferentes.

Entretanto, como podemos padronizar essa confecção de mapas? Como salvamos nosso projeto no ArcGIS e no QGIS? Como podemos utilizar inúmeras vezes o mesmo template?

Confira nossa postagem e descubra.

Já comentamos em um dos nossos vídeos como realizar a organização das suas pastas do seu SIG. Nesta postagem, iremos mostrar um passo a passo para você criar seu template de mapa, de modo que você possa sempre aproveitá-lo em outros mapas e projetos.
Criando seu mapa base no ArcGIS
No ArcGIS, logo após abrir o software, normalmente, ele abre uma janela com alguns layouts de mapas já estabelecidos. Porém, no nosso caso, queremos criar um mapa novo, então, iremos clicar em New Maps (Novos Mapas) > Meus Templates (My Templates)(1) > Mapa Vazio (Blank Map)(2) e depois em OK (3).
Criando um mapa em branco no ArcGIS.
Defina o sistema de coordenadas do seu projeto e vá no Layout View.

Para definir o sistema de coordenadas do seu projeto, clique com o botão direito sobre o Data Frame na tabela de conteúdos (Table of Contents), vá em propriedades e em seguida, na aba Sistema de Coordenadas. Nesta aba, você escolhe qual sistema irá adotar.

No Layout View, insira e organize todos os elementos do seu mapa, tais como:
- Título;
- Legenda;
- Norte;
- Escalas gráfica e numérica;
- Grade de coordenadas;
- Outros elementos (por exemplo, autor, empresa, data e sistema de coordenadas adotada).
Não esqueça também de modificar as configurações do papel, que pode ser modificada clicando sobre o fundo com o botão direito e selecionando Configuração de Página e Impressão (Page and Print Setup).
Janela para configuração da página / folha de impressão no ArcGIS.
Após inserir os elementos do seu mapa, vá em Arquivo (File) e clique em Salvar (Save) ou Salvar Como (Save As). Uma nova janela será aberta solicitando onde você quer salvar o arquivo no formato MXD, que representa um documento do ArcMap.

Ao compartilhar esse documento, você estará enviando somente o layout do mapa e não estará incluso os arquivos que você abriu nele, tais como shapefiles e rasters.

Um exemplo de mapa criado é apresentado na figura abaixo.
Exemplo de Layout criado no ArcGIS.
Agora, com este arquivo no formato MXD é possível, ao abrir um novo projeto no ArcGIS, selecioná-lo para que você abra as mesmas configurações de layout que você já havia criado.

Para replicar esse layout para outros mapas, crie cópias deste arquivo MXD e renomeie elas com o nome do novo mapa (e.g. “mapaSolos.mxd”, “Hidrografia.mxd”).
Criando seu mapa base no QGIS 2.18
Na versão 2.18 do QGIS, é possível criar vários templates utilizando o gerenciador de layouts. Este pode ser acessado clicando-se em Projeto (Project), Novo Compositor de Impressão (New Print Composer). Após clicar neste item, o QGIS irá solicitar o nome do layout. Podemos chamá-lo de “Mapa A4 Retrato”.

Após o compositor de impressão ser aberto, é possível modificar as configurações da página na aba do lado direito (1) e em seguida, adicionar os item como Norte, Legenda e Escala pelos menu lateral (2) ou pelo menu superior Layout (3).
Compositor de mapas do QGIS 2.18.
Agora, após adicionar e organizar todos os elementos do seu mapa, clique em Compositor (Composer) e em seguida em Salvar Projeto (Save Project). Isso irá salvar o projeto e os layouts do compositor de layouts.

Lembre-se que ao clicar sobre um elemento, suas propriedades podem ser mudadas na aba do lado direito do compositor.

A figura abaixo mostra um exemplo de layout criado no QGIS 2.18.
Exemplo de mapa base criado no QGIS 2.18.
O arquivo salvo estará no formato QGS, formato dos projetos salvos no QGIS 2.18. Cabe salientar que no QGIS, por meio do compositor de impressões, é possível duplicar o layout que você já criou e aproveitá-lo para outro mapa com dados similares.

Leia também: Layout de Mapas no QGIS (II) e Layout de Mapas no QGIS (III)

Com esse arquivo QGS, você poderá duplicá-lo (criar cópias) para desenvolver outros mapas utilizando sempre o mesmo template.
Criando mapa base no QGIS 3
O procedimento para criação de templates no QGIS 3 é semelhante à versão 2.18, entretanto, ao salvar seu projeto, você verá que ele estará no formato QGZ.

Lembrando também que alguns itens (botões) mudaram de posição no compositor de mapas do QGIS 3, mas nada que interfira no desenvolvimento dos seus templates de mapas.

Agora você já conhece os caminhos para desenvolver seus templates de mapas e agilizar o desenvolvimento de novos mapas, sem precisar ter que criá-los manualmente.
07 janvier 2019
22:03
QGIS Blog: User question of the Month – Jan19 & answers from Dec
sur Planet OSGeo

In December, we wanted to know what QGIS.ORG should focus on in 2019.

Portuguese

Based on these results, in today’s PSC meeting, we’ve decided that the 2019 grant programme will be focusing on bug fixing and polishing existing features. So thanks to everyone who provided feedback!
New question
This month, we’d like to know if you have ever contributed to improving QGIS and – if yes – how. As you’ll see, there are many different ways to contribute to QGIS, so please go ahead and take the survey.

The survey is available in English, Spanish, Portuguese, French, Italian, Ukrainian, and Danish. If you want to help us translate user questions in more languages, please get in touch!
05 janvier 2019
15:31
Free and Open Source GIS Ramblings: PyQGIS101 part 10 published!
sur Planet OSGeo
PyQGIS 101: Introduction to QGIS Python programming for non-programmers has now reached the part 10 milestone!

Beyond the obligatory Hello world! example, the contents so far include:
If you’ve been thinking about learning Python programming, but never got around to actually start doing it, give PyQGIS101 a try.

I’d like to thank everyone who has already provided feedback to the exercises. Every comment is important to help me understand the pain points of learning Python for QGIS.

I recently read an article – unfortunately I forgot to bookmark it and cannot locate it anymore – that described the problems with learning to program very well: in the beginning, it’s rather slow going, you don’t know the right terminology and therefore don’t know what to google for when you run into issues. But there comes this point, when you finally get it, when the terminology becomes clearer, when you start thinking “that might work” and it actually does! I hope that PyQGIS101 will be a help along the way.
04 janvier 2019
21:00
GRASS GIS: GRASS GIS 7.4.4 released
sur Planet OSGeo

We are pleased to announce the GRASS GIS 7.4.4 release
13:52
Fernando Quadro: O GeoServer e o Java 11
sur Planet OSGeo

O Java Code Sprint 2018 foi um evento ambicioso que reuniu projetos Java da OSGeo para analisar a compatibilidade com Java 11. Com uma ótima resposta da comunidade, o GeoServer agora pode rodar no Java 11!

Esse sprint foi uma resposta importante da comunidade às mudanças do Java. A motivação imediata é dar aos usuários a opção de usar o Java 11 na próxima versão do GeoServer (2.15.x), agora que o OpenJDK está assumindo como projeto principal. Também temos o prazer de informar que o Java 8 continuará sendo suportado (graças ao compromisso da RedHat e da Adopt OpenJDK), dando a todos a oportunidade de migrar quando estiverem prontos.

É importante observar que o Java 11 não suporta mais o mecanismo de extensão usado para as bibliotecas JAI e ImageIO nativas. Onde as bibliotecas podiam ser instaladas diretamente no Sistema Operacional para aumento de performance das mesmas.

Como experimentar:

1. Faça o download e instale o Java 11 (jdk.java.net/11)

2. Faça o download do GeoServer 2.15-M0

3. O GeoServer 2.15-M0 será executado sob o Java 11 sem configuração adicional no Tomcat 9 ou mais recente e no Jetty 9.4.12 ou mais recente.

4. A distribuição binária inclui uma versão compatível do Jetty.

Fonte: GeoServer Blog
03 janvier 2019
11:45
gvSIG Team: Digitalizing prospection areas. GIS applied to Archeology
sur Planet OSGeo

We publish the second module (of seven) of the free ‘GIS applied to archaeology’ course, where you will learn to create geographic information through the editing tools offered by the open source software called ‘gvSIG Desktop’.

We will start by creating a new layer of information that will represent our prospection areas. In that layer we will create new geometries and we will fill its attribute table in with alphanumeric information.

We will also see how to add information to an existing layer using the ‘table join’ tool.

To follow these exercises you will start with the project created in the first module of the course, so if you have not done already we recommend you start there.

The data to follow the course can be downloaded compressed in a .zip file on this link.

Remember that for doubts about software use, carrying out the exercises, and problems that you find when carrying out the course, you can use the gvSIG users mailing list.

Module 2: Creating geographic information

Previous modules:

Module 1: Preparation of an archaeological project
11:12
gvSIG Team: Digitalizar zonas de prospecciones. SIG aplicados a Arqueología
sur Planet OSGeo

Os traemos el segundo módulo (de siete) del curso gratuito de SIG aplicados a arqueología, donde aprenderemos a crear información geográfica mediante las herramientas de edición que ofrece el software libre y gratuito gvSIG Desktop.

Comenzaremos creando una nueva capa de información que representará nuestra zona de prospecciones. En esa capa iremos creando nuevas geometrías y rellenando su información alfanumérica.

También veremos como agregar información a una capa existente mediante la herramienta de unión de tablas.

Para realizar estos ejercicios partiremos del proyecto creado en el primer módulo del curso, por lo que si aún no lo has hecho te recomendamos que empieces por ahí.

Los datos para poder seguir el curso los podéis descargar comprimidos en un fichero .zip en este enlace.

Recuerda que para dudas en el manejo del software, realización de los ejercicios y problemas que puedan surgir en la realización del curso, podéis utilizar la lista de usuarios de gvSIG.
Tema 2: Crear información geográfica
Temas anteriores
Tema 1: Preparación de un proyecto arqueológico: vistas, capas y tablas
2:30
GeoTools Team: GeoTools Java 11 Code Sprint Results
sur Planet OSGeo

With the recent GeoTools 21-M0 Milestone release available for everyone to try out ... we can finally report back on our success at the OSGeo Java 2018 Code Sprint. References: Java 11 Quickstart shows how to use the GeoTools with the new module system Java is still free for background on Java 8 and Java 11 support options (or our Java Install page summary table) Technical Challenge The
2:23
GeoServer Team: Java 2018 Code Sprint Results
sur Planet OSGeo
The Java 2018 code sprint was an ambitious event gathering up OSGeo Java projects to look at Java 11 compatibility. With a great response from the GeoServer community we are pleased to announce the GeoServer can now run in Java 11!

This sprint was an important response by our community to changes in the Java roadmap. The immediate motivation is to give our users the option of using Java 11 in the coming year now that OpenJDK is taking over as lead project. We are also pleased to report that Java 8 will continue to be supported (thanks to the commitment of RedHat and Adopt OpenJDK) giving everyone a chance to migrate when ready.

We had 11 participants for the sprint – thanks to everyone who attended!
- David Vick (Boundless), Devon Tucker (Boundless), Jim Hughes (CCRi), Jody Garnett (Boundless) Kevin Smith (Boundless), and Torben Barsballe (Boundless) gathered at the Boundless office in Victoria, Canada.
- Andrea Aime (GeoSolutions) and Antonello Andrea (Hydrologis) participated from Andrea’s home office in Italy.
- María Arias de Reyna (GeoCat), Ian Turton (Astun), Landon Blake (BKF Engineers), and Brad Hards participated remotely.
With support from OSGeo and event sponsors many individuals were able to meet up in person, with those hosed at the Boundless Victoria office enjoying the occasional break outside.

How to try it out:
To get started:
1. Download and install Java 11 from jdk.java.net/11
2. Download the GeoServer 2.15-M0 milestone release of GeoServer
  - The GeoServer 2.15-M0 will run under Java 11 with no additional configuration on Tomcat 9 or newer and Jetty 9.4.12 or newer.
  - The binary distribution includes a compatible version of Jetty.
We need your help:
- pulling together running on java 11 instructions for additional application servers, please download the WAR bundle to help out
- testing the windows installer
Sponsors
Thanks to Gaia3D for leading with a silver sponsorship, and ASTUN Technology, OSGeo:UK, and ATOL for their bronze sponsorships.

These sprints also require people to function, and we appreciate Boundless, GeoCat, ASTUN Technology, GeoSolutions and CCRi for their in-kind participation.

Key accomplishments:
- GeoServer was upgraded to Spring 5. Thanks to David Vick and everyone for this (especially the Spring Security upgrade with James).
- Kevin went through the same steps for GeoWebCache for a difficult couple of days.
- The EMF models that drive much of our parsing technology were upgraded. Thanks to Ian who worked through this, putting us in a much better spot to make changes to our codebase.
- Torben helped remove countless references to internal “com.sun” classes.
- A big thanks to Andrea for stepping up where needed (i.e. everywhere) and encouraging everyone through out the week
- Thanks to Andrea and Brad for preliminary work that made this sprint possible.
- Thanks to Jody for doing the milestone release.
Outstanding issues:
Certain components and dependencies could not be upgraded as no Java 11-compatible versions have been released. These have been removed or replaced by different libraries where possible, but some incompatible libraries remain.

The remaining issues are primarily either Internal API usage or Illegal access. These will only cause problems when running on Java 11:
Internal API usage
Some components and dependencies reference internal Java API. Depending on the nature of this API, calling this code may cause a runtime error or a warning when running on Java 11. If you run into any unusual failures when running under Java 11, please report them.

Examples of libraries producing these warnings:
- freemarker template library (used for GetFeatureInfo responses)
- ehcache library used for image mosaic performance
- Java Advanced Imaging used for image processing
Illegal Access
Java 11 introduces stricter access constraints on reflective operations. For now, such issues will only cause a warning in the logs, of the form:

WARNING: Illegal reflective access by …

In future versions of Java, this will instead cause a runtime error.

Known warnings you can expect when running in Java 11:
- org.hsqldb.persist.RAFileNIO
- org.parboiled.transform.AsmUtils
- org.fest.reflect.util.Accessibles
- com.google.protobuf.UnsafeUtil
- net.sf.cglib.core.ReflecUtils
If you run into a reflective access warning about a package not in this list, please report it.
Community Modules
The following community modules do not yet support Java 11:
- printing
- script
- spatiallite
- monitor-hibernate
Further reading:
02 janvier 2019
10:00
Gary Sherman: Plugin Builder 3.1
sur Planet OSGeo
We've released version 3.1 of the Plugin Builder for QGIS 3.x. This version contains a number of bug fixes and performance enhancements.

Here are some of the changes included since version 3.0.3:
- Fix issue with reload on generated plugins
- Move dialog creation to run method to improve startup performance
- Move help file generation files to proper method
- Include missing tags file
- Attempt to compile resources.qrc when plugin is generated (requires pyrcc5 in path)
- Set deployment directory in Makefile based on user OS (pb_tool is recommended over make)
- Check for valid URL format for tracker and repository
Compiling Resource File
If you have the resource compiler pyrcc5 in your path, the resource file will be compiled automatically when you generate your new plugin. This means it's ready to deploy and test without any further steps.
Using pb_tool
Although a Makefile is created when you generate a new plugin, you are encouraged to use pb_tool for compiling, deploying, and managing your plugins. It runs everywhere and simplifies your plugin development.

pb_tool works with both Python 2.x and 3.x.

See [g-sherman.github.io] for information.
Installing Plugin Builder
You can install Plugin Builder 3.1 from the Plugins -> Manage and Install Plugins... menu. Version 3.1 works on QGIS versions 3.0 and up.
Plugin Builder Links
- Home page: [g-sherman.github.io]
- Issues: [https:]
- Code repository: [https:]
  - Note: Plugin Builder for QGIS 2.x is in the master branch. Code for QGIS 3.x version is in the qgis3_version branch
Pull requests and code contributions are encouraged.
01 janvier 2019
9:00
From GIS to Remote Sensing: Very Happy New Year by From GIS to Remote Sensing
sur Planet OSGeo
This post is to wish you all a very happy new year!

2018 has been a great year for the Semi-Automatic Classification Plugin (SCP), with the release of the version 6 for QGIS 3.
I'm very grateful to all the people who have contributed to this project, through their work translating the interface and the user manual, fixing bugs and reporting issues, and donating.

In 2019 I'll continue the development of the SCP, in particular designing the new version 7 (which should be released in the early 2020), and of course maintaining the version 6.
In the new version 7 I'll focus on the development of new features for automating and improving the processing of remote sensing images, in particular:
- development of parallel processing in order to reduce the computational time;
- integration with ESA SNAP software for the preprocessing of Sentinel-2 images;
- development of standalone APIs to use the tools through the command line or Python;
- improvement of the user interface.
I'll consider updating the version 6 with some of these new features during 2019. And of course I'll write new tutorials about applications and case studies.
Again, I wish you a very happy new year!

For any comment or question, join the Facebook group about the Semi-Automatic Classification Plugin.
7:48
BostonGIS: Using pg_upgrade to upgrade PostgreSQL 9.3 PostGIS 2.1 to PostgreSQL 11 2.5 on Yum
sur Planet OSGeo

In a previous article Using pg upgrade to upgrade PostGIS without installing older version I demonstrated a trick for upgrading to a newer PostgreSQL instance from PostGIS 2.2 - 2.whatever without having to install the older version of PostGIS in your new PostgreSQL service. Unfortunately that trick does not work if coming from PostGIS 2.1 because in PostGIS 2.2 we renamed a c lib function that backed sql functions in 2.1.

Fear not. There is still a way to upgrade from 2.1 to 2.5 without installing an older version of PostGIS in your new PostgreSQL instance. To do so, you need to add a step and that is to remove the functions in 2.1 that are backed by this renamed lib function. In upcoming PostGIS 3.0, we've added this function back and have it throw an error so that even coming from PostGIS 2.1, you can upgrade just the same as you do from later versions.
Continue reading "Using pg_upgrade to upgrade PostgreSQL 9.3 PostGIS 2.1 to PostgreSQL 11 2.5 on Yum"
6:56
GeoServer Team: GeoServer 2.15-M0 Milestone Release
sur Planet OSGeo
We are happy to share a GeoServer 2.15-M0 milestone release with downloads (zip|war|exe), documentation (html|pdf) and extensions.

This milestone release is provided for everyone considering a Java 11 operational environment in 2019. A milestone release provides a technology preview and a chance to support the development team with early feedback: your assistance and participation is requested!

This release is a result of participation in the OSGeo Java 2018 Code Sprint and is made in conjunction with GeoTools 21-M0 and GeoWebCache 1.15-M0. We would like to thank organizations participating in the code sprint (Boundless, GeoSolutions, GeoCat, Astun Technology, CCRi) along with sprint sponsors (Gaia3D, atol, osgeo:uk, Astun Technology). Our gratitude goes out to the individuals who worked so hard to bring everything together.
Java 11 Support
The provided binary download works with either Java 8 or Java 11. Tomcat 9 or newer is required for the WAR install.

We will the update the user guide compatibility list based on your feedback and testing of this 2.15-M0 milestone release. Please note that Java 11 no longer supports the Java 2 extension mechanism used for native JAI and native ImageIO.

This is a good opportunity to review your operational environment and help test GeoServer with the environment you intend to use in 2019.
Java Roadmap Considerations
Our initial concerns with respect to continued availability of the Java 8 platform have been alleviated by industry action and commitment. Extensive Java 8 support options are now available, with RedHat making a public commitment to contribute fixes to the OpenJDK 8 codebase, and a range of organizations committed to making OpenJDK 8 builds available on a range of platforms.

Oracle has changed to a six month release schedule, donating additional components to OpenJDK to make it lead platform. Long term support (beyond six months) is being offered from a range of organizations notably RedHat OpenJDK and Adopt OpenJDK.

The net effect of these changes:
- If you have been using Oracle JDK up until now it is time to review your options
- Java 8 will continue to be available
- The Java ecosystem is now led by the open-source Open JDK project
See the GeoTools user guide for a table outlining the Java 8 and Java 11 alternatives to consider in the year ahead.
About GeoServer 2.15 Series
Additional information on the GeoServer 2.15 series:
- Release Notes (2.15-M0)
- Running on Java 11 (User Manual)
6:36
GeoTools Team: GeoTools 21-M0 Milestone Released
sur Planet OSGeo

GeoTools project is closing out 2018 with a GeoTools 21-M0 milestone release featuring Java 11 compatibility: geotools-21-M0-bin.zip geotools-21-M0-doc.zip geotools-21-M0-userguide.zip geotools-21-M0-project.zip maven repository This milestone release is a preview of GeoTools 21 giving everyone a chance to help test Java 11 compatibility. This release works with both Java 8 and Java 9,
29 décembre 2018
20:44
Cameron Shorter: Community Inside
sur Planet OSGeo

Parody of "intel inside®" logo,
a trademark of Intel Corporation.
Image SVG.
"Community inside" is the "intel inside®" of Open Source Software. It's an indicator of quality and long term sustainability.

Why? Because in the digital economy, collaboration out-competes competition!
So prioritise software developed by healthy communities. Look for diverse, supportive, meritocratic, productive, welcoming and inspiring teams. Then engage productively with them.

More ...
28 décembre 2018
18:59
Tom Kralidis: Cheers to 2018
sur Planet OSGeo

What a year! Following on from 2017, in no particular order, here goes: Bonn Code Sprint: talk about a worthwhile event! It was great to meet old and new OSGeo hackers. Serious progress made on many projects (keep reading) OWSLib: thanks to Carsten Ebrecht, we finally have a much better handling on tests (moving to […]
14:12
Even Rouault: SRS barn raising: 7th report
sur Planet OSGeo
This is the 7th progress report of the GDAL SRS barn effort.

On the PROJ side, things are consolidating up. Tens of "random" fixes have been pushed due to the GDAL autotest suite triggering a number of interesting cases. The C API has also been enhanced to accommodate for the needs of GDAL and libgeotiff. We also have received feedback from an early adopter (a developer of a pyproj binding based on PROJ master). The major development items have been the move of the WKT 1 syntax validation from GDAL to PROJ, as well as the development of an equivalent WKT 2 syntax validator in PROJ (this task has been useful to uncover a few minor issues in the draft of the future WKT2:2018 standard). A reorganization of the PROJ source tree (with a conversion of most C files to C++ files) has also been done, as a preliminary step, for a pull request to better integrate the new ISO-19111 functionality I developed those last months with the existing C API.
Regarding libgeotiff, a v1.4.3 maintenance version has been released with the fixes of the last two years, before the new works for the integration of PROJ master are done. libgeotiff development has been moved from Subversion to [https:]] . As a preliminary step, continuous integration capability has been added to test compilation under Linux/GCC and Windows/Visual Studo, with a few runtime tests.A pull request is ready with the integration of PROJ master with libgeotiff. It features:
- PROJ master / PROJ 6 as a required depedency of libgeotiff
- Use of the proj.db database to resolve the various coded values, mostly in the GTIFGetDefn() "normalization" function of libgeotiff, instead of using the .csv files previously generated from the EPSG database. Typically an EPSG code identifying a projected CRS is resolved them into its base elements: projection method, projection parameters, base geodetic CRS, etc..
- Complete removal of those .csv files and associated functionality
This work will be merged once the above mention PROJ pull request that affects naming of the new functions of the C API has been merged and taken into account. This gdalbarn branch of libgeotiff has also been succesfully integrated in the internal libgeotiff copy of the GDAL gdalbarn branch.
Regarding GDAL, a maintenance v2.3.3 and feature v2.4.0 versions have been released for the same reasons as above.Most methods of the OGRSpatialReference class have now been re-implemented to rely on the PROJ C API to do queries and state changes, which avoids potentially lossy import / export to WKT 1. Similarly to the libgeotiff work, the ImportFromEPSG() functionality now relies on proj.db, and consequently all EPSG or ESRI related .csv files have been removed from the GDAL data directory. I've also drafted a plan regarding on how to be able to take into account WKT 2 by GDAL raster drivers, and proposed changes regarding how to better handle the gap between the axis order as mandated by the CRS authority (for example latitude first, longitude second for geographic CRS in the EPSG dataset) and the actual order of the values in raster metadata or vector geometries. The first part (use of OGRSpatialReference in raster driver) has been implemented in the gdalbarn branch and the second part is in good progress, with the drivers now advertizing their data axis to CRS axis mapping. The ongoing work is to make OGRCoordinateTransformation use the PROJ API to automatically compute the best transformation pipeline, enabling late-binding capabilities.
23 décembre 2018
10:59
GeoTools Team: GeoTools 19.4 released
sur Planet OSGeo

The GeoTools team is pleased to announce the release of GeoTools 19.4: geotools-19.4-bin.zip geotools-19.4-doc.zip geotools-19.4-userguide.zip geotools-19.4-project.zip maven repository This release is the last community sponsored maintenance release for the 19.x series and as such users and downstream projects should consider moving to the 20.x series. This release is made in conjunction
10:42
GeoServer Team: GeoServer 2.13.4 released
sur Planet OSGeo
We are happy to announce the release of GeoServer 2.13.4. Downloads are available (zip, war, and exe) along with docs and extensions.

This is a maintenance release recommended for production use (for newer projects please use the 2.14.x series, as this is last community sponsored 2.13.x release).
This release is made in conjunction with GeoTools 19.4 and GeoWebCache 1.13.4.

Highlights of this release are featured below, for more information please see the release notes (2.13.4 | 2.13.3 | 2.13.2 | 2.13.1 | 2.13.0 | 2.13-RC1 | 2.13-beta).
Improvements and Fixes
- Significant speedup in WMTS capabilities document generation when source data is not in WGS84
- GeoPackage improvements, speed up raster reading and reprojection, make sure global coverage GeoPackage can be displayed
- Removed gml:id attribute from GML 3.1 encoding (was added in previous releases by mistake)
- NetCDF output format improvement, fixing projection coefficients in some cases, make sure data packing takes into account all data slices when calculating linear transformation coefficients
- WCS 2.0 compliance improvements, DescribeCoverage output schema compliance fixes when using GeoServer extensions, returning correct exception type when subsetting request is out of range on time/elevation/custom dimension axis
- KML request handling fixes, support for sortBy and CQL filtering on layer groups
- GWC related fixes, could not modify blobstore configuration without changing its name too, wrong axis order in configuration when saving EPSG:4326 gridset configurations
- Styling related improvements, cannot read dynamic SLD served by an apache server, SLD 1.1 posted to the REST API resulted in garbled content
- And various others, please see the release notes for details
About GeoServer 2.13 Series
Additional information on the 2.13 series:
- Isolated workspaces (User Guide)
- Coverage views from heterogeneous bands (User Guide)
- State of GeoServer 2.13 (slideshare)
- See the GeoServer 2.13.0 released announcement for visual guide to new features
20 décembre 2018
1:00
PyWPS team: PyWPS-4.2.0 Released
sur Planet OSGeo
PyWPS-4.2.0 Released
At long last, a new PyWPS release is out. A new chapter in the PyWPS 4 story is now complete and as 2018 comes to an end, it is time to look back at what happened since the major release in December of 2016.

PyWPS remains fully reliant on pro bono contributions without any sort of commercial organisation around it. This means the project largely depends on how the community is able to contribute at any given time. However, it also means there are no clients to attend or shareholders to appease, no hard deadlines or sales demos. PyWPS evolves organically, to the tune of the community’s wishes.

Relying on voluntary contributions does not mean progress is slow, much to the contrary. In recent times PyWPS has had more simultaneous contributors than at any moment before in its history. This is largely due to the major step taken with PyWPS 4 towards modern Python technologies, but also to regular contributions from the Google Summer of Code (GSoC) programme and various students thesis. During these two years the Project Steering Committe (PSC) was expanded to six members, and last May PyWPS finally graduated as an OSGeo project.

In order to manage development, the Project Steering Committe (PSC) meets regularly (usually at FOSS4G events) and performs basic housekeeping: issues are reviewed and split among a series of milestones that set out a development plan. As development progresses, new issues naturally pop up, and there is always a patch or an improvement that developers wish to include in the next release. Like the Red Queen from Alice in Wonderland, the backlog for the next release seems to always have as many issues resolved as new ones to resolve. Time drags on and the perfect tidy release remains ellusive. So in the end PyWPS 4.2 is mostly a commemorative release, celebrating two years since the major release and making official all the improvements introduced since.

What are then the main improvements introduced with PyWPS 4.2? Foremost is the progressive improvement of the class model, step by step making the library more solid, agile and easier to maintain. Among other things, the input/output handling classes were refactored, the test battery was expanded and streamlined with continuous integration tools. Other relevant improvements include:
- Introduction of Jinja2 templates to automate the generation of outputs.
- Support for process execution in HPC environments, using Slurm or GridEngine.
- Support for streamed URL-based data inputs (e.g. OpenDAP).
- Automated process documentation with Sphinx.
- New validators for JSON, DODS links and netCDF formats.
And along the way, many bugs were fixed.

The organic nature of the project has meant that the goals set out by the PSC are not always met. For instance, PyWPS 4.2 still does not fully support version 2 of the WPS standard. But such is the nature of community driven projects: evolution goes along the needs of the community and not necessarily according to grand plans. There are instead other great improvements in the pipeline, like the containerisation of processes, or the automatic storage of outputs in a database. Please get involved if these developments interest you, or if you would like other improvements to be introduced.

A word also regarding the WSGI applications provided as example services using PyWPS. pywps-flask remains the most popular but a backlog of small issues has built up. During the holiday season developers will be lending attention to these issues. Usability issues with the containers of both pywps-flask and pywps-django will also be addressed during this time.

Finally the PSC would like to thank everyone that has contributed to PyWPS these past two years: from those doing the code heavy lifting, to testers, bug reporters, documentation contributors. Every bit of contribution is precious. Thank you all for making PyWPS what it is today!

Merry Christmas and a Happy New Year.
19 décembre 2018
13:10
gvSIG Team: Free Course: GIS for Archaeologists
sur Planet OSGeo
In recent years, Geographic Information Systems (GIS) have radically evolved the way in which spatial data are collected, analyzed and interpreted. Because of their characteristics, GIS have become useful applications for archaeologists. The appearance of gvSIG, with a complete suite of GIS solutions in open source software, has meant that more and more professionals in the sector are adopting gvSIG as a working tool.

However, there is still a wide ignorance of the potential of GIS in archaeology, and there are only a few free teaching materials that allows training in its use. For these reasons, from the gvSIG Association, we launch a course that will allow to know of these technologies. A free course, with open source software (gvSIG Desktop).

In this first post we want to present the main characteristics of the course, together with the first video-tutorial. As of January 2019, we will publish all the modules weekly until the agenda is completed. To follow the course, it is recommended that participants have prior knowledge about how archaeological projects are carried out.
Registration of the course and how it works
It is not necessary to register anywhere. The course is online and can be followed from anywhere in the world. The course will be available in English and Spanish.

Each week we will publish a post on this blog that will contain a video-tutorial with exercises and access to the course data. Therefore, you only need to follow the video-tutorials to complete the course.

For doubts or problems in the use of software, or about the exercises, you can use the gvSIG users mailing list:

[www.gvsig.com]
List of topics
The list of topics is composed of 7 modules. Each module will contain a video-tutorial with practical exercises. The data to follow the course can be downloaded in a .zip file from this link.
1. Preparation of an archaeological project: views, layers and tables
2. Create geographic information: Alphanumeric information and digitalization
3. Prospecting results analysis: Vector geoprocesses
4. Managing old maps: Georeferencing of images
5. Reservoir territorial analysis: Raster geoprocesses and digital terrain models (MDT)
6. Hydrological analysis to locate materials or archaeological sites
7. Graphical output of an archaeological project: Maps
And if you want to get a certificate?
Such as in other previous courses that we have published, many of you asked us about the possibility to obtain a certificate. To cover this possibility, at the end of the course we will publish a practical exercise to be solved, and to demonstrate that you have acquired the knowledge imparted during the course. This exercise is completely optional, only for those who want to get the certification.

To pay for the evaluation time of the exercises and the issuance of the certificate, apart from delivering the correctly solved exercise, 25 € will have to be paid. The way of payment will be indicated at the end of the course.
Module 1: Preparation of an archaeological project: Views, layer and tables
Ready? Well, let’s start with the first topic, where besides learning where to download gvSIG Desktop – the software that we will use for the course – we will see how to prepare an archaeological project working with the main components of a GIS: Views, layers and tables.
12:56
gvSIG Team: Curso Gratuito: SIG aplicados a Arqueología
sur Planet OSGeo
En los últimos años, los Sistemas de Información Geográfica (SIG) han evolucionado de forma radical la forma en que los datos espaciales se recopilan, analizan e interpretan. Por sus características los SIG se han convertido en aplicaciones de gran utilidad para los arqueólogos. La aparición de gvSIG, con una suite completa de soluciones SIG en software libre, ha hecho que en la actualidad sean cada vez más los profesionales del sector que adoptan gvSIG como herramienta de trabajo.

Sin embargo aún hay un amplio desconocimiento del potencial de los SIG en arqueología, y poco acceso a material educativo y gratuito que permita formarse en su uso. Por estos motivos, desde la Asociación gvSIG, lanzamos un curso que pretendemos que ayude a dar a conocer estas tecnologías. Un curso gratuito, con software libre y gratuito (gvSIG Desktop).

En este primer post queremos presentaros las características principales del curso, junto al primer vídeo-tutorial. A partir de enero de 2019 iremos publicando semanalmente todos los módulos hasta completar el temario. Para aprovechar al máximo el curso se recomienda que los participantes tengan un conocimiento previo de cómo se llevan a cabo los proyectos arqueológicos.
Inscripción y funcionamiento del curso
No es necesario inscribirse en ningún sitio. El curso es online y se podrá realizar desde cualquier parte del mundo. El curso estará disponible en inglés y castellano.

Cada semana publicaremos un post en este blog que tendrá un vídeo-tutorial con ejercicios y acceso a los datos del curso. Por tanto, sólo necesitáis ir siguiendo los vídeo-tutoriales para completar el curso.

Para dudas en el manejo del software, realización de los ejercicios y problemas que puedan surgir en la realización del curso, podéis utilizar la lista de usuarios de gvSIG:

[www.gvsig.com]
Temario
El temario se compone de 7 módulos. Cada módulo presentará un vídeo-tutorial con ejercicios prácticos. Los datos para poder seguir el curso los podéis descargar comprimidos en un fichero .zip de este enlace.
1. Preparación de un proyecto arqueológico: vistas, capas y tablas
2. Crear información geográfica: Información alfanumérica y digitalización
3. Análisis de resultados de prospecciones: Geoprocesos vectoriales
4. Gestionando mapas antiguos: Georreferenciación de imágenes
5. Análisis territorial de yacimientos: Geoprocesos ráster y modelos digitales de terreno (MDT)
6. Análisis hidrológico para localizar materiales o yacimientos arqueológicos
7. Salida gráfica de un proyecto arqueológico: Mapas
¿Y si quiere obtener un certificado?
Como en otros cursos anteriores que hemos publicado, muchos de vosotros nos preguntáis por la posibilidad de obtener un certificado. Para cubrir esta posibilidad, al final del curso publicaremos un ejercicio práctico a resolver que permitirá demostrar que habéis adquirido los conocimientos impartidos durante el curso. La realización de este ejercicio es totalmente opcional. Sólo para aquellos que quieran tener la certificación.

Para costear el tiempo de evaluación de los ejercicios y la emisión del certificado, además de entregar el ejercicio correctamente resuelto habrá que abonar 25 €. La forma de pago se indicará junto al ejercicio final.
Tema 1: Preparación de un proyecto arqueológico: vistas, capas y tablas
¿Preparados? Pues vamos a comenzar con el primer tema, donde además de aprender de dónde descargar gvSIG Desktop – el software que utilizaremos para el curso – veremos como preparar un proyecto arqueológico trabajando con los principales componentes de un SIG: Vistas, capas y tablas.
12:55
Fernando Quadro: Planet Labs adquire Boundless
sur Planet OSGeo

A Planet Labs anunciou ontem a aquisição da Boundless Spatial, Inc., empresa de soluções de software geoespacial com sede em St. Louis, que é a responsável pelo desenvolvimento do GeoServer. A aquisição expande os negócios comerciais da Planet com o governo dos EUA e clientes da agricultura comercial. A Boundless é líder em software de gerenciamento de dados geoespaciais e está alinhada com o objetivo de fornecer novos serviços de assinatura de dados geoespaciais e acelerar a adoção em clientes corporativos.

Com a aquisição da Boundless, a Planet pretende alavancar a tecnologia e o conhecimento da Boundless implantados dentro do governo dos EUA para facilitar ainda mais a utilização dos seus produtos. Além disso, este acordo também trará novos clientes do governo dos EUA para Planet.

A Planet manterá o escritório da Boundless em St Louis (MO) e em conjunto com a equipe da Planet que trabalha diretamente com o governo dos EUA formará uma subsidiária chamada Planet Federal.

Sobre os outros escritórios da Boundless nos EUA, e também sobre os funcionários que trabalham remotamente nada foi dito!

Fonte: Planet
18 décembre 2018
12:00
MapTiler: OpenMapTiles 3.9
sur Planet OSGeo
The new release of OpenMapTiles comes with many smaller changes. Version 3.9 brings rendering of docks and piers, support for multilingual street names and other improvements.
Docks and Piers
The biggest visual change is the newly added support for rendering docks and piers. Piers are now rendered in a similar way like bridges while docks were added to the water layer. This modification will have the biggest impact on complicated harbors structures like Liverpool’s harbor or spectacular structures like the famous Brighton Palace Pier.

Multilingual street names
To solve a pain in bilingual countries we added support for multilingual streets. It follows all other map elements like place names, country names or named points of interest. Street name labels are now rendered in all 58 languages and the selection of displayed language is done automatically based on end-user preferences or by an app developer decision.

Many changes for multilingual street names has been made by the Helsinki Regional Transport Authority. They use our maps on their websites, in the application as well as printed on the public transport stops. Another user of OpenMapTiles and contributor to the project is the Qwant map team.

Talking about languages, improvement have been made for Vietnamese and Azerbaijani nonlatin labels and there is brand new support for the Malayalam language.
Other changes in OpenMapTiles 3.9
- National parks in the Park layer have their names.
- Indoor POI layer contains layer and level attributes.
- The type of sport for a pitch POI can be distinguished via subclass attribute.
- Improved rendering of 3D buildings.
- Fixed issues with small islands in Scandinavia.
For detailed information, read the full changelog.

View all improvements on our map and start using it via MapTiler Cloud or self-hosted OpenMapTiles.
15 décembre 2018
16:31
Free and Open Source GIS Ramblings: Movement data in GIS #18: creating evaluation data for trajectory predictions
sur Planet OSGeo
We’ve seen a lot of explorative movement data analysis in the Movement data in GIS series so far. Beyond exploration, predictive analysis is another major topic in movement data analysis. One of the most obvious movement prediction use cases is trajectory prediction, i.e. trying to predict where a moving object will be in the future. The two main categories of trajectory prediction methods I see are those that try to predict the actual path that a moving object will take versus those that only try to predict the next destination.

Today, I want to focus on prediction methods that predict the path that a moving object is going to take. There are many different approaches from simple linear prediction to very sophisticated application-dependent methods. Regardless of the prediction method though, there is the question of how to evaluate the prediction results when these methods are applied to real-life data.

As long as we work with nice, densely, and regularly updated movement data, extracting evaluation samples is rather straightforward. To predict future movement, we need some information about past movement. Based on that past movement, we can then try to predict future positions. For example, given a trajectory that is twenty minutes long, we can extract a sample that provides five minutes of past movement, as well as the actually observed position five minutes into the future:

But what if the trajectory is irregularly updated? Do we interpolate the positions at the desired five minute timestamps? Do we try to shift the sample until – by chance – we find a section along the trajectory where the updates match our desired pattern? What if location timestamps include seconds or milliseconds and we therefore cannot find exact matches? Should we introduce a tolerance parameter that would allow us to match locations with approximately the same timestamp?

Depending on the duration of observation gaps in our trajectory, it might not be a good idea to simply interpolate locations since these interpolated locations could systematically bias our evaluation. Therefore, the safest approach may be to shift the sample pattern along the trajectory until a close match (within the specified tolerance) is found. This approach is now implemented in MovingPandas’ TrajectorySampler.
```
def test_sample_irregular_updates(self):
    df = pd.DataFrame([
        {'geometry':Point(0,0), 't':datetime(2018,1,1,12,0,1)},
        {'geometry':Point(0,3), 't':datetime(2018,1,1,12,3,2)},
        {'geometry':Point(0,6), 't':datetime(2018,1,1,12,6,1)},
        {'geometry':Point(0,9), 't':datetime(2018,1,1,12,9,2)},
        {'geometry':Point(0,10), 't':datetime(2018,1,1,12,10,2)},
        {'geometry':Point(0,14), 't':datetime(2018,1,1,12,14,3)},
        {'geometry':Point(0,19), 't':datetime(2018,1,1,12,19,4)},
        {'geometry':Point(0,20), 't':datetime(2018,1,1,12,20,0)}
        ]).set_index('t')
    geo_df = GeoDataFrame(df, crs={'init': '4326'})
    traj = Trajectory(1,geo_df)
    sampler = TrajectorySampler(traj, timedelta(seconds=5))
    past_timedelta = timedelta(minutes=5)
    future_timedelta = timedelta(minutes=5)
    sample = sampler.get_sample(past_timedelta, future_timedelta)
    result = sample.future_pos.wkt
    expected_result = "POINT (0 19)"
    self.assertEqual(result, expected_result)
    result = sample.past_traj.to_linestring().wkt
    expected_result = "LINESTRING (0 9, 0 10, 0 14)"
    self.assertEqual(result, expected_result)
```
The repository also includes a demo that illustrates how to split trajectories using a grid and finally extract samples:
14 décembre 2018
23:04
Cameron Shorter: Catching the elusive Episodic Volunteer
sur Planet OSGeo
Understanding the science of volunteer contributions is key to the success of Open Source projects. A comprehensive study of episodic volunteering, which drew upon our OSGeoLive experiences, provides insightful advice for anyone wishing to build sustained and successful Open Source Software:

Contributor MotivationFindings: Episodic volunteers with intrinsic motives are more likely to intend to remain, compared to episodic volunteers with extrinsic motives.

Recommendations:
1. Lower barriers to entry through:
2. Offer guided introductory events to help newcomers get started and to introduce the social element.
3. Provide opportunities for social interactions, such as:
Social NormsFindings: Although Open Source episodic volunteers were unlikely to see their participation as influenced by social norms, personal invitation was a common form of recruitment, especially among non-code contributors.

Recommendations:
1. Encourage existing volunteers to talk about their Open Source involvement by:
  - Highlighting the benefits of advocating broadly,
  - Providing digestible information for sharing.
Psychological Sense of CommunityFindings: Psychological sense of community is more common among long-term participants;
A policy of inclusion was a commonly mentioned reason for feeling welcomed in the community.

Recommendations:
1. Use a code of conduct to express the community’s intentions, allowing potential episodic volunteers to determine their similarity to the community.
2. Give potential episodic volunteers the opportunity to identify alignment with the community through awareness of non-coding activities:
  - Collaborate with organisations with a different focus but shared values,
  - Recognise all forms of contribution.
3. Re-enforce the psychological sense of community by:
SatisfactionFindings: Satisfaction was most commonly cited as a reason to remain;
Episodic volunteers derive satisfaction from knowing that their work is used, enjoying the work itself, and feeling appreciated.
Recommendations:
1. Encourage satisfaction by increasing feelings of appreciation, by recognising all contributors and their areas of expertise.
2. Being aware of episodic volunteers’ areas of expertise and requesting their assistance, sparingly, can:
Community CommitmentFindings: Episodic volunteers who talk about their involvement are more inclined to continue participating;
Long-term episodic volunteers often have community commitment;
Community commitment is less common among episodic volunteers with extrinsic motives.

Recommendations:
1. Encourage long-term episodic volunteers to talk about the community to strengthen their commitment to the community and:
2. Consider time-based releases for large projects to allow episodic volunteers to plan their return.
3. Use opt-in platforms to broadcast calls for participation for specific tasks to encourage episodic volunteers to return.
Episodic volunteeringFindings: Episodic volunteering is widespread in Open Source communities, but Open Source communities are often not strategically engaging with episodic volunteers;
Open Source episodic volunteers are often habitual volunteers in other communities.

Recommendations:
1. Evaluate volunteer assets, volunteer availability, and potential assignments to find opportunities for episodic volunteers.
13 décembre 2018
21:45
Stefano Costa: I libri che ho letto nel 2017
sur Planet OSGeo

Quest’anno mi porto avanti di qualche settimana rispetto al ritardo astronomico maturato negli anni e pubblico la lista dei pochi libri che ho letto nel 2017.

Il 2017 è stato l’anno in cui sono diventato papà, in cui ho traslocato dalla città dove è nato mio figlio a quella dove sono nato io ma dove non avevo quasi mai abitato, in cui ho cambiato lavoro per poter finalmente continuare a fare il mio mestiere. Quindi, scusate se ho letto troppo poco. Anche nel 2017 non ho letto nessun libro che non fosse scritto in italiano (mica intenzionalmente).

Ho anche iniziato a tenere traccia dei libri su inventaire.io.
Piero Colaprico, Trilogia della città di M.
Questo libro mi è stato donato da Stefano R., nativo della città di M., quando ho iniziato a lavorare nella città di M. il 13 dicembre 2017. L’ho consumato nei vagoni ferroviari prima dell’alba e dopo il tramonto viaggiando quotidianamente da Genova.

Le tre storie in cui si muove l’ispettore Bagni sono ambientate in una città di mezzo, tra quella vecchia che sta scomparendo e lascia il tempo ad una nuova. Muovendomi nella zona del Ticinese più volte ho provato la sensazione di essere stato in quei luoghi. Non sono un assiduo frequentatore della letteratura di genere, ma se leggo un poliziottesco, ora preferisco Bagni ai suoi colleghi più televisivi.
Wu Ming 1, Un viaggio che non promettiamo breve
Ho iniziato a leggere questo libro verso fine anno, proprio come la Trilogia ad inizio articolo, viaggiando in treno lungo uno dei percorsi che dovrebbero essere toccati dalle Grandi Opere. E poiché il libro racconta con estrema precisione la realtà, ho faticato ad ogni pagina a reggere il peso di tutte le nefandezze, e per adesso non ho letto che un terzo del libro, forse un po’ meno.

È un libro che deve essere letto, perché è una storia che non è ancora conclusa e da queste parti oltre alla inutile e dannosa TAV abbiamo un altrettanto inutile terzo valico ferroviario, mentre le linee ferroviarie già esistenti del pendolarismo quotidiano sono in frantumi.
Edmund De Waal, La strada bianca
Questo è un libro di cui mi sono innamorato vedendo la copertina in vetrina nella mia libreria. A me piacciono le copertine che raccontano qualcosa del libro, e questa era già di per sé meritevole di essere portata a casa, anche se non ci fosse stato il libro dentro.

Ho letto il primo capitolo con l’emozione di chi scopre un tesoro. L’ho dovuto persino rileggere a voce alta, per poterlo capire in tutta la sua poesia. Questo libro è stato compagno di viaggio tra la fine del 2016 e l’inizio del 2017.}

Non ho mai letto Un’eredità di avorio e ambra, che pure è sulla libreria di casa da anni. Non avevo capito bene di cosa si trattasse, potrei dire.

Per me la ceramica, e quindi anche la porcellana, è un elemento primordiale a cui ho dedicato moltissima parte della mia vita adulta, soprattutto manipolando, guardando, osservando migliaia di vasi rotti. Leggerne una storia raccontata in prima persona è stato un percorso di ricucitura, una forma di kintsugi letterario.
Marco Danielli, Uscita di sicurezza
Marco è prima di tutto un collega, che potrebbe sembrare molto versatile, vista la varietà di ruoli che ha rivestito nella vita in situazioni molto diverse. È stata una grande sorpresa scoprire che è anche un abile scrittore.

Un romanzo d’esordio finemente costruito che combina realtà autobiografica e finzione letteraria.
Igiaba Scego, Adua
Non so se i quattro lettori assidui di questa rubrica ci hanno fatto caso, ma da qualche anno leggo sempre almeno un libro di Africa. Quest’anno, addirittura due.

Leggere di Africa è inevitabilmente una seduta di colonialismo, di post-colonialismo, oltre che di un continente non misurabile, privo di significato a meno di non essere distorto da una proiezione geografica sbagliata. Nel caso clinico italiano credo si possa parlare a buon diritto di post-colonialismo assente – autrici come Igiaba Scego sono ogni giorno lì a smontare il buonismo, soprattutto con storie come quella di Adua, ben più di un romanzo storico.
Chinua Achebe, Le cose crollano

Ricevuto in dono e letto con gusto, questo primo libro della trilogia è un gioiello di grande profondità storica e soprattutto epica. Achebe muove e colpisce i suoi personaggi con l’eleganza di un vero classico. Leggendo è stato impossibile non rispecchiarsi nelle pagine di Congo, così totalmente diverso eppure intrecciato nelle stesse vicende storiche.
11:30
Fernando Quadro: Vídeos do FOSS4G Oceania
sur Planet OSGeo
Entre os dias 20 e 23 de novembro aconteceu em Melbourne (Austrália) o FOSS4G SoTM Oceania um evento regional que cria um link focado na Oceania em uma comunidade global, apoiada pela OSGeo Foundation e pela OpenStreetMap Foundation (OSMF).

No evento foram abordados diversos assuntos interessantes, e vale a pena dar uma olhada nas apresentações, que estão disponíveis gratuitamente no YouTube. Veja:

1. Apresentação (Boas Vindas)
- María Arias de Reyna (OSGeo)
- Kate Chapman (OSMF)
2. Keynotes
3. Sessões especiais
- Lightning talks
- Panel Session: Human Evolution, Ethics and Open Source Mapping
4. Palestras
12 décembre 2018
19:45
gvSIG Team: gvSIG Batoví, más que un concurso una metodología para llevar el SIG a las aulas
sur Planet OSGeo

Hace unas semanas finalizó la segunda edición del concurso denominado ‘Proyectos con estudiantes y gvSIG Batoví’. A modo de titular se podría decir de manera muy simplista que se trata de un concurso que premia los mejores trabajos realizados con gvSIG Batoví (una adaptación de gvSIG Desktop) por estudiantes de secundaria de Uruguay.

Pero este concurso es mucho más, es un ejemplo de como llevar con inteligencia y trabajo los Sistemas de Información Geográfica, el software libre gvSIG, a las aulas. El concurso, a la postre, es el incentivo que permite enganchar a los alumnos, pero en este post quería resaltar las tareas principales que hay detrás.

Fase 0. Antes de empezar
No se puede llevar una nueva tecnología a las aulas sin capacitar a los profesores, y teniendo esto en mente se consideró fundamental el realizar como primer paso unas jornadas de capacitación de docentes. Al finalizar la actividad se conformó un equipo de seguimiento de gvSIG Batoví integrado por los asistentes al taller con el objetivo de monitorear el uso y aplicación de gvSIG Batoví en tareas de docencia y en los diferentes sub-sistemas educativos.

Además, para facilitar el manejo de la herramienta, se realizaron distribuciones de gvSIG Desktop adaptadas al entorno educativo para los equipos del Plan Ceibal (Ubuntu) y particulares (disponible también para Windows) y una serie de manuales y tutoriales relacionados con la instalación y ejecución de gvSIG Batoví.

Y reincidiendo en la importancia de la capacitación a docentes, se diseño una propuesta de b-learning, que es el nombre que se conoce para definir una combinación de formación presencial y online. Más de 150 profesores de todo el país se inscribieron a esta formación. En este enlace podéis consultar la documentación que se generó para servir como guía de los talleres.

Fase 1. Desarrollando los proyectos
Pasando ya al concurso en sí, además de definir las bases y premios, se realizó un convenio con la Facultad de Ingeniería (UdelaR) para la designación de tutores estudiantes de la Tecnicatura de Cartografía para cada equipo participante. Conectando así a los estudiantes universitarios con los de los liceos.

Ya estaba todo listo para empezar a trabajar en los proyectos…y para ello se crearon espacios en la plataforma CREA (de aprendizaje en línea), para cada proyecto, para el seguimiento de los mismos, evacuación de dudas, intercambio de información, etc. El seguimiento se complementó con intercambios entre el equipo coordinador y los equipos mediante videoconferencias.

Fase 2. El concurso
Llegado el momento se procedió por cada uno de los equipos participantes a la defensa de los proyectos mediante videoconferencias. Y por último, una actividad presencial para compartir, de socialización de los resultados, entrega de premios, diplomas y certificados

El grupo que realizó el proyecto galardonado con el primer premio estudió el vínculo existente entre el rendimiento de sus compañeros y compañeras con la distancia de la que viven respecto al liceo. Para ello, referenciaron los lugares donde viven los estudiantes, el medio en el que se trasladan, el rendimiento (materias con nota insuficiente y llegadas tarde) y nivel de repetición. ¡Bravo!
Fase 3. Después de…
No todo acaba con la entrega de premios. Es el momento de la tarea de evaluación y valoración con profesores referentes y tutores.

Por cierto, aquí tenéis el acceso a los proyectos finalistas, aunque en mi opinión con este tipo de iniciativas sólo hay ganadores. Ojalá y veamos replicada esta iniciativa en otras partes del planeta. A mi me encantaría ver algo así en mi tierra ¿y a vosotros?
2:27
Stefano Costa: IOSACal on the web: quick calibration of radiocarbon dates
sur Planet OSGeo

The IOSA Radiocarbon Calibration Library (IOSACal) is an open source calibration software. IOSACal is meant to be used from the command line and installation, while straightforward for GNU/Linux users, is certainly not as easy as common desktop apps. To overcome this inconvenience, I dedicated some efforts to develop a version that is immediately usable.

The IOSACal web app is online at [https:]] .

This is a demo service, so it runs on the free tier of the commercial Heroku platform and it may take some time to load the first time you visit the website. It is updated to run with the latest version of the software (at this time, IOSACal 0.4.1, released in May).

Since it may be interesting to try the app even if you don’t have a radiocarbon date at hand, at the click of a button you can randomly pick one from the open data Mediterranean Radiocarbon dates database, and the form will be filled for you.
The random date picker in action
Unfortunately, at this time it is not possible to calibrate or plot multiple dates in the web interface (but the command-line program is perfectly capable of that).

IOSACal Web is made with Flask and the Bootstrap framework, and the app itself is of course open source.

IOSACal is written in the Python programming language and is based on Numpy, Scipy and Matplotlib. This work wouldn’t be possible without the availability of such high quality programming libraries.
11 décembre 2018
9:00
GIS for Thought: PostGIS Tutorial Data
sur Planet OSGeo

The OpenGeo (now called Boundless) PostGIS tutorial has long been my recommendation for people just getting started with PostGIS. It provides a broad intro to PostGIS, great for users new to SQL.

Unfortunately they recently removed it from their website. The guide was a fairly old now, so I don’t blame them. But it was a great free resource.

I have created a version of the data used, which can be accesses at (choose the green Clone or download button and Download ZIP):

[https:]]

The data is also not 100% the same as the original. So some of the queries will have to be modified. Particularly the neighborhoods are different. So they may have slightly different names/boundaries/capitalization. So keep that in mind when creating your queries, and interpreting the results. But overall the data is pretty similar

The course has been compiled and can be found at: [https:]]
10 décembre 2018

11:30

Fernando Quadro: Desenvolvimento de plugins do QGIS: Testando seu código

sur Planet OSGeo

Neste post (o último da série) veremos o teste do plugin que criamos para o QGIS.

1. MOCKING

Para testar um plugin do QGIS você precisa simular o ambiente que ele deve rodar. E esse ambiente é obviamente o próprio QGIS, mas não é possível iniciar o QGIS toda vez que você executa um teste. Felizmente, há o QGIS Mocking para lhe auxiliar.

2. VOCÊ PRECISA DE DADOS

Cada teste é executado repetidas vezes, o que significa que ele precisa redefinir os dados que estão sendo usados para seu estado padrão. Isso pode ser um “PIDA” se o teste alterar os dados de maneira imprevisível.

Usando as camadas do QGIS, você pode preparar novos dados para cada um dos seus testes, colocando de forma eficaz todo o processo de manipulação de dados.

3. ESCREVENDO TESTES

Cada um dos testes do plugin AttributeTransfer herda da classe unittest.TestCase, que vem com vários métodos que você pode estar familiarizado com outros idiomas: setUp() é executado antes para cada método de teste, enquanto tearDown() é executado após cada um deles. Os testes são definidos como métodos cujos nomes começam com a palavra test.

Cada teste deve chamar algum método assertWhatever que verifica se o teste foi aprovado ou não. Aqui está um exemplo de tal teste cobrindo camadas não pontuais.

#!/usr/bin/env python
# -*- coding: utf-8 -*-
# @Date    : 2017-11-18 18:40:50
# @Author  : Michal Zimmermann <zimmicz@gmail.com>

import os
import sip
import sys
import unittest
from qgis.core import QgsMapLayerRegistry, QgsVectorLayer, QgsFeature, QgsGeometry, QgsPoint
from utilities import get_qgis_app

sys.path.append(os.path.dirname(os.path.abspath(__file__)) + "/..")
from attribute_transfer import AttributeTransfer
from create_dummy_data import create_dummy_data_polygon_or_line

sip.setapi('QtCore', 2)
sip.setapi('QString', 2)
sip.setapi('QDate', 2)
sip.setapi('QDateTime', 2)
sip.setapi('QTextStream', 2)
sip.setapi('QTime', 2)
sip.setapi('QUrl', 2)
sip.setapi('QVariant', 2)

QGIS_APP = get_qgis_app()
IFACE = QGIS_APP[2]


class AttributeTransferTestPolygonOrLine(unittest.TestCase):

    def setUp(self):
        self.source_layer = QgsVectorLayer(
            "Polygon?crs=epsg:4326&field=id:integer&field=textAttr:string&field=intAttr:integer&field=decAttr:double&field=dateAttr:date&index=yes", "source layer", "memory")
        self.target_layer = QgsVectorLayer(
            "Linestring?crs=epsg:4326&field=id:integer&index=yes", "target layer", "memory")
        self.widget = AttributeTransfer(IFACE)

        registry = QgsMapLayerRegistry.instance()
        registry.removeAllMapLayers()
        registry.addMapLayers([self.source_layer, self.target_layer])
        create_dummy_data_polygon_or_line(self.source_layer, self.target_layer)
        self.widget.initGui()
        self.widget.vectors = [self.source_layer, self.target_layer]
        self.widget.editable_vectors = [self.source_layer, self.target_layer]
        self.widget.dlg.sourceLayer.addItems(["source layer", "target layer"])

    def test_text_attr(self):
        ATTRIBUTE_NAME = "textAttr"
        ATTRIBUTE_INDEX = 1

        self._test_attr(ATTRIBUTE_NAME, ATTRIBUTE_INDEX)

    def test_int_attr(self):
        ATTRIBUTE_NAME = "intAttr"
        ATTRIBUTE_INDEX = 2

        self._test_attr(ATTRIBUTE_NAME, ATTRIBUTE_INDEX)

    def test_dec_attr(self):
        ATTRIBUTE_NAME = "decAttr"
        ATTRIBUTE_INDEX = 3

        self._test_attr(ATTRIBUTE_NAME, ATTRIBUTE_INDEX)

    def test_date_attr(self):
        ATTRIBUTE_NAME = "dateAttr"
        ATTRIBUTE_INDEX = 4

        self._test_attr(ATTRIBUTE_NAME, ATTRIBUTE_INDEX)

    def test_existing_attr(self):
        ATTRIBUTE_NAME = "id"
        ATTRIBUTE_INDEX = 0

        self.widget.dlg.sourceAttribute.setCurrentIndex(ATTRIBUTE_INDEX)
        self.widget.dlg.targetAttribute.setText(ATTRIBUTE_NAME)

        self.assertEqual(
            self.widget.dlg.sourceAttribute.currentText(), ATTRIBUTE_NAME)
        self.assertFalse(self.widget.transfer())

    def _test_attr(self, attr_name, attr_index):
        self.widget.dlg.sourceAttribute.setCurrentIndex(attr_index)
        self.widget.dlg.targetAttribute.setText(attr_name)

        self.assertEqual(
            self.widget.dlg.sourceAttribute.currentText(), attr_name)

        self.widget.transfer()

        target_fields = [f.name()
                         for f in self.target_layer.dataProvider().fields()]
        self.assertIn(attr_name, target_fields)

        source_features = [f for f in self.source_layer.getFeatures()]
        target_features = [f for f in self.target_layer.getFeatures()]

        for idx, f in enumerate(source_features):
            self.assertEqual(f.attribute(attr_name), target_features[
                             idx].attribute(attr_name))


if __name__ == "__main__":
    unittest.main()

Com esse teste encerro a série de posts sobre como desenvolver um plugin para o QGIS. Espero que tenha ajudado a você que está querendo criar ou está criando o seu plugin para o QGIS.

Esta série de posts é tradução e uma adaptação livre do post escrito originalmente por Michal Zimmermann do blog “Pieces of knowledge from the world of GIS“.

Fonte: Pieces of knowledge from the world of GIS

06 décembre 2018
11:30
Fernando Quadro: Desenvolvimento de plugins do QGIS: Plugin AttributeTransfer
sur Planet OSGeo

Nesse link você pode acessar todo código-fonte do plugin QGIS AttributeTransfer.

O plugin em si se encontra no arquivo attribute_transfer.py. Quando o método run() é invocado, o formulário QT aparece com combos pré-preenchidos com camadas vetoriais disponíveis que suportam a edição de atributos.

Os combos de camada de origem e destino são mutuamente exclusivos, portanto, não é possível transferir o atributo dentro da mesma camada.

Codificando o plugin, é possível se deparar com pequenos problemas relacionados principalmente à implementação da QgsSpatialIndex. Na parte de análise do vizinho mais próximo, o método QgsSpatialIndex.nearestNeighbor foi mencionado. No entanto, esse método funciona apenas com geometrias do tipo QgsPoint. Com ele é impossível de obter geometrias do tipo QgsPolygon ou QgsPolyline, no entanto. O que podemos fazer, para lidar com esse problema? Bem… desenhe uma matriz de solução.

Usar o índice espacial traz mais um problema que foi percebido logo após a implementação dos fluxos de trabalho de comparação espacial para diferentes tipos de geometria. Há uma chance de encontrar o recurso mais próximo usando a caixa delimitadora que, na verdade, não é o recurso mais próximo. Nesse caso, optou-se por encontrar o vértice mais distante de tal recurso e utilizá-lo para construir o retângulo ao redor do recurso de destino. Se houver algum recurso de origem em tal retângulo, é muito provável que um deles seja o recurso real mais próximo.

Você pode ter pensado que esta seria a última parte da série. Mas como alguém poderia reivindicar o projeto de codificação sem escrever testes? Fique ligado no próximo post.
05 décembre 2018
11:30
Fernando Quadro: Desenvolvimento de plugins do QGIS: Criando interface gráfica
sur Planet OSGeo

Depois de mexer no console do QGIS Python e implementar a análise de vizinho mais próximo, vamos criar uma GUI muito simples para o plug-in.

Embora os documentos da API do QGIS levassem algumas horas para serem entendidos, o ecossistema do PyQGIS me ajudou muito. Aqui vem a lista de ferramentas que você deve incorporar ao seu processo de desenvolvimento o mais rápido possível.

1. CONSTRUTOR DE PLUG-INS

O QGIS Plugin Builder é um plugin criado para criar outros plugins. Ele permite você criá-lo em minutos e que você codifique em vez de configurar coisas que você não quer configurar. Note que você deve colocar o plugin dentro da pasta de plugins do QGIS (o padrão é ~/.qgis2/python/plugins) no Linux.

Lembre-se de rodar pyrcc4 -o resources.py resources.qrc dentro da pasta do seu plugin antes de adicioná-lo ao QGIS.

2. PLUGIN RELOADER

O QGIS Plugin Reloader é um plugin (possivelmente criado com o QGIS Plugin Builder) que permite que você atualize (reload) o seu plugin enquanto você codifica. Nenhuma reinicialização do QGIS é necessária.

3. QT DESIGNER

O Qt Designer vem com o pacote qt4-designer no Ubuntu. Ele é feito sob medida para projetar e construir GUIs a partir de componentes do Qt que podem ser usados dentro do QGIS. Sua interface de arrastar e soltar permite prototipar rapidamente.

Graças ao Plugin Builder, você pode carregar o arquivo attribute_transfer_dialog_base.ui diretamente no Qt Designer e ajustá-lo às suas necessidades.

Não é preciso muito, apenas um QLineEdit e alguns QComboBoxwidgets. Eles estarão disponíveis no arquivo attribute_transfer.py como propriedades da classe AttributeTransferDialog.

O nome do widget pode ser personalizado na barra lateral direita e eu aconselho você a fazer isso. Eu escolhi o seguinte:

Uma vez carregado com Plugins -> Gerenciar e Instalar Plugins -> AttributeTransfer, o plugin está disponível diretamente na barra de ferramentas ou no menu Vetor. Está faltando a lógica de negócios, mas nós já fizemos isso no post anterior. Tudo o que deve ser fazer é unir essas duas partes.
04 décembre 2018
11:30
Fernando Quadro: Desenvolvimento de plugins do QGIS: Vizinho mais próximo
sur Planet OSGeo
No último post foi descrito o básico da manipulação de camadas vetoriais. Como o nosso objetivo é criar um plug-in personalizado totalmente funcional com base em uma distância, gostaria de discutir a indexação espacial e a análise do vizinho mais próximo.

A figura acima ilustra a tarefa que pode ser resolvida apenas usando a API do QGIS. Imagine que você tenha uma camada de origem com um atributo preenchido com valores. Você recebe uma camada de destino também, infelizmente, os valores nesta camada estão faltando (não tão raros no mundo GIS, certo?). No entanto, você sabe que o valor ausente de cada recurso na camada de destino pode ser preenchido pelo valor de seu vizinho mais próximo da camada de origem. Como você faz isso?

1. GERANDO DADOS FICTÍCIOS

Vamos criar dois conjuntos de dados na memória com atributos id e value. Ambos terão dez recursos.
```
from qgis.core import QgsMapLayerRegistry, QgsVectorLayer, QgsFeature, QgsGeometry, QgsPoint, QgsSpatialIndex
from qgis.utils import iface

source_layer = QgsVectorLayer("point?crs=epsg:4326&field=id:integer&field=value:integer", "Source layer", "memory")
target_layer = QgsVectorLayer("point?crs=epsg:4326&field=id:integer&field=value:integer", "Target layer", "memory")

def create_dummy_data():

    source_layer.startEditing()
    target_layer.startEditing()

    feature = QgsFeature(source_layer.pendingFields())

    for i in range(10):
        feature.setGeometry(QgsGeometry.fromPoint(QgsPoint(i, i)))
        feature.setAttribute("id", i)
        feature.setAttribute("value", i)
        source_layer.addFeature(feature)

    feature = QgsFeature(source_layer.pendingFields())

    for i in range(10):
        feature.setGeometry(QgsGeometry.fromPoint(QgsPoint(i + i, i)))
        feature.setAttribute("id", i)
        target_layer.addFeature(feature)

    source_layer.commitChanges()
    target_layer.commitChanges()

    QgsMapLayerRegistry.instance().addMapLayer(source_layer)
    QgsMapLayerRegistry.instance().addMapLayer(target_layer)

create_dummy_data()
```
2. ESCREVENDO VALORES DO VIZINHO MAIS PRÓXIMO

A análise real do vizinho mais próximo pode ser feita em dez linhas de código! Primeiro, o qgis.core.QgsSpatialIndex é construído a partir de todos os recursos source_layer. Então, você itera sobre os recursos target_layer e para cada um deles, obtém apenas um ( nearestNeighbor(f.geometry().asPoint(), 1)[0]) vizinho mais próximo. Por fim, você apenas grava o valor do atributo do vizinho mais próximo na camada de destino e confirma as alterações.
```
def write_values_from_nn():
    source_layer_index = QgsSpatialIndex(source_layer.getFeatures())
    source_layer_features = {feature.id(): feature for (feature) in source_layer.getFeatures()}
    target_layer_features = target_layer.getFeatures()

    target_layer.startEditing()

    for f in target_layer_features:
        nearest = source_layer_index.nearestNeighbor(f.geometry().asPoint(), 1)[0]
        value = source_layer_features[nearest].attribute("value")
        target_layer.changeAttributeValue(f.id(), 1, value)

    target_layer.commitChanges()

write_values_from_nn()
```
3. O QUE VEM POR AI

Estamos um pouco mais perto do nosso objetivo. O que esta faltando?
- Verificações de capacidades: a camada de destino suporta edições? Verifique os recursos do provedor de dados da camada para descobrir.
- Registro do usuário: avisos ou erros estão completamente ausentes. Será ótimo tê-los apresentados dentro do qgis.gui.QgsMessageBar.
- Atributos customizados: essa versão espera que as duas camadas tenham o mesmo atributo com o mesmo tipo de dados.
- GUI : um widget PyQt transformará esse script baseado em console em um plug-in personalizado. Isso é o que vamos ver no próximo post.
10:00
GIS for Thought: Copying Rasters in PostGIS
sur Planet OSGeo
I ran into a process where I wanted to create copies of rasters in PostgreSQL. While seemingly a simple process this took me a bit of work to figure out.

For my workflow I had three rasters, which all have the same size, and I want to load them into the same PostGIS table with three raster geometry columns. I don’t think this will work for different sized rasters since the rid’s will not match.

Three rasters:
raster1 raster2 raster3

Which I want to copy into:
merged_raster

First to create the merged raster table:
```
CREATE TABLE merged_raster
(
  rid serial NOT NULL,
  raster1 raster,
  raster2 raster,
  raster3 raster
);
```
Then to add the rid’s. These are the id’s of the tiles that the raster was split into when loading. If your tile size is large enough then you may only have one.
```
INSERT INTO merged_raster(rid)
(SELECT rid FROM raster1);
```
Then copying the actual data is straighforward (this assumes the raster column in the raster1 datasets is called rast):
```
UPDATE merged_raster m
SET raster1 = r.rast
FROM raster1 r
WHERE r.rid = m.rid;

UPDATE merged_raster m
SET raster2 = r.rast
FROM raster2 r
WHERE r.rid = m.rid;

UPDATE merged_raster m
SET raster3 = r.rast
FROM raster3 r
WHERE r.rid = m.rid;
```
Now I still have an issue that QGIS will not load these layers. It will always load the initial raster column no matter what is chosen.
03 décembre 2018
21:59
QGIS Blog: User question of the Month – Dec 18 & answers from Nov
sur Planet OSGeo

It’s December and that means it is time to plan for the next year. Planning also means preparing a budget and to do so, we would like to learn more about what you think QGIS.ORG should focus on: features or bug fixing and polishing? Therefore, we invite you to our QGIS user question of December 2018.

We also have localized translated versions of this questionnaire for our French-speaking and Portuguese-speaking users.
Your answers in November
In November, we wanted to know which version of QGIS you use.
19:07
QGIS Blog: Call for presentations: QGIS User Conference and Developer Meeting, 2019
sur Planet OSGeo

The next International QGIS User Conference and Developer Meeting will take place in the week from 4 to 10 March 2019 in A Coruña (Spain).

The call for presentations and workshop registration is out and you can apply using the online registration form. Note that the deadline for presenting proposals is 21 Dec 2018. If you need any info please email your queries to: userconf2019@qgis.es

The International QGIS User and Developer Conference wants to be the reference conference for the QGIS community, a meeting point for the family of users and developers associated with the QGIS project. Attending the conference is an opportunity to gather experience and share knowledge about QGIS. The language of the conference is English.

The event is organized by the Spanish QGIS Association [1], the Spanish user group, and the Galician Xeoinquedos community [2] with the help of A Coruña municipality [3]. The event is under the QGIS.org umbrella. We look forward to seeing you there!

[1] [qgis.es] [2] [https:]] [3] [https:]]

11:30

Fernando Quadro: Desenvolvimento de plugins do QGIS: Usando o console

sur Planet OSGeo

Como mencionado na parte anterior da série, o console Python é um ponto de entrada para a automação do fluxo de trabalho GIS dentro do QGIS. Lembre-se de que existe um objeto iface representando a instância qgis.gui.QgisInterface dentro do console que lhe dá acesso a toda a GUI do QGIS. Vamos ver o que podemos fazer dentro do console.

1. CARREGANDO A PASTA DE CAMADAS VETORIAIS

import glob
from qgis.core import QgsMapLayerRegistry, QgsVectorLayer

def load_folder(folder):
    VALID_EXTENSIONS = ('.geojson', '.gpkg', '.shp')
    files = [f for f in glob.glob("{}/*".format(folder)) if f.endswith(VALID_EXTENSIONS)]

    for f in files:
        layer = QgsVectorLayer(f, f.split('/')[-1], 'ogr')

        if not layer.isValid():
            iface.messageBar().pushCritical("Failed to load:", f)
            continue

        QgsMapLayerRegistry.instance().addMapLayer(layer)

load_folder("path/to/your/vector/files/folder")

QgsMapLayerRegistry representa o painel de camadas presente na GUI do QGIS
iface.messageBar() retorna a barra de mensagens do aplicativo principal e permite notificar o usuário sobre o que está acontecendo
QgsVectorLayer representa uma camada vetorial com seus conjuntos de dados vetoriais subjacentes

2. EDIÇÃO DA TABELA DE ATRIBUTOS DA CAMADA ATIVA

O código a seguir demonstra a possibilidade de editar a tabela de atributos da camada vetorial por meio do console.

Qualquer atributo a ser escrito tem que vir na forma de um qgis.core.QgsField – isso é mais ou menos um encapsulamento de um nome de atributo e seu tipo (PyQt4.QtCore.QVariant)
O provedor de dados deve ser capaz de realizar a atribuição ( caps & QgsVectorDataProvider.AddAttributes)
QgsVectorLayer.addAttribute – método que retorna um booleano ao invés de lançar uma exceção

from qgis.core import QgsField
from qgis.gui import QgsMessageBar
from PyQt4.QtCore import QVariant


def edit_active_layer(attr_name, attr_type):
    layer = iface.activeLayer()
    caps = layer.dataProvider().capabilities()

    if caps & QgsVectorDataProvider.AddAttributes:
        layer.startEditing()
        if layer.addAttribute(QgsField(attr_name, attr_type)):
            iface.messageBar().pushMessage("Attribute {0} was successfully added to the active layer.".format(attr_name), QgsMessageBar.SUCCESS)
            layer.commitChanges()
        else:
            iface.messageBar().pushMessage("Attribute {0} was not added. Does it already exist?".format(attr_name), QgsMessageBar.CRITICAL)
            layer.rollBack()

edit_active_layer("new_string_attribute", QVariant.String)

3. CRIANDO UMA NOVA CAMADA VETORIAL

É possível criar uma camada vetorial totalmente nova com o console Python. Abaixo, é apresentada uma função muito simples (create_new_layer), mas espero que você possa imaginar as maneiras como ela pode ser ajustada.

from qgis.core import QgsField, QgsFields, QgsVectorLayer, QgsFeature, QgsGeometry, QgsPoint
from PyQt4.QtCore import QVariant

def create_new_layer():
    filename = "/path/to/your/vector/file.gpkg"

    fields = QgsFields()
    fields.append(QgsField("attr1", QVariant.String))
    fields.append(QgsField("attr2", QVariant.Int))

    file = QgsVectorFileWriter(
        filename,
        "UTF8",
        fields,
        QGis.WKBPoint,
        QgsCoordinateReferenceSystem(4326),
        "GPKG"
    )

    layer = QgsVectorLayer(filename, filename.split("/")[-1], "ogr")
    QgsMapLayerRegistry.instance().addMapLayer(layer)

    if not layer.dataProvider().capabilities() & QgsVectorDataProvider.AddAttributes:
        pass

    feature = QgsFeature(layer.pendingFields())
    feature.setGeometry(QgsGeometry().fromPoint(QgsPoint(0, 0)))
    feature.setAttribute("attr1", "attr1")
    feature.setAttribute("attr2", 2)

    layer.startEditing()

    if layer.addFeature(feature, True):
        layer.commitChanges()
    else:
        layer.rollBack()
        iface.messageBar().pushMessage("Feature addition failed.", QgsMessageBar.CRITICAL)

create_new_layer()

Esses foram apenas alguns exemplos do que pode ser feito com a API do QGIS e o console Python. No próximo post, vamos nos concentrar nas associações espaciais dentro do QGIS.

02 décembre 2018
17:08
Free and Open Source GIS Ramblings: Movement data in GIS #17: Spatial analysis of GeoPandas trajectories
sur Planet OSGeo
In Movement data in GIS #16, I presented a new way to deal with trajectory data using GeoPandas and how to load the trajectory GeoDataframes as a QGIS layer. Following up on this initial experiment, I’ve now implemented a first version of an algorithm that performs a spatial analysis on my GeoPandas trajectories.

The first spatial analysis algorithm I’ve implemented is Clip trajectories by extent. Implementing this algorithm revealed a couple of pitfalls:
- To achieve correct results, we need to compute spatial intersections between linear trajectory segments and the extent. Therefore, we need to convert our point GeoDataframe to a line GeoDataframe.
- Based on the spatial intersection, we need to take care of computing the corresponding timestamps of the events when trajectories enter or leave the extent.
- A trajectory can intersect the extent multiple times. Therefore, we cannot simply use the global minimum and maximum timestamp of intersecting segments.
- GeoPandas provides spatial intersection functionality but if the trajectory contains consecutive rows without location change, these will result in zero length lines and those cause an empty intersection result.
So far, the clip result only contains the trajectory id plus a suffix indicating the sequence of the intersection segments for a specific trajectory (because one trajectory can intersect the extent multiple times). The following screenshot shows one highlighted trajectory that intersects the extent three times and the resulting clipped trajectories:

This algorithm together with the basic trajectory from points algorithm is now available in a Processing algorithm provider plugin called Processing Trajectory.

Note: This plugin depends on GeoPandas.

Note for Windows users: GeoPandas is not a standard package that is available in OSGeo4W, so you’ll have to install it manually. (For the necessary steps, see this answer on gis.stackexchange.com)

The implemented tests show how to use the Trajectory class independently of QGIS. So far, I’m only testing the spatial properties though:
```
def test_two_intersections_with_same_polygon(self):
    polygon = Polygon([(5,-5),(7,-5),(7,12),(5,12),(5,-5)])
    data = [{'id':1, 'geometry':Point(0,0), 't':datetime(2018,1,1,12,0,0)},
        {'id':1, 'geometry':Point(6,0), 't':datetime(2018,1,1,12,10,0)},
        {'id':1, 'geometry':Point(10,0), 't':datetime(2018,1,1,12,15,0)},
        {'id':1, 'geometry':Point(10,10), 't':datetime(2018,1,1,12,30,0)},
        {'id':1, 'geometry':Point(0,10), 't':datetime(2018,1,1,13,0,0)}]
    df = pd.DataFrame(data).set_index('t')
    geo_df = GeoDataFrame(df, crs={'init': '31256'})
    traj = Trajectory(1, geo_df)
    intersections = traj.intersection(polygon)
    result = []
    for x in intersections:
        result.append(x.to_linestring())
    expected_result = [LineString([(5,0),(6,0),(7,0)]), LineString([(7,10),(5,10)])]
    self.assertEqual(result, expected_result) 
```
One issue with implementing the algorithms as QGIS Processing tools in this way is that the tools are independent of one another. That means that each tool has to repeat the expensive step of creating the trajectory objects in memory. I’m not sure this can be solved.
30 novembre 2018
22:16
Cameron Shorter: OSGeoLive in 4 minutes
sur Planet OSGeo

Here is our OSGeoLive lightning talk, stripped back to 4 minutes, and presented at the FOSS4G-Oceania conference. OSGeoLive talk starts at 8mins 15secs.
The shortened slide deck and notes are here.
18:05
GIScussions: Sh*thole Geography Standup Quiz
sur Planet OSGeo

Here is a little bit of fun for you to enjoy at a staff social, after your Xmas lunch or in the pub one evening

This is how the Sh*thole Geography Standup Quiz works.

“At the beginning everyone stands up.
Quizmaster asks a question with 3 suggested answers.
If you think the answer is the first one put your left hand up, if you think it is the second one both hands up, if you think it is the third one put your right hand up.
QM tells you the correct answer, everyone who got it wrong sit down. If you didn’t make a choice then you also sit down. There are some other interesting stats in the speaker notes for each slide (you need to download the presentation)
Then we go again until one person is left standing – if everyone gets knocked out you can all stand up again and keep answering questions.
There is a tie breaker question at the end”

And you could ignore all my instructions and just have a fun time after Xmas lunch or even do the quiz on your own and see how many questions you can answer, you could tweet your score with the hashtag #ShitholeQuiz

What’s the point of the maps at the end of each slide/question? I’m not sure! The quiz was first presented at FOSS4G with Ken and we had a cartographic spin to the whole talk, I changed it to become a bit more philosophical – maybe the maps disprove some of our preconceptions.

Don’t press the autoplay button on the quiz or it will go into auto advance mode which is not the ideal way to run the quiz, use the > button to advance and the < to go backwards if needed. If you want to download the quiz or access directly you can find the quiz here
15:05
Even Rouault: SRS barn raising: 6th report
sur Planet OSGeo
This is the sixth progress report of the GDAL SRS barn effort. The pace of changes has not yet slow down, with still a significant part of the work being in PROJ, and an initial integration in GDAL.
The major news item is that RFC2, implementing the new capabilities (WKT-2 support, late-binding approach, SQLite database), has now been merged into PROJ master.
An initial integration of PROJ master into GDAL has been started in a custom GDAL branch . This includes:
- PROJ master, which will be released as 6.0, is now a required dependency of GDAL. It actually becomes the only required external third-party dependency (if we except the copies of a few libraries, such as libtiff, libgeotiff, etc. that have been traditionaly included into the GDAL source tree)
- The dozen of continuous integration configurations have been modified to build PROJ master as a preliminary step.
- Related to the above, we have including into PROJ a way to "version" its symbols. If PROJ is built with -DPROJ_RENAME_SYMBOLS in CFLAGS and CXXFLAGS, all its exported symbols are prefixed with "internal_". This enables GDAL to link against PROJ master, while still using pre-compiled dependencies (such as libspatialite) that link against the system PROJ version, without a risk of symbol clash. This is particularly useful to be able to run GDAL autotests on continuous integration environments that use pre-packaged dependencies (or if you want to test the new GDAL without rebuilding all reverse dependencies of GDAL). This however remains a hack, and ultimately when PROJ 6 has been released, all reverse dependencies should be built against it. (this solution has been successfully tested in the past where GDAL had a libtiff 4.0 internal copy, whereas external libtiff used by some GDAL dependencies relied on the system libtiff 3.X)
- Compatibility mechanisms which were required to support older PROJ versions have been removed. In particular, the runtime loading (using dlopen() / LoadLibrary() mechanism) has been removed. It proved to cause code complication, and users frequently ran into headaches with different PROJ versions being loaded and clashing/crashing at runtime.
- The OGRSpatialReference class which implements CRS manipulation in GDAL has been modified to use the new PROJ functions to import and export between WKT and PROJ strings. Previously GDAL had such code, which is now redundant with what PROJ offers. This preliminary integration caused a number of fixes to be made on PROJ to have compatibility with the input and output of GDAL for WKT.1 and PROJ strings. Besides "moving" code from GDAL to PROJ, a practical consequence is that the addition of a new projection method into PROJ will no longer require changes to be made to GDAL for it to be usable for reprojection purposes.
There have been reflections on how to use the new code developped in PROJ by the existing PROJ code. A pull request is currently under review and implements:
- changes needed to remove from the data/ directory the now obsolete EPSG, IGNF, esri and esri.extra files to rely instead of the proj.db dataase
- making the proj_create_crs_to_crs() API use the new late-binding approach to create transformation pipelines
- updating cs2cs to use that new API.
- list and address backward compatibility issues related to honouring official axis order
Integration phase in GDAL continus with the aim of using more of the new PROJ code. Typically the OGRSpatialReference class that models in GDAL the CRS/SRS was up to now mostly a hierarchy of WKT nodes, where setters methods of OGRSpatialReference would directly create/modify/delete nodes, and getter methods query them. This approach was fine when you had to manage just one WKT version (with the caveat that it was also easy to produce invalid WKT representions, lacking mandatory nodes). However, this is no longer appropriate now that we want to support multiple WKT versions. Our goal is to make OGRSpatialReference act rather on osgeo::proj::CRS objects (and its derived classes). Switching between the two abstractions is a non-trivial task and doing it in a bing-bang approach seemed risky, so we are progressively doing it by using a dual internal modelling. A OGRSpatialReference instance will maintain as a primary source a osgeo::proj::CRS object, and for operations not yet converted to the new approach, will fallback to translating it internally to WKT.1 to allow direct manipulation of the nodes, and then retranslate that updated WKT.1 representation back to a osgeo::proj::CRS object. Ultimately the proportion of methods using the fallback way should decrease (it is not completely clear we can remove all of them since direct node manipulation is spread in a significant number of GDAL drivers). The task is slowly progressing, because each change can subtely modify the final WKT.1 representation (nodes being added, number of significant digits changing) and cause a number of unit tests to break (GDAL autotest suite is made of 280 000 lines of Python code) and be analyzed to see if there was a bug and or just an expected result to be slightly altered.Because of all the above impacts, we have decided to do an early release in December of GDAL master as GDAL 2.4.0 with all the new features since GDAL 2.3, in order to be able to land this PROJ integration afterwards afterwards. A GDAL 2.5.0 release will hopefully follow around May 2019 with the result of the gdalbarn work.

Other side activities regarding collecting transformation grids:
- Following a clarification from IGN France on the open licensing of their geodesy related resources, their CRS and transformation XML registry is now processed to populate the IGNF objects in the proj.db database (the previous import used the already processed IGNF file containing PROJ string, which caused information losses). The associated vertical shift grids have also been converted from their text-based format to the PROJ digestable .gtx format, integrated in the proj-datumgrid-europe package, and they have been referenced in the database for transformations that use them.
- The NGS GEOID 2012B vertical grids to convert between NAD83 ellipsoidal heights and NAVD88 heights have also been integrated in the proj-datumgrid-north-america package
29 novembre 2018
18:34
Free and Open Source GIS Ramblings: TimeManager 3.0.2 released!
sur Planet OSGeo

Bugfix release 3.0.2 fixes an issue where “accumulate features” was broken for timestamps with milliseconds.

If you like TimeManager, know your way around setting up Travis for testing QGIS plugins, and want to help improve TimeManager stability, please get in touch!
12:33
Fernando Quadro: Desenvolvimento de plugins do QGIS: Primeiros passos
sur Planet OSGeo
O QGIS é uma ferramenta brilhante para automação baseada em Python em forma de scripts personalizados ou até mesmo plugins. Os primeiros passos para escrever o código personalizado podem ser um pouco difíceis, já que você precisa entender a API do Python que é um pouco complexa. A série de Desenvolvimento de Plugins do QGIS que inicia hoje tem como objetivo o desenvolvimento de um plug-in personalizado totalmente funcional capaz de gravar valores de atributos de uma camada de origem para uma camada de destino com base em sua proximidade espacial.

Nesta parte, mencionarei o básico, o que é bom saber antes de começar.

1. Documentação

Diferentes versões do QGIS vêm com diferentes APIs do Python. A documentação deve ser encontrada em [https:]] , sendo a mais recente a versão 3.2. Observe que, se você acessar diretamente o site [qgis.org] , verá os documentos atuais.

Alternativamente, você pode o “apt install qgis-api-doc” em seu sistema Ubuntu e rodar “python -m Simpl [HTTPServer[port]”] dentro de “/usr/share/qgis/doc/api”. Você encontrará a documentação em � [localhost:8000”] (se você não fornecer o número da porta) e estará disponível mesmo quando estiver off-line.

2. Estrutura Básica da API

Abaixo uma visão resumida do que está disponível dentro da API:
- Pacote qgis.core traz todos os objetos básicos como QgsMapLayer, QgsDataSourceURI, QgsFeature, etc
- O pacote qgis.gui traz elementos GUI que podem ser usados dentro do QGIS como QgsMessageBar ou QgsInterface
- qgis.analysis, qgis.networkanalysis, qgis.server e qgis.testing pacotes que não serão abordados na série
- Módulo qgis.utils que vem com iface (muito útil no console do QGIS Python)
3. Console Python do QGIS

Usar o console Python é a maneira mais fácil de automatizar o fluxo de trabalho do QGIS. Pode ser acessado pressionando Ctrl + Alt + P ou navegando pelo menu em Plugins -> Python Console. Como mencionado acima o módulo iface do qgis.utils é exposto por padrão dentro do console, permitindo-lhe interagir com o QGIS GUI. Experimente os exemplos a seguir, para fazer um teste inicial:
```
iface.mapCanvas().scale() # retorna a escala atual do mapa
iface.mapCanvas().zoomScale(100) # zoom para escala 1:100 
iface.activeLayer().name() # obtém o nome da camada ativa 
iface.activeLayer().startEditing() # realizar edição
```
Essa foi uma breve introdução à API do QGIS, no próximo post falaremos mais profundamente sobre o console do QGIS.
7:27
GRASS GIS: GRASS GIS 7.4.3 released
sur Planet OSGeo

We are pleased to announce the GRASS GIS 7.4.3 release
28 novembre 2018
15:05
Fernando Quadro: Alteração do JTS no GeoServer
sur Planet OSGeo
A partir do GeoServer 2.14, a saída produzida pelo recurso REST para as solicitações de camadas estão usando um nome de pacote diferente (org.locationtech em vez de com.vividsolutions) para os tipos de geometria devido ao upgrade para o JTS (Java Topology Suite) 1.16.0. Por exemplo:

Antes:
```
...
<attribute>
  <name>geom</name>
  <minOccurs>0</minOccurs>
  <maxOccurs>1</maxOccurs>
  <nillable>true</nillable>
  <binding><strong>com.vividsolutions.jts.geom.Point</strong></binding>
</attribute>
...
```
Depois:
```
...
<attribute>
  <name>geom</name>
  <minOccurs>0</minOccurs>
  <maxOccurs>1</maxOccurs>
  <nillable>true</nillable>
  <binding>org.locationtech.jts.geom.Point</binding>
</attribute>
...
```
Qualquer cliente REST que dependa dessa informação de ligação deve atualizadar para suportar os novos nomes.

Fonte: GeoServer Documentation
26 novembre 2018
14:00
Paul Ramsey: Esri and Winning
sur Planet OSGeo
How much winning is enough? Have you been winning so much that you’re tired of winning now?

I ask because last month I gave a talk (PDF Download) to the regional GIS association in Manitoba, about open source and open data. My talk included a few points that warned about the downsides of being beholden to a single software vendor, and it included some information about the software available in the open source geospatial ecosystem.

In a testament to the full-spectrum dominance of Esri, the audience was almost entirely made up of Esri customers. The event was sponsored by Esri. Esri had a table at the back of the room. Before giving my talk, I made a little joke that my talk would in fact include some digs at Esri (though I left the most pointed ones on the cutting room floor).

People seemed to like the talk. They laughed at my jokes. They nodded in the right places.

Among the points I made:
- Single-vendor dominance in our field is narrowing the understanding of “what is possible” amongst practitioners in that ecosystem.
- Maintaining a single-vendor policy dramatically reduces negotiating power with that vendor.
- Maintaining a single-vendor policy progressively de-skills your staff, as they become dependant on a single set of tooling.
- Practitioners have higher market value when they learn more than just the tools of one vendor, so self-interest dictates learning tools outside the single-vendor ecosystem.
- Point’n’click GIS tools from Esri have widened access to GIS, which is a good thing, but driven down the market value of practitioners who limit themselves to those tools.
None of these points is unique to Esri – they are true of any situation where a single tool has driven competitors off the field, whether it be Adobe graphics tools or Autodesk CAD tools or Microsoft office automation tools.

Nor are any of these points indicative of any sort of ill will or malign intent on the part of Esri – they are just the systemic effects of market dominance. It is not contingent on Esri to change their behaviour or limit their success; it’s contingent on practitioners and managers to recognize the negative aspects of the situation and react accordingly.

And yet.

Despite the fact that almost all the people in the room were already their customers, that no new business would be endangered by my message, that all the students would still be taught their tools, that all the employers would still include them in job requirements, that people would continue to use the very words they choose to describe basic functions of our profession …

Despite all that, the Esri representative still went to the president of the association, complained to her about the content of my talk, and asked her to ensure that nothing I would say in my afternoon technical talk would be objectionable to himself. (In the event, I made some nasty jokes about Oracle, nobody complained.)

For some of these people, no amount of winning is enough, no position of dominance is safe, no amount of market leverage is sufficient.

It’s sad and it’s dangerous.

I was reminded of this last week, meeting an old friend in Australia and learning that he’d been blackballed out of a job for recommending software that wasn’t Esri software. Esri took away his livelihood for insufficient fealty.

This is the danger of dominance.

When the local Esri rep has a better relationship with your boss than you do, do you advocate for using alternative tools? You could be limiting or even potentially jeapardizing your career.

When Esri has locked up the local geospatial software market, do you bid an RFP with an alternative open tool set? You could lose your Esri partnership agreement and with it your ability to bid any other local contracts. Esri will make your situation clear to you.

This is the danger of dominance.

A market with only one vendor is not a market. There’s a name for it, and there’s laws against it. And yet, our profession glories in it. We celebrate “GIS day”, a marketing creation of our dominant vendor. Our publicly funded colleges and universities teach whole curricula using only Esri tools.

And we, as a profession, do not protest. We smile and nod. We accept our “free” or “discounted” trainings from Esri (comes with our site license!) and our “free” or “discounted” tickets to the Esri user conference. If we are particularly oblivious, we wonder why those open source folks never come around to market their tools to us.

We have met the enemy, and he is us.
25 novembre 2018
1:05
Nyall Dawson: Thoughts on “FOSS4G/SOTM Oceania 2018”, and the PyQGIS API improvements which it caused
sur Planet OSGeo
Last week the first official “FOSS4G/SOTM Oceania” conference was held at Melbourne University. This was a fantastic event, and there’s simply no way I can extend sufficient thanks to all the organisers and volunteers who put this event together. They did a brilliant job, and their efforts are even more impressive considering it was the inaugural event!

Upfront — this is not a recap of the conference (I’m sure someone else is working on a much more detailed write up of the event!), just some musings I’ve had following my experiences assisting Nathan Woodrow deliver an introductory Python for QGIS workshop he put together for the conference. In short, we both found that delivering this workshop to a group of PyQGIS newcomers was a great way for us to identify “pain points” in the PyQGIS API and areas where we need to improve. The good news is that as a direct result of the experiences during this workshop the API has been improved and streamlined! Let’s explore how:

Part of Nathan’s workshop (notes are available here) focused on a hands-on example of creating a custom QGIS “Processing” script. I’ve found that preparing workshops is guaranteed to expose a bunch of rare and tricky software bugs, and this was no exception! Unfortunately the workshop was scheduled just before the QGIS 3.4.2 patch release which fixed these bugs, but at least they’re fixed now and we can move on…

The bulk of Nathan’s example algorithm is contained within the following block (where “distance” is the length of line segments we want to chop our features up into):
```
for input_feature in enumerate(features):
    geom = feature.geometry().constGet()
    if isinstance(geom, QgsLineString):
        continue
    first_part = geom.geometryN(0)
    start = 0
    end = distance
    length = first_part.length()

    while start < length:
        new_geom = first_part.curveSubstring(start,end)

        output_feature = input_feature
        output_feature.setGeometry(QgsGeometry(new_geom))
        sink.addFeature(output_feature)

        start += distance
        end += distance
```
There’s a lot here, but really the guts of this algorithm breaks down to one line:
```
new_geom = first_part.curveSubstring(start,end)
```
Basically, a new geometry is created for each trimmed section in the output layer by calling the “curveSubstring” method on the input geometry and passing it a start and end distance along the input line. This returns the portion of that input LineString (or CircularString, or CompoundCurve) between those distances. The PyQGIS API nicely hides the details here – you can safely call this one method and be confident that regardless of the input geometry type the result will be correct.

Unfortunately, while calling the “curveSubstring” method is elegant, all the code surrounding this call is not so elegant. As a (mostly) full-time QGIS developer myself, I tend to look over oddities in the API. It’s easy to justify ugly API as just “how it’s always been”, and over time it’s natural to develop a type of blind spot to these issues.

Let’s start with the first ugly part of this code:
```
geom = input_feature.geometry().constGet()
if isinstance(geom, QgsLineString):
    continue
first_part = geom.geometryN(0)
# chop first_part into sections of desired length
...
```
This is rather… confusing… logic to follow. Here the script is fetching the geometry of the input feature, checking if it’s a LineString, and if it IS, then it skips that feature and continues to the next. Wait… what? It’s skipping features with LineString geometries?

Well, yes. The algorithm was written specifically for one workshop, which was using a MultiLineString layer as the demo layer. The script takes a huge shortcut here and says “if the input feature isn’t a MultiLineString, ignore it — we only know how to deal with multi-part geometries”. Immediately following this logic there’s a call to geometryN( 0 ), which returns just the first part of the MultiLineString geometry.

There’s two issues here — one is that the script just plain won’t work for LineString inputs, and the second is that it ignores everything BUT the first part in the geometry. While it would be possible to fix the script and add a check for the input geometry type, put in logic to loop over all the parts of a multi-part input, etc, that’s instantly going to add a LOT of complexity or duplicate code here.

Fortunately, this was the perfect excuse to improve the PyQGIS API itself so that this kind of operation is simpler in future! Nathan and I had a debrief/brainstorm after the workshop, and as a result a new “parts iterator” has been implemented and merged to QGIS master. It’ll be available from version 3.6 on. Using the new iterator, we can simplify the script:
```
geom = input_feature.geometry()
for part in geom.parts():
    # chop part into sections of desired length
    ...
```
Win! This is simultaneously more readable, more Pythonic, and automatically works for both LineString and MultiLineString inputs (and in the case of MultiLineStrings, we now correctly handle all parts).

Here’s another pain-point. Looking at the block:
```
new_geom = part.curveSubstring(start,end)
output_feature = input_feature
output_feature.setGeometry(QgsGeometry(new_geom))
```
At first glance this looks reasonable – we use curveSubstring to get the portion of the curve, then make a copy of the input_feature as output_feature (this ensures that the features output by the algorithm maintain all the attributes from the input features), and finally set the geometry of the output_feature to be the newly calculated curve portion. The ugliness here comes in this line:
```
output_feature.setGeometry(QgsGeometry(new_geom))
```
What’s that extra QgsGeometry(…) call doing here? Without getting too sidetracked into the QGIS geometry API internals, QgsFeature.setGeometry requires a QgsGeometry argument, not the QgsAbstractGeometry subclass which is returned by curveSubstring.

This is a prime example of a “paper-cut” style issue in the PyQGIS API. Experienced developers know and understand the reasons behind this, but for newcomers to PyQGIS, it’s an obscure complexity. Fortunately the solution here was simple — and after the workshop Nathan and I added a new overload to QgsFeature.setGeometry which accepts a QgsAbstractGeometry argument. So in QGIS 3.6 this line can be simplified to:
```
output_feature.setGeometry(new_geom)
```
Or, if you wanted to make things more concise, you could put the curveSubstring call directly in here:
```
output_feature = input_feature
output_feature.setGeometry(part.curveSubstring(start,end))
```
Let’s have a look at the simplified script for QGIS 3.6:
```
for input_feature in enumerate(features):
    geom = feature.geometry()
    for part in geom.parts():
        start = 0
        end = distance
        length = part.length()

        while start < length:
            output_feature = input_feature
            output_feature.setGeometry(part.curveSubstring(start,end))
            sink.addFeature(output_feature)

            start += distance
            end += distance
```
This is MUCH nicer, and will be much easier to explain in the next workshop! The good news is that Nathan has more niceness on the way which will further improve the process of writing QGIS Processing script algorithms. You can see some early prototypes of this work here:

I couldn’t be at the community day but managed to knock out some of the new API on the plane on the way home. **API subject to change #FOSS4G_SotM_Oceania

cc @nyalldawson @underdarkGIS this will create a alg under the hood and check for double inputs, etc pic.twitter.com/49VpyuukGU

— Nathan Woodrow (@madmanwoo) November 23, 2018

So there we go. The process of writing and delivering a workshop helps to look past “API blind spots” and identify the ugly points and traps for those new to the API. As a direct result of this FOSS4G/SOTM Oceania 2018 Workshop, the QGIS 3.6 PyQGIS API will be easier to use, more readable, and less buggy! That’s a win all round!
24 novembre 2018
1:00
PostGIS Development: PostGIS 2.3.8, 2.4.6
sur Planet OSGeo

The PostGIS development team is pleased to provide bug fix 2.3.8 and 2.4.6 for the 2.3 and 2.4 stable branches.
Continue Reading by clicking title hyperlink ..
23 novembre 2018
11:30
Fernando Quadro: Editor de estilos do GeoServer agora com modo “Tela cheia”
sur Planet OSGeo

A versão 2.14.1 do GeoServer veio com uma novidade: A página do editor de estilo agora tem um botão de “Tela cheia” no canto superior direito da janela:

Se pressionado, o editor e a visualização serão exibidos lado a lado e usarão todo o espaço da janela do navegador:

A cada atualização a seção de estilos do GeoServer tem melhorado. A equipe de desenvolvimento já a algum tempo tem dado uma atenção especial a essa seção, e já está muito mais fácil lidar com a edição e visualização dos estilo no GeoServer.
5:39
XYCarto: Wellington Elevations: Interpolating the Bathymetry
sur Planet OSGeo
It is important to note something from the very beginning. The interpolated bathymetry developed in this project does not reflect the actual bathymetry of the Wellington Harbour. It is my best guess based on the tools I had and the data I worked with. Furthermore, this interpolation is NOT the official product of any institution. It is an interpolation created by me only for the purposes of visualization.

Part of the goal when visualizing the Wellington landscape was to incorporate a better idea about what may be happening below the surface of the harbor. Various bathymetric scans in the past have gathered much of the information and institutions like NIWA have done the work visualizing that data. As for myself, I did not have access to those bathymetries; however, I did have a sounding point data set to work with, so I set about interpolating those points.

The data set, in CSV format, was over a million points; too dense for a single interpolation. I worked out a basic plan for the interpolation based on splitting the points into a grid, interpolate the smaller bits, then reassemble the grid tiles into a uniform bathymetry.

Conversion from CSV to shp
Using the open option (-oo) switch, OGR will convert CSV to shp seamlessly
```
ogr2ogr -s_srs EPSG:4167 -t_srs EPSG:4167 -oo X_POSSIBLE_NAMES=$xname* -oo Y_POSSIBLE_NAMES=$yname*  -f "ESRI Shapefile" $outputshapepath/$basenme.shp $i
```
Gridding the Shapefile
With the shapefile in place, I next needed to break it into smaller pieces for interpolation. For now, I create the grid by hand in QGIS using the ‘Create Grid’ function. This is found under Vector>Reasearch Tools>Create Grid. Determining a grid size that works best for the interpolation is a bit of trial and error. You want the largest size your interpolation can manage without crashing. Using the grid tool from QGIS in very convenient, in that it creates an attribute table of the xmin, xmax, ymin, ymax corrodinates for each tile in the grid. These attributes become very helpful during the interpolation process.

Interpolating the Points
I switched things up in the interpolation methods this time and tried out SAGA GIS. I have been looking for a while now for a fast and efficient method of interpolation that I could easily build into a scripted process. SAGA seemed like a good tool for this. The only drawback, I had a very hard time finding examples online about how to use this tool. My work around to was to test the tool in QGIS first. I noticed when the command would run, QGIS saved the last command in a log file. I found that log, copied out the command line function, and began to build my SAGA command for my script from there.

Here is look at the command I used:
```
saga_cmd grid_spline "Multilevel B-Spline Interpolation" -TARGET_DEFINITION 0 -SHAPES "$inputpoints" -FIELD "depth" -METHOD 0 -EPSILON 0.0001 -TARGET_USER_XMIN $xmin -TARGET_USER_XMAX $xmax -TARGET_USER_YMIN $ymin -TARGET_USER_YMAX $ymax -TARGET_USER_SIZE $reso -TARGET_USER_FITS 0 -TARGET_OUT_GRID "$rasteroutput/sdat/spline_${i}"
```
I tested a number of methods and landed on ‘grid_spline’ as producing the best results for the project. It was useful because it did a smooth interpolation across the large ‘nodata’ spaces.

Once the initial interpolation was complete, I needed to convert the output to GeoTIFF since SAGA exports in an .sdat format. Easy enough since GDAL_TRANSLATE recognizes the .sdat format. I then did my standard prepping and formatting for visualization:
```
gdal_translate "$iupput_sdat/IDW_${i}.sdat" "$output_tif/IDW_${i}.tif"
gdaldem hillshade -multidirectional -compute_edges "$output_tif/IDW_${i}.tif" "$ouput_hs/IDW_${i}.tif"
gdaladdo -ro "$output_tif/IDW_${i}.tif" 2 4 8 16 32 64 128
gdaladdo -ro "$ouput_hs/IDW_${i}.tif"2 4 8 16 32 64 128
```
Here is look at the interpolated harbour bathymetry, hillshaded, with Wellington 1m DEM hillshade added over top

And here is a look at the same bathy hillshade with coloring

Visualizing the Bathymetry
With the bathymetry, complete it was simply a matter of building it into the existing visualization I built for the Wellington Region. Learn more about the project here. The visualization was four steps:

Hillshade

Color

Aerial Imagery

Then merge the models together

Easy as, eh? Let me know what you think!

Note: All imagery was produced during my time at Land Information New Zealand. Imagery licensing can be found here:
“Source: Land Information New Zealand (LINZ) and licensed by LINZ for re-use under the Creative Commons Attribution 4.0 International licence."
22 novembre 2018
23:21
XYCarto: Building the Wellington Model with 1m DEM and DSM
sur Planet OSGeo

As interest in LiDAR derived elevation increases, so grows the interest in the capabilities. LiDAR derived elevation data has been great for my visualization game and in helping me communicate the story out about what LiDAR can do. It all starts with a picture to get the imagination going.

The Wellington model derived for this project is part of an ongoing project to help increase the exposure of the Wellington 1m DEM/DSM elevation data derived from LiDAR. Step one for me is getting a working model built in QGIS, capturing still images, and increasing interest in the data.

I’ve talked about the processing of the elevation data for Wellington visualizations in the past, so for this post I’m only focusing on the blending of the data sets in building the model. This project is a good model since it encompasses a number of subtle techniques to get the model to stand out. This post is one of a two part series; the second post discusses the techniques used to derive and visualize the bathymetry for the surrounding harbor.

Let’s start with the base, Aerial Imagery.

Blended with a hillshade

DSM added for texture and context

Slope added to define some edges

Some darker shading added to the bathymetry to frame the elevation data

And finally some added bathymetry to lighten the edges at the shoreline enhancing the frame a bit more.

In the end there is some post-processing in Photoshop to lighten up the image. Honestly, this could have been done in QGIS, but I was being lazy. For the images produced, there was no need to retain the georeferencing, and when that is the case, I rely on Photoshop for color and light balancing.

The greatest difficultly in this project so far has been trying to create a universal model for the data set. I’m finding that as I visualize different regions using this model, I need to adjust the hillshading quite significantly to draw out different features. Take a look at the images here. It is the same model, but with the noticeably different gradients used in the hillshades. The techniques used for the images in this post worked well for the urban region shown, but fall apart as you move further out into the more mountainous regions. Much of the blending is too harsh and turns the mountains into a black muddled mess. I am almost there, but like any project, it takes a good bit of subtle tweaking of the blending to get a universal image to work.

The entire base mapping work is completed in QGIS. The elevation data was processed using GDAL and the bathymetric interpolations were produced SAGA GIS. There are no color palettes for this project. The aerial imagery does all the work in that department.

Base data can be found here:
DEM: [https:]]
DSM: [https:]]
Aerial Imagery: [https:]]

The next post covers the development of the bathymetry for the surrounding harbor. Thanks for having a look and let me know what you think.

Note: All imagery was produced during my time at Land Information New Zealand. Imagery licensing can be found here:
“Source: Land Information New Zealand (LINZ) and licensed by LINZ for re-use under the Creative Commons Attribution 4.0 International licence.”
15:07
Fernando Quadro: Camadas WMS agora com suporte a Dimensões
sur Planet OSGeo
O GeoServer adicionou suporte a adição de dimensões específicas para camadas WMS, conforme especificado nos padrões WMS 1.1.1 e WMS 1.3.0. Existem duas dimensões pré-definidas nos padrões WMS: TIME e ELEVATION. A ativação das dimensões para uma camada permite que os usuários as especifiquem como parâmetros extras em solicitações GetMap, úteis para criar mapas ou animações a partir de dados multidimensionais.

Essas configurações podem ser ativadas e configuradas na Aba “Dimensões”:

Para cada dimensão ativa, as seguintes opções de configuração estão disponíveis:
- Atributo – Atribua o nome para escolher o valor para esta dimensão (somente vetor). Isso é tratado no início do intervalo se o atributo final também for fornecido.
- Atributo final – Atribua o nome para escolher o final do intervalo de valores para essa dimensão (opcional, somente para vetores).
- Apresentação – Tipo de apresentação para os valores disponíveis no documento de recursos.
- Valor padrão – Valor padrão a ser usado para essa dimensão, se nenhuma for fornecida com a solicitação. Selecione uma das quatro estratégias:
- Valor de referência – O especificador de valor padrão. Exibido apenas para as estratégias de valor padrão onde é usado.
- Correspondente mais próximo – Permitir, ou não, suporte ao correspondente mais próximo do WMS nesta dimensão. Atualmente suportado apenas na dimensão de tempo.
- Intervalo aceitável – uma distância máxima de pesquisa do valor especificado (disponível somente quando o correspondente mais próximo estiver ativado). Pode ser vazio (sem limite), um único valor (pesquisa simétrica) ou usado para especificar um intervalo de pesquisa assimétrica. As distâncias de tempo devem ser especificadas usando a sintaxe do ISO. Por exemplo, PT1H/PT0H permite pesquisar até uma hora antes do valor especificado pelo usuário, mas não depois.
Para dimensão de tempo, o valor deve estar no formato ISO 8601 DateTime (yyyy-MM-ddThh:mm:ss.SSSZ). Para dimensão de elevação, o valor deve ser e inteiro (ponto flutuante).

Somente para a estratégia “Valor de referência”, é possível usar intervalos de tempos e intervalos de elevação, no formulário fromValue/toValue. Somente para a estratégia “valor de referência”
também é possível usar tempos relativos como P1M/PRESENT, mas é preciso ter cuidado para que o valor de referência seja copiado integralmente no documento de recursos e, como resultado, nem todos os clientes estejam reconhecendo essa sintaxe.

Fonte: Blog do GeoServer
4:46
XYCarto: The Rejects
sur Planet OSGeo

Sometimes there is simply not enough room for all the ideas. Sometimes you want all the images to make it to the final round.

In a recent project to promote some of our elevation data, I was asked to present a number of ideas for a 2000mm x 900mm wall hanging. The piece was to act as a conversation starter and demonstrate some of the finer details elevation from LiDAR possesses.

In the end, the image above was the chosen candidate. Below are the drafts I initially presented for review. You can see the difference in treatment from the original ideas to the final product. Personally, I really enjoyed the images developed for the draft series, I liked the silvery undertones, and I thought it was a shame to merely let these images sit on my hard drive.
Below, you’ll find a brief description about a few challenges faced in the image development.

Artifacts and Finer Details
The hardest part of this job was drawing out the finer details of the chosen location. There was a strong interest in showing the ancient river bed; however, without a good bit of tweaking in the hillshades, the image is quite flat. After some trial and error, I found I could get a good contrast by limiting the hillshade values range to 170-190. That’s it, but the readability of the project really hinged on the simple tweak. It really made the details stand out.
That said, the gain in detail also revealed a significant artifact in the data. If you go back up and have a closer look, you will find diagonal depressions running across the images in equal intervals. These are lines from where the LiDAR scans overlap. I haven’t quite had the time to figure out how to remove these from the original data source, so for now I leave them in as conversational piece around improving LiDAR capture practices.
As usual, all map layout work was completed on QGIS, with the bulk of the data processing done using GDAL. The ‘Reject’ images for this post are direct exports from QGIS, with no manipulation apart from some down-sampling and cropping in Photoshop.

Base data can be found here:
DEM: [https:]]
DSM: [https:]]
Aerial Imagery: [https:]]

I produced a public repository for some of the scripting work. This repository is not specific to the above project but does contain some of the base processing I did on the Wellington elevation data: [https:]]
Hope you like and thanks for checking in!

Note: All imagery was produced during my time at Land Information New Zealand. Imagery licensing can be found here:
“Source: Land Information New Zealand (LINZ) and licensed by LINZ for re-use under the Creative Commons Attribution 4.0 International licence.”
1:00
PostGIS Development: PostGIS 2.2.8 EOL
sur Planet OSGeo
The PostGIS development team is pleased to provide bug fix 2.2.8 for the 2.2 stable branch.

This is the End-Of-Life and final release for PostGIS 2.2 series.

We encourage you to upgrade to a newer minor PostGIS version. Refer to our Version compatibility and EOL Policy for details on versions you can upgrade to.

This release supports PostgreSQL 9.1-9.6.

2.2.8
- source download
- NEWS
- PDF docs en
- html doc download
- html online en
- epub doc download
- Changelog
Continue Reading by clicking title hyperlink ..
21 novembre 2018
0:32
GeoServer Team: GeoSever 2.14.1 released
sur Planet OSGeo
We are happy to announce the release of GeoServer 2.14.1. Downloads are provided (zip|war|exe) along with docs (html|pdf) and extensions.

This is a stable release of the GeoServer 2.14 series and is recommended for all production systems. Users of prior releases of GeoServer are encouraged to upgrade.

This release is made in conjunction with GeoTools 20.1 and GeoWebCache 1.14.1. Thanks to all who contributed to this release.

For more information please see our release notes (2.14.1|2.14.0|2.14-RC).
Improvements and Fixes
This release includes a number of new features and improvements:
- New coordinate formatting options for WFS layers
- REST API granule management remove of metadata or all
- Coverage view support for indexed color model
- Check handling of WCS 2.0 time/elevation/range subsetting
- WFS 1.1.0 handling of gml:id
- WFS Shapefile and GeoJSON output now support geometry with measures
- Fixes to the start.bat and shutdown.bat scripts on windows
About GeoServer 2.14 Series
Additional information on the GeoServer 2.14 series:
- New MongoDB extension added
- Style editor improvements including side by side editing
- Nearest match support for WMS dimension handling
- Upgrade notes documenting REST API feature type definition change
- State of GeoServer 2.14 (SlideShare)
- GeoServer Ecosystem (SlideShare)
- GeoServer Developers Workshop (SlideShare)
0:30
GeoTools Team: GeoTools 20.1 Released
sur Planet OSGeo

The GeoTools team is happy to announce the release of GeoTools 20.1: geotools-20.1-bin.zip geotools-201-doc.zip geotools-20.1-project.zip geotools-20.1-userguide.zip maven repository This release is a stable release and is recommend for new development and production systems. This release is made in conjunction with GeoServer 2.14.1. Improvements GML respects formatting options for
20 novembre 2018
21:18
GIScussions: Shithole Geography re-interpreted at KortDage
sur Planet OSGeo
I was invited to deliver a keynote last week at KortDage to an audience of nearly 800 Danish (and Scandinavian) GI professionals. There were just a few snags:
- The person who invited me said “Last year Jack Dangermond gave a keynote, this year we wanted something different”. No challenge there!
- The slot was at 9.00am on the final morning after the delegates had been up till 1.30am boogying
- I don’t speak Danish and I wasn’t sure how my dry (I’d like to think) Brit humour would go down
So they wanted something different, they liked my #FAKEMAPS talk, I didn’t have a lot of time to come up with something completely new so I decided to re-invent Sh*thole Geography.

Ken Field and I had done a double act at FOSS4G in Dar es Salaam around Donald Trump’s opinion that some countries were “shithole countries” whose citizens were not welcome in the US. We had riffed on the “Orange Clown” as Ken likes to call 45, we had played with UN Sustainable Development Goal stats, Ken gave a brief cartography master class on making global maps and I had created a pub style geography quiz. I could work with that and switch the focus to make a slightly different message for this audience (with apologies to Ken) – a map may not always be the best way to present global statistics.

They didn’t record the whole talk so I can’t share that with you but you can get the drift of the points I was making from this short interview that the organisers did with me after my talk.

Kortdage 2018 – Kort fortæller ikke altid sandheden from Geoforum on Vimeo.

At the end of the talk I offered Sh*thole Geography badges and stickers to the attendees in return for donations to the FOSS4G Travel Grant Programme – the audience were amazingly generous and we raised almost $500

Apart from my talk I also got to meet Anna Webber, an artist who had made an “earth map” of Denmark

You can read more about her map and why she made it on Ken’s and my Mappery site.

Thanks to the organisers for inviting me to speak, thanks to all the people who laughed at my jokes and made me feel welcome. KortDage was a great event.
17:28
gvSIG Team: I Congreso de la Red JUST-SIDE: Infraestructuras de Datos Espaciales y Justicia Territorial
sur Planet OSGeo

Os hacemos llegar la invitación al primer evento organizado por la Red JUST-SIDE, de la que forma parte la Asociación gvSIG. Esperamos sea de vuestro interés (si podéis asistir…allí nos veremos!):

El Departamento de Geografía de la Facultad de Ciencias (UdelaR) tiene el agrado de invitarle al primer congreso de la red JUST-SIDE (Justicia y Sostenibilidad en el Territorio a través de las Infraestructuras de Datos Espaciales) a realizarse en el día 22 de noviembre de 2018 entre las 14 y 18 horas; el mismo abordará la temática de las Tecnologías Libres y contará con expositores nacionales y extranjeros.

Nuestra institución integra esta red de universidades y empresas de Argentina, Brasil, Chile, Costa Rica, España, México, Portugal y Uruguay que, coordinados por la Universidad de Coimbra (Portugal), está llevando adelante un proyecto que vincula el derecho ambiental y territorial con el acceso a la información geográfica y el desarrollo de las Infraestructuras de Datos Espaciales.

La red JUST-SIDE integra investigadores con experiencia en las áreas del derecho, ciencias sociales y tecnologías de información geográfica y busca crear una metodología basada en una Infraestructura de Datos Espaciales como apoyo a la toma de decisiones, promoviendo y fortaleciendo las políticas públicas para hacer frente a los desafíos sociales, ambientales, económicos, jurídicos y democráticos. Está iniciativa está financiada por CYTED (Programa Iberoamericano de Ciencia y Tecnología para el Desarrollo), y estudiará 16 casos concretos en territorios para la mejoría de la eficacia y justicia de las políticas públicas con incidencia territorial.

Quienes estén interesados en asistir les solicitamos realizar su registro previo a través del siguiente enlace, dado que el cupo es limitado; la actividad se desarrollará en el Hotel Dazzler (Salón Laureles), calle 21 de setiembre 2752 esquina Luis de la Torre.

Por favor realizar su registro

Se adjunta programa de la jornada.

– DEPARTAMENTO DE GEOGRAFÍA
13:16
gvSIG Team: Disponibles las ponencias y talleres de las 14as Jornadas Internacionales gvSIG
sur Planet OSGeo

Ya se han publicado las presentaciones realizadas en las 14as Jornadas Internacionales gvSIG, que se celebraron del 24 al 26 de octubre en Valencia (España).

También está disponible la grabación de las mismas y de algunos de los talleres que se impartieron. Los vídeos se encuentran disponibles en su idioma original, habiendo presentaciones y talleres en español e inglés.

Los talleres de los cuales está disponible la grabación para poder seguirlos fueron los de “gvSIG aplicado a Geología” y “Generación de informes con gvSIG Desktop“, ambos en español, y el de “Environmental modelling using gvSIG and the HortonMachine” en inglés.

Cabe destacar que el taller de Informes es la primera documentación disponible del plugin publicado hace unos días sobre la nueva versión gvSIG 2.5 que permite crear informes de forma automática en gvSIG Desktop. ¡No os lo perdáis!

Si no pudiste asistir a las jornadas ahora puedes disponer de la grabación de las distintas sesiones, y aparte de un buen material didáctico gratuito con los talleres.
11:30
Fernando Quadro: Acessando o pycsw com o QGIS
sur Planet OSGeo

Hoje vamos acessar e consultar o catálogo pycsw do GeoNode a partir de um cliente desktop, neste caso será o QGIS. Os printscreens foram gerados usando o QGIS 2.14 LTR, mas as versões seguintes devem ser semelhantes.

Em primeiro lugar, clique no menu Web e, em seguida, em MetaSearch > MetaSearch

Agora clique no botão Novo e preencha o formulário da seguinte forma:

Nome: GeoNode CSW
URL: [localhost:8000]

Clique no botão Informações do serviço e a identificação do serviço CSW deve aparecer na janela principal do MetaSearch:

Agora, clique na guia Pesquisar e tente uma pesquisa básica clicando no botão Pesquisar. Todos os registros são apresentados:

Se você clicar duas vezes em uma das camadas, verá ver as informações do registro:

Você também pode ver o pedido GetRecords realizado pelo pycsw clicando no botão “Ver pesquisa como XML”

Experimente refazer inserindo parâmetros diferentes para realizar consultas mais avançadas:

Hoje ficamos por aqui! Bons estudos!

Fonte: Paolo Corti
7:07
Blog 2 Engenheiros: Como selecionar e exportar Polígonos Sobrepostos no ArcGIS e QGIS (com e sem Python)?
sur Planet OSGeo
Quando temos nossa área de estudo, muitas vezes realizamos procedimentos de recorte para que todos elementos fiquem dentro desse limite.

Entretanto, em alguns momentos, é interessante que esse recorte não seja realizado exatamente nos limites da área de estudo. Queremos manter a forma do polígono ou linha que cruza nossa área de estudo.

Como podemos realizar essa cruzamento entre polígonos ou linhas com a nossa área de estudo?

Neste tutorial, iremos utilizar os shapefiles dos títulos minerários do estado de Santa Catarina (disponível no site da Agência Nacional de Mineração) e o shapefile das rodovias federais (disponível no site do Departamento Nacional de Infraestrutura e Transporte) para realizar esse procedimento.

Ambos podem ser baixados nos links abaixo:
- Títulos Minerais: Para baixar, vá no mapa no fim da página e selecione Santa Catarina, em seguida, clique em SHP;
- Rodovias Estaduais: Procure pelo shapefile do mapa rodiviário.
As rodovias federais podem ser baixadas no site do DNIT.

E vamos verificar quais títulos e rodovias cruzam a nossa área de estudo fictícia: Uma área circular com ponto central nas coordenadas SIRGAS UTM 700.731 / 7.058.005 e raio de 500 metros.

Leia também: Como gerar buffers no ArcGIS.

Lembre-se que para gerar buffers no QGIS, vá em Vetor (Vector) > Ferramentas de Geoprocessamento (Geoprocessing Tools) > Buffer.
Extraindo Shapefiles Sobrepostos no QGIS
Após baixar os shapefiles e criar a área de estudo fictícia no QGIS 3.2, iremos utilizar a ferramenta Seleção por Localização (Select by Location), a qual esta disponível em Vetor (Vector) > Ferramentas de Pesquisa (Research Tools).

Ao clicar nela, a janela abaixo será aberta.
Janela da Ferramenta Select by Location e Área de Estudo ao Lado.
Note que a nossa área de estudo (preenchida em vermelho) cruza com vários polígonos dos títulos minerários (linhas em preto) e com essa ferramenta, iremos separar esses polígonos (sem cortá-los, deixando-os inteiros).

Na ferramenta Select by Location, o primeiro item a ser inserido é o shapefile contendo os polígonos a serem selecionados (neste caso, os títulos minerários), em seguida, selecionamos qual tipo de operação (interseção) e baseada em quem essa operação irá ocorrer (área de estudo).

O ultimo item indica que o resultado deverá ser uma nova seleção. Veja que o resultado, apresentado abaixo, são todos os polígonos que “tocam” nossa área de estudo.
Polígonos selecionados encontram-se em amarelo.
Cabe lembrar que o processo não termina aqui. Agora precisamos salvar o resultado como shapefile, pois o que temos é apenas uma seleção.

Para salvar a seleção, clique sobre o shapefile com os itens selecionados (i.e. títulos minerários), vá em exportar (export), e clique em salvar feição como (save feature as).

Na nova janela, preencha os dados relacionados à localização do arquivo final, sistema de coordenadas e não se esqueça de marcar o item Salvar Apenas as Feições Selecionadas (Save only selected features), conforme figura abaixo.
Exportando polígonos selecionados como shapefile no QGIS.
Pronto! Você agora tem um novo shapefile com todos os títulos minerários “completos” que cruzam com a nossa área de estudo.
Extraindo Shapefiles Sobrepostos no ArcGIS
Após baixar e abrir os arquivos dos títulos minerários e da área de estudo, no ArcToolbox, vá em Data Management Tools > Layers and Table Views e dê dois cliques em Select Layer by Location.

Na janela que abrir, no campo “Input Feature Layer”, selecione o item que terá os polígonos extraídos (neste caso, os títulos minerários), em “Relationship (optional)”, você escolherá qual o processo a ser adotado (aqui, usaremos a interseção), “Selecting Features (optional)” será nossa área de estudo e por fim, em “Selection type (optional)” indica qual é o tipo de seleção que irá ocorrer (manteremos como nova seleção – New_Selection).
Ferramenta para Selecionar Feições pela localização no ArcGIS (e polígonos utilizados).
O resultado desta operação será a seleção dos títulos minerários que cruzam com a nossa área de estudo. A figura abaixo mostra os polígonos selecionados.
Polígonos que tocam nossa área de estudo foram selecionados (em azul claro).
Lembre-se que isto é apenas uma seleção. Precisamos agora exportar o resultado para termos um shapefile novo.

Para realizar esse processo, clique sobre o shapefile com os polígonos selecionados (i.e. Títulos Minerários) e vá em Data > Export Data. Na nova janela, você manterá o primeiro item (“Export:”) como “Selected features” e irá definir o local onde o novo shapefile será salvo.

O menu “Export:” define quem irá ser salvo, neste caso, ao marcar “Selected Features”, iremos exportar somente os arquivos selecionados.
Janela onde marcamos quais itens exportar e onde o shapefile será salvo.
E dessa forma, você terá os títulos minerários que cruzam a sua área de estudo.
Usando PyQGIS para Selecionar e Exportar Shapefiles
Antes de começar a escrever, abra o console Python no QGIS 3.x clicando em Plugins > Python Console. Após abrir o console, clique em “Show Editor” para abrir o editor.

Neste exemplo, vamos utilizar um shapefile diferente, utilizaremos as rodovias estaduais ao invés dos títulos minerários.

O código abaixo representa todo o procedimento que fizemos antes manualmente e após o código, explicamos os comandos utilizados.
```
# -*- coding: utf-8 -*-

# Carregando pacotes no Python
import processing
import sys
import osgeo.ogr as ogr
import osgeo.osr as osr

# Carregando nossos arquivos shapefile e criando saída
rod_est = "C:/Users/ferna/Desktop/B2E/Rodovias_SC.shp"
area_estudo = "C:/Users/ferna/Desktop/B2E/area_estudo.shp"
saida = "C:/Users/ferna/Desktop/B2E/rod_areaEst.shp"

# Abringo arquivos no QGIS
iface.addVectorLayer(rod_est, "Rodovia Estadual", "ogr")
iface.addVectorLayer(area_estudo, "Area de Estudo", "ogr")

# Selecionando rodovias que cruzam a área de estudo
# Documentação: processing.algorithmHelp('qgis:selectbylocation')
param_Sel = {'INPUT': rod_est, 'PREDICATE': 0,'INTERSECT': area_estudo, 'METHOD': 0}
processing.run('qgis:selectbylocation', param_Sel)
param_Save = {'INPUT': rod_est, 'OUTPUT': saida}
processing.run('qgis:saveselectedfeatures', param_Save)

# Carregando resultado no QGIS
iface.addVectorLayer(saida, "Rodovia", "ogr")
```
Na lista abaixo, explicamos as principais funções que utilizamos, de modo que você possa aplicar o código aos seus arquivos.
- iface.addVectorLayer( [Local do Shapefile] , [Nome] , “ogr”): Adiciona e abre o shapefile ao QGIS;
- processing.run( [Algoritmo para ser Rodado] , [Dicionário com os Parâmetros]).
Cabe lembrar que o comando processing.run() é diferente no QGIS 2.18, sendo definido como processing.runalg() e recebe argumentos diferentes.

Após rodar o nosso código, você obterá quais rodovias estaduais cruzam a nossa área de estudo.
Usando ArcPy para Selecionar e Exportar Shapefiles
No ArcGIS 10.x, clique em “Geoprocessing” > “Python”, isso irá abrir o terminal Python do ArcGIS.

O ArcGIS tem pequeno editor de python integrado ao software, nossa sugestão é que você trabalhe inicialmente em um editor de texto externo (e.g. NotePad++) e depois cole o código no terminal.

Nesta parte, iremos utilizar as rodovias estaduais que mostramos lá no começo, ao invés dos títulos minerários.

O código para executar a extração das rodovias que cruzam nossa área de estudo é apresentado abaixo.
```
# -*- coding: utf-8 -*-

# Carregando pacotes no Python
import arcpy

# Nossas variáveis
rod_est = "C:\\Users\\ferna\\Desktop\\B2E\\Rodovias_SC.shp"
area_estudo = "C:\\Users\\ferna\\Desktop\\B2E\\area_estudo.shp"
saida = "C:\\Users\\ferna\\Desktop\\B2E\\rod_areaEstArc.shp"

# Abringo arquivos no QGIS
mxd = arcpy.mapping.MapDocument("CURRENT") # Define qual projeto estamos trabalhando
dataFrame = mxd.activeDataFrame # Define o dataframe que estamos trabalhando
camada_ae = arcpy.mapping.Layer(area_estudo)

arcpy.mapping.AddLayer(dataFrame, camada_ae, "AUTO_ARRANGE") # Adiciona Camada

# Selecionando rodovias que cruzam a área de estudo
arcpy.MakeFeatureLayer_management(rod_est, "Rodovias1")
arcpy.SelectLayerByLocation_management("Rodovias1", "intersect", area_estudo)
arcpy.CopyFeatures_management("Rodovias1", saida)
```
As funções que utilizamos são explicados na lista abaixo:
- arcpy.mapping.MapDocument(“CURRENT”): Define como projeto padrão o projeto aberto;
- mxd.activeDataFrame: Define como data frame a ser utilizado o data frame aberto;
- arcpy.mapping.Layer(): Transforma o caminho com o shapefile em uma camada;
- arcpy.mapping.AddLayer( [Data Frame] , [Camada] , “AUTO_ARRANGE”): Adiciona camada à tabela de conteúdo (abre o arquivo);
- arcpy.MakeFeatureLayer_management( [Shapefile], [Arquivo Temporário]): Torna a camada selecionada como principal;
- arcpy.SelectLayerByLocation_management( [Arquivo Temporário], [Método], [Shapefile para Cruzamento]): Realiza procedimento para selecionar arquivos que se cruzam;
- arcpy.CopyFeatures_management([Arquivo Temporário], [Saída]): Salva os itens selecionados no caminho especificado.
E com isso, você terá as rodovias que cruzam nossa área de estudo.

Ficou com alguma dúvida? Utilize a seção de comentários e caso esse tutorial tenha te ajudado, compartilhe ele na sua rede social.
19 novembre 2018
21:24
QGIS Blog: Results of the MacOS bug fixing initiative
sur Planet OSGeo
Thanks to your donations, we were able to hire core developers to focus on solving Mac OS specific issues for QGIS. More than 30 MacOS QGIS users donated a little more than 3000 € for this bug fixing round.

After an effort of triage and testing, here is what has been achieved:
- The map canvas is now in HiDPI ( [https:]] ).
  The issue was that Mac OS handles differently high DPI scaling by providing a scaling ration to scale graphics output. For more information, read [doc.qt.io] and for even more info [doc.qt.io] . You may not have noticed, but the scale shown under the canvas was wrong in former QGIS 3.x versions. It has been fixed too.
- Zooming with trackpad is now fluent ( [https:]] ) The track pad is apparently sending more events than the scrolling and are now discarded to avoid scale jumps.
- Size of mouse cursors with map tools has been fixed ( [https:]] ). The scaling previously introduced could be removed thanks to the HiDPI fixes, and this has fixed by the same occasion the active point of the cursor being misplaced.
- The size of the map canvas at startup was not an issue in QGIS in itself and has been fixed by upgrading to a newer Qt version ( [https:]] )
- Rendering slowness issue *should* also have been fixed by an upgrade of Qt version ( [https:]] )
- At the moment, the “default” theme of QGIS (as opposed to the “Night Mapping” theme), has a few glitches when used in Mojave with “Dark” theme activated. Mainly text in widgets cannot be read since they are white on white. This is a Qt issue ( [https:]] ) and the upcoming Qt 5.12 released by the end of november should fix it. In the mean time, the QGIS “Night Mapping” theme is automatically applied when running on Mojave with Dark theme. While the theme is probably not perfect, it allows to work properly and the situation will be evaluated again when Qt 5.12 is released.
Unfortunately, some issues remain. Mainly, the text being rendered as outlines in PDF export ( [https:]] ) remains for now. It might be fixed in a following effort.

Thanks to all donors who helped in this effort and to Denis Rouzaud as a core developer who spent a lot of time investigating and fixing these issues!
12:10
gvSIG Team: La Suite gvSIG en Gestión Municipal: grabaciones de las ponencias de las Jornadas gvSIG disponibles
sur Planet OSGeo

En los últimos años la integración de las tecnologías de la información y la comunicación (TIC) en la gestión urbana han llevado a popularizar términos como el de ‘Smart City’ o ‘Ciudad Inteligente’. Un término en muchas ocasiones utilizado meramente como marketing y vacío de contenido. Desde la Asociación gvSIG siempre hemos defendido e impulsado la idea de que las ciudades sólo puede ser ‘inteligentes’ si se tienen políticas definidas en relación a la soberanía tecnológica y la sostenibilidad, sólo alcanzables con la adopción de tecnologías libres.

Hoy os traemos unos ejemplos de cómo las soluciones en software libre de la Asociación gvSIG mejoran la gestión municipal y pueden llegar a convertirse en herramientas imprescindibles en el día a día de un Ayuntamiento, independientemente de su tamaño.

En primer lugar el arquitecto Nacho Muñoz-Motilla nos cuenta la experiencia del Ayuntamiento de La Pobla de Vallbona, en ‘Creación de una IDE de gestión municipal con gvSIG Online en tiempo récord’

¿Pueden los pequeños ayuntamiento sumarse a la modernización que proponen las Smart Cities o han de quedar fuera? Pepe casado nos habla de la experiencia del Ayuntamiento de Villarcayo de Merindad de Castilla la Vieja en ‘gvSIG Online, como crear una IDE en un pequeño ayuntamiento’

Por último Antonio García Benlloch nos cuenta como el Ayuntamiento de Bétera está combinando el uso de distintos productos de la Suite gvSIG en ‘Captura de datos de campo con gvSIG Mobile, edición con gvSIG Desktop y publicación web en el Geoportal del Ayuntamiento de Bétera’
11:30
Fernando Quadro: A sua cidade é uma Smart City?
sur Planet OSGeo
Não há como negar que cidades inteligentes (Smart Cities) é um assunto da moda, e que muito tem se falado nisso no Brasil! Mas a pergunta que eu gostaria de lhe fazer é: Sua cidade pode ser considerada uma Smart City?

Segundo uma publicação da Juniper Research em 2017 “Uma cidade inteligente é caracterizada pela integração da tecnologia em uma abordagem estratégica para a sustentabilidade, o bem-estar dos cidadãos e o desenvolvimento econômico.”

Neste post, explorarei preceitos e visões comuns sobre o que é uma cidade inteligente, e no decorrer do post, você terá condições de responder minha pergunta inicial.

Uma cidade muitas vezes ganha seu status de “inteligente” quando processos específicos, como gerenciamento de informações e comunicações, se interconectam e usam tecnologias avançadas, como IoT (Internet of Things: dispositivos físicos conectados uns com os outros e com a internet) para melhorar suas operações e serviços. Esses serviços geralmente são baseados em automação e tecnologias inteligentes para gerenciar serviços públicos, como energia, abastecimento de água e coleta de lixo.

1. Como medimos se uma cidade é inteligente?

O tópico de como identificar ou medir como uma cidade inteligente é amplamente debatido, com muitas organizações fornecendo uma lista definitiva de critérios ou métricas “inteligentes” para as cidades. Aqui estão alguns exemplos de CITYkeys, European Smart Cities e Lendlease.

Neste blog, usaremos os critérios de um paper publicado pela Juniper Research em 2017 para medir o quão inteligente é uma cidade. No paper, a Juniper Research estabeleceu um sistema único de pontuação para analisar e avaliar a “inteligência” das cidades. A equipe mede a atividade de cada cidade observando os seguintes pontos de dados:
- Transportes: Quão conectados e eficientes são os serviços? A cidade coleta dados do monitoramento de tráfego em tempo real e usa esses dados para informar os ajustes de fluxo de tráfego e sinal, com base nos requisitos de resposta a emergências?
- Saúde: Quão inovadora e impactante é a abordagem dos serviços? A cidade tem regulamentos de emissão de ar e soluções de transporte saudáveis?
- Segurança pública: Como a tecnologia impacta a cidade para a segurança pública? A iluminação pública inteligente foi implantada? A análise inteligente de vigilância por vídeo ou software de crime preditivo e de risco de incêndio estão sendo usados?
- Produtividade: Quão colaborativo e co-criativo é a estratégia inteligente da cidade? Eles estão organizando conferências e hackathons para reunir diferentes partes interessadas da cidade? A cidade incentiva talentos que contribuam para a economia?
- Setores de energia: Quão avançado é o setor de serviços urbanos? A cidade está aplicando tecnologias, como sistemas de semáforo inteligente controlados por Inteligência Artificial (AI), soluções de armazenamento de energia em casa, painéis solares e infraestrutura para carregamento de EV?
A partir da avaliação acima, uma cidade é inteligente com base no fato de implementarem tecnologias avançadas para otimizar seus processos, custos e ofertas de serviços. Eu não preciso saber onde você mora aqui no Brasil, e nem ser vidente pra saber que sua cidade hoje (em 2018) provavelmente não é uma cidade inteligente.

2. Por que nossas cidades não são inteligentes?

Para responder essa pergunta, irei fazer algumas outras perguntas (e respostas) para você poder entender melhor os motivos.

a. A inovação é sempre visível?
A inovação pode assumir muitas formas diferentes. Muitas vezes ouvimos falar das tecnologias empolgantes e chamativas que estão em destaque nas manchetes do meios de comunicação, mas menos empolgantes, mas igualmente impactantes, são inovações como uma nova maneira de organizar uma equipe ou um novo sistema de gerenciamento.

b. As nossas expectativas sobre o que uma cidade inteligente é alta demais?
Todos os dias, ouvimos falar de uma nova tecnologia que melhorará drasticamente nossas vidas – da realidade virtual à tecnologia de reconhecimento facial e serviços robóticos. Essas tecnologias podem estar se tornando mais comuns no setor privado, mas temos expectativas reais de quão rapidamente essas tecnologias podem ser integradas em nossos serviços públicos e infraestrutura?

c. Temos observado alguma melhoria em nossos serviços de coleta de resíduos?
Por exemplo, canais de comunicação mais inovadores, como aplicativos e chatbots? Podemos considerar estes avanços em tecnologia no setor público, com a mesma velocidade com que novas tecnologias vem sendo aplicadas no setor privado?

d. As medições de uma cidade inteligente consideram a experiência do cidadão?
Os cidadãos têm diferentes prioridades para as autoridades municipais e a agenda da cidade inteligente – eles se preocupam em viver suas vidas – sobre suas famílias, suas carreiras e seu bem-estar. Talvez tudo o que eles esperam é que os serviços pelos quais pagam (através de impostos) atendam às suas necessidades e trabalhem perfeitamente ao lado de suas vidas ocupadas. Talvez esses serviços pudessem:
- Incorporar a IA (Inteligencia Artificial) nos smartphones para fornecer orientação se eles se perderem ou sofrerem atrasos durante seus deslocamentos?
- Usar VR (realidade virtual) para adicionar um valor experiencial ao sistema de transporte da cidade?
- Identificar um médico adequado para fornecer uma consulta, sem longos períodos de espera ou vários encaminhamentos?
- Serviços interconectados, para que eles não precisem fornecer repetidamente as mesmas informações em todos os lugares (ou seja, um sistema “único” em vez de um sistema desconectado)?
- Fornecer interfaces e dispositivos digitais acessíveis e intuitivos para pessoas de todas as idades para gerenciar suas utilidades e obter insights sobre seu uso de energia?
3. Como as nossas cidades podem ser inteligentes para os cidadãos?

As cidades existem para atender às necessidades das pessoas – sejam elas de trabalho, educação, sistema de saúde ou cultura. As necessidades das pessoas devem estar no centro dos nossos serviços públicos e privados. Apesar de uma tendência crescente para o design centrado no usuário, as empresas ainda tendem a ditar a agenda da cidade. Como podemos reorientar sua abordagem para garantir que as necessidades dos cidadãos sustentem o design e a execução dos serviços da cidade?

E se as autoridades e empresas da cidade utilizassem uma abordagem de baixo para cima para a inovação – garantindo que as ideias, os processos e as tecnologias se baseassem na contribuição dos cidadãos? Isso pode não apenas resultar em serviços mais inteligentes, mas também produzir uma melhor experiência do usuário para os cidadãos, que se sentiriam mais capacitados em como as cidades são definidas – em vez dos grandes provedores de serviços que atualmente dominam o mercado. Um passo para isso é uma estratégia de implementação de longo prazo que instrui as pessoas que executam esses serviços e estão em contato direto com o usuário. Ao incorporar sua opinião e experiência, poderíamos oferecer uma abordagem mais focada no ser humano, refinada e fundamentada em pesquisa?

4. A solução pode ser encontrada no design de serviços?

O que as “cidades inteligentes” podem fazer para implementar e melhorar sua experiência de serviço? E se os seguintes princípios sustentassem nossas cidades?
- Foco no Usuário: O cidadão deve estar no centro de todo serviço.
- Co-Criatividade: Cidadãos, setores público e privado devem estar envolvidos no processo.
- Sequenciamento: Os serviços devem ser visualizados por seqüências ou momentos-chave do cliente, é importante entender o que e quando as coisas acontecem. Antes, durante e depois do serviço.
- Evidenciamento: Todos os envolvidos no fornecimento do serviço precisam estar cientes da experiência do serviço.
- Holística: O design de um serviço deve levar em conta toda a experiência do serviço. O contexto é importante.
No Reino Unido, o Governo (Gov.uk) já implementou com sucesso uma estrutura de design de serviço, que segue uma abordagem de design focado no ser humano. Isso é demonstrado pelos Serviços Digitais do Governo (GDS) que estão liderando a transformação dos serviços do governo do Reino Unido de físicos para digitais. A GDS tem seu próprio conjunto de princípios de design, mas, em geral, sua estratégia de processo e serviço reflete os 5 princípios acima mencionados por Marc Stickdorn. É geralmente aceito na indústria do design que você deve “sempre começar com o usuário”.

Uma solução é que as cidades inteligentes aspirantes considerem seu próprio conjunto de princípios. Isso permitiria aos provedores de serviços do setor público introduzir uma abordagem focado no usuário, mas também ofereceria um conjunto de princípios a serem seguidos por empresas privadas. Uma questão a considerar é como poderíamos aplicar essa abordagem a todos os serviços da cidade?

Embora o governo seja um grande exemplo para o setor público e privado, a maioria dos serviços municipais não é apenas pública e, portanto, não há padrões para garantir a qualidade nos serviços urbanos. Isso parece especialmente importante para grandes cidades inteligentes como Londres ou Nova York, que lutam para competir com cidades inteligentes menores, como Copenhague ou Amsterdã, já que o tamanho populacional maior interfere nos serviços e nas experiências que uma cidade inteligente oferece.

Serviços como energia, educação e transporte são parcialmente substituídos por organizações privadas e muitas vezes não colocam o usuário em primeiro lugar ou pensam em uma experiência de serviço holística antes de criar seus serviços. Organizações privadas e públicas têm uma agenda diferente; a questão é como eles podem se tornar mais alinhados?

5. Conclusões

Como podemos ver no artigo, se tornar uma cidade realmente inteligente não é algo fácil, e não acontece da noite pro dia. É um trabalho que precisa de tempo e muito esforço para se materializar. Eu não conheço todas as cidades do Brasil, mas as que conheço não consigo ver nenhuma como um cidade inteligente.

Meu objetivo neste post, foi mostrar pra você que pra ser uma cidade inteligente é muito mais complexo do que apenas liberar acesso “Wi-fi” para a população, é um caminho duro e tortuoso até poder chegar ao nível de se intitular uma cidade de “inteligente”.

Fonte: Organicity
18 novembre 2018

16:48

Free and Open Source GIS Ramblings: Movement data in GIS #16: towards pure Python trajectories using GeoPandas

sur Planet OSGeo

Many of my previous posts in this series [1][2][3] have relied on PostGIS for trajectory data handling. While I love PostGIS, it feels like overkill to require a database to analyze smaller movement datasets. Wouldn’t it be great to have a pure Python solution?

If we look into moving object data literature, beyond the “trajectories are points with timestamps” perspective, which is common in GIS, we also encounter the “trajectories are time series with coordinates” perspective. I don’t know about you, but if I hear “time series” and Python, I think Pandas! In the Python Data Science Handbook, Jake VanderPlas writes:

Pandas was developed in the context of financial modeling, so as you might expect, it contains a fairly extensive set of tools for working with dates, times, and time-indexed data.

Of course, time series are one thing, but spatial data handling is another. Lucky for us, this is where GeoPandas comes in. GeoPandas has been around for a while and version 0.4 has been released in June 2018. So far, I haven’t found examples that use GeoPandas to manage movement data, so I’ve set out to give it a shot. My trajectory class uses a GeoDataFrame df for data storage. For visualization purposes, it can be converted to a LineString:

import pandas as pd 
from geopandas import GeoDataFrame
from shapely.geometry import Point, LineString

class Trajectory():
    def __init__(self, id, df, id_col):
        self.id = id
        self.df = df    
        self.id_col = id_col
        
    def __str__(self):
        return "Trajectory {1} ({2} to {3}) | Size: {0}".format(
            self.df.geometry.count(), self.id, self.get_start_time(), 
            self.get_end_time())
        
    def get_start_time(self):
        return self.df.index.min()
        
    def get_end_time(self):
        return self.df.index.max()
        
    def to_linestring(self):
        return self.make_line(self.df)
        
    def make_line(self, df):
        if df.size > 1:
            return df.groupby(self.id_col)['geometry'].apply(
                lambda x: LineString(x.tolist())).values[0]
        else:
            raise RuntimeError('Dataframe needs at least two points to make line!')

    def get_position_at(self, t):
        try:
            return self.df.loc[t]['geometry'][0]
        except:
            return self.df.iloc[self.df.index.drop_duplicates().get_loc(
                t, method='nearest')]['geometry']

Of course, this class can be used in stand-alone Python scripts, but it can also be used in QGIS. The following script takes data from a QGIS point layer, creates a GeoDataFrame, and finally generates trajectories. These trajectories can then be added to the map as a line layer.

All we need to do to ensure that our data is ordered by time is to set the GeoDataFrame’s index to the time field. From then on, Pandas takes care of the time series aspects and we can access the index as shown in the Trajectory.get_position_at() function above.

# Get data from a point layer
l = iface.activeLayer()
time_field_name = 't'
trajectory_id_field = 'trajectory_id' 
names = [field.name() for field in l.fields()]
data = []
for feature in l.getFeatures():
    my_dict = {}
    for i, a in enumerate(feature.attributes()):
        my_dict[names[i]] = a
    x = feature.geometry().asPoint().x()
    y = feature.geometry().asPoint().y()
    my_dict['geometry']=Point((x,y))
    data.append(my_dict)

# Create a GeoDataFrame
df = pd.DataFrame(data).set_index(time_field_name)
crs = {'init': l.crs().geographicCrsAuthId()} 
geo_df = GeoDataFrame(df, crs=crs)
print(geo_df)

# Test if spatial functions work
print(geo_df.dissolve([True]*len(geo_df)).centroid)

# Create a QGIS layer for trajectory lines
vl = QgsVectorLayer("LineString", "trajectories", "memory")
vl.setCrs(l.crs()) # doesn't stop popup :(
pr = vl.dataProvider()
pr.addAttributes([QgsField("id", QVariant.String)])
vl.updateFields() 

df_by_id = dict(tuple(geo_df.groupby(trajectory_id_field)))
trajectories = {}
for key, value in df_by_id.items():
    traj = Trajectory(key, value, trajectory_id_field)
    trajectories[key] = traj
    line = QgsGeometry.fromWkt(traj.to_linestring().wkt)
    f = QgsFeature()
    f.setGeometry(line)
    f.setAttributes([key])
    pr.addFeature(f) 
print(trajectories[1])

vl.updateExtents() 
QgsProject.instance().addMapLayer(vl)

The following screenshot shows this script applied to a sample of the Geolife datasets containing 100 trajectories with a total of 236,776 points. On my notebook, the runtime is approx. 20 seconds.

So far, GeoPandas has proven to be a convenient way to handle time series with coordinates. Trying to implement some trajectory analysis tools will show if it is indeed a promising data structure for trajectories.

1:00
PostGIS Development: PostGIS 2.5.1
sur Planet OSGeo
The PostGIS development team is pleased to provide bug fix 2.5.1 for the 2.5 stable branch.

Although this release will work for PostgreSQL 9.4 thru PostgreSQL 11, to take full advantage of what PostGIS 2.5 offers, you should be running PostgreSQL 11 and GEOS 3.7.0.

Best served with PostgreSQL 11.1 and pgRouting 2.6.1.

WARNING: If compiling with PostgreSQL+JIT, LLVM >= 6 is required Supported PostgreSQL versions for this release are: PostgreSQL 9.4 - PostgreSQL 11 GEOS >= 3.5

2.5.1
- source download
- NEWS
- PDF docs en
- html doc download
- html online en, de, ko_KR, ja
- epub doc download
- pgsql help files for non-english languages: de, ja,kr, br, es
- Changelog
Continue Reading by clicking title hyperlink ..
16 novembre 2018
14:41
gvSIG Team: Disponible el material del taller de ‘Generación de informes con gvSIG Desktop’ de las 14as Jornadas gvSIG
sur Planet OSGeo

Ya está disponible la grabación del taller sobre ‘Generación de informes con gvSIG Desktop’ impartido durante las 14as Jornadas Internacionales gvSIG.

El taller, enfocado a preparar capas y mapas en gvSIG para ser integrados en el diseñador de informes JasperSoft Studio y a crear documentos de tipo ‘Informe’ en gvSIG que incluyesen los informes creados con JasperSoft Studio, se realizó sobre la última versión de gvSIG (gvSIG 2.5), que aún está en desarrollo, para el cual se ha desarrollado este nuevo plugin. No obstante, ya hay una distribución para descargar y probar, incluyendo el creador de informes.

Para iniciar el taller es necesario disponer de:
– gvSIG Desktop 2.5.
– JasperSoft Studio: Versión oficial.
– Plugin para JasperSoft Studio de Web Service Data Source (descarga mediante registro), o enlace alternativo.

Se deberán de buscar los tres ficheros .jar del plugin que se encuentran en la carpeta “jss” y copiarlos en la carpeta de JasperSoft/plugins.

Los datos para seguir el taller pueden descargarse de la web oficial (o descarga alternativa).

El vídeo con la grabación del taller es el siguiente:
9:39
gvSIG Team: Material of the workshop about environmental modelling using gvSIG and the HortonMachine available
sur Planet OSGeo

The material of the workshop about environmental modelling using gvSIG and the HortonMachine, given at 14th International gvSIG Conference some days ago, is now available.

Although the screen image is not very clear it can be followed with the teacher instructions, and using the cartography and tutorial available to download at this link.

The HortonMachine library is an open geospatial library integrated in gvSIG containing tools for the management of environmental data and hydro-geomorphological analyses. The algorithms contained in the library are the result of more than 10 years of research, development and real application by different researchers and professionals working in the field of natural hazards.

Recording of this workshop is available here:
15 novembre 2018
15:13
gvSIG Team: Vídeos de sesión SIG aplicado a Agricultura de las 14as Jornadas gvSIG
sur Planet OSGeo

Uno de los sectores en los que estamos realizando cada vez más proyectos desde la Asociación gvSIG es el de la Agricultura. Un área donde la aplicación de la geomática es amplia y que engloba proyectos tanto de gestión como de explotación o investigación.

Durante las pasadas 14as Jornadas Internacionales de gvSIG hubo una sesión dedicada a este tipo de proyectos, de la que ahora os compartimos los vídeos de las presentaciones realizadas.

En primer lugar queremos destacar la ponencia ‘Soluciones gvSIG para Agricultura en la Generalitat Valenciana.’. Muestra como el software libre se está aplicando para poner en marcha aplicaciones tan importantes como las que conllevan la gestión de las ayudas de la PAC (la Política Agraria Común) y la gestión del Registro general de la producción agrícola (REGEPA). Aplicaciones con un alto nivel de complejidad técnica ya que deben implementar todo tipo de reglas topológicas.

En segundo lugar un proyecto de investigación muy interesante llevado a cabo por Ángel Marqués-Mateu con título ‘Cartografía de precisión de la salinidad del suelo en cultivos de arroz.’

Para finalizar la presentación de GODAN, ‘The role of open data and innovation in ensuring global food security’ que nos remarca la importancia que pueden tener los datos abiertos y el software libre.

¡Qué los disfrutéis!
14 novembre 2018
16:28
gvSIG Team: El Circuit de Cheste estrenará un sistema pionero para dar respuesta a emergencias durante el Gran Premio de Motociclismo
sur Planet OSGeo

Bajo este titular ha sido referenciado hoy en la práctica totalidad de la prensa el sistema impulsado por Diputació de València, a través del área de Información Territorial de Divalterra y el Consorcio Provincial de Bomberos y desarrollado en colaboración con la Asociación gvSIG. La plataforma, basada en gvSIG Online, fue presentada en las pasadas 14as Jornadas Internacionales de gvSIG. Ver como la Suite gvSIG está formando parte de importantes proyectos, por dimensión y por su influencia en mejorar la vida de los ciudadanos, es algo que nos produce una enorme satisfacción.

Os dejamos con la información extraída de la nota de prensa oficial y con un vídeo de presentación de la solución:

El Circuit Ricardo Tormo de Cheste estrenará el fin de semana del 16 al 18 de noviembre, durante la celebración del Gran Premio de Motociclismo, un sistema pionero de gestión de emergencias para prevenir y dar la respuesta más ágil a los incidentes que puedan registrarse tanto en el recinto como en el entorno del Circuit, incluidas las carreteras de acceso y los núcleos de población cercanos.

Se trata de un servidor cartográfico que utiliza tecnología GIS e incluye los datos procedentes del Plan de Autoprotección del Circuito, y se presenta de manera digital a través de cualquier equipo informático. De esta manera, se agilizan las consultas y se dispone al momento de todo tipo de información de utilidad para dar respuesta a las emergencias. El servidor cartográfico permitirá agilizar la gestión de cualquier incidente de seguridad. De esta manera, se agilizan las consultas y se dispone al momento de todo tipo de información de utilidad para dar respuesta a las emergencias.

Además, desde el Consorcio Provincial de Bomberos se dará acceso a esta herramienta informática a todos los organismos implicados en garantizar la seguridad en la concentración deportiva, como Delegación de Gobierno, Guardia Civil, Generalitat Valenciana y la propia seguridad del Circuit, con el fin de asegurar una máxima coordinación entre todas las partes.

El oficial jefe de Prevención del Consorcio Provincial de Bomberos de Valencia, Jorge Sánchez, explica que “la cantidad de información que nos aporta, agiliza de manera muy importante la toma de decisiones en momentos en los que cada segundo que se gane puede ser vital”.

Sánchez indica que “en España no hay ningún organismo con una herramienta GIS tan potente como esta para la gestión de emergencias” y avanza que “ya se está trabajando para aplicarla a otras infraestructuras como los edificios de la Universidad, Metrovalencia, las Fallas, etc”.

El director del Proyecto, Antonio Mas, explica que “el proceso se ha desarrollado en completa coordinación con el Consorcio de Bomberos, para asegurar que las funcionalidades que se iban desarrollando desde Divalterra eran las que realmente resultaban de utilidad para ellos”.

Para ello, se recopiló toda la documentación existente, que estaba en formato papel, y tras la definición de los términos y cartografía a implementar, se procedió a elaborar todos los mapas en formato digital. Una vez implementada la parte gráfica, se procedió al volcado de toda la documentación recopilada.

Así, recoge imágenes y datos detallados sobre los accesos y salidas del recinto, elementos de protección, recorrido hasta los accesos, rutas de evacuación, gradas, locales, recintos y edificios, y también todas las actividades de riesgo de empresas próximas, así como el posible transporte de mercancías peligrosas.

De esta manera, este servidor cartográfico permite, ante una emergencia, saber en qué punto exacto se ha producido el incidente y a qué distancia se encuentran los accesos y salidas más cercanas, cuál es el recorrido, los hidrantes más próximos, las características de la infraestructura donde se ha producido, si existen elementos colindantes o las zonas industriales próximas y sus planes de evacuación, entre otra información de utilidad.
11:30
Fernando Quadro: Usando pycsw com o GeoNode
sur Planet OSGeo
O CSW é um padrão para expor um catálogo de entidades geoespaciais através do procotolo [HTTP.] Em um portal GeoNode, os pontos de extremidade (endpoints) do CSW são fornecidos pelo pycsw, que é um componente subjacente do GeoNode. Alternativamente, se necessário, é possível substituir o pycsw pelo GeoNetwork.

No GeoNode, você pode acessar facilmente o registro CSW de uma camada, clicando no botão Baixar Metadados na página da camada. Um formulário aparecerá e você poderá acessar os metadados fornecidos pelo pycsw em uma série de diferentes formatos (Atom, Dublin, FGDC, Text, HTML e muitos outros).

Por exemplo, ao clicar no link ISO, você acessará os metadados da camada no formato ISO, que corresponde a essa solicitação GetRecordById no pycsw:
```
http://localhost:8000/catalogue/csw?outputschema=http%3A%2F%2Fwww.isotc211.org%2F2005%2Fgmd&service=CSW&request=GetRecordById&version=2.0.2&elementsetname=full&id=8bcf5bfc-5cfc-11e7-8103-02d8e4477a33
```
Você também pode notar outras informações que foram geradas pelo GeoNode nos bastidores quando a camada foi carregada:
- Identificador da camada, que identifica exclusivamente a camada no catálogo (observe que a solicitação GetRecordById usa esse identificador para acessar o registro)
- Data de criação
- Sistema de referência espacial e caixa delimitadora (BBOX)
- URL da miniatura
- Formato do recurso
- Vários endpoints do OGC
Se você quiser adicionar metadados ausentes, visite a página de metadados da camada e pressione em Editar Camada > Editar Metadados.

1. Operações pycsw

O pycsw implementa todas as operações do padrão CSW, incluindo as opcionais:
- GetCapabilities: recupera metadados de serviço do servidor
- DescribeRecord: permite que um cliente descubra elementos do modelo de informações suportado pelo serviço de catálogo de destino
- GetRecords: procura registros usando uma série de critérios
- GetRecordById: recupera metadados para um registro (camada) do catálogo por seu id
- GetDomain (opcional): recupera informações de tempo de execução sobre o intervalo de valores de um elemento de registro de metadados ou um parâmetro de solicitação
- Harvest (opcional): cria / atualiza metadados pedindo ao servidor para “puxar” metadados de algum lugar
- Transaction (opcional): criar / editar metadados “empurrando” os metadados para o servidor
O pycsw é uma implementação do serviço OGC CSW escrita em Python. Iniciado em 2010 (mais formalmente anunciado em 2011), o pycsw permite a publicação e descoberta de metadados geoespaciais por meio de várias APIs (CSW 2 / CSW 3, OpenSearch, OAI-PMH, SRU), fornecendo um componente de metadados e catálogo baseado em padrões para infraestruturas de dados espaciais. O pycsw é Open Source, lançado sob uma licença MIT e executado em todas as principais plataformas (Windows, Linux e Mac OS X).

Fonte: Paolo Corti
13 novembre 2018
11:30
Fernando Quadro: Mover o PostGIS de esquema
sur Planet OSGeo
A partir do PostGIS 2.3, a extensão do postgis foi alterada para não permitir mais a realocação. Todas as chamadas de função agora são qualificadas pelo esquema.

Embora essa alteração tenha corrigido alguns problemas com a restauração do banco de dados, ela criou o problema, pois se você instalou o PostGIS em um esquema diferente daquele que você desejava, não é intuitivo como movê-lo para um esquema diferente. Felizmente há uma maneira de fazer isso.

Para este exemplo, o PostGIS foi instalado no esquema padrão para demonstrar como movê-lo para outro esquema. Você pode executar estas etapas usando o psql, pgAdmin ou qualquer outra ferramenta desejada.

A maioria das pessoas tem seu esquema padrão configurado como public, caso você não especifique explicitamente um esquema na instalação, geralmente o PostGIS será instalado no esquema public.
```
CREATE EXTENSION postgis;
```
Agora vou criar um novo esquema para movê-lo e adicionar este esquema ao search_path:
```
CREATE SCHEMA postgis;
 
ALTER DATABASE mydb 
SET search_path = public, postgis;
```
Se você estiver executando o PostGIS 2.3 ou superior, tente mudar para um esquema diferente usando o comando abaixo:
```
ALTER EXTENSION postgis 
  SET SCHEMA postgis;
```
Isso irá falhar e o banco irá plhe apresentar o seguinte erro “ERRO: a extensão “postgis” não suporta SET SCHEMA”.

Para permitir o movimento, siga estes passos:
```
UPDATE pg_extension 
  SET extrelocatable = TRUE 
    WHERE extname = 'postgis';
 
ALTER EXTENSION postgis 
  SET SCHEMA postgis;
 
ALTER EXTENSION postgis 
  UPDATE TO "2.4.1next";
 
ALTER EXTENSION postgis 
  UPDATE TO "2.4.1";
```
Observe o uso da expressão next. É necessário utilizar a expressão next nesta etapa para atualizar todas as referências das funções do esquema atual para o novo do esquema. Next é projetado para permitir a atualização de uma extensão do postgis para uma versão já existente. Tentar executar UPDATE TO “2.4.1” quando você já está no 2.4.1 acionaria um erro, pois você já está nessa versão.

Fonte: PostGIS

9:00

Marco Bernasocchi: Visualize Postgres JSON data in QML widgets

sur Planet OSGeo

As promised some time ago in “The new QML widgets in QGIS – When widgets get unbridled” we still owe you some fancy unicorns, but first let’s have a look at another nice feature that has been introduced in QGIS 3.4 LTR, the reading of PostgreSQL JSON and JSONB types.
With JSON you have a lot of possibilities for storing unstructured data. In our case, it’s mainly interesting when the data are stored as an array or a JSON object. Let’s have a look at two examples.

Visualize Postgres JSON data with common widgets

With the usual QGIS widgets “List” and “Key/Value” you are able to display JSON arrays and simple JSON objects.

JSON array as List

[
    "European dark bee",
    "Carniolan honey bee",
    "Buckfast bee"
]

Simple JSON object as Key/Value

{
    "nomenclatura":"Apis mellifera mellifera",
    "name":"European dark bee",
    "link":"https://en.wikipedia.org/wiki/European_dark_bee"
}

Or of course both as plain text in the “Text Edit” widget:

Say hi to Postgres JSON in QML widget

Probably, your JSON data does not look really nice with the aforementioned widgets, luckily since QGIS 3.4, you are free to create your own QML widget. Since QGIS already loads the JSON data into structures that are supported by QML, we can use all the JSON data within the QML code.
Let’s assume you have the JSON array from above and you like the elegance of the blue of Jacques Majorelle. You create your personal list widget by adding the JSON field as an expression:

import QtQuick 2.0
Rectangle {
    width: 310; height: 250; color: "grey"
    Column {
        anchors.horizontalCenter: parent.horizontalCenter
        anchors.verticalCenter: parent.verticalCenter
        spacing: 5
        Repeater {
            model:expression.evaluate("\"jvalue\"")
            Rectangle {
                color: "#6050dc"
                width: 300; height: 50; radius: 10.0
                Text {
                    anchors.centerIn: parent
                    font.pointSize: 24
                    text: modelData
                }
            }
        }
    }
}

You will have your very personal list:

JSON also allows storing more complex data, like for example a list of objects. In that case, you will reach the limits of the common QGIS widgets.
Let’s assume you have a table looking like this:

nomenclatura	name	link
Apis mellifera mellifera	European dark bee	https://en.wikipedia.org/wiki/European_dark_bee
Apis mellifera carnica	Carniolan honey bee	https://en.wikipedia.org/wiki/Carniolan_honey_bee
Apis mellifera	Buckfast bee	https://en.wikipedia.org/wiki/Buckfast_bee

In JSON it would be stored like this:

[
    {"nomenclatura":"Apis mellifera mellifera","name":"European dark bee","link":"https://en.wikipedia.org/wiki/European_dark_bee"},
    {"nomenclatura":"Apis mellifera carnica","name":"Carniolan honey bee","link":"https://en.wikipedia.org/wiki/Carniolan_honey_bee"},
    {"nomenclatura":"Apis mellifera","name":"Buckfast bee","link":"https://en.wikipedia.org/wiki/Buckfast_bee"}
]

With the QML Widget you can use the QML TableView to visualize:

import QtQuick 2.0
import QtQuick.Controls 1.4
TableView {
    width: 600
    model: expression.evaluate("\"jvalue\"")
    TableViewColumn {
        role: "nomenclatura"
        title: "Nomenclature"
        width: 200
    }
    TableViewColumn {
        role: "name"
        title: "Name"
        width: 200
    }
    TableViewColumn {
        role: "link"
        title: "Wikipedia"
        width: 200
    }
}

Or, even more powerful, you can create your super individual table using the model and create each row by using a QML Repeater.
Additionally, you can use a lot of fancy stuff like:

mouse interaction
animation
opening an external link
… and so on

The QML code for that looks like this.

import QtQuick 2.0
Rectangle {
    width: 610; height: 500
    Column {
        anchors.horizontalCenter: parent.horizontalCenter
        anchors.verticalCenter: parent.verticalCenter
        Repeater {
            model: expression.evaluate("\"jvalue\"")
            Row {
                id: theRow
                height: mouseArea1.containsMouse || mouseArea2.containsMouse || mouseArea3.containsMouse ? 70 : 50;
                Rectangle { color: "lightblue";
                            border.width: 1
                            width: 150; height: parent.height
                            Text { anchors.centerIn: parent
                                   font.pointSize: 10; text: modelData.nomenclatura }
                            MouseArea {
                                id: mouseArea1
                                anchors.fill: parent
                                hoverEnabled: true
                            }
                }
                Rectangle { color: "lightgreen";
                            border.width: 1
                            width: 150; height: parent.height
                            Text { anchors.centerIn: parent
                                   font.pointSize: 10; text: modelData.name }
                            MouseArea {
                                id: mouseArea2
                                anchors.fill: parent
                                hoverEnabled: true
                            }
                }
                Rectangle {
                            id: linkField
                            color: "lightyellow";
                            border.width: 1
                            width: 300; height: parent.height
                            Text { anchors.centerIn: parent
                                   font.pointSize: 10; text: modelData.link }
                            MouseArea {
                                id: mouseArea3
                                anchors.fill: parent
                                hoverEnabled: true
                                onPressed: linkField.state = "PRESSED"
                                onReleased: linkField.state = "RELEASED"
                                onClicked: Qt.openUrlExternally(modelData.link)
                            }
                            states: [
                                State {
                                    name: "PRESSED"
                                    PropertyChanges { target: linkField; color: "green"}
                                },
                                State {
                                    name: "RELEASED"
                                    PropertyChanges { target: linkField; color: "lightyellow"}
                                }
                            ]
                            transitions: [
                                Transition {
                                    from: "PRESSED"
                                    to: "RELEASED"
                                    ColorAnimation { target: linkField; duration: 1000}
                                },
                                Transition {
                                    from: "RELEASED"
                                    to: "PRESSED"
                                    ColorAnimation { target: linkField; duration: 1000}
                                }
                            ]
                }
            }
        }
    }
}

And that’s it

I hope you liked reading and you will enjoy using it to make beautiful widgets and forms. If you have questions or inputs, feel free to add a comment.
… and in case you still asking where the promised unicorns are. Here’s is a super-fancy implementation

12 novembre 2018
11:30
Fernando Quadro: Tradução do GeoServer para o Português
sur Planet OSGeo

Prezado leitor,

Se você usa o GeoServer, já teve ter percebido que desde a versão 2.8.x existe um problema na tradução para o português brasileiro (pt-BR).

O que parece ter ocorrido, é que quando se foi criar os arquivos properties usaram como base os arquivos do idioma espanhol que ainda não estavam finalizados, por isso o GeoServer hoje tem termos em inglês, português e espanhol.

Eu gostaria então de convidar você, que utiliza o GeoServer e tem uma certa fluência em inglês, para que nos ajude a finalizar essa tradução para a próxima versão do GeoServer, que será a 2.15.x.

Quem se interessar, peço que se cadastre no Site Oficial da tradução:

– Página Oficial de tradução do GeoServer
– Página Oficial de tradução do GeoServer (Equipe pr-BR)

Quem se registrar já pode iniciar iniciar a tradução!
11 novembre 2018
20:11
QGIS Blog: User question of the Month – Nov 18
sur Planet OSGeo

QGIS 2.18 is the third LTR since we started this effort back in 2015 and next year will see the first LTR of QGIS 3. On this occasion, we want to learn more about our users and which versions of QGIS they use. Therefore, we invite you to our QGIS user question of the month.
10 novembre 2018
15:58
Free and Open Source GIS Ramblings: My favorite new recipe in QGIS Map Design 2nd ed
sur Planet OSGeo
If you follow me on Twitter, you have probably already heard that the ebook of “QGIS Map Design 2nd Edition” has now been published and we are expecting the print version to be up for sale later this month. Gretchen Peterson and I – together with our editor Gary Sherman (yes, that Gary Sherman!) – have been working hard to provide you with tons of new and improved map design workflows and many many completely new maps. By Gretchen’s count, this edition contains 23 new maps, so it’s very hard to pick a favorite!

Like the 1st edition, we provide increasingly advanced recipes in three chapters, each focusing on either layer styling, labeling, or creating print layouts. If I had to pick a favorite, I’d have to go with “Mastering Rotated Maps”, one of the advanced recipes in the print layouts chapter. It looks deceptively simple but it combines a variety of great QGIS features and clever ideas to design a map that provides information on multiple levels of detail. Besides the name inspiring rotated map items, this design combines
- map overviews
- map themes
- graduated lines and polygons
- a rotated north arrow
- fancy leader lines
all in one:

“QGIS Map Design 2nd Edition” provides how-to instructions, as well as data and project files for each recipe. So you can jump right into it and work with the provided materials or apply the techniques to your own data.

The ebook is available at LocatePress.
10:45
gvSIG Team: Exprime tu gvSIG mediante Scripting: Disponible el material del taller de las 5as Jornadas gvSIG Uruguay
sur Planet OSGeo

Con motivo de las 5as Jornadas gvSIG Uruguay (Montevideo, 18-19 octubre 2018) se impartió un taller de “Scripting: Exprimiendo / Extendiendo gvSIG ” por parte de Carlos Colombana.

Ahora se ha publicado el material de dicho taller, en el que se muestran las principales funcionalidades del módulo de Scripting, con el fin de obtener el máximo rendimiento y de expandir las capacidades de gvSIG. En él se ven las principales etapas del proceso de desarrollo de add-ons, para el caso particular de un geo-codificador de direcciones geográficas.

Podéis descargar el material para poder seguir el taller desde el siguiente enlace.
08 novembre 2018
23:34
gvSIG Team: Material del Taller de gvSIG aplicado a Geología de las 14as Jornadas gvSIG disponible
sur Planet OSGeo

Ya está disponible la grabación del Taller de gvSIG aplicado a Geología impartido durante las 14as Jornadas Internacionales gvSIG.

En este taller se aborda un análisis integrado que permite la toma de decisiones relacionadas con cuestiones tan dispares como la investigación de recursos minerales, la protección del patrimonio paleontológico, la rentabilidad económica de una explotación minera o su viabilidad ambiental. Todo ello a través de distintas herramientas de edición y geoprocesos de gvSIG desde una perspectiva geológica, realizando una cartografía a partir de datos tomados en campo y la planificación de una investigación geológica. Además se utilizan fuentes de información pública procedentes de diferentes administraciones.

Para realizar este taller podéis descargar la cartografía desde el siguiente enlace.

Si tenéis alguna duda sobre el manejo de gvSIG en la realización de este taller podéis utilizar la lista de usuarios del proyecto.

Y aunque los últimos minutos no estén disponibles aquí tenéis el vídeo con más de dos horas de taller para practicar:

Sauter à la page: 1 2 3 ... 12

← Blogs GeoRezo.net

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Carte des membres Sites amis	Annuaire Agenda Biblio Netvibes Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Carte des CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Canaux

Autres planets