Monsters nemen van rastergegevens met behulp van punten of polygonen (QGIS3)

Vele wetenschappelijke en milieu-gegevenssets zijn gerasterde rasterafbeeldingen. Hoogtegegevens (DEM) worden ook gedistribueerd als rasterbestanden. In deze rasterbestanden wordt de parameter die wordt weergegeven gecodeerd als de pixelwaarden van het raster. Vaak moet men de pixelwaarden van bepaalde locaties verzamelen of ze samenvoegen over enkele gebieden. Deze functionaliteit is beschikbaar in QGIS via algoritmes van Processing - Monster rasterwaarden voor puntlagen en Gebiedsstatistieken voor polygoonlagen.

Overzicht van de taak

Gegeven een raster van dagelijkse maximum temperaturen op het vasteland van de VS, moeten we de temperatuur uit een puntenlaag voor alle stedelijke gebieden uitnemen en de gemiddelde temperatuur voor een polygoonlaag berekenen voor elke county in de VS.

Andere vaardigheden die u zult leren

  • Meerdere lagen uit de inhoudsopgave van QGIS selecteren en verwijderen.

De gegevens ophalen

NOAA’s Climate Prediction Center verschaft GIS-gegevens gerelateerd aan temperatuur en neerslag in de VS. Download het laatste rasterbestand voor de maximum temperaturen. Het bestand zal zijn genaamd us.tmax_nohads_ll_{YYYYMMDD}_float.tif

We zullen een CSV-bestand gebruiken vanaf 2018 US Gazetteer dat de stedelijke gebieden in de VS weergeeft. Download het bestand Urban Areas Gazetteer File.

US Census Bureau verschaft TIGER/Line Shapefiles. U kunt de FTP site bezoeken en het census tracts shapefile voor Californië downloaden. Download het bestand `https://www2.census.gov/geo/tiger/TIGER2018/COUNTY/tl_2018_us_county.zip>`_.

Voor het gemak kunt u direct een kopie van de gegevensset downloaden vanaf de links hieronder:

us.tmax_nohads_ll_20190501_float.tif

2018_Gaz_ua_national.zip

tl_2018_us_county.zip

Gegevensbronnen: [NOAACPC], [USGAZETTEER] [TIGER]

Procedure

  1. Pak zowel 2018_Gaz_ua_national.zip als tl_2018_us_county.zip uit in een map op uw computer. Open QGIS en zoek naar het bestand us.tmax_nohads_ll_20190501_float.tif in de QGIS Browser en sleep het naar het kaartvenster.

../../_images/1108.png
  1. U zult een nieuwe rasterlaag us.tmax_nohads_ll_20190501_float zien geladen in het paneel Lagen. Deze rasterlaag bevat de maximum temperaturen die zijn opgenomen op elke pixel in graden Celsius. Vervolgens zullen we het bestand met punten van de stedelijke gebieden laden. Dit bestand komt als een tekstbestand in de indeling Tab Separated Values (TSV). Klik op de knop Databronnen beheren openen op de werkbalk Databronnen beheren.

../../_images/268.png
  1. Schakel naar de tab Tekengescheiden tekst. Klik op de knop naast Bestandsnaam en specificeer het pad naar het tekstbestand dat u hebt gedownload. Selecteer, in het gedeelte Bestandsindeling, Zelfgekozen tekstscheiders en selecteer Tab. Selecteer INTPTLONG als het X-veld en INTPTLAT als het Y-veld. Klik op Toevoegen en dan op Close.

../../_images/340.png
  1. Een nieuwe puntenlaag 2018_Gaz_ua_national zal worden geladen in het paneel Lagen. Nu zijn we klaar om de waarden uit de rasterlaag uit te nemen voor deze punten. Ga naar Processing ‣ Toolbox.

../../_images/425.png
  1. Zoek en lokaliseer het algoritme Rasteranalyse ‣ Monster rasterwaarden. Dubbelklik om het te starten.

