By Members of ISAAKiel
Moritz Mennenga (Lower Saxony Institute for Historical Coastal Research Wilhelmshaven)
Kay Schmütz (Graduate School - Human Developent in Landscapes / CAU Kiel)
Christoph Rinne (Institute for Pre- and Protohistoric Archaeology / CAU Kiel)
Please cite: Moritz Mennenga, Kay Schmütz, & Christoph Rinne. (2023). ISAAKiel/profileAAR: profileAAR (2.0.3). Zenodo. https://doi.org/10.5281/zenodo.8228250
Addtional information: Mennenga M, Schmütz K. profileAAR – a Plugin to transform section control points for rectifying archaeological section images in QGIS. Siedlungs- und Küstenforschung im südlichen Nordseegebiet. 2022(45):227–34.
ProfileAAR is a plugin to transform measured photogrammetry points to allow equalization of sections in GIS.
The points are rotated and shifted according to different specifications (see below),
the determination of the rotation parameters is based on a linear regression.
You require a shapefile with points. Several sections can be processed at the same time
The shapefile must have at least the following attributes
If a section is NE - SW orientated it doesn't matter if you choose E or S.
Example for an attribute list
z-Value | section_no | view | use |
---|---|---|---|
0.5 | 1 | N | 1 |
0.5 | 1 | N | 1 |
0.25 | 1 | N | 1 |
0 | 1 | N | 1 |
0 | 1 | N | 1 |
1.8 | 2 | E | 1 |
1.5 | 2 | E | 1 |
1 | 2 | E | 0 |
0.75 | 2 | E | 1 |
0.8 | 2 | E | 1 |
First, the layer and then the associated attributes in the dropdown menus must be selected. z coord, the pr number, the use and the view column boxes.
In this case use means whether the point will be used to calculate the rotation angle (coloumn with '1'). If there is an outlier in the data, mark it with a '0'.
Please select the matching attributes.
In the next two boxes, you have to choose:
In the Method box:
Sideview of a section (0). The result after choosing projected (1) or surface (2).
In the Direction box:
Before using the third option original height always check whether the given coordinate system is still constant in scale at the origin
View on a section(0). The result after choosing original (1) or horizontal (2).
Export elevation points will export an additional shapefile containing the upper right photogrammetric control point, its height, and the section number.
Export section will export an additional shapefile containing the section based on the highest outer points.
Thanks to Nils Müller-Scheeßel for the idea.
This can be used to show the line with the point of view.
Thanks to Nils Müller-Scheeßel for the idea.
Open the layer settings of the section -> Lables ->
Text -> Lable with '↓'
Placement -> [x] Line orientation dependent on position
Rendering -> Show upside-down lables -> [x] always
This can be very helpful while processing the drawings.
Some hints for a good result:
Getting values for the accuracy of the points.
When you open the attribute table of the rotated points, you will find the column distance.
This value represents the distance of each point to the calculated section between the points. The higher the value is, the worse the result will be.
If there is one point not in the line, this one will have a high distance. If you choose 0 in the original column used Points, this point will be excluded
in the calculation of the regression line, but the point will be rotated. This can be used to correct measurement errors.
An error handler will check the data during the calculation
There are a few errors and warnings possible. All Errors will be shown in a red bar and in the ErrorLOG protocol in the protocol window (View -> Panels -> Log Messages Panel)
Section XXX has 4 or less points.
The view of the sections contains different or no values for the points.
The view column holds another direction than N, E, S or W.
In the use column only 1 or 0 are allowed as values.
The coordinate system of the layer is not a projected one. Save the layer in a projected system.
Save the column containing the z-value as float
You will get the best results if your work carefully in advance on the excavation and prepare a perfectly straight section!