graph models have been used in image analysis for a long time .
the edited book @xcite gives an overview of methods in this field .
however , approaches from quantitative graph theory such as graph indices have not played a significant role in these applications so far .
this is to some extent surprising as it is not a far - fetched idea to model information contained in small patches of a textured image by graphs , and once this has been done , graph indices with their ability to extract in a quantitative form structural information from large collections of graphs lend themselves as a promising tool specifically for texture analysis .
a first step in this direction has been made in @xcite where a set of texture descriptors was introduced that combines a construction of graphs from image patches with well - known graph indices .
this set of texture descriptors was evaluated in @xcite in the context of a texture discrimination task . in @xcite , an example for texture - based image segmentation
was presented based on this work .
the present paper continues the work begun in @xcite and @xcite .
its purpose is twofold . on one hand
, it restates and slightly extends the experimental work from @xcite on texture segmentation , focussing on those descriptors that are based on entropy measures , which turned out particularly useful in the previous contributions . on the other hand ,
it undertakes a first attempt to analyse the graph - index based texture descriptors with regard to what kind of information they actually extract from a textured image .
[ [ structure - of - the - paper . ] ] structure of the paper .
+ + + + + + + + + + + + + + + + + + + + + + +    the remaining part of section  [ sec - intro ] briefly outlines the fields of research that are combined in this work , namely quantitative graph theory , see section  [ ssec - introqgt ] , graph models in image analysis with emphasis on the pixel graph and its edge - weighted variant , see section  [ ssec - introgmia ] , and finally texture analysis , section  [ ssec - introtex ] . in section  [ sec - gi ]
the construction of graph - entropy - based texture descriptors from @xcite is detailed .
section  [ sec - gac ] gives a brief account of the geodesic active contour method , a well - established approach for image segmentation that is based on the numerical solution of a partial differential equation .
texture segmentation combining the graph - entropy - based texture descriptors with the geodesic active contour method is demonstrated on two synthetic examples that represent typical realistic texture segmentation tasks , and a real - world example in section  [ sec - segex ] .
some theoretical analysis is presented in section  [ sec - tdana ] where ( one setup of ) graph - entropy - based texture descriptors is put into relation with fractal dimension measurements on a metric space derived from the pixel graph , and thus a connection is established between graph entropy methods and fractal - based texture descriptors . a short conclusion
, section  [ sec - conc ] , ends the paper .
quantitative measures for graphs have been developed for almost sixty years in mathematical chemistry as a means to analyse molecular graphs @xcite .
in the course of time , numerous _
graph indices _ have been derived based on edge connectivity , vertex degrees , distances , and information - theoretic concepts , see e.g.  the work @xcite that classifies over @xmath1 descriptors from literature and subjects them to a large - scale statistical evaluation on several test data sets .
recently , interesting new graph indices based on the so - called hosoya polynomial have been proposed @xcite .
fields of application have diversified in the last decades to include e.g.  biological and social networks , and other structures that can mathematically be modelled as graphs , see @xcite .
efforts to apply statistical methods to graph indices across these areas have been bundled in the emerging field of _ quantitative graph theory _ @xcite .
many contributions in this field group around the tasks of distinguishing and classifying graphs , and quantifying the differences between graphs .
the first task focusses on the ability of indices to _ uniquely _ distinguish large sets of individual graphs , termed _ discrimination power _
@xcite . for the latter task ,
_ inexact graph matching , _ the _ graph edit distance _ @xcite or other measures quantifying the size of substructures that are shared or not shared between two graphs are of particular importance , see also @xcite .
the concept of discrimination power has to be complemented for this purpose by the principle of _ smoothness _ of measures , see @xcite , which describes how similar the values of a measure are when it is applied to similar graphs . in @xcite , the quantitative measures of _ structure sensitivity _ and _
abruptness _ have been introduced in order to precisely analyse the discrimination power and smoothness of graph indices .
these measures are based on the average and maximum , respectively , of the changes of graph index values when the underlying graph is modified by one elementary edit step of the graph edit distance .
discrimination of graph structures by graph indices is also a crucial part of the texture analysis approach discussed here .
thus , discrimination power and the related notions of high structure sensitivity and low abruptness matter also in the present context .
however , unique identification of individual graphs is somewhat less important in texture analysis than when examining single graphs as in @xcite , as in texture analysis one is confronted with large collections of graphs associated with image pixels , and one is interested in separating these into a small number of classes representing regions . not only
will each class contain numerous graphs , but also the spatial arrangement of the associated pixels is to be taken into account as an additional source of information , as segments are expected to be connected .      as can be seen in @xcite , there are several ways in which image analysis can be linked to graph concepts .
a large class of approaches is based on graphs in which the pixels of a digital image take the role of vertices , and the set of edges is based on neighbourhood relations , with @xmath2- and @xmath3-neighbourhoods as most popular choices in 2d , and similar constructions in 3d , see ( * ? ? ?
* section  1.5.1 ) . to imprint actual image information on such a graph , one can furnish it with edge weights that are based on image contrasts . among others ,
the _ graph cut _
methods @xcite that have recently received much attention for applications such as image segmentation @xcite and correspondence problems @xcite make use of this concept .
this setup is also central for the work presented here , see the more detailed account of the pixel graph and edge - weighted pixel graph of an image in section  [ ssec - graphcons ] of the present paper .
generalising the pixel - graph framework , the graph perspective allows to transfer image processing methods from the regular mesh underlying standard digital images to non - regular meshes that can be related to scanned surface data @xcite but arise also from standard images when considering non - local models @xcite that have recently received great attention in image enhancement .
graph morphology , see e.g.@xcite , is one of these generalisations of image processing methods to non - regular meshes , but also variational and pde frameworks have been generalised in this way @xcite .
we briefly mention that graphs can also be constructed , after suitable preprocessing , from vertices representing image regions , see ( * ? ? ?
* section  1.5.2 ) , opening avenues to high - level semantic image interpretation by means of partition trees .
comparison of hierarchies of semantically meaningful partitions can then be achieved e.g.  using graph edit distance or related concepts @xcite .
returning to the pixel - graph setup which we will also use in this work , see section  [ ssec - graphcons ] , we point out a difference of our approach to those that represent the entire image in a single pixel graph .
we focus here on subgraphs related to small image patches , thus generating large sets of graphs whose vertex sets , connectivity and/or edge weights encode local image information . to extract meaningful information from such collections , statistical methods such as entropy - based graph indices are particularly suited .      in image analysis
, the term _ texture _ refers to the small - scale structure of image regions , and as such it has been an object of intensive investigation since the beginnings of digital image analysis .
for example , @xcite undertook approaches to define and analyse textures .    [ [ complementarity - of - texture - and - shape . ] ] complementarity of texture and shape .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    real - world scenes often consist of collections of distinct objects which in the process of imaging are mapped to regions in an image delineated by more or less sharp boundaries . while the geometric description of region boundaries is afforded by the concept of _ shape _ and typically involves large - scale structures , texture represents the appearance of the individual objects , either their surfaces if images in the case of reflection imaging ( such as photography of opaque objects ) , or their interior if transmission - based imaging modalities ( such as transmission microscopy , x - ray , magnetic resonance ) are being considered .
texture is then expressed in the distribution of intensities and their short - scale correlations within a region .
a frequently used mathematical formulation of this distinction is the _ cartoon - texture model _ that underlies many works on image restoration and enhancement , see e.g.  @xcite . in this approach , ( space - continuous )
images are described as the sum of two functions : a _ cartoon _ component from the space @xmath4 of functions of bounded variation , and a _ texture _ component from a suitable sobolev space . in a refined version of this decomposition @xcite
, noise is modelled as a third component assigned to a different function space .
note that also in image synthesis ( computer graphics ) the complementarity of shape and texture is used : here , textures are understood as intensity maps that are mapped on the surfaces of geometrically described objects .    the exact frontier between shape and texture information in a scene or image , however , is model - dependent .
the intensity variation of a surface is partly caused by geometric details of that surface . with a coarse - scale modelling of shape ,
small - scale variations are included in the texture description , whereas with a refined modelling of shape , some of these variations become part of the shape information .
for example , in the texture samples shown in figure  [ fi - e3](a ) and ( b ) , a geometric description with sufficiently fine granularity could capture individual leaves or blossoms as shapes , whereas the large - scale viewpoint treats the entire ensemble of leaves or blossoms as texture .
[ [ texture - models . ] ] texture models .
+ + + + + + + + + + + + + + +    capturing texture is virtually never possible on the basis of a single pixel . only the simplest of all textures , homogeneous intensity , can be described by a single intensity . for all other textures , intensities within neighbourhoods of suitable size ( that differs from texture to texture )
need to be considered to detect and classify textures .
moreover , there is a large variety of structures that can be constitutive of textures , ranging from periodic patterns in which the arrangement of intensity values follows strict geometric rules , via near - periodic and quasi - periodic structures to irregular patterns where just statistical properties of intensities within a neighbourhood are characteristic of the texture .
the texture samples in figure  [ fi - e3](a ) and ( b ) are located rather in the middle of the scale where both geometric relations and statistics of the intensities are characteristic of the texture ; near - periodic stripe patterns as in the later examples , figures  [ fi - e - stripes ] and [ fi - zebra ] , are more geometrically dominated .    with emphasis on different categories of textures within this continuum , numerous geometric and statistic approaches have been made over the decades to describe textures .
for example , frequency - based models @xcite emphasise the periodic or quasi - periodic aspect of textures .
statistics on intensities such as @xcite mark the opposite end of the scale , whereas models based on statistics of image derivative quantities like gradients @xcite or structure tensor entries @xcite attempt to combine statistical with geometrical information . a concept that differs significantly from both approaches has been proposed in @xcite where textures are described generatively via grammars .
also fractals @xcite have been proposed as a means to describe , distinguish and classify textures .
remember that a fractal is a geometric object , in fact a topological space , for which it is possible to determine , at least locally , a minkowski dimension ( or , almost identical , hausdorff dimension ) which differs from its topological dimension .
assume that the fractal is embedded in a surrounding euclidean space , and it is compact .
then it can be covered by a finite number of boxes , or balls , of prescribed size . when the size of the boxes or balls is sent to zero ,
the number of them which is needed to cover the structure grows with some power of the inverse box or ball size .
the minkowski dimension of the fractal is essentially the exponent in this power law .
the minkowski dimension of a fractal is often non - integer ( which is the reason for its name ) , however , a more precise definition is that minkowski and topological dimensions differ , which also includes cases like the peano curve whose minkowski dimension is an integer ( @xmath5 ) but still different from the topological dimension ( @xmath6 ) .
see also @xcite for different concepts of fractal dimension .
textured images can be associated with fractals by considering the image manifold , i.e.  the function graph if the image is modelled as a function over the image plane , which is naturally embedded in a product space of image domain and the range of intensity values .
for example , a planar grey - value image @xmath7 has the image manifold @xmath8 .
the dimension of this structure can be considered as a fractal - based texture descriptor .
this approach has been stated in @xcite where fractal dimension was put into relation with image roughness and the roughness of physical structures depicted in the image .
many works followed this approach , particularly in the 1980s and beginning 1990s when fractals were under particularly intensive investigation in theoretical and applied mathematics . in @xcite
several of these approaches are reviewed .
an attempt to analyse fractal dimension concepts for texture analysis is found in @xcite .
the concept has also been transferred to the analysis of 1d signals , see @xcite . during the last two decades
the interest in fractal methods has somewhat reduced but research in the field remains ongoing as can be seen from more recent publications , see e.g.  @xcite for signal analysis , @xcite for image analysis with application in material science . with regard to our analysis in section  [ sec - tdana ] that leads to a relation between graph methods and fractals , it is worth mentioning that already @xcite linked graph and fractal methods in image analysis , albeit not considering texture but shape description .
[ [ texture - segmentation . ] ] texture segmentation .
+ + + + + + + + + + + + + + + + + + + + +    the task of texture segmentation , i.e.  decomposing an image into several segments based on texture differences , has been studied for more than forty years , see @xcite .
a great variety of different approaches to the problem have been proposed since then .
many of these combine generic segmentation approaches , that could also be implemented for merely intensity - based segmentation , with sets of quantitative texture descriptors that are used as inputs to the segmentation .
for example , @xcite are based on active contour or active region models for segmentation , whereas @xcite is an example of a clustering - based method .
nevertheless , texture segmentation continues to challenge researchers ; in particular , improvements on the side of texture descriptors are still desirable .
note that the task of texture segmentation involves a conflict : on one hand , textures can not be detected on single - pixel level , necessitating the inclusion of neighbourhoods in texture descriptor computation .
on the other hand , the intended output of a segmentation is a unique assignment of each pixel to a segment , which means to fix the segment boundaries at pixel resolution . to allow sufficiently precise location of boundaries , texture descriptors should therefore not use larger patches than necessary to distinguish the textures present in an image .
the content of this paper is centred around a set of texture descriptors that have been introduced in @xcite based on graph representations of local image structure .
this model seems to be the first that exploited graph models in discrimination of textures .
note that even texture segmentation approaches in literature that use graph cuts for the segmentation task use non - graph - based texture descriptors to bring the texture information into the graph - cut formulation , see e.g.  @xcite .
our texture segmentation approach that was already shortly demonstrated in @xcite integrates instead graph - based texture descriptors into a non - graph - based segmentation framework , compare section  [ sec - gac ] .
throughout the history of texture processing , quantitative texture descriptors have played an important role .
assigning a tuple of numerical values to a texture provides an interface to established image processing algorithms that were originally designed to act on intensities , and thereby to devise modular frameworks for image processing tasks that involve texture information .    following this modular approach
, we will attack the texture segmentation task by combining a set of texture descriptors with a well - established image segmentation algorithm . in this section
we will introduce the texture descriptors whereas the following section will be devoted to describing the segmentation method .
given the variety of different textures that exist in natural images , it can not be expected that one single texture descriptor will be suitable to discriminate arbitrary textures . instead
, it will be sensible to come up with a set of descriptors that complement each other well in distinguishing different kinds of textures . to keep the set of descriptors at a manageable size
, the individual descriptors should nevertheless be able to discriminate substantial classes of textures .
on the other hand , it will be useful both for theoretical analysis and for practical computation if the set of descriptors is not entirely disparate but based on some powerful common concept .    in @xcite a family of texture descriptors
was introduced based on the application of several graph indices to graphs representing local image information . in combining six sets of graphs derived from an image ,
whose computation is based on common principles , with a number of different but related graph indices , this family of descriptors is indeed built on a common concepts .
the descriptors were evaluated in @xcite in a simple texture discrimination task , and turned out to yield results competitive with haralick features @xcite , a well - established concept in texture analysis . in this comparison ,
graph indices based on entropy measures stood out by their texture discrimination rate .    in the following , we recall the construction of texture descriptors from @xcite , focussing on a subset of the descriptors discussed there
the first step is the construction of graphs from image patches . in the second step ,
graph indices are computed from these graphs .
a discrete grey - scale image is given as an array of real intensity values sampled at the nodes of a regular grid .
the nodes are points @xmath9 in the plane where @xmath10 , @xmath11 .
the spatial mesh sizes @xmath12 and @xmath13 are often assumed to be @xmath6 in image processing , which we will do also here for simplicity .
denoting the intensity values by @xmath14 and assuming that @xmath15 , @xmath16 , the image is then described as the array @xmath17 .
the nodes of the grid , thus the pixels of the image , can naturally be considered as vertices of a graph in which neighbouring pixels are connected by edges .
we will call this graph the _ pixel graph _
@xmath18 of the image .
two common choices for what pixels are considered as neighbours are based on @xmath2-neighbourhoods , in which pixel @xmath19 has the two horizontal neighbours @xmath20 and the two vertical neighbours @xmath21 , or @xmath3-neighbourhoods , in which also the four diagonal neighbours @xmath22 are included in the neighbourhood . whereas the @xmath2-neighbourhood setting leads to a somewhat simpler pixel graph (
particularly , it is planar ) , whereas @xmath3-neighbourhoods are better suited to capture the geometry of the underlying ( euclidean ) plane . in this work
, we will mostly use @xmath3-neighbourhoods . see
1.5 ) for more variants of graphs assigned to images .
we upgrade the pixel graph to an _ edge - weighted pixel graph _ @xmath23 by defining edge weights @xmath24 for neighbouring pixels @xmath25 via @xmath26 i.e.  an @xmath27 sum of the spatial distance of grid nodes @xmath28 ( where @xmath29 denotes the euclidean norm ) , and the contrast @xmath30 of their corresponding intensity values , weighted by a positive _ contrast scale _ @xmath31 .
this construction can of course be generalised by replacing the euclidean norm in the image plane , and the @xmath27 sum by other norms . with various settings for these norms
, it has been used e.g.  in @xcite to construct larger spatially adaptive neighbourhoods in images , so - called _ morphological amoebas . _
see also @xcite for a more detailed description of the amoeba framework in a graph - based terminology .
all graphs that will enter the texture descriptor construction are derived from the pixel graph or the edge - weighted pixel graph of the image .
first , given a pixel @xmath32 and a radius @xmath33 , we define the euclidean patch graph @xmath34 as the subgraph of @xmath23 which includes all nodes @xmath35 with euclidean distance @xmath36 . in this graph ,
image information is encoded solely in the edge weights .
second , we define the adaptive patch graph @xmath37 as the subgraph of @xmath23 which includes all nodes @xmath35 for which @xmath23 contains a path from @xmath32 to @xmath35 with total weight less or equal @xmath38 . in the terminology of @xcite ,
the node set of @xmath37 is a morphological amoeba of amoeba radius @xmath38 around @xmath32 , which we will denote by @xmath39 .
note that the graph @xmath37 encodes image information not only in its edge weights , but also in its node set @xmath39 .
one obvious way to compute @xmath39 is by dijkstra s shortest path algorithm @xcite with @xmath32 as starting point .
a natural by - product of this algorithm , which is not used in amoeba - based image filtering as in @xcite etc . , is the dijkstra search tree , which we denote as @xmath40 .
this is the third candidate graph for our texture description .
image information is encoded in this graph in three ways : in the edge weights , the node set , and the connectivity of the tree .
dropping the edge weights from @xmath40 , we obtain an unweighted tree @xmath41 which still encodes image information in its node set and connectivity .
finally , a dijkstra search tree @xmath42 and its unweighted counterpart @xmath43 can be obtained by applying dijkstra s shortest path algorithm within the euclidean patch graph @xmath34 .
whereas @xmath42 encodes image information in the edge weights and connectivity , @xmath43 does so only in the connectivity .
applying these procedures to all pixels @xmath44 of a discrete image @xmath45 , we have therefore six collections of graphs which represent different combinations of three cues to local image information ( edge weights , node sets , connectivity ) and can therefore be expected to be suitable for texture discrimination . in the following
we will drop the arguments @xmath32 , @xmath38 and use simply @xmath46 etc .  to refer to the collections of graphs .
in order to turn the collections of graphs into quantitative texture descriptors suitable for texture analysis tasks , the realm of graph indices developed in quantitative graph theory lends itself as a powerful tool .    in @xcite , a selection of graph indices
was considered for this purpose , including on one hand distance - based graph indices ( the wiener index @xcite , the harary index @xcite and the balaban index @xcite ) and on the other hand entropy - based indices ( bonchev - trinajsti indices @xcite and dehmer entropies @xcite ) .
the so - obtained set of @xmath47 texture descriptors was evaluated in @xcite with respect to their discrimination power and diversity . using nine textures from a database @xcite ,
texture discrimination power was quantified based on simple statistics ( mean value and standard deviation ) of the descriptor values within single - texture patches , calibrating thresholds for certain and uncertain discrimination of textures within the set of textures .
diversity of descriptors was measured based on the overlap in the sets of texture pairs discriminated by different descriptors . despite the somewhat ad - hoc character of the threshold calibration , the study gives valuable hints for the selection of powerful subsets of the @xmath47 texture descriptors .    among the descriptors being analysed , particularly the entropy - based descriptors ranked medium to high regarding discrimination power for the sample set of textures . for the present work ,
we focus therefore on three entropy measures which we will recall in the following , namely the dehmer entropies @xmath48 and @xmath49 as well as bonchev and trinajsti s mean information on distances @xmath50 .
the latter is restricted by its construction to unweighted graphs , and is therefore used with the unweighted dijkstra trees @xmath51 and @xmath52 .
the dehmer entropies can be combined with all six graph collections .    in @xcite , the dehmer entropies on the patch graphs @xmath46 and @xmath53 achieved the highest rates of certain discrimination of textures , and outperformed the haralick features included in the study . some of the other descriptors based on dehmer entropies as well as the bonchev - trinajsti information measures achieved middle ranks , combining medium rates of certain discrimination with uncertain discrimination of almost all other texture pairs ; thereby , they were still comparable to haralick features and distance - based graph indices .
the measures considered here are based on shannon s entropy @xcite @xmath54 that measures the information content of a discrete probability measure @xmath55 , @xmath56 , @xmath57 .
( note that for @xmath58 one has to set in @xmath59 by limit . )
following @xcite , a discrete probability measure can be assigned to an arbitrary nonnegative - valued function @xmath60 , an _ information functional _ , via @xmath61 an entropy measure on an arbitrary information functional @xmath62 is then obtained by applying to
.      introduced in @xcite and further investigated in @xcite , the _
mean information on distances _ is the entropy measure resulting from an information functional on the path lengths in a graph .
let a graph @xmath63 with @xmath64 vertices @xmath65 be given , and let @xmath66 denote the length of a shortest path from @xmath67 to @xmath68 in @xmath63 ( unweighted , i.e.  each edge counting @xmath6 ) .
let @xmath69 be the diameter of @xmath63 .
for each @xmath70 , let @xmath71 the mean information on distances then is the entropy measure based on the information functional @xmath72 , i.e.@xmath73 let us shortly mention that @xcite also introduces the _ mean information on realised distances _
@xmath74 , which we will not further consider here . as an entropy measure
, @xmath74 can be derived from the information functional @xmath66 on the set of all vertex pairs @xmath75 , @xmath76 , of @xmath63 .
as pointed out in @xcite , @xmath74 can be generalised straightforwardly to edge - weighted graphs by measuring distances @xmath66 with edge weights , but a similar generalisation of @xmath77 would be degenerated because in generic cases all edge - weighted distances in a graph will be distinct , leading to @xmath78 for all realised values @xmath79 .
therefore we will use the mean information on distances @xmath77 only with the unweighted graphs @xmath52 and @xmath51 .
the two entropy measures @xmath80 and @xmath81 for unweighted graphs @xmath63 were introduced in @xcite .
their high discriminative power for large sets of graphs was impressively demonstrated in @xcite .
both measures rely on information functionals on the vertex set @xmath82 of @xmath63 whose construction involves spheres @xmath83 of varying radius @xmath79 around @xmath67 . note that the sphere @xmath83 in @xmath63 is the set of vertices @xmath68 with @xmath84 .    for @xmath48 the information functional @xmath85 on vertices @xmath67 of an unweighted graph @xmath63 is defined as @xcite @xmath86 where @xmath87 is the cardinality of the @xmath79-sphere around @xmath68 , with positive parameters @xmath88 .
( note that @xcite used a general exponential with base @xmath89 .
for the purpose of the present paper , however , this additional parameter is easily eliminated by multiplying the coefficients @xmath90 with @xmath91 . )    for @xmath49 the information functional @xmath92 relies on the quantities @xmath93 i.e.  @xmath94 is the sum of distances from @xmath67 to all points in its @xmath79-sphere . with similar parameters @xmath88 as before , one defines then @xmath95 as pointed out in @xcite , both information functionals , and thus the resulting entropy measures @xmath48 , @xmath49 , can be adapted to edge - weighted graphs @xmath63 via @xmath96 where distances @xmath66 are now measured using the edge weights , and @xmath97\to\mathbb{r}^+_0 $ ] is a decreasing function interpolating a reverse partial sum series of the original @xmath98 coefficients .    further following @xcite
, we focus on the specific choice @xmath99 ( an instance of the _ exponential weighting scheme _ from @xcite ) and obtain accordingly @xmath100 with a positive constant @xmath101 , which yields @xmath102 with a positive constant @xmath101 , as the final form of the information functionals for our construction of texture descriptors .
we use for our experiments a well - established segmentation method based on partial differential equations ( pde ) . introduced in @xcite , _ geodesic active contours _ ( gac ) perform a contrast - based segmentation of a ( grey - scale ) input image @xmath62 .
of course , other contrast - based segmentation methods could be chosen , including clustering @xcite or graph - cut methods @xcite .
advantages or disadvantages of these methods in connection with graph - entropy - based texture descriptors may be studied in future work . for the time being we focus on the texture descriptors themselves , thus it matters to use just one well - established standard method .      from the input image @xmath62 ,
a gaussian - smoothed image @xmath103 is computed , where @xmath104 is a gaussian kernel of standard deviation @xmath105 . from @xmath106 ,
one computes an edge map @xmath107 with the help of a decreasing and bounded function @xmath108 with @xmath109 .
a popular choice for @xmath110 is @xmath111 which has originally been introduced by perona and malik @xcite as a diffusivity function for nonlinear diffusion filtering of images .
herein , @xmath112 is a contrast parameter that acts as a threshold distinguishing high gradients ( indicating probable edges ) from small ones .
in addition to the input image , gac require an initial contour @xmath113 ( a regular closed curve ) specified e.g.  by user input .
this contour is embedded into a _ level set function _
@xmath114 , i.e.  @xmath114 is a sufficiently smooth function in the image plane whose zero level - set ( set of all points @xmath115 in the image plane for which @xmath116 ) is the given contour . for example , @xmath114 can be introduced as a signed distance function : @xmath117 is zero if @xmath115 lies on @xmath113 ; it is minus the distance of @xmath115 to @xmath113 if @xmath115 lies in the region enclosed by @xmath113 , and plus the same distance if @xmath115 lies in the outer region .
one takes then @xmath114 as initial condition at time @xmath118 for the parabolic pde @xmath119 for a time - dependent level - set function @xmath120 . at each time @xmath121 ,
an evolved contour can be extracted from @xmath122 as zero level - set . for suitable input images and initialisations and with appropriate parameters ,
the contours lock in at a steady state that provides a contrast - based segmentation .
to understand equation one can compare it to the curvature motion equation @xmath123 that would evolve all level curves of @xmath124 by an inward movement proportional to their curvature .
in , this inward movement of level curves is modulated by the edge map @xmath107 , which slows down the curve displacement at high - contrast locations , such that contours stick there .
the name _ geodesic active contours _ is due to the fact that the contour evolution associated to can be understood as gradient descent for the curve length of the contour in an image - adaptive metric ( a riemannian metric whose metric tensor is @xmath107 times the unit matrix ) , thus yielding a final contour that is a geodesic with respect to this metric .      in its pure form
, geodesic active contours require the initial contour ( at least most of it ) to be placed outside the object to be segmented . in some situations
, however , it is easier to specify an initial contour inside an object , particularly if the intensities within the object are fairly homogeneous but many spurious edges irritating the segmentation exist in the background region .
moreover , despite being able to handle also topology changes such as a splitting from one to several level curves encircling distinct objects to some extent , it has limitations when the topology of the segmentation becomes too complex .    as a remedy to both difficulties
, one can modify by adding a _ force term _
@xmath125 to its right - hand side .
adding a force term was proposed first in @xcite ( by the name of _ balloon force _ ) whereas the specific form of the force term weighted with @xmath110 was proposed in @xcite .
depending on the sign of @xmath126 , this force exercises an inward ( @xmath127 ) or outward ( @xmath128 ) pressure on the contour , which ( i ) speeds up the curve evolution , ( ii ) supports the handling of complex segmentation topologies , and ( iii ) enables for @xmath128 also segmentation of objects from initial contours placed inside .    the modified gac evolution with force term
, @xmath129 will be our segmentation method when performing texture segmentation based on only one quantitative texture descriptor . in this case
, the texture descriptor will be used as input image @xmath62 from which the edge map @xmath110 is computed .
it is straightforward to extend the gac method , including its modified version with force term to multi - channel input images @xmath130 where each location @xmath115 in the image plane is assigned an @xmath131-tuple @xmath132 of intensities . a common case , with @xmath133 ,
are rgb colour images .
in fact , equations and incur almost no change as even for multi - channel input images , one computes the evolution of a simple real - valued level - set function @xmath124 .
what is changed is the computation of the edge map @xmath110 : instead of the gradient norm @xmath134 one uses the frobenius norm @xmath135 of the jacobian @xmath136 where @xmath137 is the gaussian - smoothed input image , @xmath138 with @xmath139 , yielding @xmath140 as edge map .    equation with this edge map will be our segmentation method when performing texture segmentation with multiple texture descriptors .
the input image @xmath130 will have the individual texture descriptors as channels . to weight the influence of texture descriptors , the channels
may be multiplied by scalar factors .      for numerical computation
we rewrite pde as @xmath141 ( where we have omitted the argument of @xmath110 which is a fixed input function to the pde anyway ) .    following established practice ,
we use then an explicit ( euler forward ) numerical scheme where the right - hand side is spatially discretised as follows
. the first term , @xmath142 , is discretised using central differences . for the second term , @xmath143 , an upwind discretisation @xcite
is used in which the upwind direction for @xmath124 is determined based on the central - difference approximations of @xmath144 . the third term , @xmath145 , is discretised with an upwind discretisation , too . here
, the upwind direction depends on the components of @xmath146 and the sign of @xmath126 .
although a detailed stability analysis for this widely used type of explicit scheme for the gac equation seems to be missing , the scheme works for time step sizes @xmath147 up to ca .
@xmath148 ( for spatial mesh sizes of @xmath149 ) for @xmath150 , which needs to be reduced somewhat for non - zero @xmath126 . in our experiments in section  [ sec - segex ]
we use consistently @xmath151 .    for the level - set function @xmath124
, we use the signed distance function of the initial contour as initialisation .
since during the evolution the shape of the level - set function changes , creating steeper ascents in some regions but flattening slopes elsewhere , we re - initialise @xmath124 to the signed distance function of its current zero level set every @xmath152 iterations .
in this section , we present experiments on two synthetic and one real - world test image that demonstrate that graph - entropy - based texture descriptors can be used for texture - based segmentation .
an experiment similar to our second synthetic example was already presented in @xcite .
0.001    ( 1000 , 238 ) ( 127 , 0 ) pixels .
 * ( b ) * texture patch _ leaves _ , same size .  * ( c ) * test image composed from ( a ) and ( b ) , @xmath153 pixels .
 both texture patches are converted to greyscale , downscaled and clipped from the _ vistex _ database @xcite .
