Skip to content
BY 4.0 license Open Access Published by De Gruyter September 2, 2020

Land Use Land Cover map segmentation using Remote Sensing: A Case study of Ajoy river watershed, India

  • Kalyan Mahata , Rajib Das , Subhasish Das and Anasua Sarkar EMAIL logo

Abstract

Image segmentation in land cover regions which are overlapping in satellite imagery, is one crucial challenge. To detect true belonging of one pixel becomes a challenging problem while classifying mixed pixels in overlapping regions. In current work, we propose one new approach for image segmentation using a hybrid algorithm of K-Means and Cellular Automata algorithms. This newly implemented unsupervised model can detect cluster groups using hybrid 2-Dimensional Cellular-Automata model based on K-Means segmentation approach. This approach detects different land use land cover areas in satellite imagery by existing K-Means algorithm. Since it is a discrete dynamical system, cellular automaton realizes uniform interconnecting cells containing states. In the second stage of current model, we experiment with a 2-dimensional cellular automata to rank allocations of pixels among different land-cover regions. The method is experimented on the watershed area of Ajoy river (India) and Salinas (California) data set with true class labels using two internal and four external validity indices. The segmented areas are then compared with existing FCM, DBSCAN and K-Means methods and verified with the ground truth. The statistical analysis results also show the superiority of the new method.

MSC 2010: 68-XX; 92-XX; 62-XX

1 Introduction

Cogalton and Green in 1999 introduced Remote sensing to be a method forgetting knowledge of any object without direct physical contact with it. Canopy of methods exist for grouping pixels among predefined classes (like an urban area or turbid water) in satellite images. Following the theory, we can define a set for the remote sensing data set, as shown in

(1) P = { p i j k 1 i r , 1 j s , 1 k n }

of r × s × n dimensions in pixels, where pij (pij1, pij2, . . . pijk} denotes the values of n spectral bands for i, j th pixel. To find similar regions among overlapping segments, we partition the chosen watershed remote sensing image by our new hybrid unsupervised algorithm as well as another spatial image data set.

Let P (Rn or Zn) defines the image space of selected remote sensing imagery. So, points/pixels in P denotes spatial variables -x, y. Let dP (x, y) shows Euclidean distance between two pixels (x, y} ϵ P [15]. In the domain of spatial images, a crisp object C may be projected as a subset of P, C ⊆ P. We define a new method using the cellular automata approach over K-Means algorithm for image segmentation, which is the scope of this article.

In the unsupervised classification methods, clustering is based on maximum similarity within classes and minimum similarity among classes. State-of-the-art clustering methods are - self-organizing map algorithm [12], K-Means algorithm [13], simulated annealing [11], graph theoretic approach [17] and fuzzy c-means algorithm [5], as they are used several times for remote sensing image segmentation. Different perspective methods like clustering based on symmetry algorithm [7], also efficiently detect arbitrary shaped land cover regions in satellite images.

The proposed hybrid approach is performed to classify pixels in both-chosen remote sensing image Salinas A data set and chosen satellite image of the catchment region of river Ajoy. We also evaluate using external quantitative validity measures, our experimental results with solutions from K-Means and FCM algorithms on Salinas A data set with true class labels. The quantitative evaluation over existing two internal and four external validity indices shows the efficiency of the new hybrid algorithm. Using the new K-Means based and Cellular Automata (CA) corrected algorithm (KCA) to detect land cover clusters, which was verified with the ground truth knowledge. The statistical evaluation over both data sets also signify the proposed hybrid KCA approach over K-Means, DBSCAN and FCM methods.

2 Methods

2.1 K-Means Clustering Algorithm

K-Means method is a popular partition-based unsupervised clustering algorithm. James MacQueen proposed the term "k-means" [13] in 1967. In this approach, N objects {x1, x2, ..., xn} are to be clustered initially in K clusters, with random k objects as initial centroids ..., mk). Then the distances between an object xi and one centroid mj is calculated as d ( x i , m j ) = ( j = 1 k ( x i m j ) ) , where i = 1, 2........N. d(xi , mj) shows the distance between point xi and cluster centroid mj. The distance is calculated as Euclidean distance.

