Exercise : Measuring heights of trees outside the forest (TOF) in digital orthophotos (DOP)
(→Calculating individual tree heights) |
(→Calculating individual tree heights) |
||
Line 46: | Line 46: | ||
# Digitze additional tree shadows. | # Digitze additional tree shadows. | ||
# Right click {{button|text=Open Atrribute table}}. {{button|text=Toggle editing}} mode [[file:QGIS_2.0_Edit.png|25px]]. Add a new column [[file:QGIS_2.0_AddColumn.png|25px]] Name: height, Type: Decimal number (real), Width: 10, Precision: 1. {{button|text=OK}}. Save the edits. | # Right click {{button|text=Open Atrribute table}}. {{button|text=Toggle editing}} mode [[file:QGIS_2.0_Edit.png|25px]]. Add a new column [[file:QGIS_2.0_AddColumn.png|25px]] Name: height, Type: Decimal number (real), Width: 10, Precision: 1. {{button|text=OK}}. Save the edits. | ||
− | # {{button|text=Open Field Calculator}} [[file:QGIS_2.0_FieldCalcualtor.png|25px]] and click the checkbox "Update existing field" on. Select height from the drop list and type an expression to | + | # {{button|text=Open Field Calculator}} [[file:QGIS_2.0_FieldCalcualtor.png|25px]] and click the checkbox "Update existing field" on. Select height from the drop list and type an expression to calculate the tree heights.{{button|text=OK}}. |
Revision as of 11:05, 16 July 2014
Contents |
Shadows in orthorectified remote sensing images
A simple application of the intercept theorem in elementary geometry is to determine the height of a tree by measuring the shadow length. According to Fig. A the shadow length \(s_1\) and the height \(h_1\) of a man-made object such as an overhead line tower were measured on the ground. If the lenght of a tree's shadow \(h_2\) is determined at the same time of the day we are able to compute \(h_2=s_2*\frac{h_1}{s_1}\).
We may also use trigonometry if we know the sun elevation \(\alpha\). The height of a tree is then computed \(h_2=s_2*\tan \alpha\).
Digitizing shadow length and measuring sun azimuth in digital orthophotos
Load and display German Geoadata using QGIS 2.4 following the Exercise: Displaying German Geobasis data (GBD) or load the saved project file .\GBData\display_gbd.qgs where the project coordinate reference system (CRS) is set to ETRS89/UTM32N (EPSG:25832).
- Create a new shapefile layer
- Select Layer --> New --> New Shapefile Layer.
- As layer type, select Line. Click the Specify CRS button and select ETRS89/UTM32N (EPSG:25832).
- To add an attibute:
- For the attribute's name type Class into the Name field of the New attribute section.
- Select Whole number as data type.
- Confirm with Add to attributes list.
- Confirm with OK and enter path and name (e.g. measure_shadow_length.shp) in the following menu.
- Zoom in to the overhead line tower 220kV shown in Fig. A. Click Add feature
to start digitizing. Simply click at first on the top of the shadow of the tower second add another node in the center of the base of the tower (Fig. B). Finish the line geometry by right-clicking and entering the attributes in the appearing window (just use an increasing number for the ID, and Class = 1 in case of a man-made object).
- Change the display style of the new line. Select the shapefile in the Layers window. Layer --> Properties --> Style --> Single Symbol. Click on "Simple line" and change the Color according to your preferences and increae the Pen width to 0.5.
- Do not forget to save your changes by clicking again the Toggle editing button
! Confirm Save to stop editing and save the changes permanently.
- Select Layer --> New --> New Shapefile Layer.
- Measure the sun azimuth at the time of image acquisition
- Install the CadTools plugin using Plugins --> Manage and install plugins.
- We use the "Show azimuth" function of CadTools to calculate the azimuth between two points which have to be selected first. Zoom in to the overhead line tower and the digitized shadow lenght. Select the shapefile in the Layers window. Activate Select 2 Vertex Points
and click at first on the node on top of the shadow of the line tower and second on the node near the base of the line tower. Calculate the azimuth by clicking Show Azimuth
. Make a note of the result in unit degrees.
Estimating sun elevation
If the position of the area of interest (WGS84 / Long/Lat) and date and exact time of the digital image acquisition is given we are able to compute the sun elevation \(\alpha\) using the Solar Geometry Calculator of the U.S. National Oceanic and Atmospheric Administration (NOAA). We assume that the date of the image flight is correctly indicated in the metadata as 2010-04-17 but the exact acquisition time is not known. We may construct a table for the day of the image aqcuision listing sun elevation and azimuth with a time intervall of one second.
- Find the position of the overhead line tower changing the project coordinate system from UTM to Long/Lat. Open Project Properties --> Project Properties --> CRC. Type into the text field "Filter": 4326 and than click on WGS84 to select geographic coordinates. Place the cursor on the line tower and make a note of the coordinates shown below the viewer in unit decimal degrees. Do not forget to switch the project CRS back to ETRS89 / UTM zone 32N, EPSG:25832 before we continue.
- Fill out the fields of the Solar Geometry Calculator and click Submit Query.
- Search in the resulting list for the sun azimuth calculated using the CadTools plugin (Fig. C) and report the corresponding sun elevation and acquisition time.
Calculating individual tree heights
- Vector --> Geometry Tools --> Export/Add geometry columns and select the vector layer measure_shadow_length. click OK and Yes. Close the dialog with OK and Close.
- Select measure_shadow_length in the layers window. Right click Open Atrribute table. Report the shadow length of the overhead line tower. The height of the tower was measured as 32.1m by a clinometer device. Compute the sun elevation \( \alpha \) and compare the result of the solar geometry calculator.
- Digitze additional tree shadows.
- Right click Open Atrribute table. Toggle editing mode
. Add a new column
Name: height, Type: Decimal number (real), Width: 10, Precision: 1. OK. Save the edits.
- Open Field Calculator
and click the checkbox "Update existing field" on. Select height from the drop list and type an expression to calculate the tree heights.OK.