Classification of Sweden's Forest and Alpine Vegetation Using Optical Satellite and Inventory Data

Creation of accurate vegetation maps from optical satellite data requires use of reference data to aid in interpretation or to verify map results. Reference data may be taken, for example, from field visits, aerial photo-interpretation, or ground-based inventories. National inventories are a potential source of reference data useful in land cover mapping projects. This thesis addresses aspects of mapping forest and alpine vegetation in Sweden through combined use of optical satellite data and inventory data. Issues such as reference and satellite data pre-processing, spatial scale, quantity and quality of reference data, and classification methods have been examined. Optical satellite data with pixel sizes ranging from 10 to 300 m have been used together with reference data from the Swedish National Forest Inventory (NFI), National Inventory of Landscapes in Sweden (NILS), a point sample based on the Terrestrial Habitat Monitoring program (THUF), and a forest stand database. Results include modifications to common remote sensing methods, such as introducing iterative adjustment of prior probabilities in Maximum Likelihood classification, and improved topographic normalization (C-correction) of satellite data. Probability-based samples such as NFI, NILS and THUF provide data necessary for assignment of prior probabilities, estimation of continuous values, and are useful as training and validation data. For managed boreal forest stands, coarser pixel (60 m) AWiFS data were nearly as effective for stem volume estimation as SPOT 5 data (10 m). On the other hand, the most accurate classification of detailed alpine vegetation types (72.9% overall accuracy) was from SPOT 5 data combined with elevation derivatives, while classifications of Landsat TM (25 m), AWiFS, and MERIS (300 m) were less accurate. Non-parametric methods (e.g., random forests, decision/regression trees) produced higher classification accuracies than traditional parametric methods for alpine vegetation. The quantity of reference data affected classification accuracy, as more reference data produced higher map accuracy, although other factors such as distribution and quality of the reference data should be considered. As seen in this thesis, the characteristics of the landscape exert an influence on satellite and training data requirements, classification methods and resulting map accuracy.

[1]  Eric D. Kolaczyk,et al.  On the choice of spatial and categorical scale in remote sensing land cover classification , 2005 .

[2]  Alejandro C. Frery,et al.  The influence of training errors, context and number of bands in the accuracy of image classification , 2009 .

[3]  Leo Breiman,et al.  Bagging Predictors , 1996, Machine Learning.

[4]  C. Brodley,et al.  Decision tree classification of land cover from remotely sensed data , 1997 .

[5]  Michael A. Wulder,et al.  Validation of a large area land cover product using purpose-acquired airborne video , 2007 .

[6]  T. Remmel,et al.  Effect of configuration on spectral signatures , 2010 .

[7]  M. Nilsson,et al.  Combining national forest inventory field plots and remote sensing data for forest databases , 2008 .

[8]  Giles M. Foody,et al.  Mapping a specific class with an ensemble of classifiers , 2007 .

[9]  V. Lakshmi Observation of the Earth and Its Environment—Survey of Missions and Sensors , 2003 .

[10]  Qihao Weng,et al.  A survey of image classification methods and techniques for improving classification performance , 2007 .

[11]  Jonas Ardö,et al.  Volume quantification of coniferous forest compartments using spectral radiance recorded by Landsat Thematic Mapper , 1992 .

[12]  Eric W. T. Ngai,et al.  Customer churn prediction using improved balanced random forests , 2009, Expert Syst. Appl..

[13]  Steven W. Running,et al.  Remote sensing of temperate coniferous forest leaf area index The influence of canopy closure, understory vegetation and background reflectance , 1990 .

[14]  Ulf. Söderberg,et al.  Funktioner för skogliga produktionsprognoser , 1986 .

[15]  J. Townshend,et al.  Global land cover classi(cid:142) cation at 1 km spatial resolution using a classi(cid:142) cation tree approach , 2004 .

[16]  BIRGER,et al.  Comparison of satellite imagery and infrared aerial photography as vegetation mapping methods in an arctic study area; Jameson Land, East Greenland , 2007 .

[17]  F. J. Ahern,et al.  Forestry information content of Thematic Mapper data , 1986 .

[18]  R. D. Ramsey,et al.  Landscape Cover-Type Modeling Using a Multi-Scene Thematic Mapper Mosaic , 1997 .

