Journal of Sensors

Volume 2015, Article ID 327123, 8 pages

http://dx.doi.org/10.1155/2015/327123

## Automatic Change Detection Method of Multitemporal Remote Sensing Images Based on 2D-Otsu Algorithm Improved by Firefly Algorithm

^{1}Faculty of Land Resource Engineering, Kunming University of Science and Technology, Kunming 650093, China^{2}Kunming Surveying and Mapping Institute, Kunming 650051, China

Received 30 October 2014; Accepted 3 February 2015

Academic Editor: Xue Cheng Tai

Copyright © 2015 Liang Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper presents a new automatic change detection method of multitemporal remote sensing images based on 2D-Otsu algorithm improved by Firefly algorithm. The proposed method is designed to automatically extract the changing area between two temporal remote sensing images. First, two different temporal remote sensing images were acquired through difference value method of remote sensing images; then, the 2D-Otsu threshold segmentation principles are analyzed and the optimal threshold of 2D-Otsu threshold segmentation method is figured out by using the Firefly algorithm, where the difference images are conducted with binary classification to obtain the changing category and the nonchanging category; finally, the proposed method is used to carry out change detection experiments on the two selected areas, where a variety of methods are compared. Experimental results show that the proposed method can effectively and quickly extract the changing area between the two temporal remote sensing images; thus, it is an effective method of change detection for remote sensing images.

#### 1. Introduction

With the development of society and technology, human activities are changing the landscape of ground surface and its utilization mode every day. The rapid population growth and frequent natural disasters have also accelerated the speed of such changes. Therefore, quickly and effectively monitoring the information of such changes and analyzing the characteristics, causes, impacts and results of the changes have great significance to the global sustainable development [1]. Emergence and development of remote sensing technology have provided technical support for quickly and effectively monitoring the information of these changes. Detection of remote sensing image changes aims to acquire the changing information of desired surface features by analyzing and processing two or multiple remote sensing images of the same area at different times [2, 3]. Currently, the change detection technology of remote sensing image has been widely used in many fields, such as disaster assessment [4], land use/coverage monitoring [5], environmental change monitoring [6], agricultural survey [7], and urban planning [8].

Currently, a number of methods have been proposed for the change detection of remote sensing image. At present, change detection mainly includes two strategies: direct comparison and postclassification comparison, of which the postclassification comparison strategy is an approach to classifying (supervised or unsupervised method) different temporal images of the same area to compare and analyze the classified results in order to obtain the location and type of changing information. The advantage of this strategy is to minimize the impact of nonfeature change factors, but the change detection accuracy through this strategy is substantially equal to the product of the two classification accuracy values, and the classification error of each temporal classification result will be amplified in the process of comparison, thus inevitably exaggerating the extent of changes. Meanwhile, it is often more complex and difficult to get high-accuracy classification results, resulting in low-accuracy and uncertain change detection results [9]. Therefore, the strategy of direct comparison has always been highlighted in researches on change detection. Direct comparison is a strategy to obtain difference images by directly operating and alternating the pixel values of different temporal remote sensing images of the same area that have been calibrated so as to figure out the changing area by analyzing the difference images. In analyzing these difference images, it would be the most simple thing to directly use one-dimensional threshold processing methods, but the processing accuracy of threshold is not high [10, 11]. In order to improve the processing accuracy of threshold, many scholars have proposed to extend the one-dimensional threshold selection method to two dimensions (2D). Currently, a number of two-dimensional threshold selection methods have been raised, such as 2D entropy (2D maximum entropy, 2D minimum cross entropy) [12], 2D-Otsu [13], and 2D maximum fuzzy entropy [14], but it is still difficult to select the threshold values of the two-dimensional histogram, for which unremitting exploration efforts have been made by many scholars. To solve this problem, many scholars attempt to achieve fast optimization for the 2D threshold by combining 2D threshold segmentation method and optimization algorithm. Shen et al. [15], Zheng et al. [16], and Alim et al. [17] sought to figure out the optimal threshold of 2D maximum entropy, respectively, through genetic algorithm, ant colony algorithm, PSO (Particle Swarm Optimization) algorithm, and ABCO (Artificial Bee Colony Optimization). Qian used PSO algorithm to find the optimal threshold of 2D-Otsu [18]. Tian and Zeng used QPSO (Quantum-behaved Particle Swam Optimization) algorithm to carry out image threshold segmentation in combination with 2D maximum fuzzy entropy [19]. The applications of PSO algorithm, QPSO algorithm, bee colony algorithm, ant colony algorithm, and genetic algorithm have improved the speed of figuring out thresholds, but the optimization results may be inaccurate as these algorithms are prone to falling into local extremum.

