Remote Sensing Detecting Mountain Peaks and Delineating Their Shapes Using Digital Elevation Models, Remote Sensing and Geographic Information Systems Using Autometric Methodological Procedures

The detection of peaks (summits) as the upper parts of mountains and the delineation of their shape is commonly confirmed by inspections carried out by mountaineers. In this study the complex task of peak detection and shape delineation is solved by autometric methodological procedures, more precisely, by developing relatively simple but innovative image-processing and spatial-analysis techniques (e.g., developing inventive variables using an annular moving window) in remote sensing and GIS domains. The techniques have been integrated into automated morphometric methodological procedures. The concepts of peaks and their shapes (sharp, blunt, oblong, circular and conical) were parameterized based on topographic and morphologic criteria. A geomorphologically high quality DEM was used as a fundamental dataset. The results, detected peaks with delineated shapes, have been integratively enriched with numerous independent datasets (e.g., with triangulated spot heights) and information (e.g., etymological information), and mountaineering criteria have been implemented to improve the judgments. This holistic approach has proved the applicability of both highly standardized and universal parameters for the geomorphologically diverse Kamnik Alps case study area. Possible applications of this research are numerous, e.g., a comprehensive quality control of DEM or significantly improved models for the spatial planning proposes.


Introduction
The aim of this study is the detection of mountain (called also topographic) peak points and the analysis of their surroundings in order to determine their shape.Automated morphometric techniques and procedures are required to solve this complex task.The methods developed are intended to be studied holistically according to a number of other points of view in the context of spatial studies, in order to arrive at a reliable conception of a peak.The following three associated introductory challenges are relevant here: (1) semantics, definitions, conception and standardizations (2) appropriate data sources (3) morphometric algorithms based on a digital elevation model (DEM)