1995 massachusetts institute of technology .
developed by rosalind picard , chris graczyk , steve mann , josh wachman , len picard , and lee campbell at the media laboratory , mit , cambridge , massachusetts .
under general permission for scholarly use .
, title="fig : " ] ( 128 , 10 ) ( 381 , 0 ) pixels .  * ( b ) * texture patch _ leaves _ , same size .  * ( c ) * test image composed from ( a ) and ( b ) , @xmath153 pixels .
 both texture patches are converted to greyscale , downscaled and clipped from the _ vistex _ database @xcite .
1995 massachusetts institute of technology . developed by rosalind picard , chris graczyk , steve mann , josh wachman , len picard , and lee campbell at the media laboratory , mit , cambridge , massachusetts . under general permission for scholarly use .
, title="fig : " ] ( 382 , 10 ) ( 635 , 0 ) pixels .  * ( b ) * texture patch _ leaves _ , same size .  * ( c ) * test image composed from ( a ) and ( b ) , @xmath153 pixels .
 both texture patches are converted to greyscale , downscaled and clipped from the _ vistex _ database @xcite .
1995 massachusetts institute of technology .
developed by rosalind picard , chris graczyk , steve mann , josh wachman , len picard , and lee campbell at the media laboratory , mit , cambridge , massachusetts . under general permission for scholarly use .
, title="fig : " ] ( 636 , 10 )    in our first example we use a synthetic image , shown in figure  [ fi - e3](c ) , which is composed from two textures , see figure  [ fi - e3](a ) an ( b ) , with a simple shape ( the letter ` e ' ) switching between the two . note that the two textures were also among the nine textures studied in @xcite for the texture discrimination task . with its use of real - world textures ,
this synthetic example mimicks a realistic segmentation task .
its synthetic construction warrants at the same time a ground truth to compare segmentation results with .    0.001    ( 1000 , 492 ) ( 0 , 254)(c ) ( values rescaled to @xmath154 $ ] ) .
patch radii were fixed to @xmath155 , contrast scale to @xmath156 , weighting parameter to @xmath157 .
* top row , left to right : * * ( a ) * @xmath158 .
 * ( b ) * @xmath159 .  * ( c ) * @xmath160 .  * ( d ) * @xmath161 .  * bottom row , left to right : * * ( e ) * @xmath162 .  * ( f ) * @xmath163 .  * ( g ) * @xmath164 .  * ( h ) * @xmath165 .
, title="fig : " ] ( 1 , 264 ) ( 254 , 254)(c ) ( values rescaled to @xmath154 $ ] ) .
patch radii were fixed to @xmath155 , contrast scale to @xmath156 , weighting parameter to @xmath157 .
* top row , left to right : * * ( a ) * @xmath158 .  * ( b ) * @xmath159 .
* ( c ) * @xmath160 .  * ( d ) * @xmath161 .  * bottom row , left to right : * * ( e ) * @xmath162 .
 * ( f ) * @xmath163 .
 * ( g ) * @xmath164 .  * ( h ) * @xmath165 .
, title="fig : " ] ( 255 , 264 ) ( 508 , 254)(c ) ( values rescaled to @xmath154 $ ] ) .
patch radii were fixed to @xmath155 , contrast scale to @xmath156 , weighting parameter to @xmath157 .
* top row , left to right : * * ( a ) * @xmath158 .
 * ( b ) * @xmath159 .
 * ( c ) * @xmath160 .
 * ( d ) * @xmath161 .  * bottom row , left to right : * * ( e ) * @xmath162 .
* ( f ) * @xmath163 .  * ( g ) * @xmath164 .  * ( h ) * @xmath165 . , title="fig : " ] ( 509 , 264 ) ( 762 , 254)(c ) ( values rescaled to @xmath154 $ ] ) .
patch radii were fixed to @xmath155 , contrast scale to @xmath156 , weighting parameter to @xmath157 .
* top row , left to right : * * ( a ) * @xmath158 .
 * ( b ) * @xmath159 .  * ( c ) * @xmath160 .  * ( d ) * @xmath161 .  * bottom row , left to right : * * ( e ) * @xmath162 .
* ( f ) * @xmath163 .  * ( g ) * @xmath164 .  * ( h ) * @xmath165 . ,
title="fig : " ] ( 763 , 264 ) ( 0 , 0)(c ) ( values rescaled to @xmath154 $ ] ) .
patch radii were fixed to @xmath155 , contrast scale to @xmath156 , weighting parameter to @xmath157 .
* top row , left to right : * * ( a ) * @xmath158 .  * ( b ) * @xmath159 .  * ( c ) * @xmath160 .  * ( d ) * @xmath161 .  * bottom row , left to right : * * ( e ) * @xmath162 .
 * ( f ) * @xmath163 .  * ( g ) * @xmath164 .  * ( h ) * @xmath165 . , title="fig : " ] ( 1 , 10 ) ( 254 , 0)(c ) ( values rescaled to @xmath154 $ ] ) .
patch radii were fixed to @xmath155 , contrast scale to @xmath156 , weighting parameter to @xmath157 .
* top row , left to right : * * ( a ) * @xmath158 .
 * ( b ) * @xmath159 .  * ( c ) * @xmath160 .  * ( d ) * @xmath161 .  * bottom row , left to right : * * ( e ) * @xmath162 .  * ( f ) * @xmath163 .  * ( g ) * @xmath164 .  * ( h ) * @xmath165 .
, title="fig : " ] ( 255 , 10 ) ( 508 , 0)(c ) ( values rescaled to @xmath154 $ ] ) .
patch radii were fixed to @xmath155 , contrast scale to @xmath156 , weighting parameter to @xmath157 .
* top row , left to right : * * ( a ) * @xmath158 .  * ( b ) * @xmath159 .  * ( c ) * @xmath160 .  * ( d ) * @xmath161 .  * bottom row , left to right : * * ( e ) * @xmath162 .
 * ( f ) * @xmath163 .  * ( g ) * @xmath164 .  * ( h ) * @xmath165 . , title="fig : " ] ( 509 , 10 ) ( 762 , 0)(c ) ( values rescaled to @xmath154 $ ] ) .
patch radii were fixed to @xmath155 , contrast scale to @xmath156 , weighting parameter to @xmath157 .
* top row , left to right : * * ( a ) * @xmath158 .
 * ( b ) * @xmath159 .  * ( c ) * @xmath160 .  * ( d ) * @xmath161 .
 * bottom row , left to right : * * ( e ) * @xmath162 .
 * ( f ) * @xmath163 .  * ( g ) * @xmath164 .  * ( h ) * @xmath165 . , title="fig : " ] ( 763 , 10 )    figure  [ fi - e3-gi ] shows the results of eight graph - entropy - based texture descriptors for the test image .
in particular , the combination of the dehmer entropy @xmath48 with all six graph variants from section  [ ssec - graphcons ] is shown as well as @xmath49 on the weighted dijkstra trees in non - adaptive and adaptive patches .
patch radii were fixed to @xmath155 for both non - adaptive and adaptive patches , whereas the contrast scale was chosen as @xmath156 .
these parameter settings have already been used in @xcite ; they are based on values that work across various test images in the context of morphological amoeba image filtering .
further investigation of variation of these parameters is left for future work .    visual inspection of figure  [ fi - e3-gi ] indicates that for this specific textured image , the entropy measure @xmath48 separates the two textures well in combination in particular with the weighted dijkstra tree settings , in both adaptive and non - adaptive patches , see frames ( b ) and ( f ) .
the other @xmath48 results in frames ( c ) , ( e ) and ( g ) show insufficient contrast along some parts of the contour of the letter ` e ' .
the index @xmath166 in frame ( a ) , which was identified in @xcite as a descriptor with high texture discrimination power , does not distinguish these two textures clearly but creates massive over - segmentation within each of them . in a sense , this over - segmentation is the downside of the high texture discrimination power of the descriptor .
note , however , that also other @xmath53-based descriptors tend to this kind of over - segmentation .    regarding the @xmath49 index , figure  [ fi - e3-gi ] ( d ) and ( h ) ,
there is a huge difference between the adaptive and non - adaptive patch setting .
distinction of the two textures is much better when using non - adaptive patches .    0.001    ( 1000 , 238 ) ( 0 , 0)(c ) .
the edge map is computed from the graph - entropy - based texture descriptor @xmath163 from figure  [ fi - e3-gi](f ) using pre - smoothing with @xmath167 , perona - malik edge - stopping function with @xmath168 , with expansion force @xmath169 . *
left to right : * * ( a ) * initial contour ( @xmath118 ) .  * ( b ) * @xmath170 .  * ( c ) * @xmath171 .  * ( d ) * @xmath172 ( steady state ) . , title="fig : " ] ( 1 , 10 ) ( 254 , 0)(c ) .
the edge map is computed from the graph - entropy - based texture descriptor @xmath163 from figure  [ fi - e3-gi](f ) using pre - smoothing with @xmath167 , perona - malik edge - stopping function with @xmath168 , with expansion force @xmath169 .
* left to right : * * ( a ) * initial contour ( @xmath118 ) .  * ( b ) * @xmath170 .  * ( c ) * @xmath171 .
 * ( d ) * @xmath172 ( steady state ) . , title="fig : " ] ( 255 , 10 ) ( 508 , 0)(c ) .