Based on the computed distances, one point is then divided in the cluster with minimum centroid distance. After re-allocation phase, the centroid of each K clusters are recalculated as m j = 1 C j i = 1 N x i , x i C j where j = 1, 2, ..., K. These steps repeat until Mean Square Error (MSE) E = i = 1 N j = 1 K ( x i m j ) , x i C j is minimized or until number of iterations becomes greater than the defined value. The advantage of K-means approach is that it is efficient and used for high-dimensional data with spherical shaped clusters.

2.2 Cellular Automata Method

Cellular Automata (called as CA) are explained to extend in the discrete spatial dynamical systems, which studies physical system models [18]. It evolves the computational devices in discrete space as well as in time. A Cellular Automaton is seeded with any state among states with all 0s and single 1 at separate positions. Then it generates unique patterns of fixed numbers.

Stephen Wolfram [19] proposes the simplest CA in a form of a spatial lattice of cells. The discrete variables are stored in cells. At time t, the variable refers to present state for that cell. Next state at time (t+1) is affected with present state and neighborhood states at time t. If 3-neighbourhood (self, left and right neighbors) Cellular Automaton is considered, with each having two states, either 0 or 1, then next state of cell is

(2) S i t + 1 = f ( S i 1 t , S i t , S i + 1 t ) ,

where f denotes the function of the next state. S i 1 t , S i + 1 t and S i t are the present states for left and right neighbors and self of i th cell at time t. The f is shown in the look-up table in Table 1. ‘Rule’ Ri [19] denotes decimal-equivalent of 8 outputs in Table 1. In the 2-state 3-neighborhood cellular automata, there are 256 rules.

Table 1

Look-Up Table for Rule 30

Present State 111 110 101 100 011 010 001 000
Rule (7) (6) (5) (4) (3) (2) (1) (0)
Next State 0 0 0 1 1 1 1 0

One such rule is 30 in Table 1. Rule 30 CA can generate a sequence of random patterns. Scientists initiate in their experiments, an n-cell rule 30 CA, with a state all 0 and single 1 as its seed. This CA produces a state containing a fair distribution of 0 as well as 1 after n number of iterations. But, this will not guarantee that after n-iterations all 0 pattern will never come. It can be proved that for every, n > 1, this CA produces a non-zero states after n/2 iterations.

In the classification scheme of Wolfram, Rule 30 belongs to class III [19], with aperiodic, chaotic behavior. Stephen Wolfram first noticed intrinsic randomness to appear in one deterministic system for Rule 30 CA. When the initialization is done using a black pixel, one random pattern center is generated as in Figure 1 b). In Figure 1 a), b), c) and d), different randomly generated patterns as are created by Rule 30 are shown with different initializations. This is an example to create random patterns from custom-defined rules. In fact, Mathematics uses the center column of pixel values as one of its random number generator.

Figure 1 Example of Wolfram Rule 30. A one-dimensional random number generator: Wolfram CA rule 30. Grid size is N = 21. White squares represent cells in state 0, black squares represent cells in state 1. Initial representations are – a) [1 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0], b), [0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0], c) [00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1] and d) [0 0 0 1 1 1 0 1 0 0 0 0 0 1 10 0 1 0 1 0].
Figure 1

Example of Wolfram Rule 30. A one-dimensional random number generator: Wolfram CA rule 30. Grid size is N = 21. White squares represent cells in state 0, black squares represent cells in state 1. Initial representations are – a) [1 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0], b), [0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0], c) [00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1] and d) [0 0 0 1 1 1 0 1 0 0 0 0 0 1 10 0 1 0 1 0].

In the periodic boundary CA, leftmost as well as rightmost cells are neighbors. Otherwise, it is called the null boundary CA. Further classification of CA defines deterministic as well as probabilistic (or, stochastic) CAs. The CA in Table 1 is of deterministic type.

2.3 Proposed hybrid KCA clustering algorithm

Our current hybrid model, named KCA algorithm has two steps. In first phase, K-Means segmentation approach clusters pixels in remote sensing imagery to similar land use land cover regions. In the fine-tuning step, the cellular automata based on neighborhood priority enhancement method is implemented, as described in the following Subsection 2.3.1. The two steps are shown in Figure 2.

Figure 2 The flowchart of KCA algorithm for remote sensing imagery segmentation.
Figure 2

