Configureren externe toepassingen¶

Een opmerking voor gebruikers van Windows¶

If you are not an advanced user and you are running QGIS on Windows, you might not be interested in reading the rest of this chapter. Make sure you install QGIS in your system using the standalone installer. That will automatically install SAGA, GRASS and OTB in your system and configure them so they can be run from QGIS. All the algorithms from these providers will be ready to be run without needing any further configuration. If installing through OSGeo4W application, make sure you select for installation SAGA, GRASS and OTB as well.

If you want to know more about how these providers work, or if you want to use some algorithms not included in the simplified toolbox (such as R scripts), keep on reading.

Een opmerking met betrekking tot bestandsindelingen¶

Bij het gebruiken van externe software, betekent het openen van een bestand in QGIS niet dat het kan worden geopend en ook verwerkt kan worden in die andere software. In de meeste gevallen kan andere software lezen wat u hebt geopend in QGIS , maar in sommige gevallen hoeft dat niet zo te zijn. Bij het gebruiken van databases of ongebruikelijke bestandsindelingen, voor raster- of vectorlagen, zouden problemen kunnen optreden. Als dat gebeurd probeer dan goed bekende bestandsindelingen te gebruiken waarvan u weet dat zij door beide programma’s worden begrepen, en controleer de uitvoer in de console (in de dialoogvensters Geschiedenis en Log) om meer te weten te komen over wat er fout gaat.

Bijvoorbeeld het gebruiken van GRASS rasterlagen is één geval waarbij u problemen kunt ondervinden en niet in staat zijn uw werk te voltooien als u een extern algoritme aanroept met een dergelijke laag als invoer. Deze lagen zullen, om deze reden, niet verschijnen als beschikbaar voor algoritmen.

You should, however, find no problems at all with vector layers, since QGIS automatically converts from the original file format to one accepted by the external application before passing the layer to it. This adds extra processing time, which might be significant if the layer has a large size, so do not be surprised if it takes more time to process a layer from a DB connection than it does to process one of a similar size stored in a shapefile.

Providers die geen externe toepassingen gebruiken kunnen elke laag verwerken die u kunt openen in QGIS, omdat zij het voor analyse openen via QGIS.

Met betrekking tot uitvoerindelingen kunnen alle indelingen die worden ondersteund door QGIS als uitvoer worden gebruikt, zowel voor raster- als voor vectorlagen. Sommige providers ondersteunen bepaalde indelingen niet, maar zij kunnen allemaal worden geëxporteerd naar veelvoorkomende indelingen voor rasterlagen die later automatisch kunnen worden getransformeerd door QGIS. Net als in het geval van invoerlagen, als deze conversie nodig is, zou dat de verwerkingstijd kunnen verhogen.

Een opmerking over selecties van vectorlagen¶

Externe toepassingen kunnen ook bewust worden gemaakt van de selecties die bestaan in vectorlagen binnen QGIS. Dat vereist echter het opnieuw schrijven van alle vectorlagen voor de invoer, net als wanneer zij origineel in een indeling waren die niet wordt ondersteund door de externe toepassing. Alleen wanneer er geen selectie bestaat, of de optie Alleen geselecteerde objecten gebruiken niet is ingeschakeld in de algemene configuratie van Processing, kan een laag direct worden doorgegeven aan een externe toepassing.

In andere gevallen is slechts het exporteren van de geselecteerde objecten nodig, wat er voor zorgt dat de benodigde tijd voor uitvoering langer wordt.

SAGA¶

Algoritmen voor SAGA kunnen worden uitgevoerd vanuit QGIS als u SAGA hebt geïnstalleerd op uw systeem en u het framework Processing juist heeft geconfigureerd, zodat het de uit te voeren bestanden van SAGA kan vinden. In het bijzonder is het uit te voeren bestand voor de opdrachtregel in SAGA nodig om algoritmen van SAGA uit te voeren.

Als u werkt op Windows bevatten zowel het zelfstandige installatieprogramma als het installatieprogramma OSGeo4W SAGA tezamen met QGIS, en het pad wordt automatisch geconfigureerd, dus hoeft er verder niets te worden gedaan.