the edge map is computed from the graph - entropy - based texture descriptor @xmath163 from figure  [ fi - e3-gi](f ) using pre - smoothing with @xmath167 , perona - malik edge - stopping function with @xmath168 , with expansion force @xmath169 .
* left to right : * * ( a ) * initial contour ( @xmath118 ) .  * ( b ) * @xmath170 .  * ( c ) * @xmath171 .  * ( d ) * @xmath172 ( steady state ) .
, title="fig : " ] ( 509 , 10 ) ( 762 , 0)(c ) .
the edge map is computed from the graph - entropy - based texture descriptor @xmath163 from figure  [ fi - e3-gi](f ) using pre - smoothing with @xmath167 , perona - malik edge - stopping function with @xmath168 , with expansion force @xmath169 . *
left to right : * * ( a ) * initial contour ( @xmath118 ) .  * ( b ) * @xmath170 .  * ( c ) * @xmath171 .  * ( d ) * @xmath172 ( steady state ) .
, title="fig : " ] ( 763 , 10 )    finally , we show in figure  [ fi - e3-gac ] geodesic active contour segmentation of the test image with the descriptor @xmath163 .
we start from an initial contour inside the ` e ' shape , see figure  [ fi - e3-gac](a ) , and use an expansion force ( @xmath169 ) to drive the contour evolution in an outward direction . frames ( b ) and ( c ) show two intermediary stages of the evolution , where it is evident that the contour starts to align with the boundary between the two textures .
frame ( d ) shows the steady state reached after @xmath173 iterations ( @xmath172 ) . here
, the overall shape of the letter ` e ' is reasonably approximated , with deviations coming from small - scale texture details .
precision of the segmentation could be increased slightly by combining more than one texture descriptor .
we do not follow this direction at this point .
0.001    ( 1000 , 492 ) ( 0 , 254 ) pixels ) showing a stripe - textured shape in front of a noise background .  * ( b  d ) * texture descriptors based on graph entropies applied in adaptive patches , @xmath155 , @xmath156 , @xmath157 ; values rescaled to @xmath154 $ ] . *
( b ) * @xmath162 .
 * ( c ) * @xmath163 .
 * ( d ) * @xmath164 .  * bottom row : * geodesic active contour segmentation of ( a ) using the texture descriptor ( @xmath48 on @xmath174 ) from ( b ) , same parameters as in figure  [ fi - e3-gi ] except for @xmath175 .