[19]  Xiaofeng Li,et al.  Improved land cover mapping using random forests combined with landsat thematic mapper imagery and ancillary geographic data. , 2010 .

[20]  Sarah Parks,et al.  An effective assessment protocol for continuous geospatial datasets of forest characteristics using USFS Forest Inventory and Analysis (FIA) data , 2010 .

[21]  H. Olsson Changes in satellite-measured reflectances caused by thinning cuttings in Boreal forest , 1994 .

[22]  Leo Breiman,et al.  Classification and Regression Trees , 1984 .

[23]  Laurence Hubert-Moy,et al.  A Comparison of Parametric Classification Procedures of Remotely Sensed Data Applied on Different Landscape Units , 2001 .

[24]  Joanne C. White,et al.  Monitoring Canada’s forests. Part 1: Completion of the EOSD land cover project , 2008 .

[25]  David M. Johnson A Comparison of Coincident Landsat-5 TM and Resourcesat-1 AWiFS Imagery for Classifying Croplands , 2008 .

[26]  A. Gillespie,et al.  Topographic Normalization of Landsat TM Images of Forest Based on Subpixel Sun–Canopy–Sensor Geometry , 1998 .

[27]  G. Dedieu,et al.  SMAC: a simplified method for the atmospheric correction of satellite measurements in the solar spectrum , 1994 .

[28]  C. Woodcock,et al.  The use of variograms in remote sensing. I - Scene models and simulated images. II - Real digital images , 1988 .

[29]  J. Chen Spatial Scaling of a Remotely Sensed Surface Parameter by Contexture , 1999 .

[30]  Robert H. Fraser,et al.  Mapping northern land cover fractions using Landsat ETM , 2007 .

[31]  A. Winsor Sampling techniques. , 2000, Nursing times.

[32]  Richard A. Wadsworth,et al.  Final Report for LCM2007 - the new UK land cover map. Countryside Survey Technical Report No 11/07 , 2011 .

[33]  Hermann Kaufmann,et al.  Comparison of Topographic Correction Methods , 2009, Remote. Sens..

[34]  F. M. Danson,et al.  Preliminary evaluation of the relationships between SPOT-1 HRV data and forest stand parameters , 1987 .

[35]  M. Wijk,et al.  Using Information Theory to Determine Optimum Pixel Size and Shape for Ecological Studies: Aggregating Land Surface Characteristics in Arctic Ecosystems , 2009, Ecosystems.

[36]  Michael A. Wulder,et al.  Remote sensing methods in medium spatial resolution satellite data land cover classification of large areas , 2002 .

[37]  Leo Breiman,et al.  Random Forests , 2001, Machine Learning.

[38]  Keiji Kushida,et al.  Spectral vegetation indices for estimating shrub cover, green phytomass and leaf turnover in a sedge‐shrub tundra , 2008 .

[39]  Ramakrishna R. Nemani,et al.  A remote sensing based vegetation classification logic for global land cover analysis , 1995 .

[40]  E. Small,et al.  The impact of soil reflectance on the quantification of the green vegetation fraction from NDVI , 2005 .

[41]  R. Latifovic,et al.  Development of a circa 2000 land cover map of northern Canada at 30 m resolution from Landsat , 2009 .

[42]  D. Civco Topographic normalization of landsat thematic mapper digital imagery , 1989 .

[43]  Philip J. Howarth,et al.  High Spatial Resolution Remote Sensing Data for Forest Ecosystem Classification: An Examination of Spatial Scale , 2000 .

[44]  K. Davis,et al.  The MODIS (Collection V005) BRDF/albedo product: Assessment of spatial representativeness over forested landscapes , 2009 .

[45]  C. Woodcock,et al.  The factor of scale in remote sensing , 1987 .

[46]  Erik Næsset,et al.  Using remotely sensed data to construct and assess forest attribute maps and related spatial products , 2010 .

[47]  C. Goulding,et al.  Estimation of timber volume in a coniferous plantation forest using Landsat TM , 1997 .

[48]  Damiano Gianelle,et al.  Ecosystem carbon fluxes and canopy spectral reflectance of a mountain meadow , 2009 .

[49]  B. Datt,et al.  On the relationship between training sample size and data dimensionality: Monte Carlo analysis of broadband multi-temporal classification , 2005 .