The Firefly algorithm [20] is a global optimization algorithm proposed by Yang, which can overcome the problem of easily falling into local optimum. Chen et al. [21] and Alomoush et al. [22], respectively, applied the Firefly algorithm to find optimal threshold and then used for image segmentation and achieved good segmentation results. Given the fast global search capability of the Firefly algorithm and the good segmentation result of 2D-Otsu threshold segmentation method, this paper segments the thresholds of difference images by combing the Firefly algorithm and 2D-Otsu threshold method to obtain the binary change detection map.

#### 2. Architecture of the Proposed Method and Problem Formulation

Assume and are the remote sensing images of the same area at different times and that underwent standard product preprocessing and coregistration; the remote sensing image sizes are . The ultimate aim of the proposed method is to generate the binary change map. The proposed method consists of three steps: construction of difference image; threshold optimization based on Firefly algorithm; and generation of final binary change map.

First, it is construction of difference images. Currently, there are two methods, namely, image algebraic operation and image transformation, to construct difference images. The method based on image algebraic operation includes difference method, ratio method, and the combination of difference method and ratio method; the method based on image transformation includes principal component analysis, change vector analysis and correlation analysis. Since the image algebraic operation method requires simple algorithm but can obtain high accuracy, this paper applied the image difference method of the image algebraic operation, which can be expressed below: where is the difference image constructed.

Second, it is optimization of Firefly threshold. After acquiring the difference images, the difference images need to be analyzed so as to obtain the binary change map. Binary classification through threshold segmentation is one of the most frequently used methods. Selecting thresholds via one-dimensional histogram is the easiest method, of which the main method is OTSU algorithm. However, as OTSU algorithm does not use local spatial data of the image, the segmentation will be ineffective when the image is subject to noise disturbance or gray scale distribution intersection. Therefore, many scholars have expanded it to two-dimensional approach, and better segmentation results have been achieved [23]. To this end, this paper used the 2D-Otsu algorithm to analyze difference images acquired. However, compared with the thresholds selected in one-dimensional histogram, 2D-Otsu algorithm requires a large amount and a long time of computation. To solve the above problem, the Firefly algorithm is introduced to optimize the threshold of 2D-Otsu algorithm.

Third, it is generation of the final binary change map. After figuring out the optimal threshold of 2D-Otsu algorithm, threshold segmentation is given to difference images to obtain the binary change detection map of the two temporal remote sensing images.

#### 3. Change Detection Based on 2D-Otsu Algorithm Improved by Firefly Algorithm

##### 3.1. 2D-Otsu Threshold Segmentation Method

The size of the difference image is , where , , the gray scale of the image is , and the average gray-scale of pixel neighborhood is also divided into level. A binary group is obtained by calculating the average gray scale of its neighborhood at each pixel point, that is, the gray-scale value of the pixel point and the average gray-scale value of its neighborhood. Assume the probability of the binary group is ; the corresponding joint probability density can be defined as [13, 24]where is the number of pixels of the difference image, , and

Assume that there are two categories, namely, (changing category) and (nonchanging category) in the 2D histogram, as well as two different probability density distribution functions. Set the threshold value as ; then, the probabilities of and are shown below:

The corresponding mean vectors of and are

The total mean vector of 2D histogram is

Usually, the probability away from the histogram diagonal can be negligible; then, it can be assumed in the two areas: ; and ; have , where

Define a between-class dispersion matrix:

Use the trace of as a measure of the between-class dispersion matrix; then, the 2D-Otsu function of the threshold value corresponding to the difference image can be defined asBy using Formula (8), it can be simplified intowhere , .

The segmentation criteria corresponding to the 2D-Otsu threshold segmentation method has to maximize the 2D-Otsu function for the purpose of obtaining the optimal threshold .