* left to right : * * ( e ) * initial contour ( @xmath118 ) .  * ( f ) * @xmath176 .
 * ( g ) * @xmath177 .  * ( h ) * @xmath178 ( steady state ) .
, title="fig : " ] ( 1 , 264 ) ( 254 , 254 ) pixels ) showing a stripe - textured shape in front of a noise background .
 * ( b  d ) * texture descriptors based on graph entropies applied in adaptive patches , @xmath155 , @xmath156 , @xmath157 ; values rescaled to @xmath154 $ ] . * ( b ) * @xmath162 .  * ( c ) * @xmath163 .  * ( d ) * @xmath164 .  * bottom row : * geodesic active contour segmentation of ( a ) using the texture descriptor ( @xmath48 on @xmath174 ) from ( b ) , same parameters as in figure  [ fi - e3-gi ] except for @xmath175 . * left to right : * * ( e ) * initial contour ( @xmath118 ) .  * ( f ) * @xmath176 .  * ( g ) * @xmath177 .  * ( h ) * @xmath178 ( steady state ) . , title="fig : " ] ( 255 , 264 ) ( 508 , 254 ) pixels ) showing a stripe - textured shape in front of a noise background .  * ( b  d ) * texture descriptors based on graph entropies applied in adaptive patches , @xmath155 , @xmath156 , @xmath157 ; values rescaled to @xmath154 $ ] . * ( b ) * @xmath162 .  * ( c ) * @xmath163 .  * ( d ) * @xmath164 .  * bottom row : * geodesic active contour segmentation of ( a ) using the texture descriptor ( @xmath48 on @xmath174 ) from ( b ) , same parameters as in figure  [ fi - e3-gi ] except for @xmath175 . * left to right : * * ( e ) * initial contour ( @xmath118 ) .  * ( f ) * @xmath176 .  * ( g ) * @xmath177 .  * ( h ) * @xmath178 ( steady state ) . , title="fig : " ] ( 509 , 264 ) ( 762 , 254 ) pixels ) showing a stripe - textured shape in front of a noise background .  * ( b  d ) * texture descriptors based on graph entropies applied in adaptive patches , @xmath155 , @xmath156 , @xmath157 ; values rescaled to @xmath154 $ ] . * ( b ) * @xmath162 .  * ( c ) * @xmath163 .  * ( d ) * @xmath164 .  * bottom row : * geodesic active contour segmentation of ( a ) using the texture descriptor ( @xmath48 on @xmath174 ) from ( b ) , same parameters as in figure  [ fi - e3-gi ] except for @xmath175 . * left to right : * * ( e ) * initial contour ( @xmath118 ) .  * ( f ) * @xmath176 .  * ( g ) * @xmath177 .  * ( h ) * @xmath178 ( steady state ) . , title="fig : " ] ( 763 , 264 ) ( 0 , 0 ) pixels ) showing a stripe - textured shape in front of a noise background .  * ( b  d ) * texture descriptors based on graph entropies applied in adaptive patches , @xmath155 , @xmath156 , @xmath157 ; values rescaled to @xmath154 $ ] . * ( b ) * @xmath162 .  * ( c ) * @xmath163 .  * ( d ) * @xmath164 .  * bottom row : * geodesic active contour segmentation of ( a ) using the texture descriptor ( @xmath48 on @xmath174 ) from ( b ) , same parameters as in figure  [ fi - e3-gi ] except for @xmath175 . * left to right : * * ( e ) * initial contour ( @xmath118 ) .  * ( f ) * @xmath176 .  * ( g ) * @xmath177 .  * ( h ) * @xmath178 ( steady state ) .
