Exercise : Measuring heights of trees outside the forest (TOF) in digital orthophotos (DOP)
(→Measuring shadow length and sun azimut in digital orthophotos) |
(→Measuring shadow length and sun azimut in digital orthophotos) |
||
Line 18: | Line 18: | ||
##** Confirm with {{button|text=Add to attributes list}}. | ##** Confirm with {{button|text=Add to attributes list}}. | ||
##* Confirm with {{button|text=OK}} and enter path and name (e.g. {{typed|text=measure_shadow_length.shp}}) in the following menu. | ##* Confirm with {{button|text=OK}} and enter path and name (e.g. {{typed|text=measure_shadow_length.shp}}) in the following menu. | ||
− | ## Zoom in to the | + | ## Zoom in to the overhead line tower 220kV shown in Fig. '''A'''). Click {{button|text=Add feature}} [[file:QGIS_2.0_AddLine.png|25px]] 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. 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). |
==Estimating sun elevation== | ==Estimating sun elevation== | ||
==Calculating individual tree heights== | ==Calculating individual tree heights== |
Revision as of 15:19, 15 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\).
Measuring shadow length and sun azimut 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. 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).
- Select Layer --> New --> New Shapefile Layer.