Semantics, Definitions, Conception and Standardizations
A landform peak can be defined as a point higher in elevation than the adjacent area.The shape is a geometric property of the peak point on an arbitrarily selected surrounding area (depending on the scale or level of detail) that is morphologically expressed as being sharp, blunt, oblong, circular, conical, or other.Both definitions are vague and subject to further discussion.
Professional mountaineers have perhaps the best knowledge of the peaks and other characteristics of the mountains from different regions worldwide, but it partly depends on their subjective perception and experiences.Researchers have been trying to measure and standardize mountain peaks for at least the past 120 years [1][2][3][4][5][6].Defining a peak is a conceptual problem that requires a certain level of abstraction which is an essential part of the model, and which should have subsequent feasibility for the applicable automation process.The concept of a peak is contextually associated with several other concepts, such as the concept of the whole mountain, its roughness, locations of prominence and others, all of which are important for an enriched understanding of the problem in question.Size and shape (in context) are important criteria for the definition of peaks as a landform category or even for the construction of the landform taxonomy [7].Until now, definitions of peaks have been based mostly on tradition, experience and visual appreciations-on a subjective perception of the natural appearance of the peaks [1].The different views of the problem studied in this research can serve as sources of a comprehensive definition of a mountain peak, which can be further used as a source for an enhanced ontology of mountain peaks.
The concepts and definitions of the term "peak" must be considered in association with the wider topography (regional and global) and its relation to the term "mountain", considering the entire mountain or range.The concept of a peak as a topographic and morphologic extreme in a specific geographical space has several other meanings in a semantic sense.However, our aim here is not to find absolute definitions, but merely to discuss the complexity of the subject of our study.
Just as the most significant peaks are visible from a distance, in the immediate neighborhood the peaks on the lower rising ground are also visible.The peaks may appear different when viewed from various distances (e.g., the peak appears oblong from 100 m, but circular from 1 km), from different directions-oblique or vertical views (e.g., sharp from one valley and blunt from another), when view at different times of the day (Figure 1), and also in different weather conditions (including the visibility phenomenon) with different snow cover or during different seasons.Furthermore, peaks are conceived in different ways worldwide, which impacts on the terminology.For instance, the characteristics of the Alps, Carpathians, Himalayas or even the Pannonian Plain are quite distinctive.These factors depend on physical and environmental properties, such as the elevation above sea level, geological properties, the shapes of peaks, the illumination (Figure 1), vegetation, and on the particular perception of an individual observer.These findings may lead to mentally various perceptions and consequently to different identification, conceptualization and classification by the observer.In order to provide a proper description of the peaks, the landscape in question and its properties has to be considered.An image processing technique of multidirectional visibility index (MVI) [8] shows another aspect in appearance of peaks and other terrain features (a composition of series of 57 photograph visibility masks between 6:30 AM and 17:15 PM (GMT+1), 26 September; view from Kredarica to peak Tosc; webcam of ARSO).
The relation of the peak to the roughness concept of the relief refers to the irregularity of a terrain, which may even be defined with self-similar fractals.The terrain roughness is geomorphometrically related to the periodic functions, where a texture is used to refer to the shortest wavelength in topography, and the grain used to refer to the longest significant wavelength [5].Three types of peaks will be discussed according to their different size-scale (and level of detail): local, regional and global.The local peaks represent significant maximum points of the texture, while the regional and global peaks are related to the grain.The concept of roughness raises the many further complexities of the peaks e.g., the relation to aesthetic quantification.
Within common concepts, the peak is defined as being any point on a surface that is elevated by a certain difference in height with respect to its surroundings.A peak must be autonomous in the sense of possessing individuality, interest and other characteristics.The peaks might be considered in addition to the features that satisfy the especially subjective criterion of a "well defined morphology".These primary criteria are referred to as topographic, morphologic and mountaineering [2].
According to dictionary definitions, standards, and encyclopedias, e.g., [7,9,10], the term peak is less technically defined, i.e., mostly with reference to the morphologic properties of the surface that surrounds the higher point.These definitions may therefore be vastly dissimilar.Possible peak definitions are: "any rising ground that splits the mountain in its higher part"; "upper, usually folded part of the mountain"; "the highest point of the mountain"; "the point on a surface that is higher in elevation than all points immediately adjacent to it".More precise definitions of this phenomenon concern size and shape and distinguish between the terms "peak" and "summit".The summit means the highest elevation on a profile, whilst the concept of a peak depicts part of a summit with moderate to very steep sides [11].However, this last definition already reflects a corresponding ontology and cultural semantic aspects (e.g., concepts are not the same in English and the Slovenian language) [1,12].The term "peak" is used in the study to correspond to a point.
The concepts of regional and global peaks should also be discussed.Regional peaks are complex features which occur in the most common definitions of peaks that are found in basic dictionaries, described by the UIAA, and defined by people who live in the mountains.This term is defined by a broader range of topographic, morphologic and other so-called mountaineering properties on a regional scale [2].Global peaks are a subset of regional peaks.Global peaks are limited to the most prominent or salient peaks of a given mountain (range), e.g., the entire Kamnik Alps.Global peaks have acquired cultural importance and have been of immense interest to humans over centuries or millennia.Therefore, they are defined according to the mountain range and/or according to their prominent location.
Peaks are also defined by referencing topographically related mountain ranges (i.e., groups of mountains), mountains, massifs, hills, or hill-like features [7].Similarly, the terms "mountain" and "hill" also vary in definition worldwide.Thus, only the term "mountain" is discussed here in order to support its wider geo-connoted conception.The concept of a mountain is more complicated and subjective than the concept of a peak, and it involves aesthetical, ecological, chronological, and other factors.The term mountain is not a real physical (bona fide) object.It can even be interpreted as a fiat object, which only exists in the human understanding and division of the landscape [12][13][14][15].However, it is clear that mountains exist as real correlations of everyday human thought and action, and that they form the archetype for geographic objects [12].
A mountain can be defined as a natural rise in the Earth's surface that usually has a peak (or a summit) [8].Similarly, the term mountain is defined as a landform that sufficiently extends above the surrounding terrain over a limited area with comparatively steep sides and a peak.A peak (a highest point or a highest ridge) and the mountainside (part of a mountain between the peak and the foot) are considered to be the main sections of a mountain [16].
The mountain as an object comprises the following morphometrical parameters: size, perimeter length, maximum elevation, slope, mean or average slope, slope gradient, relative or local relief, relative massiveness (according to stage of erosion), shape, roughness, texture and grain, etc. [5,17,18].The following are those definitions that contain quantitative operational definitions: mountains are usually steeper and taller than hills; they are often considered to be hills that rise over 600 m (a.s.l.) in height; another study utilizes the higher vegetation zones as the important criterion for defining it as a mountain [4]; a mountain is defined with an altitude (height), slope and relative height distinction; the lowest limit is 300 m in the altitudes of 65°N and 55°S and up to 1,000 m at the equator.A team of Kapos [19] used criteria based on altitude and slope in combination with a depiction of the world's mountain environments to form their definition.They proposed seven classes based on slope and local elevation range.The most important for the definition of a mountain are: (classes 1 to 3) elevations >2,500 m; (class 4) elevations 1,500 to 2,500 m and slope ≥2°; (class 5) elevations 1,000 to 1,500 m and slope ≥ 5° or local elevation range >300 m (radius of 7 km); and (class 6) elevations 300 to 1,000 m and local elevation range >300 m outside 23°N to 19°S (radius of 7 km).

Appropriate Data Sources
The mapping of peaks and mountains has generally been performed manually through fieldwork and the visual interpretation of topographic maps, aerial photographs and satellite images.An example (Figure 2) presents an old Josephine military topography (scale 1:28,800) from the end of 18th century that was verified with an overlaid modern topographic map demonstrating that the peak positions are in this case considered to be correct (scale 1:25,000, DTK 25).On the old map, some peaks are clearly detected but with their positions have been incorrectly mapped.One of the reasons for this is that some peaks had not been climbed yet when this old map was produced.
The problem of the remote sensing techniques that use aerial photographs and satellite images, from which modern maps are produced, is that a large amount of semantic information is inherently hidden in these data sources.Moreover, the detection and interpretation of topographic peaks is problematic in forest areas.
Another potential data source for studying topographic peaks is a digital elevation model-DEM (or digital terrain model-DTM).The DEM allows a very high level of automation for the study of geomorphologic features [20].Its geomorphological quality and appropriate (high) scale/resolution (Table 1) are crucial for successful applications [21,22].However, the lack of both properties is a distinctive problem for real data sources application.Typical sources of the DEM are maps, aerial photographs and satellite images, as mentioned before, which is not encouraged.Fortunately, more independent and higher quality data sources are available nowadays based on LIDAR or SAR technologies as important remote sensing techniques for topographic data collection.
A geomorphologically higher quality DEM would be required in order to study the peaks of flatter areas, compared to the requirements in the hilly areas, such as the Pannonian Plain.The resolution of the DEM controls an optimal scale and the level of detail of the peak determination, particularly local peaks.It should reflect the discussed topographically-relevant terrain roughness.The integrative study of more alternative datasets (and techniques) in a geographical information systems (GIS) environment can improve the judgments in morphologically enhanced studies [22].Table 1.The characteristics and detection potential of local, regional and global peaks.Regional wider surroundings, area of summit regional extreme, topographic, morphologic and mountaineering criteria complex extraction problem: fulfilling possibly standardized criteria and parameters [2] Global related to entire mountain, range "globally" prominent points, many different criteria complex perception and definition problem; regional peak implementation criteria are expanded

