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/1141.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/2109.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/360.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/430.png
  1. Zoek en lokaliseer het algoritme Rasteranalyse ‣ Monster rasterwaarden. Dubbelklik om het te starten.

../../_images/528.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/627.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/727.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/826.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/927.png
  1. Een nieuwe laag tl_2018_us_county zal worden geladen in het paneel Lagen. Ga naar Processing ‣ Toolbox.

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

../../_images/1142.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/1236.png
  1. Selecteer alleen de waarde Gemiddelde en klik op OK.

../../_images/1334.png
  1. Now, click on the next to Zonal Statistics and select Save to File option to save the output.

    ../../_images/1431.png
  2. Browse to the folder where you want to save the output. Save the output layer as us_county_tmax.shp. Click Run to start the processing. The algorithm may take a few minutes to complete and you will see the output layer added as soon as the processing finishes. Click Close.

../../_images/1529.png
  1. Right-click the us_county_tmax layer, and select Open Attribute Table.

../../_images/1627.png
  1. You will see a new column called tmax_mean added in the attribute table. This contains the average temperature value extracted over the polygon for each feature. There are some null values because those counties (belonging to Alaska, Hawaii and Puerto Rico) are outside of the raster layer’s extent.

../../_images/1728.png

If you want to give feedback or share your experience with this tutorial, please comment below. (requires GitHub account)