Als u SAGA zelf hebt geïnstalleerd en uw installatieprogramma van QGIS bevatte dat niet, moet het pad naar het uit te voeren bestand van SAGA worden geconfigureerd. Open het dialoogvenster Opties om dat te doen. In het blok SAGA vindt u een instelling genaamd SAGA Folder. Voer het pad in naar de map waar SAGA is geïnstalleerd. Sluit het dialoogvenster Opties en nu bent u gereed om algoritmen van SAGA uit te voeren vanuit QGIS.

Als u werkt op Linux zijn de binaries van SAGA niet opgenomen in Processing, dus moet u de software zelf downloaden en installeren. Bekijk de website van SAGA voor meer informatie.

In dit geval is het niet nodig om het pad naar het uit te voeren bestand van SAGA te configureren, en u zult deze items voor de mappen niet zien. In plaats daarvan dient u er voor te zorgen dat SAGA juist is geïnstalleerd en dat de map ervan is toegevoegd aan de omgevingsvariabele PATH. Open eenvoudigweg een console en typ saga_cmd om te controleren of het systeem kan vinden waar de binaries van SAGA zijn opgeslagen.

R. Creating R scripts¶

Integratie van R in QGIS is anders dan die van SAGA op die manier dat er geen voorgedefinieerde set van algoritmen is die u kunt uitvoeren (uitgezonderd een aantal voorbeelden). In plaats daarvan zou u uw scripts moeten schrijven en opdrachten in R moeten aanroepen, net zoals u zou doen vanuit R, en op een hele soortgelijke manier als die welke we zagen in het gedeelte over scripts voor Processing. Dit gedeelte toont u de te gebruiken syntaxis om deze opdrachten in R te gebruiken vanuit QGIS en hoe objecten van QGIS (lagen, tabellen) er in te gebruiken.

Het eerste dat u moet doen, zoals we zagen in het geval van SAGA, is om QGIS te vertellen waar uw binaries van R zijn opgeslagen. U kunt dit doen door middel van het item map R in het dialoogvenster Opties. Als u die parameter eenmaal hebt ingesteld, kunt u beginnen met het maken en uitvoeren van uw eigen scripts in R.

Nogmaals dit is anders in Linux en u dient er voor te zorgen dat de map R is opgenomen in de omgevingsvariabele PATH. Als u R kunt starten door slechts R in een console te typen, dan bent u klaar om te beginnen.

U dient een scriptbestand te maken dat het framework Processing vertelt hoe die bewerking moet worden uitgevoerd en de corresponderende opdrachten in R om dat te doen om een nieuw algoritme toe te voegen dat een functie in R aanroept (of een meer complex script in R dat u heeft ontwikkeld en dat u beschikbaar zou willen hebben vanuit QGIS).

Scriptbestanden van R hebben de extensie .rsx, en het maken ervan is redelijk eenvoudig als u basiskennis bezit van de syntaxis en scripten van R. Zij zouden moeten worden opgeslagen in de map voor scripts van R. U kunt deze map instellen in de groep met instellingen R (beschikbaar vanuit het dialoogvenster Opties), net zoals u doet met de map voor normale scripts voor Processing.

Laten eens kijken naar een heel eenvoudig scriptbestand, dat de methode in R spsample aanroept om een willekeurig raster te maken binnen de begrenzing van de polygonen in een bepaalde polygoonlaag. Deze methode behoort tot het pakket maptools. Omdat bijna alle algoritmen die u zou willen inbedden in QGIS ruimtelijke gegevens zullen gebruiken of genereren, is kennis van ruimtelijke pakketten zoals maptools en, speciaal, sp, verplicht.

##polyg=vector
##numpoints=number 10
##output=output vector
##sp=group
pts=spsample(polyg,numpoints,type="random")
output=SpatialPointsDataFrame(pts, as.data.frame(pts))