Morphometric Algorithms Based on a Digital Elevation Model (DEM)
The extraction techniques of peaks have evolved considerably-from manual in the past, through computer assisted, to various automated methods over the last decades.The automated application of the detection of peaks and the delineation of their shape is discussed through morphometric algorithms based on a geomorphologically accurate DEM.The algorithms for local peak detection are considered as an ordinary problem, while the algorithms for topographic (regional and global) peak detection demand more complex solutions to the problem of feature extraction and recognition, all based on morphometric image processing or spatial analysis techniques.
Local peaks are computed using the most elementary numerical (computer-based) peak detection approach from the DEM.The algorithms are based on an assessment of local extremes, e.g., the Peucker and Douglas method [23] that detects every rise in ground elevation.The problem can be resolved with tools as part of both: standard remote sensing (image processing) and GIS (spatial analysis).
The existing solutions for the numerical implementation of regional peaks incorporate mostly only topographic criteria, e.g., an algorithm to extract the number of peaks based on a morphological smoothing method [24], and an algorithm based on the minimum drop surrounding peak criterion [29].
The procedures for global peak detection use a considerably larger surrounding area than the procedures for regional peak detection.Their definition is either related to (A) the entire mountain range, or (B) a prominent location.
(A) Most of the developed automated techniques concentrate on the detection or classification of whole mountains or mountain ranges, e.g., the algorithms for feature-matching [3], the componentlabeling algorithm applied to the mountain terrain class [17], a local histograms analysis [25], and mathematical-morphological based algorithms [16].
(B) A number of studies have attempted to establish prominent locations (peaks) using numerical algorithms.Their characteristic technique is to study specifically cultural responses that influence human behavior, including beliefs, taboos, and rituals, etc. Prominent peaks are explored by analyzing their properties according to their wider surroundings.This includes: an exposure index, cumulative viewsheds, defining location hierarchies (nested peaks), taking into account the fuzziness of general multi-scale landscape morphometry, a comparative evaluation according to a multi-scale prototypicality of peaks identified, a computation of horizon lines based on a viewshed analysis, and the identification of the salient points on the horizon based on a hierarchical decomposition of the projected horizon, etc. [14,[26][27][28][29][30].A comparison of the automated implementation of local, regional and global peaks can be seen in Table 1.
The comprehensive peak detection procedure should also consider the morphological shape of peaks.The topographic shape delineation of a peak is a more complex conceptual and implementational problem than its detection.The algorithms that touch on the problem of shape analysis relate to the issue of peak prominence discussed in the previous paragraph.The automated shape analysis procedures can additionally improve the strategy for detecting peaks.

The Problem
Considering the challenges introduced, the technical description of a regional peak has been established according to the rules of the International Mountaineering and Climbing Federation (UIAA) [2].The rules of the UIAA are formulated for those peaks higher than 4,000 m in the Alps according to the following three criteria: topographic, morphologic and mountaineering (alpine).While the mountaineering criterion is the most subjective and qualitative, the topographic and morphologic criteria have a certain potential to be processed by quantitative automated procedures.All criteria together can be considered as a qualitative peak analysis and reasoning [31].

The Concept of Regional Peak Parameterization
The topographic criterion defines the level of height difference between each peak, and specifies that the highest adjacent pass or notch should be at least 30 m (calculated as an average of the peaks at the limit of acceptability).An additional criterion can be the horizontal distance between a peak and the base of another adjacent peak [2].The suggested application of this rule is by a geometric analysis on the DEM: the main criterion being a definition of the relative height (dH) of the potential peak that should be split by the adjacent peaks with a certain gap (which should be 30 m lower than the peak).An additional criterion is defined as the minimum horizontal distance (dL) between the peaks, which is related to peak density [5], which in turn is related to grain as part of the roughness concept.These values help to define the dominance of the potential peak.The larger the measured values, the more extensive and significant the peak becomes.
The morphologic criterion takes into account the overall morphology and aspect of the peak, especially when applied to shoulders, secondary peaks, rocky outcrops, etc. [2].The proposed application focuses various techniques of spatial analysis for determining the various sizes of the immediate peak neighborhood (determination of mountain area G p , potential surfaces P ij , see Table 2 for appropriate indexes).The shape of the upper part of the mountain is therefore considered, and much depends upon its geological characteristics.
Table 2. Digits that determine the categories (i) and (j) of a peak scale and shape.The category (j) = 9 comprises (any) undefined shape.The following mountaineering criterion takes into account the importance of the peak from the viewpoint of mountaineers: the quality of the routes that lead to the peak, the historical significance of the peak, and the frequency it has been climbed [2].This criterion is the weakest and most difficult for standardization and numerical application.It considers the inspection of terrain as a subjective human perception.With the evolution of mountaineering, the criterion may change in the future.The mountaineering criteria are not intended for parameterization; they are simply used as means of verifying and interpreting the results.

The Concept of Peak Shape Parameterization
The morphologic criterion of the regional peak is studied in greater detail.Morphological shape of the peak can be explained through typical categories of sharp, blunt, oblong, circular, conical shapes, etc.The shape categorization is a complex semantic and conceptual problem that may be used to support the peak definition.The definition of the shape might be developed on different levels.
In this study, a two-digit numeric code system is proposed.Category i explains the properties of scale/size of the peak, whilst category j describes the (dimensionless) shape of the peak.Thus, this code system relates the determination of the peak to the shape delineation problem.For example, an oblong regional peak is, according to Table 2, assigned the number 23, which is a composition of category i and category j to ij.The characteristics of the shape of a peak are dependent on each other and are not universal.Table 3 shows the linkages between the shapes (category j; Table 2) that are applicable in the shape delineation modeling.For example, the most complex, conical shape should include sharp and circular shapes (black backgrounds of field), and should exclude blunt and oblong shapes (white backgrounds).Table 3.The linkage between the studied shapes of peaks (category j).White backgrounds of fields denote no connection, grey denote a possible connection, and black denote a definite connection between the shapes.The arrows on the left indicate the direction of reading along the rows and obtaining the attributes from the columns (e.g., a conical peak should be sharp, but a sharp peak is not necessarily conical).
Additional morphological shape categories might be studied beyond those that are listed in the peak taxonomy of Tables 2 and 3 (e.g., a pyramidal peak).The category of shape featured in Table 2 may be further (sub)classified and additionally coded to categories that describe the xy ratio (2D footprint or floor plan, e.g., eccentricity or orientation) and the vertical to horizontal ratio (cross section, e.g., slope or concave/convex) parameters of the shape of a peak.

Methodology
No single criteria or parameter can be used to determine the complex different terrain forms.A combination of significant parameters can increase the reliability of terrain measurements [11].Two procedures are proposed.The first has been developed to detect the points of topographic (regional and partly global) peaks.Topographic and some morphologic criteria based on the UIAA directive [2] are employed using a numerical procedure and algorithms, while selected mountaineering criteria are used for their verification and interpretation.The second procedure has been developed to analyze the surrounding area in order to determine their shapes.The morphologic criteria are realized according to algorithms that use classified peak shapes.
The methodology of both interconnected procedures is considered to be part of the specific (in contrast to general) (geo)morphometry that applies to and describes discrete landforms such as peaks or talus cones and can involve arbitrary decision-making in the quantification of its concepts [5,32].The expectation for this methodology is that it demonstrates a certain potential for bridging qualitative peak concepts and terminology to achieve more quantitative definitions of geographically delineated peaks.

The Detection of Peaks by Topographic and Morphologic Criteria
The following five-step procedure is based on the classification and extraction algorithms for fulfilling the basic topographic and morphologic criteria for detecting regional peaks: (1) the calculation of local peaks as a local maximum in elevation [23] (2) the selection of peaks that are morphologically not on flat areas (peaks on plains are denoted as G p ) (3) a further selection process considering a minimum horizontal distance (dL) between peaks (4) a further selection process considering a relative elevation (dH) of peaks (5) a final selection of the regional peak points considering the potential surfaces of peaks (P ij , see Table 2 for indexes) In each step of the procedure only those peaks are processed which satisfy all conditions in previous steps.Thus, in every new iteration a smaller subset of peaks is generated.The procedure for the detection of regional peaks (and global in the additional sixth step) within the five-step procedure is presented in Figure 3.The steps are described below in greater detail.
(1) In the first step the local peaks are calculated by applying a local moving window with the kernel of size 3 × 3 cells, dedicated to raster-based spatial data analysis, e.g., [33].This step is based on topographic criteria.
(2) All of the local peaks on the flat areas (G p ) are eliminated setting criteria NOT G p .The removed local peaks are part of insignificantly convex areas.The classification of G p is part of a landform regionalization process applied to plains, low hills, hills, and mountains.This process applies a combination of slope, curvature and elevation of terrain [22].With an additional threshold those peaks which do not surpass a certain absolute minimum elevation (e.g., <600 m) can be eliminated.This step employs morphologic criteria in order to optimize the procedure significantly.
(3) The minimum horizontal distance (dL) is applied.Only those peaks which represent the highest within the circle of radius dL are retained.The adjacent peaks are therefore arbitrarily removed with dL > [150 m, 200 m].This step is based on the topographic criteria and follows the grain concept.
(4) The analysis of the relative height of the peaks (dH) serves as a basic topographic criterion.The algorithm determines the areas around the individual peaks that register as being up to the dH lower than the corresponding peaks.Next, the entire area is examined and the other peaks that were identified in the third step are counted.If no peak has been identified, then the examined peak is adopted.Starting with the highest peak, the procedure is repeated until the lowest peak has been processed.The peaks eliminated in the previous iterations of this procedure are not further examined.The condition for their removal is set to dH > [25 m, 30 m].There are further options to calculate the results, e.g., with contour lines in vector form [11,34] or by filling the pits of the inverse DEM in raster form.
(5) The final step is based on the morphologic criteria.A continuous potential surface of regional peaks of any shape (P 29 ) is applied (see Table 2 and the final step of the six-step procedure in the following section).The appropriate threshold setting for P 29 determines the elimination of several peaks on morphologically less rough areas.The progressive elimination of peak points can be applied with a continual increase in the surface threshold.The remaining points are finally adopted as regional peaks with respect to the topographic and morphologic criteria.
The five-step procedure runs twice by applying two slightly different parameters in steps 3 and 4. In the first execution the parameters are set to dL = 200 m and dH = 30 m, and in the second to dL = 150 m and dH = 25 m.These two calculations yield slightly different but complementary results and the minima of peak points of both are considered.The global peaks are derived using the same procedure (Table 2 and Figure 3/6).In this case, the parameter dL in step 3 is arbitrarily set to 2,000 m, and a potential surface of global peaks (P 15 ) is applied in step 5.

The Delineation of Peak Shapes by Morphologic Criteria
The application of the morphologic criterion allows for the delineation of peak shapes.However, some of the potential surfaces can also be used for a more advanced detection of peaks that has been primarily determined by various topographic criteria (see step five in the previous section).The shape analysis requires a database of peaks, in our case due to the peak detection section.The following sixstep procedure is proposed as a means of obtaining a shape analysis of peaks: (1) generating the variables through the operations of spatial analysis based on the raster format (2) relating the variables with two-digit codes ij (see Table 2) (3) converting the variables into binary indexes (I k ) (4) the further selection of the significant (descriptive) I k (5) combining the I k to generate a continuous potential surfaces (P ij ) ( 6) assigning the values of the P ij to the points of peaks (determined in the five-step procedure) The steps are described below in greater detail.
(1) The standard and more innovatively composed variables are generated in this initial step.Specific spatial patterns related to the shape of peaks may be discovered in these variables.The variables are generated through a raster-based spatial analysis, as provided by GIS-tools.The basic variables are slope and curvature [21,35].The standard variables that are based on both, a local moving window [21] of different sizes or on a local annular moving window (AMW) are: the minimum (MIN); the maximum (MAX); the mean (M); the standard deviation (STD); and the range (R; values from-to).The following innovative variables have been derived (see Table 4, Variable): -the "relief above" (RA = [MAX − (value at local window center)]) -the "relief below" (RB = [(value at focal center)-MIN]) -the "ridge-drainage" (RD = [MAX+MIN-2*(value at local window center)] = RA-RB) -the "rim" (RI = [MAX+MIN-2M]) -the "concave-convex" (CC = [(value at local window center)-M]) (as a kind of relative or local relief [5]) One of the most applicable variables that was developed is a relative relief based on a "multidirectional visibility index" (MVI).It is based on a visibility masks calculation with regard to the variations in azimuths and zenith angles [8].The catchment analysis (CA) is based on a calculation of cost surfaces around the peaks.
The variables are generated using DEMs of different resolutions and differently smoothed surfaces (see Table 4, Scale).
(2) Each variable is related to a single or multiple two-digit codes ij to determine the peak's size and shape (12 possible codes) according to Table 2.Each variable can therefore be related to more peak categories (see Table 4, Code (ij)).
(3) The variables as continuous functions are converted into binary (0/1) indexes I 1 , I 2 … I k by applying appropriate thresholds.One variable can be converted into more indexes (Figure 4).The threshold classification of variables to their Boolean forms I k is determined by the following primary rule: variables are sought that have values which characterize those areas that are least likely to reflect the target property (according to codes ij); however, they should simultaneously exclude as extensive a portion of the area as possible, using the upper limits that allow all pixels of the target I k to be retained.It is expected that the individual index I k should cover the target areas in their entirety.The settings of the threshold parameters may be derived through various approaches or combinations (e.g., fitting the I k areas according to a histogram analysis involving a certain portions of area, or to a certain number of detected/known peaks) based on iterative adjustment.In our case, the thresholds of every index I k are adjusted on the basis of a few preliminary selected typical peaks (according to codes ij), derived from the five-step procedure.Nevertheless, the quality of such results depends greatly on the operator's experiences.
(4) The indexes I k are further selected through independence and significance tests, if necessary.These tests can be processed using different statistical approaches (e.g., through regression analysis).In our case the I k are visually assessed by (a) comparing all sets of I k that reflect the same codes ij in order to recognize their independence, and by (b) assessing the portion of the whole area coverage to recognize their significance.Smaller portions reflect a higher level of significance.From among the many indexes, k = 35 are selected as significant (see Table 4, Index).However, the independence and significance of the indexes I k do not greatly influence the result of this procedure, but rather decrease the number of indexes I k .
(5) The aim of this step is to combine more I k in order to generate continuous potential surfaces (P ij ) (see Figure 4).The combinations of the I k can be preceded by fundamental analytical operations, such as Boolean algebra, and arithmetical and relational operations [36].In our case, we have applied a so-called graphical approach to our model which is based on the Boolean algebra of the indexes I k , which have been proved in some way to influence the location patterns of the target peaks ij.A simple sum of the indexes I k equally determines the location potential of the response data (final model), where the expected range is a number of the total included indexes I k .Figure 5 presents a potential surface for the determination of global peaks independent of the shape, with code ij = 19 (P 19 ) that comprises the following indexes: I 1 , I 14-17 , (optionally additionally included, but were not: I 20 , I 24-27 , I 29 as surfaces) (Table 4).The linkages between the shapes of P ij (category j; Table 2) are presented in Table 3.
Table 4.The variables, scale as the target resolution of the DEM, two-digit codes ij, and final k = 35 significant indexes I k .The indexes that are based on the AMW (twofold variable: as a local moving window and as an AMW) are only optionally involved in this procedure (6-SP), while those based on CA (this variable is not continuous) are not involved.Codes 14 and 24 (circular) are automatically part of codes 15 and 25 (conical), see Table 3   ), where a higher potential is signified by darker areas.
(6) In the final step the values of the P ij surface are assigned to the points of peaks determined in the five-step procedure.The grater is the value of the P ij in the surrounding of the particular peak, the higher is its significance (and potential).Additionally, an appropriate threshold can be applied to P ij to assign a binary decision significance/non-significance to the particular peaks.Figure 4 schematically illustrates this procedure.
There are alternative methods of analyzing the shape of peaks using the DEM (Table 4).We will consider two alternatives involving just the first step of the six-step procedure.These are based on the AMW technique and the CA variable.The applications of the AMW are suitable for the detection of conical peaks (codes 15, 25 and consequently 14, 24; Table 2, and 1st figure in Section 4.3).An independent procedure in two steps is proposed: (a) The AMW is defined by two radii (the inner r 1 and the outer r 2 ) measured in the number of cells along the x-(=row) and y-direction (=column) of the matrix.Since the grid spacing is usually specified in real world units, the radii of the AMW may also be expressed in meters or in any other standard unit (Table 5).All cells located within the ring defined by these two radii are included in the calculation.(This step generates a continuous variable that might be used through the six-step procedure for the conical peaks delineation).(b) The results are measured in the second part of this procedure.The annular-shaped ring (ASR) is constructed with the same dimensions as the AMW used in step (a).It is centred on the points of the peaks.The average values of the variables (MIN, MAX, MEAN) are calculated for the ASR and assigned to the points of the peaks.These values are then compared to the values of the variables on the points of the peaks or compared to the heights of the peaks.The most significant results are calculated with the variable MAX.The second approach is based on the CA, particularly on cost surfaces [2].The cost surfaces are computed using the points of the peaks, where the slope is assigned as a cost.Therefore, this procedure may not be used in the six-step procedure.A geometric analysis of the shapes and sizes of the areas generated around the peaks in certain cost distances indicates the shapes of the peaks, especially the oblongness and the sharpness (Figure 6).The cost surfaces can be calculated to different scales using the DEMs of various resolutions.Utilizing different scales, the general shape of the whole mountain range can be assessed hierarchically starting with the small areas around the peaks and increasing to the areas at the foot of the mountains [29].

