6.4. Lesson: Ruimtelijke statistieken

Notitie

Les ontwikkeld door Linfiniti en S Motala (Cape Peninsula University of Technology)

Ruimtelijke statistieken stellen u in staat te analyseren en te begrijpen wat er gaande is in een bepaalde vector gegevensset. QGIS bevat vele nuttige gereedschappen voor statistische analyses.

Het doel voor deze les: Weten hoe de ruimtelijke statistische gereedschappen in de Toolbox van Processing te gebruiken.

6.4.1. basic Follow Along: Een test-gegevensset maken

We zullen een willekeurige set punten maken, om een gegevensset te krijgen waar we mee kunnen werken.

U heeft, om dat te doen, een polygoon gegevensset nodig die het bereik van het gebied, waarin u de punten wilt maken, definieert.

We zullen het gebied gebruiken dat wordt bedekt door straten.

  1. Begin een nieuw project.

  2. Voeg uw laag roads toe, als ook srtm_41_19 (hoogtegegevens) te vinden in exercise_data/raster/SRTM/.

    Notitie

    U zou kunnen merken dat de laag SRTM DEM een ander CRS heeft dan die van de laag roads. QGIS projecteert beide lagen in één enkel CRS. Voor de volgende oefeningen is dit verschil niet van belang, maar het staat u vrij een laag opnieuw te projecteren (zoals eerder weergegeven in deze module).

  3. Open de Toolbox van Processing

  4. Gebruik het gereedschap Vector geometrie ► Minimum begrenzing geometrie om een gebied te maken dat alle wegen omsluit door Convex Hull te kiezen als Type geometrie:

    ../../../_images/roads_hull_setup.png

    Zoals u weet maakt Processing , als u geen uitvoer specificeert, tijdelijke lagen. Het is aan u om de lagen onmiddellijk op te slaan of op een later moment.

Willekeurige punten genereren

  • Maak in dit gebied 100 willekeurige punten met het gereedschap Vector maken ► Willekeurige punten binnen grenzen van laag, met een minimum afstand van 0.0:

    ../../../_images/random_points_setup.png

    Notitie

    Het gele waarschuwingsteken vertelt u dat die parameter afstanden betreft. De laag Begrenzing geometrie staat in een Geografisch coördinatensysteem en het algoritme herinnert u daar slechts aan. Voor dit voorbeeld zullen we deze parameter niet gebruiken, dus kunt u het negeren.

Verplaats, indien nodig, de gemaakte willekeurige punten naar bovenaan in de legenda om hem beter te kunnen zien:

../../../_images/random_points_result.png

Een monster uit de gegevens

U dient het algoritme Rasteranalye ► Monster rasterwaarden te gebruiken om een gegevensset met monsters uit het raster te nemen. Dit gereedschap neemt monsters uit het raster op de locaties van de punten en voegt de rasterwaarden toe in nieuwe veld(en), afhankelijk van het aantal banden in het raster.

  1. Open het dialoogvenster voor het algoritme Monsters rasterwaarden

  2. Selecteer Willekeurige punten als de laag die de monsterpunten bevat en het SRTM-raster als de band waaruit de waarden moeten komen. De standaardnaam voor het nieuwe veld is rvalue_N, waar N het nummer is van de rasterband. U kunt de naam van het voorvoegsel wijzigen als u dat wilt.

    ../../../_images/sample_raster_dialog.png
  3. Klik op Uitvoeren

Nu kunt u de gegevens voor het monster uit het rasterbestand controleren in de attributentabel van de laag Monsterpunten.Zij zullen in een nieuw veld staan, met de naam die u hebt gekozen.

Een mogelijke laag voor het monster wordt hier weergegeven:

../../../_images/random_samples_result.png

De monsterpunten zijn op hun veld rvalue_1 geclassificeerd, zodat rode punten op een hogere hoogte liggen.

U zult deze monsterlaag gaan gebruiken voor de rest van de statistische oefeningen.

6.4.2. basic Follow Along: Basisstatistieken

Nu nog de basisstatistieken voor deze laag ophalen.

  1. Klik op het pictogram sum Statistisch overzicht weergeven op de werkbalk Attributen. Een nieuw paneel zal openen.

  2. Specificeer, in het dialoogvenster dat verschijnt, de laag Monsterpunten als de bron.

  3. Selecteer het veld rvalue_1 in het combinatievak voor velden. Dit is het veld waarvoor u de statistieken wilt berekenen.

  4. Het paneel Statistieken zal automatisch worden bijgewerkt met de berekende statistieken:

    ../../../_images/basic_statistics_results.png

    Notitie

    U kunt de waarden kopiëren door te klikken op de knop editCopy Statistieken naar klembord kopiëren en de resultaten in een werkblad plakken.

  5. Sluit het paneel Statistieken indien gereed