[50]  Conghe Song,et al.  Forest mapping with a generalized classifier and Landsat TM data , 2001 .

[51]  Janet Franklin,et al.  Thematic mapper analysis of coniferous forest structure and composition , 1986 .

[52]  Joanne C. White,et al.  A new data fusion model for high spatial- and temporal-resolution mapping of forest disturbance based on Landsat and MODIS , 2009 .

[53]  P. Mather,et al.  Classification Methods for Remotely Sensed Data , 2001 .

[54]  Thomas M. Lillesand,et al.  Statewide land cover derived from multiseasonal Landsat TM data: A retrospective of the WISCLAND project , 2002 .

[55]  Estimating annual cuttings using multi-temporal satellite data and field data from the Swedish NFI , 2009 .

[56]  Joanne C. White,et al.  Multiscale satellite and spatial information and analysis framework in support of a large-area forest monitoring and inventory update , 2010, Environmental monitoring and assessment.

[57]  John P. Wilson,et al.  Multi-scale linkages between topographic attributes and vegetation indices in a mountainous landscape , 2007 .

[58]  Alan H. Strahler,et al.  On the nature of models in remote sensing , 1986 .

[59]  Peng Gong,et al.  An assessment of some factors influencing multispectral land-cover classification , 1990 .

[60]  Stein Rune Karlsen,et al.  Monitoring the spatio-temporal dynamics of geometrid moth outbreaks in birch forest using MODIS-NDVI data , 2009 .

[61]  Giles M. Foody,et al.  An evaluation of some factors affecting the accuracy of classification by an artificial neural network , 1997 .

[62]  Andrew K. Skidmore,et al.  Forest mapping accuracies are improved using a supervised nonparametric classifier with SPOT data , 1988 .

[63]  T. Gregoire Design-based and model-based inference in survey sampling: appreciating the difference , 1998 .

[64]  Hans Tømmervik,et al.  Monitoring vegetation changes in Pasvik (Norway) and Pechenga in Kola Peninsula (Russia) using multitemporal Landsat MSS/TM data , 2003 .

[65]  Giles M. Foody,et al.  Applications of the self-organising feature map neural network in community data analysis , 1999 .

[66]  W. Gould,et al.  Canadian Arctic vegetation mapping , 2002 .

[67]  Pol Coppin,et al.  Satellite inventory of Minnesota forest resources , 1994 .

[68]  Håkan Olsson,et al.  Combining optical satellite data and airborne laser scanner data for vegetation classification , 2012 .

[69]  Comparison of satellite imagery and infrared aerial photography as vegetation mapping methods in an Arctic study area; Jameson Land, East Greenland , 1994 .

[70]  Jin Chen,et al.  Analysis of NDVI and scaled difference vegetation index retrievals of vegetation fraction , 2006 .

[71]  T. M. Lillesand,et al.  Remote Sensing and Image Interpretation , 1980 .

[72]  C. Woodcock,et al.  Classification and Change Detection Using Landsat TM Data: When and How to Correct Atmospheric Effects? , 2001 .

[73]  L. Pedroni Improved classification of Landsat Thematic Mapper data using modified prior probabilities in large and complex landscapes , 2003 .

[74]  A. Huete,et al.  Normalization of multidirectional red and NIR reflectances with the SAVI , 1992 .

[75]  Stephen V. Stehman,et al.  Design and Analysis for Thematic Map Accuracy Assessment: Fundamental Principles , 1998 .

[76]  D. Lu The potential and challenge of remote sensing‐based biomass estimation , 2006 .

[77]  Craig A. Coburn,et al.  SCS+C: a modified Sun-canopy-sensor topographic correction in forested terrain , 2005, IEEE Transactions on Geoscience and Remote Sensing.

[78]  P. Aplin,et al.  On scales and dynamics in observing the environment , 2006 .

[79]  Darrel L. Williams,et al.  Historical record of Landsat global coverage: mission operations, NSLRSDA, and International Cooperator stations , 2006 .

[80]  Ronald J. P. Lyon,et al.  Influence of rock-soil spectral variation on the assessment of green biomass , 1985 .

[81]  Peter M. Atkinson,et al.  Spatial variation in land cover and choice of spatial resolution for remote sensing , 2004 .

[82]  Tiit Nilson,et al.  Radiative Transfer in Nonhomogeneous Plant Canopies , 1992 .