Results and Interpretation of the Case Study
The Kamnik Alps will serve as our case study area.It is located in the northern part of Slovenia along the border with Austria (Figure 7).The area is partly karstic.The dataset selected for the modeling of the peaks is morphologically high quality DEMs of Slovenia produced from various data sources with a resolution of 12.5 m and resampled to 25 and 100 m respectively [22].Various DEM resolutions are suitable for different scale analyses of peak shapes.The reference data is from maps DTK 25, and geographical names REZI 25 [37].Additional reference data is drawn from two independent databases of 2000ers in Slovenia [19,38].

Result of Detected Peak Heights and Their Verification with Lists of Peaks
The results of the proposed automated five-step procedure (Figures 3 and 8) are that we can detect peaks, starting with local (37,218), via immediate (3493, 2529, 366), to regional (346), and additionally to global (3).
Regional peaks over 2,000 m (a.s.l.) are documented for the entire Slovenian Alps and ranges between 179 and 349 according to different lists.We found two fully documented lists [19,38].
For the evaluation we limited the problem to the case study area and those peaks over 2,000 m (a.s.l.).The list of 42 peaks was reclassified from the 346 numerically determined regional peaks, which we refined and revised using the five-step procedure (Figure 8).Let us call it List 0. We will call the two reclassified verification lists List I [38] of 41 peaks and List II [39] of 48 peaks.The third verification dataset will be called List W, which is generated with a morphometric method developed by Wood [29].Wood's method applies a topographic criterion called "minimum drop surrounding peak", which is similar to our relative elevation (dH), set to 30 m.
The List I is presented in the upper section of Table 6.Verification clearly shows that 10 of the peaks from list 0 highlighted in bold italic letters do not exist in List 1. Below, all lists are compared with a kind of confusion matrix, where the diagonal values represent the number of peaks in individual lists, and the other cell values represent the number of omitted peaks.Approximately 1/4 of all peaks from the different lists do not correlate across the lists, however the majority of prominent peaks are identifiable in all lists.
The results presented are comparable to the manually defined peaks in the whole Alpine region for peaks over 4,000 m, and in the Himalayas for peaks higher than 8,000 m.There are 128 peaks higher than 4,000 m in the entire Alpine region, but only 82 when we consider the stricter criteria [2].For comparison, the oldest list from 1923 contains only 59 peaks [40].The number of peaks in the Himalayas above 8,000 m ranges from 14 to 21 [41].
Figure 8.The number of peak points for a case study area: from local (1) to regional peaks (5), and an additional attempt to attain global peaks (6), calculated with reference to the topographic and geomorphic criteria.Table 6.Upper: List I of 41 peaks (described with names and heights in meters) over 2,000 m (a.s.l.) [38] in the case study area.Peaks highlighted in bold italic letters do not feature in our results, i.e., the List 0. Below: The correlation between the Lists I, II, W and 0. The arrows indicate the direction to read the data.