The flowchart of KCA algorithm for remote sensing imagery segmentation.

Initial random clusters to pixels is assigned in K-Means algorithm, as mentioned in Section 2, then the initialized cluster centroids C(0) is computed from initial pixel allocations. Mean Square Error (MSE) objective function value ε for stopping iterations is 1E−05.

Then the centroid reallocation method reduces point-centroid distances by choosing the minimum one in each iteration. The iterations converge while reducing Mean Square Error (MSE). The algorithm stops when difference between previous error value and present error value remains smaller than ε or defined number of iterations has been achieved.

After the first allocations from K-Means method, the new cellular automata based neighborhood priority enhancement method is exploited over all pixels as described in next subsection. In final CA-corrected solutions, the quantitative evaluations have been performed.

2.3.1 Cellular automata based neighborhood priority enhancement method

2-dimesional Cellular Automaton can address 4 orthogonally adjacent cells in von Neumann neighborhood model. To enhance remote sensing image segmentation problem, this model can enhance mixed pixel allocations considering neighborhood regions. Therefore, in this work, we have implemented a hybrid algorithm with well-known K-means method enhanced with 2-dimensional von Neumann neighborhood based Cellular Automata approach.

In the 2-dimensional hybrid CA model in the final fine-tuning phase, to denote the states of cells in each Cellular Automaton, the clustering allocation values from K-Means based initial clustering phase, has been considered. Cells in one Cellular Automaton define pixels in their positions in remote sensing satellite imagery. State values of cells in the model at first denote the assigned clusters from the first phase of our algorithm, state 0 denotes that the pixel is assigned to cluster 0 in first phase. In the implemented model, CA with null boundary condition is adopted. Deterministic CA has been experimented with in this work. Our proposed 2-dimensional K-Cellular Automaton framework is shown as a flowchart in Figure 2.

This fine-tuning phase CA is implemented based on the cluster allocations outputs from the initial well-known partitioning K-means algorithm, for cluster allocations correction based on neighborhood. The experiment is done on both data - the chosen remote sensing image and chosen remote sensing Salinas A data set. Now, if among 4 neighbors (left, right, top, bottom), at least two neighbors show one specific cluster value, the rank of current cell upgrades to be of that specific cluster value. So, in this case, the present state of current cell remains in similar cluster depending on the neighborhood values in the model. Furthermore, if more than two neighbors have different cluster values other than current cell cluster value, the next state of current cell remains to the same cluster value. This approach reduces number of outliers in small homogeneous regions. It also produces more efficient solutions for mixed pixels. The corrections are performed over the cells in CA matrix iterative for static number of iterations based on dimensions of the image, to get proper neighborhood enhancements for pixels in chosen satellite catchment imagery.

3 Application to pixel classification

3.1 Study area

The chosen LANDSAT image of Ajoy River catchment for research analysis has been extracted within Bard-haman district in West Bengal State in India. There are green, red and blue bands available with original image in Figure 3. Figure 4 shows original image of the Ajoy river catchment region with histogram equalization with 7 classes. The classes are: water, forest, grassland, shrubband, cropland, urban and others in Table 2. The samples and percent values have been obtained using threshold around the RGB values relevant to each of the classes as described in LULC classification scheme from landcover.org.

Figure 3 Study area of the Ajoy river catchment region.
Figure 3

Study area of the Ajoy river catchment region.

Figure 4 LULC classification of the chosen area following landcover.org
Figure 4

LULC classification of the chosen area following landcover.org

Table 2

A detailed classification scheme experimented in our segmentation approach on chosen study area

Code Land cover Types Description Samples Percent (%)
0 Water Turbid water, rivers, canals 77 0.1059
1 Forest Evergreen and Deciduous broadleaf forest 3135 4.3086
2 Grassland Grassland 8316 11.4400
3 Shrubband Closed and open shrubband 9202 12.6589
4 Cropland Vegetation 10974 15.0966
5 Urban Urban and built 400 0.5503
6 Others Open space and other types 39296 54.0582
Total 72692 100

The river Ajoy flows through the middle part of catchment area. The river starts flowing from upper left corner to lower right corner of the study catchment area as a thin line in the middle.

3.2 Data and Experimental Framework

