dems ( digital elevation models ) are used for the extraction of land - surface parameters and objects through geomorphometric analysis ( pike , 2000 ; pike et al . , 2009 ) . as general geomorphometry applies to continuous land surface ,  specific geomorphometry applies to discrete landforms ( evans , 1972 ) . most land - surface parameters and objects vary with spatial scale , which in the digital realm is widely understood as a function of cell size or grid resolution ( wilson and gallant , 2000 ) . however , grid resolution is not a particularly appropriate representation of scale ( gallant and hutchinson , 1996 ) . the dependence of land - surface parameters on grid resolution has been described by evans ( 1972 ) as  a basic problem in geomorphometry ( shary et al . , 2002 ) . since most geomorphometric algorithms work through a fixed neighborhood operation ( pike et al . , 2009)usually within a 3    3 kernel  the scale of analysis is tied to the resolution of the input dem and changes as the resolution changes ( zhu et al . , 2008 ) . in the absence of scale optimization techniques ( li , 2008 ) , the geomorphometric analysis is conducted at rather arbitrary scales , which rely on the user 's experience . in fact , the scale of analysis often depends on data availability as many users perform what schmidt and andrew ( 2005 ) called the  let 's take a dem  approach , without much concern for scale effects in analysis . obviously there is no sound scientific justification for a direct linkage between natural phenomena and data acquisition techniques ( strobl , 2008 ) . lately , a large amount of literature has developed particularly relevant to applications of geomorphometric analysis ( terrain - based environmental modeling ) , underpinning the impact of scale mismatches between the target variable and the explanatory ones leading to statistical bias . with the advent of increasingly high resolution dems , scale is becoming an important issue in geomorphometry ( macmillan et al . , 2003 ) . particularly when mapping landforms from gridded land - surface models such as dems , scale - related shortcomings are partly generated by the implicit assimilation of data model elements ( cells or pixels ) to geographic objects ( fisher , 1997 ) . thus , in a somehow counter - intuitive manner , most landform classification systems work through the classification of cells , which are clustered to define the extents of objects  instead of delineating the objects first and then classifying them . the approach of classifying cells directly is limited in several aspects , including the scattered aspect of classification in the so - called  salt - and - pepper effect , tying the scale of analysis to the raster resolution , difficulties in including topological relationships in classification and also in developing hierarchies of landforms . although delineation techniques have been proposed , quantitative evaluations of their performance have only recently been done ( van niekerk , 2010 ) . still , fundamental questions pertinent to the nature of objects and their delineation relative to scale have only been touched upon . minar and evans ( 2008 ) provided a useful review on object ontology , showing the limitations of cell - classification approach in mapping landforms and introducing segmentation of elementary forms as an alternative . deng ( 2007 ) and hengl and reuter ( 2009 ) reviewed contributions pertinent to scale and object representation from dems , but rather in wider contexts . this paper provides an overview and prospects for object representation from dems as a function of scale . in the next section , relationships between object delineation and classification or regionalization are explored , in the context of differences between general and specific geomorphometry . rather than being an exhaustive overview , this section will complement minar and evans ( 2008 ) with a focus on generating objects with the aid of a multi - resolution segmentation ( mrs ) algorithm . a review of scales issues in geomorphometry  ranging from scale effects to scale optimization techniques is followed by an analysis of pros and cons of using cells and objects in dem analysis , in section  3 . this section focuses on methods to establish non - arbitrary scales in the land surface . prospects for coupling multi - scale analysis and object delineation are then discussed , in section  4 . specific geomorphometry applies to discrete spatial features or geomorphic objects ( evans , 1972 ) . landforms and landform elements ( see macmillan and shary , 2009 for a detailed review ) are particular cases of geomorphic objects . extracting landform divisions ( landform elements , landforms , etc . ) from dems usually involves applying an object model on a raster data structure , which is the most popular format for spatial modeling ( pike et al . , 2009 ) . this is because raster elements themselves , i.e. cells or pixels , do not have any meaning in reality regardless of their size ( fisher , 1997 ) . cells as artificial , discrete units exist merely for the purpose of representation ( goodchild et al . , 2007 ) . clearly , footprints of cells have nothing to do with the size and shape of real - world entities such as landforms . therefore , a model to translate from the continuous land surface to discrete entities is required . this is challenging as the land surface is smoother than other terrain variables , such as vegetation ( hengl and evans , 2009 ) . consequently land - surface objects are less obvious on dems than land cover patches on satellite images , for instance . in the absence of comprehensive conceptual data models for land surfaces ( brndli , 1996 ) , a landform is rather represented as a collection of cells that exhibit similar morphometric characteristics ( schmidt and dikau , 1999 ) . although this strategy is straightforward two major shortcomings have been identified : firstly , cells are often treated as spatially independent from each other , and thus adjacent cells are frequently assigned to different classes resulting in a highly scattered spatial representation of landforms ( burrough et al . , 2001 ) , the so called  salt - and - pepper effect. secondly , since landforms are classified on a cell by cell basis after applying statistical rules , they are solely defined thematically , but not spatially ( e.g. by location , context or topology ) and hence , they do not represent spatially configured objects ( deng , 2007 ; minar and evans , 2008 ) . often , in a kind of procrustean bed approach , cells are allocated to classes following pre - defined categories of , or thresholds in , land - surface parameters ( e.g. the threshold value of slope to define a flat surface ) . landform boundaries are then given by the edges of the aggregated cells ( as resulting after some filtering , needed to reduce the  salt - and - pepper effect ) . but these boundaries may not coincide with morphologic discontinuities in a given landscape ; they are merely conceptual or fiat boundaries ( smith , 1995 ) . this problem of arbitrary incidence has been exemplified for curvatures by minar and evans ( 2008 ) . the authors observed that  isoline boundaries may create artificial areas without sufficient respect to the natural structure of landforms with various types of homogeneity ( p. 241 ) . such crude representation of landscapes is part of  the conceptual and computational gap between local geometry and meaningful landforms , which broadens paradoxically with improving quality and resolution of dems ( mark , 2009 ) . moving on from collections of geomorphometric objects ( schmidt and dikau , 1999 ) requires delineating the objects first , then classifying them . field ( cova and goodchild , 2002 ) , dems can be partitioned into discrete , spatially intact land - surface objects , following data - driven approaches rather than pre - defined classification templates . a strategy of clustering similar cells in property space by means of image analysis methods to delimit form types has already been envisioned by pike ( 1995 ) and applied by irvin et al . while this method produces less scattered objects , the problem of matching land - surface discontinuities still remains , since clusters are created using global thresholds instead of local contrasts ( van niekerk , 2010 ) . the same applies to classification methods using dynamic but global thresholds ( iwahashi and pike , 2007 ) . identification of spatial discontinuities in land - surface parameters seems to be more appropriate for object delineation . this idea was presented by minar and evans ( 2008 ) as an axiom :  at a given scale , the land surface may  exhibit discontinuities ; these may be recognized as natural boundaries of geomorphic objects  . a manual technique of mapping based on morphological discontinuities was proposed by savigear ( 1965 ) . ( 1995 ) described an algorithm for automated mapping of land components through approximation of slope breaks . recently , image segmentation techniques have increasingly been used to generate objects based on the concept of heterogeneity . the most known algorithm is mrs ( baatz and schpe , 2000 ) as implemented in the ecognition software . this is a region - merging technique to create objects from pixels through an optimization process that minimizes the internal weighted heterogeneity of each object at a given scale . these objects are then merged or split to create objects at consecutive scales , either higher , created in a bottom - up approach , or lower , created in a top - down one . therefore , each decision of merging or splitting is based on the attributes of homogeneous structures of a recent scale ( baatz and schpe , 2000 ) and on the user - defined heterogeneity threshold , called scale parameter . drgu and blaschke ( 2006 ) and van asselen and seijmonsbergen ( 2006 ) introduced this algorithm to the analysis of dems . this approach has lately been increasingly used in delineation of landforms or land entities ( drgu and blaschke , 2008 ; mller et al . , 2008 ; , 2008 ; anders et al . , 2009 ; blanco et al . , 2009 ; kringer et al . , 2009 ; while segmentation relies on local contrasts in drawing meaningful boundaries of objects , classification uses global thresholds to facilitate interpretation of landform classes . van niekerk ( 2010 ) recently found that an mrs algorithm is more sensitive to morphological discontinuities than two other alternatives , alcom and isodata , and he proposed it as the most suitable technique for delineating land components from dems . an image segmentation approach should not be confused with landform segmentation as used by pennock and corre ( 2001 ) , which is actually a classification procedure based on pre - defined morphometric categories . other segmentation algorithms have also been successfully applied for geomorphological mapping or delineation of homogeneous areas from dems . for instance , miliaresis ( 2001a , b , 2006 ) used a region - growing algorithm for extraction of bajadas from dems and satellite imagery , and for geomorphometric mapping at regional scales ; lucieer and stein ( 2005 ) proposed a region - growing segmentation procedure based on texture to extract landform objects from lidar data ; stepinski et al . ( 2006 , 2007 ) , stepinski and bagaria ( 2009 ) , and ghosh et al . ( 2010 ) developed a segmentation approach combined with artificial intelligence to automatically map planetary surfaces ; jellema et al . a particularly simple and appealing procedure for watershed segmentation of curvature was proposed by romstad and etzelmller ( 2009 ) for the purpose of geomorphological mapping . due to increasing availability and easier access to dems at a broad range of spatial resolutions ( from lidar at several centimeters up to gtopo30 at approx . 1  km ) , multi - scale analysis of the land surface is becoming more feasible . the modeling of scale effects with respect to both changing resolution and varying window size for surface calculations has been identified as a major research topic not only in geomorphometry , but rather in all disciplines dealing with dems including hydrology , soil science , and geomorphology . in a recent paper li ( 2008 ) provided a valuable review on the numerous approaches that examine scale dependencies in terrain - based modeling , and has outlined a general strategy for their analysis . as has been shown parameters such as slope ( deng et al . , 2008 ) and curvatures ( schmidt and andrew , 2005 ) , do not only change in magnitude , but may even shift their topographic meaning , as for sign of curvature . moreover , scale dependencies in terrain analysis are driven by terrain characteristics ( e.g. simple or complex , carter , 1992 ) , and also vary across different landform types ( gao , 1997 ; deng et al . , 2007 ) . an important development in examining scale effects is that scale has become an integrated part of terrain analysis ( deng , 2007 ) . for example , wood ( 1996 , 1998 ) proposed calculating surface parameters at various window sizes for a constant resolution , and for each cell recording the results as a series of values linked with scale information . his open - source software package landserf for  multi - scale surface characterization offers powerful visualization tools for exploring scale effects . following wood 's approach several researchers focused on exploring the effects of neighborhood size on computed terrain parameters as well as on the application of various window sizes for multi - scale characterization of landforms ( fisher et al . , 2004 ; schmidt and hewitt , 2004 ; schmidt and andrew , 2005 ; reuter et al . , 2006 ; lately , it has been found that parameters are less sensitive to dem resolution changes than to variations in neighborhood size ( zhu et al . pure modeling of scale effects barely gives clues on how to select non - arbitrary scales for given analyses . they conducted either experimental testing in the context of terrain - based environmental modeling , or theoretical analysis in data - driven approaches ( for a review see li , 2008 ) . in terrain - based environmental modeling it is essential to fit the spatial scale of the terrain data to the scale of the processes or features under investigation in order to obtain valid model results . in doing so , several authors compared results from multi - resolution terrain analysis with reference data such as field measurements of a landscape property ( bian and walsh , 1993 ; florinsky and kuryakova , 2000 ; hengl , 2006 ; drgu et al . , 2009a ) , or model outputs ( smith et al . , 2006 they were able to identify the grid size or resolution range with the most powerful predictions . however , the optimal grid size might be different for different target variables , so that one can hardly select a single optimal resolution ( hengl , 2006 ) . especially in soil - landscape modeling , recent efforts have been made towards optimizing neighborhood sizes for the prediction of soil classes ( smith et al . , 2006 ; zhu et al . , 2008 ; behrens et al . , 2010 ) . indeed , it is more desirable to develop methods that work data - driven , and without reference to dependent variables outside the dems . depending on the method one obtains a measure for each scale , and when plotting all these measures or simply the results from terrain analysis against scale , one finally gets a schmidt and andrew ( 2005 ) introduced a spatially adaptive scale detection technique exemplified for curvatures in order to recognize dominant scale ranges of landforms . gallant and dowling ( 2003 ) proposed an algorithm to produce a multiresolution index of valley bottoms based on their topographic signatures at multiple scales . a particularly appealing concept to describe aspects of scale dependency in spatial objects is fractals ( mandelbrot , 1975 ) . however , empirical evidence ( chase , 1992 ; evans and mcclean , 1995 ; perron et al . , 2008 ) suggests that fractal models are not appropriate to the land surface , which has a statistically multidimensional character ( evans , 1998 ) . tate and wood ( 2001 ) provide a comprehensive review on fractals and scale dependencies in the land surface . probably the most promising approach for data - driven scale detection is the method of local variance ( lv ) . this method was originally developed in image analysis for the purpose of scale detection ( woodcock and strahler , 1987 ) . the approach is based on the relationship between the size of objects in the real world and pixel resolution , as expressed by the spatial structure of images . thus , in a high - resolution scene ( strahler et al . , 1986 ) , objects in the real world are represented by multiple pixels , hence spatial autocorrelation is high ( fig .  1 , top ) . local variance , computed as average value of standard deviation measured in a small neighborhood ( 3    3 moving window ) , is therefore small . when successive coarser scales are produced from the initial dataset through resampling , local variance increases with the scale levels up to the point where pixels start approximating the representative objects in the scene ( fig .  1 , top ) . at this scale level the maximum value of local variance is recorded as the likelihood of neighbors being similar decreases . at coarser pixel sizes local variance decreases again as a consequence of including more objects within a pixel , hence the spectral difference between neighbor pixels is reduced ( fig .  1 , top ) . despite its simplicity and usefulness this method was not widely adopted in remote sensing and gis ( cao and lam , 1997 ) . li ( 2008 ) suggested the lv method could prove useful as a scale detection technique in dem analysis . two recent studies tested the suitability of the lv method for multi - scale pattern analysis in geomorphometry ( drgu et al . scale levels were simulated from the same high - resolution datasets through resampling and image segmentation respectively , in a bottom - up approach . the authors found that the lv method performed better when scale levels were created with image segmentation as compared to cell aggregation . up - scaling through resampling induces an isotropic smoothing that levels out spatial patterns at coarser resolutions . in contrast , image segmentation maintains distinct boundaries of objects ( karl and maurer , 2010 ) by steadily adding heterogeneity to objects at each coarser scale level ( fig .  1 , bottom ) . thus , at each coarser scale level contiguous objects with similar property values are merged into larger ones . talking about the specific geomorphometry , mark ( 1975 , p. 165 ) observed that  the specific approach can only be applied once an area has been identified as a drainage basin , an alluvial fan , a drumlin , etc.. this means that delineation of objects that satisfy the condition of maximizing internal homogeneity and external differences is not a specific approach , as long as the objects do not bear a meaning other than statistical . as minar and evans ( 2008 , p. 238239 ) pointed out  segmentation of the land surface can provide a transition from the field model to the object model , and from general geomorphometry to specific geomorphometry , thus connecting the continuous and atomistic hypotheses. the main issue is that this transition has not been conceptualized in geomorphometry so far . discrete geomorphometry would apply to and describe land - surface divisions defined solely by the criteria of homogeneity in respect to a given land - surface parameter or a combination of several parameters . homogeneity , in its turn , should always be relative to scale . the main aim of discrete geomorphometry is to produce morphometrically meaningful objects . this approach has a general and objective character as general geomorphometry ; however the two approaches are framed into different conceptual models field vs. object . on the other hand , discrete geomorphometry shares the object model with specific geomorphometry ; however , in the specific approach the character of objects is defined before segmentation , while in the discrete the meanings are assigned after or along with segmentation . discrete geomorphometry thus centers on objects , defined as homogeneous areas delineated by discontinuities in land - surface parameters , either on individual or combined layers , which reveal the land - surface patterns at a given scale or across scales . therefore , spatial pattern , which  has been missing from most quantitative work ( evans , in press ) would come into focus . possible candidates for such objects have been called land components by dymond et al . ( 1995 ) , terrain facets by rowbotham and dudycha ( 1998 ) , elementary forms by minar and evans ( 2008 ) , pattern elements by drgu et al . these objects can be seen as intermediate building blocks in the translation from cells to landform divisions . once the objects are delineated in the digital realm , further statistical , relational and semantic rules can be applied to map each object to the landform concept to which it comes closest ( macmillan et al . , 2004 ; minar and evans , 2008 ; bishop , 2009 ; eisank , 2010 ) , by incorporating expert knowledge ( macmillan et al . , 2005 ) , or to use terrain objects as basic areal divisions for the study of land - surface processes ( rowbotham and dudycha , 1998 ) . once the objects are assigned to classes of elementary form , or given interpretations such as  terrace , segmentation is a good candidate as the main method of discrete geomorphometry . in the second section of this paper we presented its advantages over clustering in property space and cell classification ( also see minar and evans , 2008 for a comprehensive discussion on segmentation techniques ) . two main strategies of segmentation can be followed : segmentation into pre - defined types and degrees of homogeneity ( minar and evans , 2008 ) and data - driven approach ( drgu et al . constant values of altitude or its derived morphometric properties ( minar and evans , 2008 , p. 244 ) . a unified system of elementary forms was predefined based on variation in the general fitted function . the forms are data - based and assigned to a form class after or along with delimitation . this approach is potentially likely to facilitate the reference of terrain divisions to process or genesis . although the approach should be in principle independent of scale , the authors acknowledge that many aspects of the land surface are scale dependent ( evans , 2003 , 2009 , 2010 ) , therefore further investigation is needed before applying it at broader scales ( minar and evans , 2008 ) . the data - driven approach was developed with the aid of the multiresolution segmentation algorithm ( see section  2 of this paper ) . selection of the scale parameter was turned into an objective choice , with the aid of the esp tool ( estimation of scale parameters ) , by using the concept of local variance ( drgu et al . , 2010 ) . segmentation of land - surface parameters such as slope gradient ( eisank and drgu , 2010 ) was performed with the esp tool in a bottom - up approach , where objects at finer scale were steadily merged into more heterogeneous objects at broader scales ( fig .  1 ) , by changing the scale parameter in a constant increment . was measured by first calculating the standard deviation of objects , and then averaging object values to obtain the pattern mean . values of lv were plotted against scale parameter ; breaks in the lv graph and its rate of change indicated the scales where the probability that delineated land - surface objects match a group of similar - sized real - world forms should be higher than for other scales ( drgu et al . although the two methods are in incipient stages , they offer good prospects for delineating homogeneous morphometric primitives either independent of scale , or at multiple scales , through multi - scale pattern analysis . however , much work lies ahead before they can be fully operational for what olaya ( 2009 ) called the  discrete analysis of the land surface. object ontology needs particular attention . for instance , more research on relationships between real land - surface features and segmented objects across scales is required . this is to make sure that we delineate real objects and do not create artificial constructs ; arbitrary partitions at various scales would represent a particular case of the modifiable area unit problem ( maup ) with all its acknowledged shortcomings ( openshaw and taylor , 1979 ) . a possible strategy for linking object ontologies of segmented objects with concepts of real landforms is semantic modeling ( eisank et al . , since land - surface features show smooth transitions , how can we account for fuzziness on both conceptual and spatial domains ? preliminary results ( drgu et al . , in press ) show that discontinuities in land - surface parameters seem to relate to the quality of object boundaries : sharper contrasts are expressed by smoother lines , while soft transitions give more indented boundaries . if this proves true , the degree of fuzziness could be measured via edge analysis . models like core vs. transitional areas ( drgu and blaschke , 2008 ) , attractors ( minar and evans , 2008 ) , or spatial gradation of slope positions ( qin et al . , 2009 ) would be further useful in classification of objects . here we presented methods producing non - overlapping objects , potentially amenable to nested hierarchies . however , overlapping areas might be better suited for given applications . romstad and etzelmller ( 2009 ) proposed an interesting approach based on overlapping convexities and concavities . this paper has provided a critical review on object representation from digital elevation models . translation from continuous models of elevation and its derivatives ( slope , aspect , and curvatures ) to landform divisions ( landforms and landform elements ) depending on the dem resolution and on the size of the analysis window , land - surface parameters have different values at the same location ; consequently one landscape can be represented in multiple ways . thus , scale impacts heavily on the results of geomorphometric analysis . whereas scale effects on geomorphometric analysis are now relatively well understood , scale optimization techniques remain a priority for future research . scale becomes more important with increasing dem resolution ( macmillan et al . , 2003 ) : while the number of landscape representations from a low resolution dataset is limited , many more versions of the landscape can be represented from a single very high resolution dem through up - scaling . which of those representations are best suited for a given purpose ? the answer to this question prompts a reliance either on expert knowledge ( gustavsson and kolstrup , 2009 ) , or an exploratory attitude in the use of local land - surface parameters ( deng , 2007 ) . an exploratory attitude is unfortunately hindered by poor technical implementation of scale issues in most gis software packages . landserf ( wood , 1996 ) is a remarkable exception , providing a suite of solutions for scale / multi - scale analysis ( wood , 2009 ) . although landform divisions are discrete features by definition ( evans , 1972 ; macmillan and shary , 2009 ; pike et al . , 2009 ) , most of the procedures to produce them in the digital era have essentially been tributary to a field model : individual cells are allocated to pre - defined classification schemes . objects emerge then as aggregations of cells through decisions that are not easily applicable to other landscapes . boundaries of such aggregates are likely to fail in matching land - surface discontinuities . image segmentation procedures create the technical framework for delineating homogeneous objects delimited by real discontinuities in land - surface parameters . abstracting complex surfaces as objects is innately human ( gessler et al . , 2009 ) , so that this technique can bridge the gap between a pixel - based approach and manual mapping based on visual interpretation , while preserving the advantages of speed and objectivity given by computers . discrete geomorphometry would apply to and describe land - surface divisions defined solely by the criteria of homogeneity in respect to a given land - surface parameter or a combination of several parameters . the emerging idea of smoothing local variability in land - surface parameters into homogeneous entities looks promising for future developments . delineation of such objects should produce morphometric patterns that match landscape patterns at given scales . this approach does not incorporate a priori knowledge on a specific landscape , therefore results are transferable . landform classification can further give meaning to the resulting building blocks ( elementary forms , pattern elements or morphometric primitives ) by adding semantics . minar and evans ( 2008 ) set up the stage by creating a unified system of geometric primitives . coupling multi - scale pattern analysis  with the help of local variance  with object delineation is another promising approach towards hierarchies built on such objects .
in the last decade landform classification and mapping has developed as one of the most active areas of geomorphometry . 
 however , translation from continuous models of elevation and its derivatives ( slope , aspect , and curvatures ) to landform divisions ( landforms and landform elements ) is filtered by two important concepts : scale and object ontology . 
 although acknowledged as being important , these two issues have received surprisingly little attention.this contribution provides an overview and prospects of object representation from dems as a function of scale . 
 relationships between object delineation and classification or regionalization are explored , in the context of differences between general and specific geomorphometry . a review of scales issues in geomorphometry 
 ranging from scale effects to scale optimization techniques 
 is followed by an analysis of pros and cons of using cells and objects in dem analysis . 
 prospects for coupling multi - scale analysis and object delineation are then discussed . within this context 
 , we propose discrete geomorphometry as a possible approach between general and specific geomorphometry . 
 discrete geomorphometry would apply to and describe land - surface divisions defined solely by the criteria of homogeneity in respect to a given land - surface parameter or a combination of several parameters . 
 homogeneity , in its turn , should always be relative to scale .