Veel verschillende statistieken zijn beschikbaar:

Telling

Het aantal monsters/waarden.

Som

De waarden bij elkaar opgeteld.

Gemiddelde

De gemiddelde waarde is eenvoudigweg de som van de waarden, gedeeld door het aantal waarden.

Mediaan

Als u alle waarden schikt van de kleinste tot de grootste, is de middelste waarde (of het gemiddelde van de twee middelste waarden als N een even getal is) de mediaan van de waarden.

St afw (pop)

De standaardafwijking. Geeft een indicatie over hoe dicht de waarden zijn geclusterd rondom het gemiddelde. Hoe kleiner de standaardafwijking, hoe meer waarden neigen naar het gemiddelde.

Minimum

De laagste waarde.

Maximum

De hoogste waarde.

Bereik

Het verschil tussen de laagste en de hoogste waarden.

Kw1

Eerste kwartiel van de gegevens.

Kw3

Derde kwartiel van de gegevens.

Ontbrekende (null) waarden

Het aantal ontbrekende waarden.

6.4.3. basic Follow Along: Statistieken berekenen van afstanden tussen punten

  1. Maak een nieuwe tijdelijke puntenlaag.

  2. Ga naar modus Bewerken en digitaliseer drie punten ergens tussen de andere punten.

    Gebruik, als alternatief, dezelfde methode voor het maken van willekeurige punten als hiervoor, maar specificeer slechts drie punten.

  3. Sla uw nieuwe laag op als distance_points in de indeling die u wilt.

Statistieken maken over de afstand tussen punten in de twee lagen:

  1. Open het gereedschap Vector-analyse ► Afstandsmatrix.

  2. Selecteer de laag distance_points als de invoerlaag en de laag Monsterpunten als de doellaag.

  3. Stel het als volgt in:

    ../../../_images/distance_matrix_setup.png
  4. Als u wilt, kunt u de uitvoerlaag opslaan als een bestand of alleen het algoritme uitvoeren en de tijdelijke uitvoerlaag later opslaan.

  5. Klik op Uitvoeren om de laag voor de afstandsmatrix te maken.

  6. Open de attributentabel van de gemaakte laag: waarden verwijzen naar de afstanden tussen de objecten van Afstandsmatrix en hun twee dichtstbijzijnde punten op de laag Monsterpunten:

    ../../../_images/distance_matrix_example.png

Met deze parameters berekent het gereedschap Afstandsmatrix statistieken over de afstanden voor elk punt van de invoerlaag, ten opzichte van de dichtstbijzijnde punten van de doellaag. De velden van de uitvoerlaag bevatten gemiddelde, standaardafwijking, minimum en maximum voor de afstanden tot de dichtstbij gelegen buren van de punten op de invoerlaag.

6.4.4. basic Follow Along: ‘Dichtstbijzijnde buur’-analyse (in laag)

Een ‘Dichtstbijzijnde buur’-analyse uitvoeren op een puntenlaag:

  1. Kies Vector-analyse ► ‘Dichtstbijzijnde buur’-analyse.

  2. Selecteer, in het dialoogvenster dat verschijnt, de laag Random samples en klik op Uitvoeren.

  3. De resultaten zullen verschijnen in het paneel van Processing Resultaten bekijken.

    ../../../_images/result_viewer.png
  4. Klik op de blauwe link om de HTML- pagina met de resultaten te bekijken:

    ../../../_images/nearest_neighbour_example.png

6.4.5. basic Follow Along: Gemiddelde coördinaten

De gemiddelde coördinaten van een gegevensset verkrijgen:

  1. Start Vector-analyse ► Gemiddelde coördin(a)t(en)

  2. Specificeer, in het dialoogvenster dat verschijnt, Willekeurige punten als Invoerlaag, en laat de optionele keuzen ongewijzigd.

  3. Klik op Uitvoeren.

Laten we dit eens vergelijken met de centrale coördinaat van de polygoon die werd gebruikt om het willekeurige monster te maken.

  1. Start Vector geometrie ► Zwaartepunten

  2. Selecteer, in het dialoogvenster dat verschijnt, Begrenzing geometrie als de invoerlaag.

Zoals u kunt zien vallen de gemiddelde coördinaten (roze punt) en het middelpunt van het gebruikte gebied (in groen) niet noodzakelijkerwijze samen.

Het zwaartepunt is het barycenter van de laag (het barycenter van een vierkant is het midden van het vierkant), terwijl de gemiddelde coördinaten het gemiddelde weergeven van alle coördinaten van de knopen.