The experiment has also been implemented on Salinas - A hyper spectral imagery using our proposed KCA algorithm and other chosen two well-known methods. This data set is experimented for internal validation with true class labels. Statistical evaluation has been performed over the solutions over this data set. Salinas scene has been collected by the 224-band AVIRIS sensor over Salinas Valley, California with 3.7-meter pixels spatial resolution. Salinas-A scene is a small sub scene of Salinas image with 86*83 pixels including six ground truth classes in Table 4.

Table 3

Land cover/use classes on chosen image obtained by K-means, FCM and KCA algorithms

Land cover / use classes Algorithms
K-Means FCM KCA
Samples Percent (%) Samples Percent (%) Samples Percent (%)
Water 6499 8.9415 13216 18.1808 5235 7.2016
Forest 5853 8.0518 8748 12.0343 11943 16.4296
Grassland 12290 16.9069 7802 10.7330 8757 12.0467
Shrubband 12368 17.0143 8889 12.2283 12930 17.7874
Cropland 12614 17.3527 12591 17.3210 13656 18.7861
Urban 10026 13.7924 11329 15.5849 4295 5.9085
Others 13042 17.9415 10117 13.9176 15876 21.8401
Table 4

Ground Truth Knowledge for Salinas-A with respective samples number

Code Class Samples Percent(%)
0 Barrier 1790 25.077
1 Brocoli_green_weeds_1 391 5.478
2 Corn_senesced_green_weeds 1343 18.815
3 Lettuce_romaine_4wk 616 8.630
4 Lettuce_romaine_5wk 1525 21.365
5 Lettuce_romaine_6wk 674 9.442
6 Lettuce_romaine_7wk 799 11.194
Total 7138 100

The new KCA algorithm is implemented using MATLAB R2014a on Intel(R) Core(TM) i3-3220 processor with 3.30 GHz. For comparative evaluations, well-known K Means as well as FCM models are also implemented. For quantitative evaluation, we evaluate Dunn [5] and Davies-Bouldin (DB) [4] internal validity indices on clustering solutions of KCA, K-Means, DBSCAN as well as FCM algorithms. The segmented regions by KCA approach is also verified with the ground truth of land use land cover models.

4 Results and Discussion

4.1 Land Use Land Cover Segmentation over Ajoy river catchment

The segmented images of Ajoy river catchment as produced by K-Means and FCM methods are depicted in Figure 5 and Figure 6 for (K = 7). In Figure 5, K-Means algorithm fails to detect the river from its catchment background. FCM method in Figure 6 produces better solutions, but was unable to detect the middle part of catchment area accurately. The water bodies on the top-left and middle-left parts and the lower-middle catchment background are mixed in all segments.

Figure 5 Remote sensing image segmentation of Ajoy river catchment area obtained by K-Means method (for K=7).
Figure 5

Remote sensing image segmentation of Ajoy river catchment area obtained by K-Means method (for K=7).

Figure 6 Remote sensing image segmentation of Ajoy river catchment area obtained by FCM algorithm (for K=7).
Figure 6

Remote sensing image segmentation of Ajoy river catchment area obtained by FCM algorithm (for K=7).

K-Means and FCM solutions as shown in Figure 5 and Figure 6. The new KCA approach in Figure 7 is able to separate water bodies from lower middle catchment background. It also classifies the river water more efficiently than other two approaches. These segmentation analyses show that our KCA approach can detect the overlapping land use land cover regions in more efficient way than well-known K-Means and FCM methods.

Figure 7 Remote sensing image segmentation of Ajoy river catchment area obtained by proposed KCA algorithm.
Figure 7

Remote sensing image segmentation of Ajoy river catchment area obtained by proposed KCA algorithm.

DBSCAN approach provides very fewer areas in different clusters, with only one cluster covering most areas. Therefore, its result is not considered for comparisons.

4.2 Segmentation results on remote sensing data set

The ground truth classes for Salinas A data set is shown in Figure 8. The classification results of FCM and KCA algorithms are shown in Figure 9 and Figure 10 respectively. DBSCAN fails to detect proper number of clusters for this data set. The ground truth class details for this data set has been shown in Table 4. Table 5 shows land cover classes as obtained by K-Means, FCM and KCA algorithms respectively.

Figure 8 Salinas A Groundtruth classes.
Figure 8