, title="fig : " ] ( 1 , 10 ) ( 254 , 0 ) pixels ) showing a stripe - textured shape in front of a noise background .
 * ( b  d ) * texture descriptors based on graph entropies applied in adaptive patches , @xmath155 , @xmath156 , @xmath157 ; values rescaled to @xmath154 $ ] . * ( b ) * @xmath162 .  * ( c ) * @xmath163 .  * ( d ) * @xmath164 .  * bottom row : * geodesic active contour segmentation of ( a ) using the texture descriptor ( @xmath48 on @xmath174 ) from ( b ) , same parameters as in figure  [ fi - e3-gi ] except for @xmath175 . * left to right : * * ( e ) * initial contour ( @xmath118 ) .
 * ( f ) * @xmath176 .  * ( g ) * @xmath177 .  * ( h ) * @xmath178 ( steady state ) . , title="fig : " ] ( 255 , 10 ) ( 508 , 0 ) pixels ) showing a stripe - textured shape in front of a noise background .  * ( b  d ) * texture descriptors based on graph entropies applied in adaptive patches , @xmath155 , @xmath156 , @xmath157 ; values rescaled to @xmath154 $ ] . *
( b ) * @xmath162 .
 * ( c ) * @xmath163 .
 * ( d ) * @xmath164 .  * bottom row : * geodesic active contour segmentation of ( a ) using the texture descriptor ( @xmath48 on @xmath174 ) from ( b ) , same parameters as in figure  [ fi - e3-gi ] except for @xmath175 . * left to right : * * ( e ) * initial contour ( @xmath118 ) .  * ( f ) * @xmath176 .  * ( g ) * @xmath177 .  * ( h ) * @xmath178 ( steady state ) . , title="fig : " ] ( 509 , 10 ) ( 762 , 0 ) pixels ) showing a stripe - textured shape in front of a noise background .  * ( b  d ) * texture descriptors based on graph entropies applied in adaptive patches , @xmath155 , @xmath156 , @xmath157 ; values rescaled to @xmath154 $ ] . *