Verification of the Heights and Positions of the Peaks
Lists I [38] and II [39] contain peak heights, three of which differ by up to 66 m (a possible gross error).A vectorized topographic map to scale 1:25,000 (DTK 25) was used to verify the results of automatic peak detection [22,37].The triangulated spot heights of the DTK 25 are common on the highest peaks and with the highest precision.These heights are comparable to the heights featured on Lists I and II.As can be seen in Figure 9, most heights derived from the DEM with a resolution of 12.5 m are usually a few meters lower than the triangulated spot heights of the DTK 25 (with the exception of a few gross errors).The reason for this lies in the generalization and the limited resolution of the DEM where the grid cell values are slightly biased towards the points of height of the peaks.
The positional accuracy of the determined peaks is limited by the resolution of the DEM.However, several peaks are dislocated by up to several 100 m (a gross error) comparing the spot heights of the DTK 25 and the peaks derived on the basis of the DEM.The gross errors of the heights and the positions of the peaks were not investigated.These errors might have sources in all lists.

Figure 9.
Assessment of the peaks generated by the automated method (dots with altitude) (List 0), compared to those peaks from DTK 25 (triangles with altitude in brackets) and also compared to Wood's method [29] (crosses) (List W), and labeled with geographical names from REZI 25.A zoom-in area (see Figure 7) with dimensions of 8 km × 6 km is presented.Figure 10.Regional peaks (points) and global peaks (triangular).Additional labels denote the shape of the peaks: crosslet: sharp peak; empty circle: blunt peak; rhombus: oblong peak (a full rhombus denotes a greater degree of oblongness).It can be seen that the sharp and blunt shapes of the peaks cancel each other out, while the oblong shapes can add to each other.Some peaks have not been classified, as their shape is not clear.