De eerste regels, die worden aangeduid met een Python commentaar symbool (##), geven QGIS informatie over de parameters van het algoritme en de resultaten die het zal genereren. Zij werken met exact exact dezelfde syntaxis als de scripts voor Processing die we al eerder hebben gezien en zal dus hier niet nogmaals worden beschreven.

Bekijk eens de hoofdstukken R Intro en R Syntax in de Trainingshandleiding om meer informatie te krijgen over hoe u uw eigen scripts voor R schrijft-

Wanneer u een parameter voor de invoer declareert, gebruikt QGIS die informatie voor twee dingen: het maken van de interface om de gebruiker te vragen naar de waarde voor die parameter en het maken van een overeenkomstige variabele in R die later kan worden gebruikt als invoer voor opdrachten in R.

In het bovenstaande voorbeeld declareren we een invoer van het type vector genaamd polyg. Bij het uitvoeren van het algoritme zal QGIS in R de laag openen die is geselecteerd door de gebruiker en die opslaan in een variabele die ook is genaamd polyg. Dus de naam van een parameter is ook de naam van de variabele die we in R kunnen gebruiken voor de toegang tot de waarde van die parameter (dus zou u moeten vermijden om door R gereserveerde woorden als namen voor parameters te gebruiken).

Ruimtelijke elementen zoals vector- en rasterlagen worden gelezen met behulp van de opdrachten readOGR() en brick() (u hoeft zich geen zorgen te maken over het toevoegen van deze opdrachten aan uw bestand voor de beschrijving – QGIS zal dit voor u doen), en zij worden opgeslagen als objecten Spatial*DataFrame. Tabelvelden worden opgeslagen als tekenreeksen die de naam van het geselecteerde veld bevatten.

Tabellen worden geopend met behulp van de opdracht read.csv(). Als een door de gebruiker ingevoerde tabel niet in de indeling CSV is, zal die worden geconverteerd, voorafgaande aan het importeren in R.

Aanvullend kunnen rasterbestanden worden gelezen met behulp van de opdracht readGDAL() in plaats van met brick() door de ##usereadgdal te gebruiken.

Als u een gevorderde gebruiker bent en niet wilt dat QGIS het object maakt dat de laag vertegenwoordigt, kunt u de tag ##passfilename gebruiken om aan te geven dat u in plaats daarvan een tekenreeks met de bestandsnaam prefereert. In dat geval is het aan u om het bestand te openen vóórdat een bewerking wordt uitgevoerd op de gegevens die het bevat.

Met bovenstaande informatie kunnen we nu de eerste regel van ons eerste voorbeeldscript begrijpen (de eerste regel die niet begint met een opmerking in Python).

pts=spsample(polyg,numpoints,type="random")

De variabele polygon bevat al een object SpatialPolygonsDataFrame, dus kan het worden gebruikt om de methode spsample aan te roepen, net als numpoints, die het aantal punten aangeeft die moeten worden toegevoegd aan het gemaakte voorbeeldraster.

Omdat we al een uitvoer hebben gedeclareerd van het type vector genaamd out, moeten we een variabele genaamd out maken en er een object Spatial*DataFrame in opslaan (in dit geval een SpatialPointsDataFrame). U kunt elke naam gebruiken voor uw tussentijdse variabelen. Zorg er alleen voor dat de variabele die uw uiteindelijke resultaat opslaat dezelfde naam heeft als die welke u gebruikt om het te declareren, en dat het een geschikte waarde bevat.

In dit geval zal het resultaat dat wordt verkregen uit de methode spsample expliciet moeten worden geconverteerd naar een object SpatialPointsDataFrame, omdat het zelf een object van de klasse ppp is, wat geen geschikte klasse is om te worden teruggegeven aan QGIS.

Als uw algoritme rasterlagen genereert, is de manier waarop zij worden opgeslagen afhankelijk van het feit of u al dan niet de optie ##dontuserasterpackage heeft gebruikt. Wanneer u die heeft gebruikt worden lagen opgeslagen met behulp van de methode writeGDAL(). Indien niet, zal de methode writeRaster() uit het pakket raster worden gebruikt.

Als u de optie ##passfilenames gebruikte, wordt de uitvoer gegenereerd met behulp van het pakket raster (met writeRaster()), zelfs als het niet is gebruikt voor de invoer.

Als uw algoritme geen laag genereert, maar in plaats daarvan een tekstresultaat in de console, dient u aan te geven dat u wilt dat de console wordt weergegeven als de uitvoering eenmaal is voltooid. Start eenvoudigweg de opdrachtregels die de resultaten produceren die u wilt afdrukken met het teken > (‘groter dan’) om dat te doen. De uitvoer van alle andere regels zal niet worden weergegeven. Hier is bijvoorbeeld het bestand voor de beschrijving van een algoritme dat een test voor normalen uitvoert op een bepaald veld (kolom) van de attributen van een vectorlaag:

##layer=vector
##field=field layer
##nortest=group
library(nortest)
>lillie.test(layer[[field]])

De uitvoer van de laatste regel wordt afgedrukt, maar de uitvoer van de eerste wordt dat niet (en ook de uitvoer van de andere opdrachtregels, die automatisch door QGIS werden toegevoegd, worden dat niet).

Als uw algoritme iets grafisch maakt (met behulp van de methode plot()), voeg dan de volgende regel toe:

##showplots

Dit zal er voor zorgen dat QGIS alle grafische uitvoer voor R zal omleiden naar een tijdelijk bestand, wat zal worden geopend als de uitvoering van R is voltooid.

Beide grafische en console-resultaten zullen worden weergegeven in het beheer van de resultaten van Processing.

Bekijk, voor meer informatie, de scriptbestanden die zijn opgenomen in Processing. De meeste daarvan zijn redelijk eenvoudig en zullen u enorm helpen te begrijpen hoe u uw eigen scripts kunt maken.

GRASS¶

Configuring GRASS is not much different from configuring SAGA. First, the path to the GRASS folder has to be defined, but only if you are running Windows. Additionally, a shell interpreter (usually msys.exe, which can be found in most GRASS for Windows distributions) has to be defined and its path set up as well.

By default, the processing framework tries to configure its GRASS connector to use the GRASS distribution that ships along with QGIS. This should work without problems in most systems, but if you experience problems, you might have to configure the GRASS connector manually. Also, if you want to use a different GRASS installation, you can change that setting and point to the folder where the other version is installed. GRASS 6.4 is needed for algorithms to work correctly.

Als u werkt op Linux hoeft u er slechts voor te zorgen dat GRASS correct is geïnstalleerd, en dat het zonder problemen kan worden uitgevoerd vanaf een console.

Algoritmen van GRASS gebruiken een regio voor berekeningen. Deze regio kan handmatig worden gedefinieerd met behulp van waarden die soortgelijk zijn aan die welke werden gebruikt in de configuratie van SAGA, of automatisch, met het minimum bereik dat alle gebruikte invoerlagen bedekt bij het elke keer uitvoeren van het algoritme. Als de laatste benadering het gedrag is dat u prefereert, selecteer dan de optie Use min covering region in de configuratie van de parameters in GRASS.

GDAL¶

No additional configuration is needed to run GDAL algorithms. Since they are already incorporated into QGIS, the algorithms can infer their configuration from it.

Orfeo Toolbox¶

Orfeo Toolbox (OTB) algorithms can be run from QGIS if you have OTB installed in your system and you have configured QGIS properly, so it can find all necessary files (command-line tools and libraries).

As in the case of SAGA, OTB binaries are included in the stand-alone installer for Windows, but they are not included if you are running Linux, so you have to download and install the software yourself. Please check the OTB website for more information.

Once OTB is installed, start QGIS, open the processing configuration dialog and configure the OTB algorithm provider. In the Orfeo Toolbox (image analysis) block, you will find all settings related to OTB. First, ensure that algorithms are enabled.

Then, configure the path to the folder where OTB command-line tools and libraries are installed:

• Usually OTB applications folder points to /usr/lib/otb/applications and OTB command line tools folder is /usr/bin.
• If you use any of the installers that include OTB, such as OSGeo4W, there is no need for further configuration. Processing will detect the path automatically and will not show the corresponding configuration entries. Otherwise, fill the OTB applications folder and OTB command line tools folder parameters with the to the corresponding values for your installation.

TauDEM¶

TauDEM (Terrain Analysis Using Digital Elevation Models) is a tools for the extraction and analysis of hydrological information from Digital Elevation Models (DEM). TauDEM can be used from QGIS if you have it installed in your system and configured QGIS properly, so it can find all necessary files.

There are two versions of TauDEM tools: singlefile (TauDEM 5.0.6 or 5.1.2) and multifile (TauDEM 5.2.0). The difference between these versions in the supported inputs/outputs. Single files version accepts only single raster file and write single file as output. Multifile version accepts a directory with rasters and writes directory with rasters as output. Such directory should contain rasters that will be treated as a single DEM grid.

TauDEM Processing provider supports both single- and multifile versions of TauDEM and even allows to use them simultaneously.

mkdir build
cd build

CXX=mpicxx cmake -DCMAKE_INSTALL_PREFIX=/usr/local ..
make

sudo make install

#include "mpi.h"

#include <stdint.h>

#include "mpi.h"
#include <stdint.h>

#include <stdint.h>