( b ) * @xmath162 .
 * ( c ) * @xmath163 .
 * ( d ) * @xmath164 .  * bottom row : * geodesic active contour segmentation of ( a ) using the texture descriptor ( @xmath48 on @xmath174 ) from ( b ) , same parameters as in figure  [ fi - e3-gi ] except for @xmath175 . * left to right : * * ( e ) * initial contour ( @xmath118 ) .  * ( f ) * @xmath176 .  * ( g ) * @xmath177 .  * ( h ) * @xmath178 ( steady state ) . , title="fig : " ] ( 763 , 10 )    in our second experiment , figure  [ fi - e - stripes ] , we use again a synthetic test image where foreground and background segments are defined using the ` e ' letter shape such that again the desired segmentation result is known as a ground truth .
also in this image we combine two realistic textures which can be seen as a simplified version of the foreground and background textures of the real - world test image , figure  [ fi - zebra](a ) , used in the next section .
this time , the foreground is filled with a stripe pattern whereas the background is noise with uniform distribution in the intensity range @xmath154 $ ] , see figure  [ fi - e - stripes](a ) . in frames ( b)(d ) of figure  [ fi - e - stripes ]
we show the texture descriptors based on @xmath48 with the three graph settings in adaptive patches , using again @xmath155 and @xmath156 .
the descriptor @xmath179 that was used for the segmentation in section  [ ssec - segex-1 ] visually does not distinguish foreground from background satisfactorily here , whereas @xmath180 that provided no clear distinction of the two textures in section  [ ssec - segex-1 ] clearly stands out here .
this underlines the necessity of considering multiple descriptors which complement each other in distinguishing textures .
our gac segmentation of the test image shown in frames ( e)(h ) is based on the texture descriptor @xmath180 and quickly converges to a fairly good approximation of the segment boundary .
0.001    ( 1000 , 758 ) ( 0 , 516 ) pixels ) , converted to greyscale . original image from ` http://en.wikipedia.org/wiki/file:grants_zebra.jpg ` , author ruby 1x2 , released to public domain .  * ( b ) * texture descriptor @xmath158 , @xmath181 , @xmath156 , @xmath157 .  * ( c ) * texture descriptor @xmath159 , parameters as in ( b ) .
 * middle row , left to right : * * ( d ) * texture descriptor @xmath161 , parameters as in ( b ) .  * ( e ) * texture descriptor @xmath182 , parameters as in ( b ) .  * ( f ) * initial contour for geodesic active contour segmentation ( @xmath118 ) .
 * bottom row , left to right : * * ( g ) * geodesic active contours with edge - stopping function computed from the texture descriptors @xmath161 shown in ( d ) and @xmath182 shown in ( e ) with @xmath183 , perona - malik edge - stopping function , @xmath184 , @xmath185 .  * ( h ) * same as ( g ) but @xmath186 .
 * ( i ) * same as ( g ) but @xmath187 ( steady state ) .
, title="fig : " ] ( 1 , 526 ) ( 339 , 516 ) pixels ) , converted to greyscale .
original image from ` http://en.wikipedia.org/wiki/file:grants_zebra.jpg ` , author ruby 1x2 , released to public domain .  * ( b ) * texture descriptor @xmath158 , @xmath181 , @xmath156 , @xmath157 .  * ( c ) * texture descriptor @xmath159 , parameters as in ( b ) .
 * middle row , left to right : * * ( d ) * texture descriptor @xmath161 , parameters as in ( b ) .  * ( e ) * texture descriptor @xmath182 , parameters as in ( b ) .  * ( f ) * initial contour for geodesic active contour segmentation ( @xmath118 ) .
 * bottom row , left to right : * * ( g ) * geodesic active contours with edge - stopping function computed from the texture descriptors @xmath161 shown in ( d ) and @xmath182 shown in ( e ) with @xmath183 , perona - malik edge - stopping function , @xmath184 , @xmath185 .  * ( h ) * same as ( g ) but @xmath186 .
 * ( i ) * same as ( g ) but @xmath187 ( steady state ) .
, title="fig : " ] ( 340 , 526 ) ( 678 , 516 ) pixels ) , converted to greyscale .
original image from ` http://en.wikipedia.org/wiki/file:grants_zebra.jpg ` , author ruby 1x2 , released to public domain .  * ( b ) * texture descriptor @xmath158 , @xmath181 , @xmath156 , @xmath157 .  * ( c ) * texture descriptor @xmath159 , parameters as in ( b ) .
 * middle row , left to right : * * ( d ) * texture descriptor @xmath161 , parameters as in ( b ) .
 * ( e ) * texture descriptor @xmath182 , parameters as in ( b ) .
 * ( f ) * initial contour for geodesic active contour segmentation ( @xmath118 ) .  * bottom row , left to right : * * ( g ) * geodesic active contours with edge - stopping function computed from the texture descriptors @xmath161 shown in ( d ) and @xmath182 shown in ( e ) with @xmath183 , perona - malik edge - stopping function , @xmath184 , @xmath185 .  * ( h ) * same as ( g ) but @xmath186 .  * ( i ) * same as ( g ) but @xmath187 ( steady state ) . , title="fig : " ] ( 680 , 526 ) ( 0 , 258 ) pixels ) , converted to greyscale . original image from ` http://en.wikipedia.org/wiki/file:grants_zebra.jpg ` , author ruby 1x2 , released to public domain .  * ( b ) * texture descriptor @xmath158 , @xmath181 , @xmath156 , @xmath157 .  * ( c ) * texture descriptor @xmath159 , parameters as in ( b ) .
* middle row , left to right : * * ( d ) * texture descriptor @xmath161 , parameters as in ( b ) .  * ( e ) * texture descriptor @xmath182 , parameters as in ( b ) .  * ( f ) * initial contour for geodesic active contour segmentation ( @xmath118 ) .  * bottom row , left to right : * * ( g ) * geodesic active contours with edge - stopping function computed from the texture descriptors @xmath161 shown in ( d ) and @xmath182 shown in ( e ) with @xmath183 , perona - malik edge - stopping function , @xmath184 , @xmath185 .  * ( h ) * same as ( g ) but @xmath186 .
 * ( i ) * same as ( g ) but @xmath187 ( steady state ) . ,
title="fig : " ] ( 1 , 268 ) ( 339 , 258 ) pixels ) , converted to greyscale .
original image from ` http://en.wikipedia.org/wiki/file:grants_zebra.jpg ` , author ruby 1x2 , released to public domain .  * ( b ) * texture descriptor @xmath158 , @xmath181 , @xmath156 , @xmath157 .  * ( c ) * texture descriptor @xmath159 , parameters as in ( b ) .
 * middle row , left to right : * * ( d ) * texture descriptor @xmath161 , parameters as in ( b ) .  * ( e ) * texture descriptor @xmath182 , parameters as in ( b ) .  * ( f ) * initial contour for geodesic active contour segmentation ( @xmath118 ) .
 * bottom row , left to right : * * ( g ) * geodesic active contours with edge - stopping function computed from the texture descriptors @xmath161 shown in ( d ) and @xmath182 shown in ( e ) with @xmath183 , perona - malik edge - stopping function , @xmath184 , @xmath185 .  * ( h ) * same as ( g ) but @xmath186 .
 * ( i ) * same as ( g ) but @xmath187 ( steady state ) .
, title="fig : " ] ( 340 , 268 ) ( 678 , 258 ) pixels ) , converted to greyscale .
original image from ` http://en.wikipedia.org/wiki/file:grants_zebra.jpg ` , author ruby 1x2 , released to public domain .  * ( b ) * texture descriptor @xmath158 , @xmath181 , @xmath156 , @xmath157 .  * ( c ) * texture descriptor @xmath159 , parameters as in ( b ) .
 * middle row , left to right : * * ( d ) * texture descriptor @xmath161 , parameters as in ( b ) .
 * ( e ) * texture descriptor @xmath182 , parameters as in ( b ) .
 * ( f ) * initial contour for geodesic active contour segmentation ( @xmath118 ) .  * bottom row , left to right : * * ( g ) * geodesic active contours with edge - stopping function computed from the texture descriptors @xmath161 shown in ( d ) and @xmath182 shown in ( e ) with @xmath183 , perona - malik edge - stopping function , @xmath184 , @xmath185 .
 * ( h ) * same as ( g ) but @xmath186 .
 * ( i ) * same as ( g ) but @xmath187 ( steady state ) .
, title="fig : " ] ( 680 , 268 ) ( 0 , 0 ) pixels ) , converted to greyscale .
original image from ` http://en.wikipedia.org/wiki/file:grants_zebra.jpg ` , author ruby 1x2 , released to public domain .
 * ( b ) * texture descriptor @xmath158 , @xmath181 , @xmath156 , @xmath157 .  * ( c ) * texture descriptor @xmath159 , parameters as in ( b ) .  * middle row , left to right : * * ( d ) * texture descriptor @xmath161 , parameters as in ( b ) .  * ( e ) * texture descriptor @xmath182 , parameters as in ( b ) .  * ( f ) * initial contour for geodesic active contour segmentation ( @xmath118 ) .  * bottom row , left to right : * * ( g ) * geodesic active contours with edge - stopping function computed from the texture descriptors @xmath161 shown in ( d ) and @xmath182 shown in ( e ) with @xmath183 , perona - malik edge - stopping function , @xmath184 , @xmath185 .  * ( h ) * same as ( g ) but @xmath186 .  * ( i ) * same as ( g ) but @xmath187 ( steady state ) .
, title="fig : " ] ( 1 , 10 ) ( 339 , 0 ) pixels ) , converted to greyscale .
original image from ` http://en.wikipedia.org/wiki/file:grants_zebra.jpg ` , author ruby 1x2 , released to public domain .  * ( b ) * texture descriptor @xmath158 , @xmath181 , @xmath156 , @xmath157 .  * ( c ) * texture descriptor @xmath159 , parameters as in ( b ) .
 * middle row , left to right : * * ( d ) * texture descriptor @xmath161 , parameters as in ( b ) .  * ( e ) * texture descriptor @xmath182 , parameters as in ( b ) .  * ( f ) * initial contour for geodesic active contour segmentation ( @xmath118 ) .
 * bottom row , left to right : * * ( g ) * geodesic active contours with edge - stopping function computed from the texture descriptors @xmath161 shown in ( d ) and @xmath182 shown in ( e ) with @xmath183 , perona - malik edge - stopping function , @xmath184 , @xmath185 .  * ( h ) * same as ( g ) but @xmath186 .
 * ( i ) * same as ( g ) but @xmath187 ( steady state ) .
