Monsters nemen van rastergegevens met behulp van punten of polygonen

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 twee plug-ins - Point Sampling Tool en de plug-in Gebiedssstatistieken.

Overzicht van de taak

Gegeven een raster van maximum temperaturen in de VS, moeten we de temperatuur voor alle stedelijke gebieden uitnemen en ook de gemiddelde temperatuur berekenen voor elke county in de VS.

Andere vaardigheden die u zult leren

  • Een vectorlaag opnieuw projecteren.

  • 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 2013 US Gazetteer dat de stedelijke gebieden in de VS weergeeft. Download het bestand Urban Areas Gazetteer File.

Omdat we temperaturen willen samenvoegen over counties, zullen we 2013 TIGER/Line Shapefiles gebruiken. Download het shapefile Counties (and equivalents).

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

us.tmax_nohads_ll_20140525_float.tif

2013_Gaz_ua_national.zip

tl_2013_us_county.zip

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

Procedure

  1. Ga naar Kaartlagen ‣ Laag toevoegen ‣ Rasterlaag toevoegen en blader naar het gedownloade bestand us.tmax_nohads_ll_{YYYYMMDD}_float.tif en klik op Openen.

../_images/1131.png
  1. Selecteer, als de laag eenmaal is geladen, het gereedschap Objecten identificeren en klik ergens op de laag. U zult de waarde van de temperatuur zien in Celsius als de waarde van Band 1 op die locatie.

../_images/279.png
  1. Pak nu het gedownloade bestand 2013_Gaz_ua_national.zip uit en neem daaruit het bestand 2013_Gaz_ua_national.txt op uw schijf. Ga naar Kaartlagen ‣ Laag toevoegen ‣ Tekstgescheiden kaartlaag toevoegen.

../_images/340.png
  1. klik, in het dialoogvenster maak een kaartlaag uit een tekstgescheiden bestand op Bladeren en open 2013_Gaz_ua_national.txt. Kies Tab onder Zelfgekozen tekstscheiders. De puntcoördinaten staan in Latitude en Longitude, dus selecteer INTPTLONG als X-veld en INTPTLAT als Y-veld. Selecteer het vak Ruimtelijke index gebruiken en klik op OK.

../_images/430.png
  1. Nu zijn we klaar om de waarden voorde temperatuur uit te nemen vanuit de rasterlaag. Installeer de plug-in Point Sampling Tool. Bekijk Plug-ins gebruiken voor details over hoe plug-ins te installeren.

../_images/531.png
  1. Open het dialoogvenster voor de plug-in via Plug-ins beheren en installeren ‣ Analyses ‣ Point sampling tool.

../_images/629.png
  1. Selecteer, in het dialoogvenster van Point Sampling Tool, 2013_Gaz_ua_national als de Layer containing sampling points. We moeten expliciet de velden kiezen uit de invoerlaag die we in de uitvoerlaag willen hebben. Kies de velden GEOID en NAME uit de laag 2013_Gaz_ua_national layer. We kunnen in één keer monsterwaarden uit meerdere rasterbanden halen, maar omdat ons raster slechts 1 band heeft, kies de us.tmax_nohads_ll_{YYYYMMDD}_float: Band 1. Noem de uitvoer-vectorlaag max_temparature_at_urban_locations.shp. Klik op OK om het proces van monsters nemen te beginnen. Klik op Close als het proces is voltooid.

../_images/728.png
  1. U zult een nieuwe laag zien max_temparature_at_urban_locations, die is geladen in QGIS. gebruik het gereedschap Objecten identificeren om op een willekeurig punt te klikken om de attributen te zien. U zult het veld us.tmax_no zien - dat de pixelwaarde voor het raster op de locatie van het punt bevat.

../_images/827.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 max_temparature_at_urban_locations en 2013_Gaz_ua_national layers. Klik met rechts en selecteer Verwijderen om ze uit de inhoudsopgave van QGIS te verwijderen.

../_images/925.png
  1. Ga naar Kaartl;agen ‣ Laag toevoegen ‣ Vectorlaag toevoegen. Blader naar het gedownloade bestand tl_2013_us_county.zip en klik op Open. Selecteer tl_2013_us_county.shp als de laag en klik op OK.

../_images/1026.png
  1. Het bestand tl_2013_us_county zal worden toegevoegd aan QGIS. Deze laag staat in de projectie EPSG:4269 NAD83. Dat komt niet overeen met de projectie van de rasterlaag. We zullen deze laag opnieuw projecteren in de projectie EPSG:4326 WGS84.

../_images/1132.png
  1. Klik met rechts op de laag tl_2013_us_county en selecteer Opslaan als...

../_images/1227.png
  1. Klik, in het dialoogvenster Vectorlaag opslaan as.., op Bladeren en noem het uitvoerbestand counties.shp. Kies Geselecteerd CRS uit het keuzemenu CRS. Klik op Bladeren en selecteer WGS 84 als het CRS. Selecteer Voeg opgeslagen bestand toe aan kaart en klik op OK.

../_images/1324.png
  1. Een nieuw laag, genaamd counties, zal worden toegevoegd aan QGIS.

../_images/1423.png
  1. Schakel de plug-in Gebiedsstatistieken in. Dit is een bron plug-in dus die is al geïnstalleerd. Bekijk Plug-ins gebruiken om te weten te komen hoe bron plug-ins worden ingeschakeld.

../_images/1522.png
  1. Ga naar Raster ‣ Gebiedsstatistieken ‣ Gebiedsstatistieken.

../_images/1621.png
  1. Selecteer us.tmax_nohads_ll_{YYYYMMDD}_float als de Rasterlaag en counties als de Polygoonlaag met zones. Voer ZS_ in als Voorvoegsel uitvoerkolom. Klik op OK.

../_images/1720.png
  1. De analyse kan enige tijd vergen, afhankelijk van de grootte van de gegevensset.

../_images/1820.png
  1. Selecteer de laag counties als het proces is voltooid. Gebruik het gereedschap Objecten identificeren en klik op een willekeurig polygoon van een county. U zult zien dat drie nieuwe attributen zijn toegevoegd aan de laag: ZS_count, ZS_mean en ZS_sum. Deze attributen bevatten respectievelijk het totaal aantal rasterpixels, de gemiddelde waarde van de rasterpixels, en de som van de waarden van de rasterpixels. Omdat we zijn geïnteresseerd in de gemiddelde temperatuur, is het veld ZS_mean het enige om te gebruiken.

../_images/1918.png
  1. Laten we deze laag opmaken een een kaart voor de temperatuur te maken. Klik met rechts op de laag counties en selecteer Eigenschappen.

../_images/2015.png
  1. Schakel naar de tab Stijl. Kies de stijl Gradueel en selecteer ZS_mean als de Kolom. Kies een Kleurverloop en Modus van uw keuze. Klik op Classificeren om de klassen te maken. Klik op OK. (Zie Basis opmaken van vectoren voor meer details over opmaken.)

../_images/2119.png
  1. U zult zien dat de polygonen van de counties zijn opgemaakt met behulp van de gemiddelde temperatuur die uit het raster is uitgenomen.

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

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