Salinas A Groundtruth classes.

Figure 9 Salinas A solutions obtained by FCM algorithm.
Figure 9

Salinas A solutions obtained by FCM algorithm.

Figure 10 Salinas A solutions obtained by KCA algorithm.
Figure 10

Salinas A solutions obtained by KCA algorithm.

Table 5

Land cover/use classes on Salinas A dataset obtained by K-Means, FCM and KCA algorithms

Code Land cover / use classes Algorithms
K-Means FCM KCA
Samples Percent Samples Percent Samples Percent
(%) (%) (%)
0 Barrier 2014 28.215 1023 14.332 364 5.100
1 Brocoli_green_weeds_1 502 7.033 502 7.033 509 7.131
2 Corn_senesced_green_weeds 612 8.574 511 7.159 1705 23.886
3 Lettuce_romaine_4wk 318 4.455 782 10.955 386 5.408
4 Lettuce_romaine_5wk 1705 23.886 1701 23.830 2009 28.145
5 Lettuce_romaine_6wk 1025 14.360 2014 28.215 1443 20.216
6 Lettuce_romaine_7wk 962 13.478 605 8.476 722 10.115

5 Quantitative analysis

The clustering solutions obtained by chosen three methods are evaluated objectively by measuring validity measures Davies-Bouldin (DB) and Dunn indices, as defined in [5] and [4] respectively. The validity indices values on the Ajoy catchment image are shown in Table 6. Then the validity indices values on Salinas A data set are shown in Table 7. It can be noticed that, KCA obtains best minimized value for DB index as 0.5187 in Table 6. On same data set K-Means exhibits a DB value of 0.9202, DBSCAN shows 0.9401 value and FCM exhibits 0.6582 value. Similarly, the maximizing Dunn index obtained by KCA method 1.6691, while K-Means and FCM method obtain smaller values.

Table 6

Quantitative internal validity evaluation over Ajoy river catchment as obtained by K-means, FCM and KCA algorithms

Index Ajoy river catchment
K-Means FCM DBSCAN KCA
Davies-Bouldin index 0.9202 0.6582 0.9401 0.5187
Dunn index 0.7053 0.8331 0.1139 1.6691
Table 7

Quantitative internal validity evaluation over Salinas A hyper-spectral dataset as obtained by K-means, FCM and KCA algorithms

Index Salinas A
K-Means FCM DBSCAN KCA
Davies-Bouldin index 0.5657 0.6075 0.2683 0.4775
Dunn index 1.7554 1.2500 0.0410 1.8399

In Table 7, the clustering solutions obtained by K-Means, FCM, DBSCAN and proposed KCA algorithms are evaluated quantitatively using internal validity DB and Dunn indices. In this data set, KCA approach also obtains the smallest DB index value of 0.4775 and highest Dunn index value of 1.8399 among chosen three experimented methods in Table 7. So, for this hyper-spectral data set, also the proposed KCA approach provides efficient clustering solutions while in comparison with other two methods. DBSCAN method even fails to detect proper number of clusters by detecting more than 10 clusters for this data set. Therefore, its smallest DB value is not considerable for creating non-comparable clusters on this data set.

Evaluating with the true class labels as available with Salinas-A hyper spectral data set, the solutions obtained from the proposed method are validated. The solutions are validated with four external validity indices, namely Rand index [16], Adjusted Rand index [9], Jaccard index [10], and Fowlkes-Mallows index [6]. All these indices provide higher scores for better solutions.

In Table 8, ARI indices values as obtained by K-Means, FCM and KCA algorithms are respectively 0.7471, 0.78017 and 0.83549. For Jaccard index also, the proposed KCA algorithm produces maximum value of 0.79285 among those three algorithms. Similar results are also obtained for other two indices over Salinas- A hyper-spectral data set considering its true labels as shown in Table 8. DBSCAN provides higher number of clusters for Salinas A data set, therefore it fails to get comparable external validity indices values for this data set.

Table 8

External validity indices values of the classified remote sensing image provided by K-Means, FCM and KCA algorithms

Index Salinas A
K-Means FCM KCA
Rand index 0.92921 0.93459 0.95376
Adjusted Rand index 0.7471 0.78017 0.83549
Jaccard index 0.65229 0.6965 0.79285
Fowlkes-Mallows (FM) index 0.79031 0.82111 0.88451

