Page 96 - Geochemical Anomaly and Mineral Prospectivity Mapping in GIS
P. 96

Fractal Analysis of Geochemical Anomalies                             95





















           Fig. 4-7. Log-log plot of concentration-area relationship for the soil Fe data. By careful inspection
           of the plots,  three inflection points (or thresholds) corresponding to 1.6%, 7.2% and 8.6% Fe
           (pointed by short arrows) can be defined and, thus, four straight lines can be fitted by least squares
           through the concentration-area plots. The individual straight lines fitted through the plots to the
           left of any threshold satisfy the power-law relation in equation (4.6), whilst the straight line fitted
           through the plots to the right of the rightmost threshold satisfies the power-law relation in equation
           (4.7).



           Generation and discretisation of geochemical surfaces

              Because most  exploration geochemical data are recorded as attributes of sampling
           points, point-to-surface transformations through spatial interpolation (see Chapter 2) are
           essential in analysis of geochemical thresholds via the concentration-area fractal method.
           Spatial interpolation entails analysis of spatial correlation and variability of point
           geochemical  data in  order to determine the precise  way in  which to generate  a
           geochemical surface for a certain element. This topic is, however, beyond the scope of
           this volume. The  references  to this topic cited earlier in  this chapter and many other
           relevant publications can be consulted for further details.
              Most GIS software  packages support spatial interpolation of point  data via either
           triangulation or gridding techniques. The former techniques are not appropriate but the
           latter techniques are appropriate for the application of the concentration-area  fractal
           method.  That is because  interpolation  via gridding techniques  provide surfaces
           represented as a raster of pixels (see Figs. 2-5 and 2-14), which are amenable to the box-
           counting method for estimation of areas enclosed by certain uni-element concentration
           levels (Fig. 4-1). A  geochemical surface generated  from point data  set  must then be
           discretised  or  classified according to some intervals of the data. Classification is, of
           course, a basic functionality of a GIS (see Chapter 2).
   91   92   93   94   95   96   97   98   99   100   101