Interpoleren van puntgegevens

Interpoleren is een algemeen gebruikte techniek in GIS om een doorlopend oppervlak te maken uit afzonderlijke punten. Veel van de fenomenen in de echte wereld zijn doorlopend - hoogten, bodems, temperaturen etc. Als we deze oppervlakten voor analyses willen modelleren, is het onmogelijk metingen te doen over het gehele oppervlak. Daarom worden veldmetingen genomen op verschillende punten in het oppervlak en dan worden tussenliggende waarden berekend door een proces, genaamd ‘interpolatie’. In QGIS wordt interpolatie verricht met behulp van de ingebouwde gereedschappen voor interpolatie in de toolbox van Processing.

Overzicht van de taak

We zullen veld-dieptemetingen gebruiken voor Lake Arlington in Texas en een hoogte-reliëfkaart en contouren maken uit deze metingen.

Andere vaardigheden die u zult leren

• Contouren maken uit puntgegevens.

• Maskeren van waarden Geen gegevens uit een rasterlaag.

• Labels toevoegen aan een vectorlaag.

De gegevens ophalen

Texas Water Development Board provides the shapefiles for completed lake surveys.

Download the 2007-12 survey shapefiles for Lake Arlington.

Gemakshalve kunt u de in deze handleiding gebruikte voorbeeldgegevens direct downloaden vanaf de koppelingen hieronder.

Shapefiles.zip

Gegevensbronnen: [TWDB]

Procedure

1. Open QGIS, zoek in Browser en sleep de laag Arlington_Soundings_2007_stpl83.shp naar het kaartvenster.

2. Een dialoogvenster Transformatie voor Arlington_Soundings_2007_stpl83 selecteren zal verschijnen, laat de selectie op de standaardwaarden en klik op OK.

3. The layer will be added, now locate and drag the Boundary2004_550_stpl83.shp layer to canvas.

4. The layer will be added to the canvas, now turn off this layer to visualize the Arlington_Soundings_2007_stpl83 .

5. Klik op het pictogram Inzoomen en selecteer een klein gebied op het scherm. Als u meer inzoomt zult u de punten zien. Elk punt vertegenwoordigt een meting door een Depth Sounder op die locatie opgenomen met een gereedschap DGPS.

6. Selecteer het gereedschap Objecten identificeren en klik op een punt. U zult zien dat het paneel Identificatieresultaten aan de linkerkant zal worden weergegeven met de waarden voor de attributen van dat punt. In dit geval bevat het attribuut ELEVATION de diepte van het meer op die locatie. Omdat onze taal is een diepteprofiel en hoogte-contouren te maken, zullen we deze waarden gebruiken als invoer voor de interpolatie.

7. From Browser locate and drag the Islands_2004_550_stpl83.shp layer to canvas.

8. De laag zal aan het kaartvenster worden toegevoegd, deze laag bevat de informatie over de eilanden in de regio die een constante hoogte zouden moeten hebben (zouden niet moeten worden geïnterpoleerd).

9. Zoek en lokaliseer, vanuit de Toollbox van Processing, het algoritme Interpolatie ‣ TIN-interpolatie. Dubbelklik erop om het te starten.

Notitie

Resultaten van interpolatie kunnen significant variëren, gebaseerd op de methode en de parameters die u kiest. Interpolatie in QGIS ondersteunt de methoden Triangulated Irregular Network (TIN) en Inverse Distance Weighting (IDW) voor interpolatie. De methode TIN wordt algemeen gebruikt voor hoogtegegevens waar de methode IDW wordt gebruikt voor het interpoleren van andere typen gegevens zoals minerale concentraties, bevolking etc. Bekijk de module Ruimtelijke analyse van de documentatie van QGIS voor meer details.

10. Selecteer, in het dialoogvenster TIN-interpolatie, Arlington_Soundings_2007_stpl83 als de Vectorlaag, Elevation als het Interpolatie attribuut. Klik dan op het pictogram Invoerlaag(lagen).

11. Selecteer nu Islands_2004_550_stpl83 als de Vectorlaag, Elevation als het Interpolatie attribuut. Klik dan op het pictogram Invoerlaag(lagen). Wijzig nu het Type van de laag naar Lijnen opbreken.

Notitie

Lijnen opbreken stelt ons in staat onverwachte onderbrekingen in de hoogte te modelleren tijdens het modelleren van de lagen voor het oppervlak. Specificeren van het type laag naar Lijnen opbreken zal het algoritme voor interpolatie vertellen een constante hoogte te gebruiken voor de eilanden, in plaats van geïnterpoleerde waarden uit de punten.

12. Klik, in Bereik, op de … en selecteer de laag Boundary2004_550_stpl83.

13. In Grootte uitvoerraster, stel de Grootte pixel X en Grootte pixel Y in op 5. Klik dan op de … onder Geïnterpoleerd om de laag op te slaan als elevation_tin.tif. Klik op Uitvoeren.

14. Nu zal een nieuwe laag elevation_tin worden toegevoegd aan het kaartvenster.

15. Zoek en lokaliseer in Toolbox van Processing het algoritme GDAL ‣ Raster uitnemen ‣ Raster op maskerlaag clippen. Dubbelklik om het te starten.

16. In het dialoogvenster Raster op maskerlaag clippen, selecteer elevation_tin als de Invoerlaag, Boundary2004_550_stpl83 als de Maskeerlaag. Klik dan op de … onder Geclipt (masker) om de laag op te slaan als elevation_tin_clipped.tif. Klik op Uitvoeren.

17. Nu zal een nieuwe laag elevation_tin_clipped worden toegevoegd aan het kaartvenster. Klik op het pictogram Paneel Laag opmaken openen.

18. Stel de Symbologie in op Enkelbands pseudokleur, klik op de pijl in Kleurverloop en selecteer Kleurverloop omdraaien, voer 0 in bij Precisie label. Klik op Classificeren.

19. Zoek en lokaliseer in Toolbox van Processing het algoritme GDAL ‣ Raster uitnemen ‣ Contour. Dubbelklik om het te starten.

20. In het dialoogvenster Contour, selecteer elevation_tin_clipped als Invoerlaag, voer 5.000 in bij Interval voor contourlijnen. Klik op … in Contouren en voer de naam in als contour.gkpg. Klik op Uitvoeren.

Notitie

De interval wordt gespecificeerd in de eenheid van het CRS van de laag. Onze brongegevens staan in het EPSG:2276 NAD83 / Texas North Central (ftUS) - dus de interval voor coutouren zal worden geïnterpreteerd als 5 feet.

21. Nu zal een nieuwe laag contour worden toegevoegd aan het kaartvenster. Klik op het pictogram Paneel Laag opmaken openen. Schakel naar Labels.

22. Selecteer Enkele labels, in Waarde kies ELEV.

23. Schakel nu naar Plaatsing en wijzig de Modus naar Gebogen.

comments powered by Disqus