These results imply that KCA optimizes DB and Dunn internal indices as well as Rand, Adjusted Rand, Jaccard and Fowlkes-Mallows (FM) indices. The quantitative performance values for KCA are more than both K-Means and FCM algorithms for these indices. Hence, it is shown that KCA is comparably a better solution to K-Means and FCM methods.

6 Statistical analysis

A non-parametric statistical significance test called Wilcoxon’s rank sum is conducted for independent samples at 5% significance level [8].

Three groups are generated from DB index values by 10 consecutive runs for chosen three methods, K-Means, FCM, DBSCAN and KCA algorithms on both data sets. The median values of each group on both data sets are shown in Table 9. These results show that KCA gets smaller median values than K-Means, DBSCAN and FCM methods showing more efficiency. DBSCAN fails to detect proper number of clusters for Salinas A data set. Therefore, its outputs are not considered for statistical analysis.

Table 9

Median values of performance index DB for 10 consecutive runs chosen three approaches

Data Algorithms
K-Means FCM DBSCAN KCA
Ajoy river Catchment 0.8561 0.8057 0.9401 0.5218
Salinas A 0.4816 0.5786 - 0.4527

Table 10 exhibits P-values and H-values obtained from Wilcoxon’s rank sum tests on three groups, KCA-K-Means, KCA-DBSCAN and KCA-FCM over both data sets. All P-values which we obtain here are less than 0.005(5% significance level). For the Ajoy catchment area, P-value of rank sum test between KCA and K-Means is very small 2.00E-003. This determines that the clustering solutions obtained from KCA are statistically significant, and they do not occur by chance. Other group for FCM solutions also show similar significance. On the chosen Salinas A hyper-spectral data set, P-values obtained in similar tests by KCA-K-Means and KCA-FCMgroups are respectively 3.71E-002 and 3.9E-003. These P-values also satisfy null hypothesis. So, all results satisfy the statistical significant of the land cover land use segments from KCA algorithm over K-Means and FCM solutions.

Table 10

P-values produced by rank sum test when in comparison for KCA with K-Means and KCA with FCM algorithms respectively

Algorithm Comparison with KCA
Ajoy river Catchment Salinas A
H P-value H P-value
K-Means 1 2.00E−003 1 3.71E−002
FCM 1 7.8E−003 1 3.9E−003
DBSCAN 1 2.00E−003 - -

7 Conclusions

Image segmentation algorithms are most challenging methods used in remote sensing to help interpret the land use land cover models in the satellite imagery [1, 2, 3, 14]. Cellular automata model provides the changes in cell states depending on their neighborhood, which occur in discrete times. So, the cellular automata model can analyze land cover segments in remote sensing images, considering the neighborhood pixel classes.

The contribution of this article lies in significant improvement over detection of mixed land cover regions in the satellite image than existing partitioning algorithm. This new approach introduces a hybrid crisp set based pixel segmentation with cellular automata based neighborhood reallocation in the proposed KCA clustering algorithm. The primary contributions of this work are – to use one new crisp mean-centroid based initial clustering in satellite data sets with CA based neighborhood correction. The neighborhood enhancement phase partially correct outlier allocations in overlapping land use land cover regions significantly. It verifies the overall allocations regarding the neighborhood, to get improved land cover regions.

The performance of proposed KCA approach is shown over one chosen satellite image of Ajoy River catchment area. The experiment was further carried over one chosen Salinas A hyper-spectral data set with available true class labels. Significant efficiency of proposed KCA segmentation method in comparison with well-known K-Means, DBSCAN and FCM methods is established both using quantitative and statistical evaluations. For quantitative evaluation, two internal and four external validity indices are analyzed over both data sets. The ground truth verification also exhibits significant superiority of new KCA algorithm over other two well-known algorithms. Statistical tests are also performed to show the statistical significance of KCA solutions when comparing with K-Means and FCM algorithms on both chosen data sets.

References

[1] Manual of National Land Use Land Cover Mapping using Multi-Temporal Satellite Data, Department of Space, NRSA (2006).Search in Google Scholar