[83]  Björn Nilsson,et al.  National Inventory of Landscapes in Sweden (NILS)—scope, design, and experiences from establishing a multiscale biodiversity monitoring system , 2011, Environmental monitoring and assessment.

[84]  C. Justice,et al.  The topographic effect on spectral response from nadir-pointing sensors , 1980 .

[85]  C. Justice,et al.  Characterization and classification of South American land cover types using satellite data , 1987 .

[86]  Peter M. Atkinson,et al.  Defining an optimal size of support for remote sensing investigations , 1995, IEEE Trans. Geosci. Remote. Sens..

[87]  T. Lin,et al.  The Lambertian assumption and Landsat data. , 1980 .

[88]  Paul J. Curran,et al.  Simple size for ground and remotely sensed data , 1986 .

[89]  Michael A. Wulder,et al.  An accuracy assessment framework for large‐area land cover classification products derived from medium‐resolution satellite data , 2006 .

[90]  J. Sheail,et al.  The Land of Britain, 1930-1990: a comparative study of Field Mapping and Remote Sensing Techniques , 1994 .

[91]  Paul M. Mather,et al.  Support vector machines for classification in remote sensing , 2005 .

[92]  Mark H. Hansen,et al.  Using a land cover classification based on satellite imagery to improve the precision of forest inventory area estimates , 2002 .

[93]  Russell Congalton,et al.  Assessing the Accuracy of Remotely Sensed Data: Principles and Practices, Second Edition , 1998 .

[94]  P. Atkinson,et al.  Introduction Neural networks in remote sensing , 1997 .

[95]  R. Clarke,et al.  National-scale vegetation change across Britain; an analysis of sample-based surveillance data from the Countryside Surveys of 1990 and 1998. , 2003, Journal of environmental management.

[96]  Variable Selection Strategies for Small-area Estimation Using FIA Plots and Remotely Sensed Data , 2005 .

[97]  Giles M. Foody,et al.  The use of small training sets containing mixed pixels for accurate hard image classification: Training on mixed spectral responses for classification by a SVM , 2006 .

[98]  Pierre-Yves Henry,et al.  Habitat monitoring in Europe: a description of current practices , 2008, Biodiversity and Conservation.

[99]  P. Teillet,et al.  On the Slope-Aspect Correction of Multispectral Scanner Data , 1982 .

[100]  M. Nilsson,et al.  Countrywide Estimates of Forest Variables Using Satellite Data and Field Data from the National Forest Inventory , 2003, Ambio.

[101]  Jan Åge Riseth,et al.  Above ground biomass changes in the mountain birch forests and mountain heaths of Finnmarksvidda, northern Norway, in the period 1957-2006 , 2009 .

[102]  Normalisation of within-scene optical depth levels in multispectral satellite imagery using National Forest Inventory plot data. , 2005 .

[103]  L G Firbank,et al.  Changing landscapes, habitats and vegetation diversity across Great Britain. , 2003, Journal of environmental management.

[104]  G. Staahl Optimizing the utility of forest inventory activities , 1994 .

[105]  Erik Næsset,et al.  Advances and emerging issues in national forest inventories , 2010 .

[106]  D. Stow,et al.  Reflectance characteristics of arctic tundra vegetation from airborne radiometry , 1993 .

[107]  T. W. Berge,et al.  Modelling biomass and leaf area index in a sub-arctic Scandinavian mountain area , 2004 .

[108]  Limin Yang,et al.  COMPLETION OF THE 1990S NATIONAL LAND COVER DATA SET FOR THE CONTERMINOUS UNITED STATES FROM LANDSAT THEMATIC MAPPER DATA AND ANCILLARY DATA SOURCES , 2001 .

[109]  D. Stow,et al.  THE EFFECT OF TRAINING STRATEGIES ON SUPERVISED CLASSIFICATION AT DIFFERENT SPATIAL RESOLUTIONS , 2002 .

[110]  R. McRoberts,et al.  Comparisons of National Forest Inventories , 2010 .

[111]  Ronald E. McRoberts,et al.  Probability- and model-based approaches to inference for proportion forest using satellite imagery as ancillary data , 2010 .

[112]  Stephen V. Stehman,et al.  Introduction to special issue on map accuracy , 2003, Environmental and Ecological Statistics.

