qgSurf - Plane distances from points/DEM cells - help
This tool allows for the calculation of the perpendicular distance between a plane and points
or centers of cells of a grid.
The plane is expressed via its attitude (dip direction and dip angle) and a point located on it (the "source point").
The calculated distance is positive when the point or cell center is above the plane, and negative
when below the plane.
Input parameters
There are two alternative ways to input the analysis parameters: by loading
previously defined and exported parameters via the "Load analysis parameters
... from file" or by defining the source point and the plane orientation in the
"Plane -> source point definition" and "Plane -> orientation definitions" tabs
(Figs. 1, 2).
The geological plane is defined by a "source point" laying on the considered plane (Fig. 1)
and its attitude (dip direction and dip angle, Fig. 2).
The source point position is defined via its latitude and longitude, as expressed in the
WGS-84 datum, and its elevation above a reference level (usually the geoid).
The source point can be positioned graphically in the map canvas, by pressing the
"Set source point in map and afterward clicking in the map. When the "Get elevation from DEM"
check box is on, the elevation is derived from the chosen DEM (Fig. 1).
Fig. 1. Example of source point definition in the tool GUI.
The attitude of the plane is expressed via its dip direction and dip angle (Fig. 2). The dip direction azimuth may refer
to the Geographic North or alternatively to the top direction of the current map. The top direction, depending on the chosen CRS,
can differ from the geographic North by a few degrees.
Fig. 2. Example of definition of the plane attitude.
Calculations
As previously said, two types of calculations can be performed:
distances between the plane and the cell centers of a grid, within a chosen interval;
distances between the plane and a set of points stored in a layer.
Distances between the plane and a grid
The cell center-by-cell center orthogonal distances between a grid and a plane
are calculated by defining the grid and within a range defined by a maximum and a minimum distance
thresholds from the plane (Fig. 3).
Generally the maximum distance threshold is positive (i.e., when the cell center is above the plane) while the
minimum distance threshold is negative (i.e., when the cell center is below the plane).
Fig. 3. Example of calculation of distances of grid centers within a range
comprised within 0 up to 100 meters above the defined plane. The result is a grid of the distances, that can be automatically loaded into the project (Fig. 4).
In the example illustrated in Fig. 4, the maximum orthogonal thicknesses of the Triassic meta-limestones in the
hanging-wall of the Tempa Carlone thrust (Mt. Alpi, Basilicata, Southern Italy) are 50 meters (Fig. 4).
Fig. 4. Example of resulting grid ("distances_04") of orthogonal distances
between the topographic surface ("malpi_raparo_w4u2n"),
from a geological plane with attitude 143.3/20.6 (dip direction/dip angle), the
Best-Fit-Plane to the trace of the Tempa Carlone thrust (Mt. Alpi, Basilicata, Southern Italy).
Distances between the plane and a set of points
In the case of the calculation of the distances between the plane and a set of points
stored in a layer, the user has to define two required fields in the input point layer (Fig. 5):
id: the record identifier field
z: the point elevation field
Moreover, the user defines the output point layer storing the calculated distances.
Fig. 5. Parameters for the calculation of the perpendicular distances
between the points of a layer and a plane.
The resulting variables are stored in a new shapefile, in WGS84 (EPSG:4326), where
the orthogonal distances are stored in the "distance" field (Fig. 6).
Fig. 6. The distances (in meters) between the Best-Fit-Plane to the Tempa Carlone thrust (Mt. Alpi, Basilicata, Southern Italy) trace
and the trace points.