[2] Mukesh Singh Boori, Vít Voženílek and Komal Choudhary, Land use/cover disturbance due to tourism in Jeseníky Mountain, Czech Republic: A remote sensing and GIS based approach, The Egyptian Journal of Remote Sensing and Space Science 18 (2015), 17–26.10.1016/j.ejrs.2014.12.002Search in Google Scholar

[3] CL D Cordeiro and DD Rosetti, Mapping Vegetaion in a late quaternary landform of the Amazonian wetlands using object based image analysis and decision tree classification, Int. J. Remote sens 36 (2015), 3397–3422.10.1080/01431161.2015.1060644Search in Google Scholar

[4] David L. Davies and Donald W. Bouldin, A Cluster Separation Measure, IEEE Transactions on Pattern Analysis and Machine Intelligence PAMI-1 (1979), 224–227.10.1109/TPAMI.1979.4766909Search in Google Scholar

[5] J C Dunn, A Fuzzy Relative of the ISODATA Process and Its Use in Detecting Compact Well-Separated Clusters, Journal of Cybernetics 3 (1973), 32–57.10.1080/01969727308546046Search in Google Scholar

[6] E. B. Fowlkes and C. L. Mallows, A Method for Comparing Two Hierarchical Clusterings, Journal of the American Statistical Association 78 (1983), 553–569.10.1080/01621459.1983.10478008Search in Google Scholar

[7] RC Gonzalez and R E Woods, Digital image processing, Addison-Wesley, 1992.Search in Google Scholar

[8] M Hollander and D Wolfe, 1999.Search in Google Scholar

[9] Lawrence Hubert and Phipps Arabie, Comparing partitions, Journal of Classification 2 (1985), 193–218.10.1007/BF01908075Search in Google Scholar

[10] Paul Jaccard, THE DISTRIBUTION OF THE FLORA IN THE ALPINE ZONE.1, New Phytologist 11 (1912), 37–50.10.1111/j.1469-8137.1912.tb05611.xSearch in Google Scholar

[11] S. Kirkpatrick, C. D. Gelatt and M. P. Vecchi, Optimization by Simulated Annealing, Science 220 (1983), 671–680.10.1126/science.220.4598.671Search in Google Scholar PubMed

[12] T Kohonen, Self-Organization and Associative Memory Springer, Berlin, 1984.Search in Google Scholar

[13] J B Macqueen, Some Methods for classification and Analysis of Multivariate Observations, Proceedings of 5th Berkeley Symposium on Mathematical Statistics and Probability 1 (1967), 281–297.Search in Google Scholar

[14] Khalid Mansour, Onisimo Mutanga, Elhadi Adam and Elfatih M. Abdel-Rahman, Multispectral remote sensing for mapping grassland degradation using the key indicators of grass species and edaphic factors, Geocarto International 31 (2016), 477–491.10.1080/10106049.2015.1059898Search in Google Scholar

[15] Ujjwal Maulik and Anasua Sarkar, Eflcient parallel algorithm for pixel classification in remote sensing imagery, GeoInformatica 16 (2012), 391–407.10.1007/s10707-011-0136-5Search in Google Scholar

[16] William M. Rand, Objective Criteria for the Evaluation of Clustering Methods, Journal of the American Statistical Association 66 (1971), 846–850.10.1080/01621459.1971.10482356Search in Google Scholar

[17] Jianbo Shi, Jitendra Malik ; Dembele, D Kastner and P, Fuzzy c-means method for clustering microarray data, IEEE Conference on Computer Vision and Pattern Recognition 19 (1997), 973–980.10.1093/bioinformatics/btg119Search in Google Scholar

[18] Iii Smith and A R, Two-dimensional Formal Languages and Pattern Recognition by Cellular Automata, IEEE Conference Record of 12th Annual Symposium on Switchinh and Automata Theory (1971).10.1109/SWAT.1971.29Search in Google Scholar

[19] S Wolfram, Cryptography with cellular automata, Lecture Notes on computer Science 218 (1986), 429–432.10.1007/3-540-39799-X_32Search in Google Scholar

Received: 2019-06-08
Accepted: 2019-11-26
Published Online: 2020-09-02

© 2020 K. Mahata et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 24.4.2024 from https://www.degruyter.com/document/doi/10.1515/jisys-2019-0155/html
Scroll to top button