The Result of the Delineated Peak Shapes
The result of the six-step procedure is a determination of peak shapes (Tables 2 and 3, Figures 4-6).52 sharp (code 21), 153 blunt (code 22) and dissimilar oblong peaks (code 23) were detected from among the 346 regional peaks.Three global peaks were also defined (code 19).The results are presented in Figure 10.
Figure 11 shows a conical formation of middling scale, with a radius r 2 = 500 m from the points of the peaks (see Table 5).The size of the black circles is determined by the difference between: (a) the average value of the variable MAX on annular-shaped ring (ASR); and (b) the value of the variable MAX at the peak point.The size of the white circles is defined by the standard deviation of the variable MAX on ASR.Both the sizes of the circles are proportional to the values of variables.
Figure 11.The extent of conical formation in regional peaks.Large black circles with small white dots denote evident "conicalness", whilst smaller grey circles denote some degree of conical form, and the white circles of different sizes denote the peaks that are not conical.

Verification of Delineated Shapes Using Peak Toponyms
The last is toponomastical or, to be more precise, an etymological comparison between the shapes of the derived regional peaks and the toponyms from the register of geographical names (REZI 25) [37].
Our assumption was that the shapes of peaks, coded in Table 2 and determined in the study, would be reflected in the toponyms of REZI 25.The main ontological key for semantic approval lies in dictionaries, e.g., [10] and other thematic etymological publications, e.g., [1,42] for the Slovenian language.The peak's toponyms have also been related to their visual interpretation and aesthetic judgments.This comparison and interpretation is related to the mountaineering criteria for peak delineation.
The instantly recognizable examples from the case study area are: Ojstrica (which means a sharp peak); Kopa, Tolsti vrh (both blunt peaks); Planjava (meaning a bare plain peak), Storžič (a conical or cone peak); or, more indirectly, Grintovec (a crumbly peak) (Figure 9).It can be seen from a visual interpretation of Figure 5 that the most significant peaks (the darkest areas) are: Storžič (Figure 12), Kozji vrh, Grintovec, Goli vrh, Skuta, Brana, and Ojstrica.It can be seen from Figure 10 that all of those peaks are also sharp.Furthermore, three global peaks are identified: Grintovec, Planjava, and Raduha (Figure 3(f)).A high degree of conical formation among the peaks of the entire Kamnik Alps is reflected within the toponyms Storžič and Kärntner Storschitz (Pristovški Storžič).Both toponyms can be confirmed using the automated method based on AMW.In addition, the following peaks can be identified as being conical: Kozji vrh, Goli vrh, Grintovec, and Ojstrica (Figure 11).There are some peaks whose shapes are unclassified, and their meaning does not fit with any of the coded expressions in Table 2.These are marked only with black dots in Figure 10.A good example is Kočna, whose three peaks, all within a distance of 500 m, are visible from a great distance-from the Slovenian capital Ljubljana in the south, and from Jezersko near to the Austrian border in the north.These three peaks are geographically significant, but no one criterion developed in this study (sharp, blunt, oblong, conical) can be used to describe their actual shape.The different lists also record variations in the measurements of heights and positions.Such examples present a challenge to discovering and optimizing the algorithms and procedures required for a standardized peak detection and delineation method.
A cross-analysis of a range of interpretations confirms that Grintovec is considered to be a particularly important peak, considering topographic, morphologic and mountaineering criteria.It is also the highest point in the Kamnik Alps.On the other hand, it could be considered to be less important since it is geologically crumbly and is, according to typical views, quite overshadowed by some of the other smaller mountains.This example illustrates the highly complex problem of morphometrically interpreting peaks.Other interesting examples from the results are the plateaux Velika Planina and Menina Planina, both of which are relatively significant on the potential surface of the global peaks (P 19 ; Figure 4), but the visualization in Figure 10 shows that they are in fact high plains with clusters of several small, blunt peaks.

