Exercise: Mapping storm damaged forest areas
(→Estimating area and perimeter) |
(→Estimating area and perimeter) |
||
Line 31: | Line 31: | ||
# Make sure that the edit mode of the vector polygon layer is switched off. | # Make sure that the edit mode of the vector polygon layer is switched off. | ||
− | # {{mitem|text=Vector --> Geometry Tools --> Export/Add geometry columns}}. Select the vector layer clearings and click {{button|text=OK}} and {{button|text= | + | # {{mitem|text=Vector --> Geometry Tools --> Export/Add geometry columns}}. Select the vector layer clearings and click {{button|text=OK}} and {{button|text=Yes}}. Close the dialog with {{button|text=OK}} and {{button|text=Close}}. |
# Select the clearing layer in the layers window. Right click {{button|text=Open Atrribute table}}. Report area of the clearing in hectares and perimeter in meters. | # Select the clearing layer in the layers window. Right click {{button|text=Open Atrribute table}}. Report area of the clearing in hectares and perimeter in meters. | ||
==Using Google Maps for digitizing== | ==Using Google Maps for digitizing== |
Revision as of 11:59, 11 July 2014
Contents |
Definition of storm damage
Mapping of forest damages caused by storms using visual interpretation of very high resolution digital ortho photos (DOP) need clear definitions and characteristics which are usually documented as interpretation key. Here, clearings are defined with a minimum area = 0.2ha, minimum width = 20m, minimum tree crown cover = 10%.
Characteristics of windthrow areas depend on the time between storm events and the availabiltiy of optical remote sensing images but also on the status and progress of logging and clearing operations at the time of the image aqcuisition. Typical characteristics for recognizing storm damages are: horizontal complete trees with root plates, root plates, stumps with fresh cross section, piles of logs for storage / skidding, piles of branches, skid roads and tracks, no or few ground vegetation.
Digitizing windthrow areas
Load and display German Geoadata 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 Polygon. 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 the path and name (e.g. clearing.shp) in the following menu.
- Select Layer --> New --> New Shapefile Layer.
- Digitizing areas
- Select the new shapefile in the Layers window. Start the edit mode by clicking the Toggle editing button .
- Zoom in on a cleared area in the forest (figure A). Click Add feature to start digitizing. Draw polygons according to the definition of clearings. Simply click the desired point on the canvas to start. After selecting a second point, the third point will construct a triangle and each additional point will add another vertex to a polygon. Finish the 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 windthrow area). Take your time to get used to this way of digitizing -- you can always delete the last polygon by simply clicking Cancel in the attributes dialogue, or by selecting it and deleting with the Delete selected button. Do not forget to save your changes by clicking again the Toggle editing button ! Confirm Save to stop editing and save the changes permanently.
- Change the display style of the new polygon. Select the shapefile in the Layers window. Layer --> Properties --> Style --> Single Symbol. Click on "Simple fill" and change the Fill style: diagonal, Border width: 0.5, and the colors of "Colors Fill" and "Borders" according to your preferences. An example of a finished polygon can be seen in figure B.
Estimating area and perimeter
- Make sure that the edit mode of the vector polygon layer is switched off.
- Vector --> Geometry Tools --> Export/Add geometry columns. Select the vector layer clearings and click OK and Yes. Close the dialog with OK and Close.
- Select the clearing layer in the layers window. Right click Open Atrribute table. Report area of the clearing in hectares and perimeter in meters.