-
Notifications
You must be signed in to change notification settings - Fork 1
/
index-en.html
140 lines (137 loc) · 7.61 KB
/
index-en.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
<!doctype html>
<html lang="en">
<head>
<meta charset="utf-8">
<meta name="viewport" content="width=device-width, initial-scale=1.0">
<title>profileAAR help</title>
</head>
<body>
<p>By Members of <a href="https://isaakiel.github.io">ISAAKiel</a> </p>
<p>Moritz Mennenga (Lower Saxony Institute for Historical Coastal Research Wilhelmshaven)</p>
<p>Kay Schmütz (Graduate School - Human Developent in Landscapes / CAU Kiel)</p>
<p>Christoph Rinne (Institute for Pre- and Protohistoric Archaeology / CAU Kiel)</p>
<br>
<p><b>Please cite: Moritz Mennenga, Kay Schmütz, & Christoph Rinne. (2023). ISAAKiel/profileAAR: profileAAR (2.0.3). Zenodo. https://doi.org/10.5281/zenodo.8228250</b></p>
<br>
<p>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.</p>
<br>
<b>No warranty. Check your data after conversion</b>
<br>
<br>
<h1>How to use the profileAAR-Plugin</h1>
<p>ProfileAAR is a plugin to transform measured photogrammetry points to allow equalization of sections in GIS.</p>
<p>The points are rotated and shifted according to different specifications (see below), </p>
<p>the determination of the rotation parameters is based on a linear regression.</p>
<h2>1. Input data</h2>
<p>You require a shapefile with points. Several sections can be processed at the same time
<p>The shapefile must have at least the following attributes</p>
<ul>
<li> the <b>z-value</b> [float/double/real] </li>
<li> the label of the <b>section</b> (e.g. number of the section [text or number])</li>
<li> the <b>view</b> direction (defined by the position of the observer related to the section (view from ...) as N, E, S or W [text])</li>
<li> the <b>used Points</b> (points that are used to determine the slope of the regression [0/1])</li>
<img src="./img/view.jpg" alt="View" style="width:250px;height:228px;">
<p>If a section is NE - SW orientated it doesn't matter if you choose E or S.</p>
</ul>
<p></p>
<p>Example for an attribute list</p>
<table>
<tr>
<th>z-Value</th> <th>section_no</th> <th>view</th><th>use</th>
</tr>
<tr>
<td>0.5</td> <td>1</td> <td>N</td> <td>1</td>
</tr>
<tr>
<td>0.5</td> <td>1</td> <td>N</td> <td>1</td>
</tr>
<tr>
<td>0.25</td> <td>1</td> <td>N</td> <td>1</td>
</tr>
<tr>
<td>0</td> <td>1</td> <td>N</td> <td>1</td>
</tr>
<tr>
<td>0</td> <td>1</td> <td>N</td> <td>1</td>
</tr>
</tr>
<tr>
<td>1.8</td> <td>2</td> <td>E</td> <td>1</td>
</tr>
<tr>
<td>1.5</td> <td>2</td> <td>E</td> <td>1</td>
</tr>
<tr>
<td>1</td> <td>2</td> <td>E</td> <td>0</td>
</tr>
<tr>
<td>0.75</td> <td>2</td> <td>E</td> <td>1</td>
</tr>
<tr>
<td>0.8</td> <td>2</td> <td>E</td> <td>1</td>
</tr>
</table>
<p></p>
<h2>2. Options</h2>
<p>First, the layer and then the associated attributes in the dropdown menus must be selected. <b>z coord</b>, the <b>pr number</b>, the <b>use</b> and the <b>view column</b> boxes. </p>
<p>In this case <b>use</b> 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'.</p>
<p>Please select the matching attributes.</p>
<p>In the next two boxes, you have to choose:</p>
<ul>
<p>In the <b>Method</b> box:</p>
<li>1.: <b>Projected</b>. The section will be projected on a vertical projection pane. This is like a drawing in the field.</li>
<li>2.: <b>Surface</b>. The section will be projected on a pane, that is parallel to the section. This is like drawing the surface of the section.</li>
<img src="./img/method.jpg" alt="Method" style="width:630px;height:234px;">
<p> Sideview of a section (0). The result after choosing projected (1) or surface (2).</p>
<p>In the <b>Direction</b> box:</p>
<li>1.: <b>Original</b>. The points will be parallel to the section.</li>
<li>2.: <b>Horizontal</b>. The points will be parallel to the x-axis.</li>
<li>3.: <b>original height</b>. Same as horizontal but the section will be at x = 0 and y is the measured height. This allows an export with a coordinate frame.</li>
<p><i>Before using the third option <b>original height</b> always check whether the given coordinate system is still constant in scale at the origin</i></p>
<p></p>
<img src="./img/direction.jpg" alt="Direction" style="width:630px;height:234px;">
<p> View on a section(0). The result after choosing original (1) or horizontal (2).</p>
</ul>
<p><b>Export elevation points</b> will export an additional shapefile containing the upper right photogrammetric control point, its height, and the section number.</p>
<p><b>Export section</b> will export an additional shapefile containing the section based on the highest outer points.</p>
<p>Thanks to Nils Müller-Scheeßel for the idea.</p>
<p>This can be used to show the line with the point of view. </p>
<p>Thanks to Nils Müller-Scheeßel for the idea.</p>
<p>Open the layer settings of the section -> Lables -> </p>
<p>Text -> Lable with '↓'</p>
<p>Placement -> [x] Line orientation dependent on position</p>
<p>Rendering -> Show upside-down lables -> [x] always</p>
<p></p>
<p>This can be very helpful while processing the drawings.</p>
<h2>3. Run the calculation</h2>
<p>Some hints for a good result:</p>
<ul>
<li>If you have a straight section you will have the best results.</li>
<li>If sections are very long and not very straight, sometimes it is better to split them (we plan a function for the next version to do it the easy way for the user).</li>
</ul>
<h2>4. Accuracy</h2>
<p>Getting values for the accuracy of the points.</p>
<p>When you open the attribute table of the rotated points, you will find the column <b>distance</b>.</p>
<p>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.</p>
<p>If there is one point not in the line, this one will have a high distance. If you choose 0 in the original column <b>used Points</b>, this point will be excluded</p>
<p>in the calculation of the regression line, but the point will be rotated. This can be used to correct measurement errors.</p>
<h2>5. Error handling</h2>
<p>An error handler will check the data during the calculation</p>
<p>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)</p>
<ul>
<li> <i> A section needs min. 4 points. Error on section: XXX </i> </li>
<p>Section XXX has 4 or less points.</p>
<li> <i>The view column of your data is inconsistent (either none or two different views are present). Error on section: XXX </i> </li>
<p> The view of the sections contains different or no values for the points. </p>
<li> <i>The view value is not one of the four cardinal directions. Error on section: XXX </i> </li>
<p>The view column holds another direction than N, E, S or W.</p>
<li> <i>Only 0 or 1 are allowed in the selection/use. Error on section: XXX </i> </li>
<p>In the use column only 1 or 0 are allowed as values.</p>
<li> <i>Layer XXX is not projected. Please choose a projected reference system.</i> </li>
<p>The coordinate system of the layer is not a projected one. Save the layer in a projected system.</p>
<li> <i>The z-value needs to be a float. Check the field type of the z-value</i> </li>
<p>Save the column containing the z-value as float</p>
<p><b>You will get the best results if your work carefully in advance on the excavation and prepare a perfectly straight section!</b></p>
</ul>
</body>
</html>