Discussions
The method of peak detection is based on the parameterization of topographic and morphologic criteria, and on the mountaineering criterion of the International Mountaineering and Climbing Federation (UIAA) [2].The proposed five-step procedure relied on the application of topographic and morphologic criteria, i.e., the morphometric measurements of relative heights and distances between the peaks with the support of potential surfaces.The quality of algorithms developed was confirmed by comparing them with the official lists of all peaks over 2000 m (a.s.l.) [38][39][40][41] that are human-based lists, thus less objective, and also by comparison with Wood's algorithm [29].The degree of variability within the official lists is similar to the differences between the official lists and our results.
The analysis of the shapes of the peaks is based on a method of shape categorization and the modeling of morphologic criteria that were proposed and applied in this study.A six-step procedure (with two exceptions) to generate the potential surfaces was developed for sharp, blunt, oblong, and circular peaks.The high quality of the results was proved by cross-checking against the topographic maps [10,37] and by etymologically comparing the meanings of the toponyms with their shapes [1,42].
The selected study area is relatively geomorphologically diverse with many relief types, and with the availability of different datasets and other information that supports the study.These circumstances have enabled a more integrative study of the topographic peak problem.The methods have been developed within the GIS environment and tested using morphologically high quality DEM across the case study area.The object-based image analysis is promising set of techniques for multi-scale peak detection and shape analysis, especially when using a DEM based LIDAR datasets [43].
The possible applications of this research are numerous, for example (some of which have already been applied): to improve DEM production [22]; the comprehensive quality control of DEMs; enhanced automatic generalization of DEM [44]; landscape classification; an increased knowledge of multi-scale morphology [45], applicable to geological and tectonic studies; karst research; generic prominence analysis [46,47]; silhouette generation; the improved cartographic presentation of terrain [8]; examining anthropogenic influences on the landscape morphology; and ethnophysiographic, sociolinguistic, cognitive psychologic, psycholinguistic, or ergonomic [11] studies.