##### 3.2. 2D-Otsu Threshold Segmentation Method Based on Firefly Algorithm

In order to solve the problem that 2D-Otsu threshold segmentation method requires a large amount and long time of computation, this paper applied the Firefly algorithm to the threshold optimization for 2D-Otsu threshold segmentation method and proposed the 2D-Otsu threshold segmentation method based on the firefly algorithm, so as to convert the threshold value selection of 2D-Otsu threshold segmentation method into the optimization of 2D-Otsu function based on the Firefly algorithm.

The Firefly algorithm [25], proposed by Yang in the Cambridge University, is a bionic swarm intelligent optimization algorithm by simulating the natural behaviors of fireflies. The algorithm is characterized by simple parameter setting, high-accuracy optimization, and powerful global optimization [20]. The principle of Firefly algorithm is to simulate firefly individuals in nature by searching for space dots. In the optimization process, by using the phototaxis characteristics of firefly, the searching and optimization process is simulated into the attraction and movement process of firefly individuals, where the target function of the problem is measured by the merits of the position of firefly individuals. In other words, the target function of the problem is converted into the firefly seeking maximum brightness. The process of survival of the fittest among firefly individuals is also the iterative process of feasible solutions in the course of target function optimization. Therefore, the Firefly algorithm can rapidly conduct global optimization [21].

Firefly algorithm consists of two main elements: brightness and attraction, where brightness reflects the merits of the position of fireflies and determines the movement direction of firefly individuals, while attraction determines the movement distance of firefly individuals. Through the continuous updating and iteration of brightness and attraction, the optimal solution to the target function will be achieved [21, 26, 27]. In the Firefly algorithm, the relative fluorescence brightness of firefly individuals is defined aswhere is the maximum fluorescence brightness of the firefly, is the attraction coefficient of light intensity, which is usually set to a constant; is the spatial distance between Fireflies and , that is, , where and are, respectively, the spatial position of Fireflies and .

The attraction degree of firefly is defined aswhere is the biggest attraction of the firefly.

In the case that Firefly is attracted by Firefly and moves towards the position of Firefly , it can be defined in the formula below: where is the step size factor which is a constant between ; is the uniformly distributed random factor between ; is the random disturbance to avoid untimely falling into local optimum.

By setting the 2D-Otsu function as the target function of the Firefly algorithm, the optimization result of the Firefly algorithm is the position of the firefly with maximum brightness and is the desired threshold. The implementation process is as follows:(1)initialize the basic parameters of the Firefly algorithm. Set the number of fireflies as , of which the maximum attraction is , the light intensity attraction coefficient is , the step factor is , and the maximum number of iterations is ;(2)by randomly initializing the position of fireflies, calculate the 2D-Otsu function value corresponding to each firefly and take the target function as their brightness to get the position of the firefly with maximum brightness;(3)calculate the relative brightness and attraction degree of fireflies based on (12) and (13), and determine their movement direction according to the relative brightness;(4)update the spatial position of fireflies based on Formula (14), and carry out random disturbance on fireflies in the best position;(5)recalculate the brightness of fireflies based on the updated position of fireflies;(6)when the maximum number of searches is reached, go to (6); otherwise, the number of searches should plus 1; then, go to (3) for the next search;(7)output the position and brightness of the firefly with maximum brightness, and segment the difference images by taking as a threshold.

#### 4. Experiments and Results

##### 4.1. Description of Experimental Data

In order to verify the validity and reliability of the proposed method as well as the advantages of the proposed method on the operating speed, the paper selected two groups of experimental data set, both of which were Landsat5 TM remote sensing images, with a spatial resolution of 30 m. The data set is provided by International Scientific & Technical Data Mirror Site, Computer Network Information Center, and Chinese Academy of Sciences (http://www.gscloud.cn/).

The first group of data set consists of two Landsat5 TM remote sensing images of the local area of Poyang Lake in Jiangxi Province, as shown in Figures 1(a) and 1(b). Figures 1(a) and 1(b) were, respectively, acquired on July 15, 1989, and on January 14, 2010, and the size of two remote sensing images was 510 × 510 pixels, with a gray-scale of 256; their reference changing map is shown in Figure 1(c), where the black area is the changing area, and 69,413 pixels were changed and 190,687 pixels were unchanged.