, title="fig : " ] ( 340 , 10 ) ( 678 , 0 ) pixels ) , converted to greyscale .
original image from ` http://en.wikipedia.org/wiki/file:grants_zebra.jpg ` , author ruby 1x2 , released to public domain .  * ( b ) * texture descriptor @xmath158 , @xmath181 , @xmath156 , @xmath157 .  * ( c ) * texture descriptor @xmath159 , parameters as in ( b ) .
 * middle row , left to right : * * ( d ) * texture descriptor @xmath161 , parameters as in ( b ) .  * ( e ) * texture descriptor @xmath182 , parameters as in ( b ) .  * ( f ) * initial contour for geodesic active contour segmentation ( @xmath118 ) .
 * bottom row , left to right : * * ( g ) * geodesic active contours with edge - stopping function computed from the texture descriptors @xmath161 shown in ( d ) and @xmath182 shown in ( e ) with @xmath183 , perona - malik edge - stopping function , @xmath184 , @xmath185 .  * ( h ) * same as ( g ) but @xmath186 .
 * ( i ) * same as ( g ) but @xmath187 ( steady state ) .
, title="fig : " ] ( 680 , 10 )    in our last experiment , figure  [ fi - zebra ] , we consider a real - world image showing a zebra , see frame ( a ) . in a sense , this experiment resembles the synthetic case from section  [ ssec - segex-2 ] because again a foreground dominated by a clear stripe pattern is to be distinguished from a background filled with small - scale detail . in frames ( b)(e )
four texture descriptors are shown . with regard to the higher resolution of the test image ,
the patch radius has been chosen slightly larger than in the previous examples , @xmath181 , whereas @xmath156 was retained . as can be seen in frame ( b )
, @xmath158 shows the same kind of over - segmentation behaviour as observed in section  [ ssec - segex-1 ] , however , it also separates a large part of the zebra shape well from the background .
the second descriptor , @xmath159 in frame ( c ) , appears unsuitable here because it does not yield sufficiently similar values within the black and white stripes to recognise these as a common texture .
in contrast , @xmath161 and @xmath182 in figure  [ fi - zebra](d ) and ( e ) , respectively , achieve this largely .
our gac segmentation in frames ( f)(i ) uses a larger gaussian kernel for pre - smoothing than before , @xmath183 , to flatten out small - scale inhomogeneities in the texture descriptors , and combines the two descriptors from ( d ) and ( e ) . with these data , a large part of the zebra including the head and front part of the torso is segmented in the final steady state .
not included are the rear part and the forelegs .
note that in the foreleg part the stripes are much thinner than in the segmented region , apparently preventing the recognition of this texture as a continuation of the one from the head and front torso .
in contrast , the rear part of the torso shown very thick stripes which under the patch size chosen decompose into separate ( homogeneous ) textures for black and white stripes , as is also visible in the texture descriptors ( d ) and ( e ) themselves .
further investigation of parameter variations and inclusion of more texture descriptors might improve this preliminary result in the future .
in this section , we undertake an attempt to analyse the texture descriptors based on the entropy measures @xmath48 and @xmath49 , focussing on the question what properties of textures are actually encoded in their information functionals @xmath85 and @xmath92 .
part of this analysis is on a heuristic level at the present stage of research , and future work will have to be invested to add precision to these arguments .
this applies to the limiting procedure in section  [ ssec - graphfrac ] as well as to the concept of local fractal dimension arising in section  [ ssec - fracdim ] .
we believe , however , that even in its present shape , the analysis provided in the following gives valuable intuition about the principles underlying our texture descriptors .      for the purpose of our analysis
, we generalise the information functional @xmath85 from directly to edge - weighted graphs by replacing the series @xmath188 of cardinalities from with the monotone increasing function @xmath189 , @xmath190 that measures volumes of spheres with arbitrary radius . assuming the exponential weighting scheme and large @xmath191 this yields @xmath192 an analogous generalisation of is @xmath193      we assume now that the image is sampled successively on finer grids , with decreasing @xmath194 .
note that the number of vertices of any region of the edge - weighted pixel graph , or any of the derived edge - weighted graphs introduced in section  [ ssec - graphcons ] , grows in this process with @xmath195 . by using the volumes of spheres instead of the original cardinalities ,
provides a re - normalisation that compensates this effect in .
thus , it is possible to consider the limit case @xmath196 . in this limit case ,
graphs turn into metric spaces representing the structure of a space - continuous image .
additionally , these metric spaces are endowed with a volume measure which is the limit case of the discrete measures on graphs given by vertex counting .    in simple cases , these metric spaces with volume measure
can be manifolds .
for example , for a homogeneous grey image without any contrast , the limit of the edge - weighted pixel graph is an approximation to a plane , i.e.  a @xmath5-dimensional manifold . for an image with extreme contrasts in one direction ,
e.g.  a stripe pattern , the edge - weighted pixel graphs will be path graphs , resulting in a metric space as limit which is essentially a @xmath6-dimensional manifold .
finally , in the extreme case of a noise image in which neighbouring pixels have almost nowhere similar grey - values , the graph will practically decompose into numerous isolated connected components , corresponding to a discrete space of dimension @xmath197 .    for more general textured images
, the limit spaces will possess a more complicated topological structure . at the same time
, it remains possible , of course , to measure volumes of spheres of different radii in these spaces .
clearly , sphere volumes will increase with sphere radius .
if they fulfil a power law , the ( possible non - integer ) exponent can immediately be interpreted as a dimension .
the space itself is then interpreted as some type of fractal @xcite .
the dimension concept underlying here is almost that of the minkowski dimension ( closely related to hausdorff dimension ) that is frequently used in fractal theory , with the difference that the volume measure here is inside the object being measured instead of in an embedding space .
based on the above reasoning , values of the dimension will range between @xmath197 and @xmath5 .
note that even in situations in which there is no global power law for the sphere volumes , and therefore no global dimension , power laws , possibly with varying exponents , will still be approximated for a given sphere centre in suitable ranges of the radius , thus allowing to define the fractal dimension as a quantity varying in space and resolution .
this resembles the situation with most fractal concepts being applied to real - world data : the power laws that are required to hold for an ideal fractal across all scales will be found only for certain ranges of scales in reality .
dijkstra trees , too , turn into @xmath6-dimensional manifolds in the case of sharp stripe images ; for other cases they will also yield fractals .
fractal dimensions , wherever applicable , will be below those observed with the corresponding full edge - weighted pixel graphs , thus , the relevant range of dimensions is again bounded by @xmath197 from below and @xmath5 from above
.    one word of care must be said at this point .
the fractal structures obtained here as limit cases of graphs for @xmath196 are not identical with the image manifolds whose fractal structures are studied as a means of texture analysis in @xcite and others .
in fact , fractal dimensions of the latter , measured as minkowski dimensions by means of the embedding of the image manifold of a grey - scale image in three - dimensional euclidean space , range from @xmath5 to @xmath198 with increasing roughness of the image , whereas the dimensions measured in the present work go down from @xmath5 to @xmath197 for increasing image roughness . whereas it can be conjectured that these two fractal structures are related , future work will be needed to gain clarity about this relationship .      based on the findings from the previous section ,
let us now assume that the limit @xmath196 from one of the graph structures results in a measured metric space @xmath199 of dimension @xmath200 $ ] , in which sphere volumes are given by the equation @xmath201 where @xmath202 is the volume of a unit sphere , @xmath203 denoting the gamma function .
thus , we assume that @xmath204 interpolates the sphere volumes of euclidean spaces for integer @xmath205 .
note that this assumption has indeed two parts .
the first , , means that a volume measure on the metric space @xmath199 exists that behaves homogeneously with degree @xmath205 with regard to distances . in the manifold case ( integer @xmath205 ) , this is the case of vanishing curvature ; for general manifolds of integer dimension @xmath205 , would hold as an approximation for small radii .
the second assumption , , corresponds to the euclideanness of the metric . for edge - weighted pixel graphs based on @xmath2- or @xmath3-neighbourhoods
, the volume of unit spheres actually deviates from , even in the limit .
however , with increasing neighbourhood size , is approximated better and better .
most of the following analysis does not depend specifically on ; thus we will return to only later for numerical evaluation of information functionals .
with we have @xmath206 where we have substituted @xmath207 . as a result , we obtain @xmath208 analogous considerations for @xmath92 from lead to @xmath209 as pointed out before , the metric structure of the fractal @xmath199 will in general be more complicated such that it does not possess a well - defined global dimension .
however , such a dimension can be measured at each location and scale .
the quantities @xmath85 and @xmath92 as stated in , can then be understood as functions of the local fractal dimension in a neighbourhood of vertex @xmath67 where the size of the neighbourhood  the scale  is controlled by the decay of the function @xmath210 in the integrands of and , respectively .    as a result
, we find that the information functionals @xmath85 and @xmath92 represent distribution over the input pixels of an image patch ( non - adaptive or adaptive ) in which the pixels are assigned different weights dependent on a local fractal dimension measure .
the entropies @xmath48 and @xmath49 then measure the local homogeneity or inhomogeneity of this dimension distribution : for very homogeneous dimension values within a patch , the density resulting from each of the information functionals @xmath85 , @xmath92 will be fairly homogeneous , implying high entropy .
the more the dimension values are spread out , the more will the density be dominated by a few pixels with high values of @xmath85 or @xmath92 , thus yielding low entropy .
the precise dependency of the entropy on the dimension distribution will be slightly different for @xmath85 and @xmath92 and will also depend on the choice of @xmath211 .
details of this dependency will be a topic of future work .
0.001    ( 1000 , 505 ) ( 5,500 ) ( 0,370 ) ( 190,260)@xmath205 ( 15,260)(a ) ( 333,500 ) ( 338,370 ) ( 525,260)@xmath205 ( 350,260)(b ) ( 661,500 ) ( 666,370 ) ( 853,260)@xmath205 ( 680,260)(c ) ( 5,240 ) ( 0,110 ) ( 190 , 0)@xmath205 ( 15 , 0)(d ) ( 333,240 ) ( 338,110 ) ( 525 , 0)@xmath205 ( 350 , 0)(e ) ( 661,240 ) ( 666,110 ) ( 853 , 0)@xmath205 ( 680 , 0)(f )    to give a basic intuition , we present in figure  [ fi - fvfp ] graphs of @xmath212 and @xmath213 as functions of the dimension @xmath200 $ ] for selected values of @xmath211 . in computing these values , the specific choice for @xmath214 has been used .    in the left column ,
frames ( a ) and ( b ) , we use @xmath157 as in our experiments in section  [ sec - segex ] . here , both @xmath212 and @xmath213 increase drastically by almost @xmath215 from @xmath216 to @xmath217 and even by @xmath218 from @xmath217 to @xmath219 . in the resulting information functionals @xmath85 and @xmath92 , i.e.  after applying @xmath220 to the functions shown in the figure , even pixels with only slightly higher values of the dimension strongly dominate the entire information density within the patch .    for increasing @xmath211 ,
the rate of increment in @xmath85 and @xmath92 with @xmath205 becomes lower . for @xmath221 as shown in the second column , frames ( b ) and ( e ) , of figure  [ fi - fvfp ] , the variation of @xmath212 and @xmath213 is already reduced to @xmath5 and @xmath2 , respectively , such that vertices across the entire dimension range @xmath222 $ ] will have a relevant influence on the information density .
for even larger @xmath211 the dependency of @xmath85 and @xmath92 on @xmath205 becomes non - monotonic ( as shown in ( c ) for @xmath85 with @xmath223 ) and even monotonically decreasing ( for both @xmath85 and @xmath92 at @xmath224 ; not shown ) .
it will therefore be interesting for further investigation to evaluate also the texture discrimination behaviour of the entropy measures for varying @xmath211 , as this may be a way to targeting the sensitivity of the measures specifically at certain dimension ranges .    in this context
, however , it becomes evident that the parameter @xmath211 plays two different roles at the same time .
first , it steers the approximate radius of influence for the fractal dimension estimation . here
, it is important that this radius of influence is smaller than the patch size underlying the graph construction , such that the cut - off of the graphs has no significant influence on the values of the information functional at the individual vertices .
second , @xmath211 determines the shape and steepness of the function ( compare figure  [ fi - fvfp ] ) that relates the local fractal dimension to the values of the information functionals .
this makes it desirable to refine in future work the parametrisation of the exponential weighting scheme so that the two roles of @xmath211 are distributed to two parameters .
in this paper , we have presented the framework of graph - index - based texture descriptors that has first been introduced in @xcite .
particular emphasis was put on entropy - based graph indices that have proven in @xcite to afford medium to high sensitivity for texture differences .
we have extended the work from @xcite in two directions .
firstly , we have stated an approach to texture - based image segmentation in which the texture descriptor framework was integrated with geodesic active contours @xcite , a standard method for intensity - based image segmentation .
this approach was already briefly introduced in @xcite and is demonstrated here by a larger set of experiments , including two synthetic and one real - world example .
secondly , we have analysed one representative of the graph - entropy - based texture descriptors in order to gain insight about the image properties that this descriptor relies on .
it turned out that it stands in close relation to measurements of fractal dimension of certain metric spaces that arise from the graphs in local image patches that underly our texture descriptors .
although this type of fractal dimension measurement in images differs from existing applications of fractal theory in image ( and texture ) analysis , as the latter treat the image manifold as fractal object , results indicate that the two fractal approaches are related .
our texture descriptor framework as a whole and also both novel contributions presented here require further research . to mention some topics , we start with parameter selection of the texture descriptors . in @xcite as well as in this paper
, most parameters were fixed to specific values based on heuristics
. a systematic investigation of the effect of variations of all these parameters is part of ongoing work .
inclusion of novel graph descriptors proposed in the literature , e.g.@xcite , is a further option .
the algorithms currently used for the computation of graph - entropy - based texture descriptors need computation times in the range of minutes already for small images . as the algorithms have not been designed for efficiency
, there is much room for improvement which will also be considered in future work .
both texture discrimination and texture segmentation have been demonstrated so far rather on a proof - of - concept level .
extensive evaluation on larger sets of image data is ongoing .
this is also necessary to gain more insight about the suitability of particular texture descriptors from our set for specific classes of textures .    regarding the texture segmentation framework ,
the conceptual break between the graph - based set of texture descriptors and the partial differential equation for segmentation could be reduced by using e.g.  a graph - cut segmentation method .
it can be asked whether such a combination even allows for some synergy between the computation steps .
this is not clear so far since the features used to weight graph edges are different : intensity contrasts in the texture descriptor phase ; texture descriptor differences in the graph - cut phase . in further course
, the integration of graph - entropy - based texture descriptors into more complex segmentation frameworks will be a goal .
unsupervised segmentation approaches are not capable to handle involved segmentation tasks ( like in medical diagnostics ) where highly accurate segmentation can only be achieved by including prior information on the shape and appearance of the objects to be segmented .
state - of - the - art segmentation frameworks therefore combine the mechanisms of unsupervised segmentation approaches with model - based methods as introduced e.g.  in @xcite .    on the theoretical side , the analysis of the fractal limit of the descriptor construction will have to be refined and extended to include all six graph settings from section  [ ssec - graphcons ] .
relations between the fractal structures arising from the graph construction and the more image manifold more commonly treated in fractal - based image analysis will have to be analysed .
generally , much more theoretical work deserves to be invested in understanding the connections and possible equivalences between the very disparate approaches to texture description that can be found in the literature
. a graph - based approach like ours admits different directions of such comparisons .
it can thus be speculated that it could play a pivotal role in understanding the relations between texture description methods , and create a unifying view on different methods that would also have implications for the understanding of texture itself .
y.  boykov , o.  veksler , and r.  zabih .
markov random fields with efficient approximation . in _ proc .
1998 ieee computer society conference on computer vision and pattern recognition _ , pages 648655 , santa barbara , ca , june 1998 .
a.  buades , b.  coll , and j.  morel . a non - local algorithm for image denoising . in _ proc .
2005 ieee computer society conference on computer vision and pattern recognition _ , volume  2 , pages 6065 , san diego , ca , june 2005 .
ieee computer society press .
d.  comaniciu and p.  meer . robust analysis of feature spaces : color image segmentation . in _ proc .
1997 ieee computer society conference on computer vision and pattern recognition _ , pages 750755 , san juan , puerto rico , june 1997 .
ieee computer society press .
a.  elmoataz , o.  lezoray , v .- t .
ta , and s.  bougleux .
partial difference equations on graphs for local and nonlocal image processing .
in o.  lezoray and l.  grady , editors , _ image processing and analysis with graphs : theory and practice _ , chapter  7 , pages 174206 .
crc press , boca raton , 2012 .
m.  ferrer and h.  bunke .
graph edit distance  theory , algorithms , and applications . in o.
lezoray and l.  grady , editors , _ image processing and analysis with graphs : theory and practice _ , chapter  13 , pages 383422 .
crc press , boca raton , 2012 .
b.  georgescu , i.  shimshoni , and p.  meer .
mean shift based clustering in high dimensions : a texture classification example . in _ proc .
2003 ieee international conference on computer vision _ ,
volume  1 , pages 456463 , nice , oct . 2003 .
h.  ishikawa .
graph cuts  combinatorial optimization in vision . in o.
lezoray and l.  grady , editors , _ image processing and analysis with graphs : theory and practice _ , chapter  2 , pages 2563 .
crc press , boca raton , 2012 .
s.  kichenassamy , a.  kumar , p.
olver , a.  tannenbaum , and a.  yezzi .
gradient flows and geometric active contour models . in _ proc .
fifth international conference on computer vision _ , pages 810815 , cambridge , ma , june 1995 .
ieee computer society press .
r.  lerallut ,  .
decencire , and f.  meyer .
image processing using morphological amoebas . in c.
ronse , l.  najman , and e.  decencire , editors , _ mathematical morphology : 40 years on _ , volume  30 of _ computational imaging and vision _ , pages 1322 .
springer , dordrecht , 2005 .
o.  lezoray and l.  grady .
graph theory concepts and definitions used in image processing and analysis . in o.
lezoray and l.  grady , editors , _ image processing and analysis with graphs : theory and practice _ , chapter  1 , pages 124 .
crc press , boca raton , 2012 .
l.  lucchese and s.  k. mitra .
unsupervised segmentation of color images based on k - means clustering in the chromaticity plane . in _ proc .
ieee workshop on content - based access of image and video libraries _ ,
pages 7478 , fort collins , co , 1999 .
p.  maragos .
fractal aspects of speech signals : dimension and interpolation . in _ proc .
1991 international conference on acoustics , speech , and signal processing _
, volume  1 , pages 417420 , toronto , ontario , canada , may 1991 .
l.  najman and f.  meyer . a short tour of mathematical morphology on edge and vertex weighted graphs . in o.
lezoray and l.  grady , editors , _ image processing and analysis with graphs : theory and practice _ , chapter  6 , pages 141173 .
crc press , boca raton , 2012 .
b.  sandberg , t.  chan , and l.  vese .
a level - set and gabor - based active contour algorithm for segmenting textured images .
technical report cam-02 - 39 , department of mathematics , university of california at los angeles , ca , u.s.a . ,
july 2002 .
m.  welk .
discrimination of image textures using graph indices . in m.
dehmer and f.  emmert - streib , editors , _ quantitative graph theory : mathematical foundations and applications _ , chapter  12 , pages 355386 .
crc press , 2014 .
m.  welk and m.  breu .
morphological amoebas and partial differential equations . in p.
w. hawkes , editor , _ advances in imaging and electron physics _ ,
volume 185 , pages 139212 .
elsevier academic press , 2014 .
l.  zhu , w.  ng , and s.  han . classifying graphs using theoretical metrics : a study of feasibility . in j.
xu , g.  yu , s.  zhou , and r.  unland , editors , _ database systems for advanced applications _ , volume 6637 of _ lecture notes in computer science _ , pages 5364 .
springer , berlin , 2011 .