Conclusions
This study demonstrates the strength, reliability and robustness of the relatively simple autometric procedures that have been developed and integrated into innovative procedures and algorithms in order to resolve a complex task: the task of detecting (regional) topographic (mountain) peaks with a delineation of their shape.The method of peak detection is based on the parameterization of topographic and morphologic criteria, and on the mountaineering criterion.The analysis of the shapes of the peaks is based on a method of shape categorization and the modeling of morphologic criteria that were proposed and applied in this study.
Compared to previous studies, the proposed procedures have been designed to adopt a holistic view in order to arrive at a reliable conception of topographic peak.These procedures include inventive numerical morphometric image processing or spatial analysis techniques, which have been integratively enriched with alternative algorithms, independent peak lists, topographic maps, triangulated spot heights, etymological research, visualizations of different slant views, geological interpretation, and even with aesthetic judgments.
The results confirm the efficacy of parameterization and numerical application, and represent an improvement in the standards of peak detection.They demonstrate the capability to analyze complex shapes and to perform a more complex terrain analysis (i.e., the determination of sinkholes, or entire mountains in the context of the peaks).The analysis of shapes also enables the improved detection of peaks.The advantages of the automated approaches are the highly standardized and universal parameters (at least for the topographic peak detection), and the comparative results yielding a higher overall quality.However, the automated methods will never fully satisfy the requirements that stem from the complexity of nature.Some details will have to be discovered and solved by human beings using our other senses rather than eyesight only.This problem is related to the problem of generalization in automatic cartography where a reliable solution can be solved in a semantic context and indicated by ontological principles.
A subsequent study could compare and reconsider the results of the proposed approaches with the number of available lists of peaks, e.g., >1,000 m in Quebec >2,000 m in Slovenia, or >4,000 m in the Alps.Future challenges would lie in interdisciplinary research of various peaks properties for numerous applications, using high-resolution point clouds or the DEMs based on LIDAR datasets, discovering applicable criteria and variables, and optimizing the algorithms and automated procedures.

Figure 1 .
Figure 1.Peaks at different times of day, in different weather, and differently illuminated.An image processing technique of multidirectional visibility index (MVI)[8] shows another aspect in appearance of peaks and other terrain features (a composition of series of 57 photograph visibility masks between 6:30 AM and 17:15 PM (GMT+1), 26 September; view from Kredarica to peak Tosc; webcam of ARSO).

Figure 2 .
Figure 2. The evident discrepancy in the position of peaks between old (hatched) and modern (contour lines) maps (for the plateau Pokljuka).

Figure 3 .
Figure 3.The detection (identification) of the regional peaks (black dots) according to the five-step procedure.Picture (f) presents an additional attempt at detecting global peaks.

Figure 5 .
Figure 5.The potential surface for the determination of global peaks (P 19), where a higher potential is signified by darker areas.

Figure 6 .
Figure 6.CA based on the slope of the DEM, calculated for a cost distance, equivalent to 100 m form the peaks (white crosslets).The black areas are nearly circular, the grey are oblong (1:2) and white are oblong to a greater degree (1:3) Smaller areas denote sharper or rougher peaks.A zoomed-in area is present (see Figure7) with dimensions of 13 km × 9.5 km.

Figure 7 .
Figure 7.The case study area: the Kamnik Alps in Slovenia, with a primary area of dimensions 30 km × 20 km (dark grey) that reflect the areas for the Figures 3, 4, 10 and 11.The two additional frames illustrate zoom-in areas for the Figures 6 (13 km × 9.5 km), and 9 (8 km × 6 km).The three black dots in the primary area represent the detected global peaks.

Figure 12 .
Figure12.A view of the conical peak of Storžič as seen from Kokrsko sedlo, 10 km east from the peak, with the corresponding insertion of the topographic map DTK 25 below with 1 km gridlines[37].The peak appears very conical from all sides. .

Table 5 .
Dimensions of annular-shaped moving window (AMW) according to different DEM resolutions.