../../_images/525.png
  1. Selecteer 2018_Gaz_ua_national als de Invoer puntenlaag. Selecteer us.tmax_nohads_ll_20190501_float als de Rasterlaag waarvan monsters moeten worden genomen. Vergroot de Gevorderde parameters en voer tmax in als Voorvoegsel uitvoerkolom. Klik op Uitvoeren. Klik, als de verwerking is voltooid, op Close.

../../_images/624.png
  1. Een nieuwe laag Monsterpunten zal worden geladen in het paneel Lagen. Selecteer het gereedschap Objecten identificeren op de werkbalk Attributen en klik op een willekeurig punt. U zult de attributen zien worden weergegeven in het paneel Identificatieresultaten. U zult een nieuw attribuut zien toegevoegd aan elk object, genaamd tmax_1. Dat is de pixelwaarde van de rasterlaag die is uitgenomen op de locatie van het punt. De 1 geeft het bandnummer weer van het raster. Als de rasterlaag meerder banden zou hebben gehad, zou u meerdere nieuwe kolommen in de uitvoerlaag hebben gezien.

../../_images/724.png
  1. Het eerste gedeelte van onze analyse is voltooid. Laten we de lagen die niet meer nodig zijn verwijderen. Houdt de Shift-toets ingedrukt en selecteer de lagen Monsterpunten en 2018_Gaz_ua_national. Klik met rechts en selecteer Verwijderen om ze uit QGIS te verwijderen. Selecteer OK indien gevraagd wordt 2 legenda-items verwijderen?.

../../_images/823.png
  1. Nu zullen we de laag met county’s gebruiken om monsters te nemen uit het raster en de gemiddelde temperatuur berekenen voor elke county. Zoek naar het bestand tl_2018_us_county.shp in de QGIS Browser en sleep het naar het kaartvenster.

../../_images/924.png

Notitie

De meeste algoritmes van Processing zullen de invoerlaag lezen en een nieuwe laag maken. Maar het algoritme Gebiedsstatistieken is anders. Het past de invoerlaag aan en voegt daar nieuwe attributen aan toe. Dat is waarom het belangrijk is om eerst de invoerbestanden uit te pakken. QGIS kan een laag direct vanuit een ZIP-archief laden, maar het kan geen gezipte laag aanpassen. Het algoritme van Processing zal falen en het kan de invoerlaag niet aanpassen.

  1. Een nieuwe laag tl_2018_us_county zal worden geladen in het paneel Lagen. Ga naar Processing ‣ Toolbox.

../../_images/1022.png
  1. Zoek en lokaliseer het algoritme Rasteranalyse ‣ Gebiedsstatistieken en dubbelklik om het te starten.

../../_images/1126.png
  1. Selecteer us.tmax_nohads_ll_20190501_float als de Rasterlaag en tl_2018_us_county als de Vectorlaag bevat zones. Voer tmax_ in als Voorvoegsel uitvoerkolom. Klik op de knop naast Statistieken om te berekenen.

../../_images/1224.png
  1. Selecteer alleen de waarde Gemiddelde en klik op OK.

../../_images/1324.png
  1. Klik op Uitvoeren om de verwerking te beginnen. Het algoritme kan er een paar minuten over doen om te voltooien. Klik op Close.

../../_images/1422.png
  1. Zoals eerder al vermeld maakt het algoritme GebiedssStatistieken geen nieuwe laag, maar past de laag met zones aan. Klik met rechts op de laag tl_2018_us_county en selecteer Attributentabel openen.

../../_images/1522.png
  1. U zult zien dat een nieuwe kolom, genaamd tmax_mean, is toegevoegd aan de attributentabel. Deze bevat de gemiddelde waarde voor de temperatuur die is uitgenomen uit de polygoon voor elk object. Er zijn enkele nulwaarden omdat die county’s (behorende tot Alaska, Hawaii en Puerto Rico) buiten het bereik van de rasterlaag liggen.

../../_images/1620.png
comments powered by Disqus

This work is licensed under a Creative Commons Attribution 4.0 International License