[113]  Tiit Nilson,et al.  Age dependence of forest reflectance: Analysis of main driving factors , 1994 .

[114]  Janne Heiskanen,et al.  Assessment of multispectral, -temporal and -angular MODIS data for tree cover mapping in the tundra-taiga transition zone , 2008 .

[115]  A. Strahler,et al.  Geometric-Optical Bidirectional Reflectance Modeling of a Conifer Forest Canopy , 1986, IEEE Transactions on Geoscience and Remote Sensing.

[116]  O. Hagner,et al.  A method for calibrated maximum likelihood classification of forest types , 2007 .

[117]  Mark A. Friedl,et al.  Using prior probabilities in decision-tree classification of remotely sensed data , 2002 .

[118]  Juha Hyyppä,et al.  Effects of stand size on the accuracy of remote sensing-based forest inventory , 2001, IEEE Trans. Geosci. Remote. Sens..

[119]  Taskin Kavzoglu,et al.  Increasing the accuracy of neural network classification using refined training data , 2009, Environ. Model. Softw..

[120]  Donald A. Walker,et al.  The Circumpolar Arctic vegetation map , 2005 .

[121]  Christopher Small,et al.  The Landsat ETM+ spectral mixing space , 2004 .

[122]  Rasim Latifovic,et al.  Selecting Representative High Resolution Sample Images for Land Cover Studies. Part 2 Application to Estimating Land Cover Composition , 2000, Remote Sensing of Environment.

[123]  M. Pal,et al.  Random forests for land cover classification , 2003, IGARSS 2003. 2003 IEEE International Geoscience and Remote Sensing Symposium. Proceedings (IEEE Cat. No.03CH37477).

[124]  Alan H. Strahler,et al.  The Use of Prior Probabilities in Maximum Likelihood Classification , 1980 .

[125]  Matt Aitkenhead,et al.  Automating land cover mapping of Scotland using expert system and knowledge integration methods , 2011 .

[126]  Yuhong He,et al.  Landsat-comparable land cover maps using ASTER and SPOT images: a case study for large-area mapping programmes , 2011 .

[127]  Paul Treitz,et al.  Remote Sensing of Arctic Vegetation: Relations between the NDVI, Spatial Resolution and Vegetation Cover on Boothia Peninsula, Nunavut , 2009 .

[128]  C. Woodcock,et al.  The use of variograms in remote sensing: I , 1988 .

[129]  F. Hall,et al.  Improved topographic correction of forest image data using a 3‐D canopy reflectance model in multiple forward mode , 2008 .

[130]  J. Townshend,et al.  A new global 1‐km dataset of percentage tree cover derived from remote sensing , 2000 .

[131]  Frédéric Achard,et al.  GLOBCOVER : The most detailed portrait of Earth , 2008 .

[132]  Stephen V. Stehman,et al.  Statistical Rigor and Practical Utility in Thematic Map Accuracy Assessment , 2001 .

[133]  Michael D. Bankier,et al.  Power Allocations: Determining Sample Sizes for Subnational Areas , 1988 .

[134]  G. B. Groom,et al.  A standardized procedure for surveillance and monitoring European habitats and provision of spatial data , 2007, Landscape Ecology.

[135]  M. Nilsson Estimation of forest variables using satellite image data and airborne Lidar , 1997 .

[136]  R. Fournier,et al.  Remote sensing and the measurement of geographical entities in a forested environment. 2. The optimal spatial resolution , 1994 .

[137]  J. Wickham,et al.  Impacts of Patch Size and Land-Cover Heterogeneity on Thematic Image Classification Accuracy , 2002 .

[138]  J. Dash,et al.  The MERIS terrestrial chlorophyll index , 2004 .

[139]  Scott Mitchell,et al.  Integration of forest inventory and satellite imagery:: a Canadian status assessment and research issues , 2005 .

[140]  S. Ekstrand,et al.  Assessment of forest damage with landsat TM: Correction for varying forest stand characteristics , 1994 .

[141]  J. Wickham,et al.  Completion of the 2001 National Land Cover Database for the conterminous United States , 2007 .

[142]  A. Belward,et al.  GLC2000: a new approach to global land cover mapping from Earth observation data , 2005 .

[143]  Damien Sulla-Menashe,et al.  MODIS Collection 5 global land cover: Algorithm refinements and characterization of new datasets , 2010 .