../../../_images/polygon_centroid_mean.png

6.4.6. basic Follow Along: Histogrammen van afbeeldingen

Het histogram van een gegevensset geeft de verdeling van de waarden ervan weer. De eenvoudigste manier om dit in QGIS te demonstreren is via het histogram van een afbeelding, beschikbaar in het dialoogvenster Laageigenschappen van elke afbeeldingslaag (raster gegevensset).

  1. Klik, in uw paneel Lagen, met rechts op de laag srtm_41_19.

  2. Selecteer Eigenschappen

  3. Kies de tab Histogram. U moet misschien op de knop Histogram herberekenen klikken om de grafiek te maken. U zult een grafiek zien die verdeling van de frequentie van de rasterwaarden weergeeft.

    ../../../_images/histogram_export.png
  4. De grafiek kan met de knop fileSave Plot opslaan worden geëxporteerd als een afbeelding

  5. U kunt meer gedetailleerde informatie over de laag vinden op de tab Informatie (het gemiddelde en maximum waarde zijn geschat, en zouden niet exact kunnen zijn).

De gemiddelde waarde is 332.8 (geschat op 324.3), en de maximum waarde is 1699 (geschat op 1548)! U kunt zoomen in het histogram. Het histogram ziet er verticaal gecomprimeerd uit, omdat er heel veel pixels met de waarde 0 zijn. Door in te zoomen om alles, behalve de piek op 0, te bedekken zult u meer details zien:

../../../_images/histogram_export_zoom.png

Notitie

Als de waarden gemiddelde en maximum niet hetzelfde zijn als hierboven, kan dat komen door de berekening van de waarden min/max. Open de tab Symbologie en vergroot het menu Instellingen min / max waarden. Kies radioButtonOn Min / max en klik op Apply.

Onthoud dat een histogram u de verdeling van de waarden laat zien en dat niet alle waarden noodzakelijkerwijze ook zichtbaar zijn in de grafiek.

6.4.7. basic Follow Along: Ruimtelijke interpolatie

Laten we zeggen dat u een collectie monsterpunten heeft van waaruit u gegevens wilt extrapoleren. U heeft misschien toegang tot de gegevensset Monsterpunten die we eerder gemaakt hebben en wil een idee krijgen over hoe het terrein eruit ziet.

  1. Start, om te beginnen, het gereedschap GDAL ► Rasteranalyse ► Raster (IDW met zoeken naar Nearest neighbor) in de Toolbox van Processing.

  2. Selecteer ``Monsterpunten``als Puntenlaag

  3. Stel Gewogen macht in op 5.0

  4. Stel, in Gevorderde parameters, rvalue_1 in voor Z-waarde uit veld

  5. Klik tenslotte op Uitvoeren en wacht tot de verwerking is voltooid.

  6. Sluit het dialoogvenster

Hier is een vergelijking van de originele gegevensset (links) met die welke we hebben gemaakt uit onze monsterpunten (rechts). Die van u kan er anders uitzien vanwege de willekeurige herkomst van de locatie van de monsterpunten.

../../../_images/interpolation_comparison.png

Zoals u kunt zien zijn 100 monsterpunten niet echt genoeg om een gedetailleerde impressie van het terrein te maken. Het geeft een zeer algemeen idee, maar het kan ook misleidend zijn.

6.4.8. moderate Try Yourself Verschillende methoden voor interpolatie

  1. Gebruik de hierboven weergegeven processen om een nieuwe set van 10.000 willekeurige punten te maken

    Notitie

    Als het aantal punten echt heel groot is kan de verwerkingstijd behoorlijk oplopen.

  2. Gebruik deze punten om een monster te nemen uit de originele DEM

  3. Gebruik het gereedschap Raster (IDW met zoeken naar Nearest neighbor) op deze gegevensset

  4. Stel Macht en Afvlakken respectievelijk in op 5.0 en 2.0.

De resultaten (afhankelijk van de positie van uw willekeurige punten) zal er min of meer zo uitzien:

../../../_images/interpolation_comparison_10000.png

Dit is een betere weergave van het terrein, vanwege de hogere dichtheid van monsterpunten. Onthoud: meer monsters geven betere resultaten.

6.4.9. In Conclusion

QGIS heeft een aantal gereedschappen voor het analyseren van ruimtelijke statistische eigenschappen van gegevenssets.

6.4.10. What’s Next?

Nu we vectoranalyse hebben behandeld, waarom niet eens kijken wat er met rasters gedaan kan worden? Dat is wat we zullen gaan doen in de volgende module!