[144]  J. Townshend,et al.  Towards an operational MODIS continuous field of percent tree cover algorithm: examples using AVHRR and MODIS data , 2002 .

[145]  Giles M. Foody,et al.  Status of land cover classification accuracy assessment , 2002 .

[146]  John E. Estes,et al.  Institutional issues affecting the integration and use of remotely sensed data and geographic information systems , 1991 .

[147]  J. Heiskanen Estimating aboveground tree biomass and leaf area index in a mountain birch forest using ASTER satellite data , 2006 .

[148]  R. McRoberts,et al.  Remote sensing support for national forest inventories , 2007 .

[149]  Joanne C. White,et al.  Supporting large-area, sample-based forest inventories with very high spatial resolution satellite imagery , 2009 .

[150]  R. Fernandes,et al.  Approaches to fractional land cover and continuous field mapping: A comparative assessment over the BOREAS study region , 2004 .

[151]  Mahesh Pal,et al.  Random forest classifier for remote sensing classification , 2005 .

[152]  R. Hill,et al.  The UK Land Cover Map 2000: Construction of a Parcel-Based Vector Map from Satellite Images , 2002 .

[153]  Michael J. Choate,et al.  Effects of Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper plus radiometric and geometric calibrations and corrections on landscape characterization , 2001 .

[154]  E. Tomppo National Forest Inventories : pathways for common reporting , 2010 .

[155]  Zhanqing Li,et al.  The bidirectional effects of AVHRR measurements over boreal regions , 1996, IEEE Trans. Geosci. Remote. Sens..

[156]  Håkan Olsson,et al.  A method for using Landsat time series for monitoring young plantations in boreal forests , 2009 .

[157]  Ronald E. McRoberts,et al.  Comprar National Forest Inventories · Pathways for Common Reporting | Tomppo, Erkki | 9789048132324 | Springer , 2010 .

[158]  M. S. Moran,et al.  Bidirectional reflectance factors of agricultural targets: A comparison of ground-, aircraft-, and satellite-based observations , 1990 .

[159]  Carla E. Brodley,et al.  Identifying Mislabeled Training Data , 1999, J. Artif. Intell. Res..

[160]  Zhiliang Zhu,et al.  Combining Accuracy Assessment of Land-Cover Maps with Environmental Monitoring Programs , 2000 .

[161]  Paul M. Mather,et al.  An assessment of the effectiveness of decision tree methods for land cover classification , 2003 .

[162]  David L. Verbyla,et al.  Conservative bias in classification accuracy assessment due to pixel-by-pixel comparison of classified images with reference grids , 1995 .

[163]  F. Chapin,et al.  Subgrid-scale variability in the surface energy balance , 1998 .

[164]  Tiit Nilson,et al.  Successional reflectance trajectories in northern temperate forests , 1993 .

[165]  Ranga B. Myneni,et al.  Remote sensing of vegetation and land-cover change in Arctic Tundra Ecosystems , 2004 .

[166]  C. Woodcock,et al.  Theory and methods for accuracy assessment of thematic maps using fuzzy sets , 1994 .

[167]  Antoine Guisan,et al.  Evaluating effects of spectral training data distribution on continuous field mapping performance , 2009 .

[168]  Chao Chen,et al.  Using Random Forest to Learn Imbalanced Data , 2004 .

[169]  Jesslyn F. Brown,et al.  Development of a land-cover characteristics database for the conterminous U.S. , 1991 .

[170]  Stein Rune Karlsen,et al.  Monitoring vegetation changes on Finnmarksvidda, Northern Norway, using Landsat MSS and Landsat TM/ETM+ satellite images , 2005 .

[171]  Ronald J. Hall,et al.  Operational mapping of the land cover of the forested area of Canada with Landsat data: EOSD land cover program , 2003 .

[172]  R. D. Ramsey,et al.  Mapping moderate-scale land-cover over very large geographic areas within a collaborative framework : A case study of the Southwest Regional Gap Analysis Project (SWReGAP) , 2007 .

[173]  J. Colby,et al.  Anisotropic reflectance correction of SPOT-3 HRV imagery , 2002 .

[174]  J. Fry,et al.  Updating the 2001 National Land Cover Database land cover classification to 2006 by using Landsat imagery change detection methods , 2009 .