Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
An Improved Shape Contexts Based Ship Classification in SAR Images
Previous Article in Journal
Mapping Seasonal Inundation Frequency (1985–2016) along the St-John River, New Brunswick, Canada using the Landsat Archive
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Change Vector Analysis to Monitor the Changes in Fuzzy Shorelines

1
Faculty of Geo-Information Science and Earth Observation (ITC), University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
2
Geospatial Information Agency (BIG), Jl. Raya Jakarta-Bogor Km. 46, Cibinong 16911, Indonesia
*
Author to whom correspondence should be addressed.
Remote Sens. 2017, 9(2), 147; https://doi.org/10.3390/rs9020147
Submission received: 23 November 2016 / Accepted: 6 February 2017 / Published: 10 February 2017

Abstract

:
Mapping of shorelines and monitoring of their changes is challenging due to the large variation in shoreline position related to seasonal and tidal patterns. This study focused on a flood-prone area in the north of Java. We show the possibility of using fuzzy-crisp objects to derive shoreline positions as the transition zone between the classes water and non-water. Fuzzy c-means classification (FCM) was used to estimate the membership of pixels to these classes. A transition zone between the classes represents the shoreline, and its spatial extent was estimated using fuzzy-crisp objects. In change vector analysis (CVA) applied to water membership of successive shorelines, a change category was defined if the change magnitude between two years, T1 and T2, differed from zero, while zero magnitude corresponded to no-change category. Over several years, overall change magnitude and change directions of the shoreline allowed us to identify the trend of the fluctuating shoreline and the uncertainty distribution. The fuzzy error matrix (FERM) showed overall accuracies between 0.84 and 0.91. Multi-year patterns of water membership changes could indicate coastal processes such as: (a) high change direction and high change magnitude with a consistent positive direction probably corresponding to land subsidence and coastal inundation, while a consistent negative direction probably indicates a success in a shoreline protection scheme; (b) low change direction and high change magnitude indicating an abrupt change which may result from spring tides, extreme waves and winds; (c) high change direction and low change magnitude which could be due to cyclical tides and coastal processes; and (d) low change direction and low change magnitude probably indicating an undisturbed environment, such as changes in water turbidity or changes in soil moisture. The proposed method provided a way to analyze changes of shorelines as fuzzy objects and could be well-suited to apply to coastal areas around the globe.

Graphical Abstract

1. Introduction

The study of changing shorelines is essential to assist in the design of effective coastal protection [1,2], verifying numerical models [3,4], developing hazard maps [5], formulating policies regarding coastal development [6], and for coastal research and monitoring [7]. A shoreline is defined as the intersection of coastal land and water surface indicating the water edge movements of which the position is changing through time due to different water levels during high tide and low tide [8,9,10]. Oertel [11] referred to a shoreline as the line associated with sea level rather than with high and low tides. When considering only the tide, many shorelines are due to the shifting of water with tidal differences. Tidal differences vary and are influenced by the changes in the magnitude of gravitational attractions on the water body of the Earth, winds and waves. Furthermore, shorelines have changed their dynamics at varying rates as a response to coastal processes such as sediment erosion, transportation and deposition along the shore. Rapid changes occur during an extreme event such as storms, whereas gradual changes occur during an intervening period [12]. Shoreline changes can be estimated over various time scales and result into long-term, cyclic and local random variation. Long-term variation includes variation due to the land subsidence, relative sea level rise and sediment storage. Cyclic variation is related to the tide cycles or seasons, whereas waves and storms cause random variation of a local character.
As the shoreline positions vary over time, shoreline indicators are used as proxies to represent shoreline positions, including: (a) distinguishable coastal features, for example, a previous high-tide line; (b) the intersection of coastal profile with specific vertical water elevation, e.g., the mean sea level; and (c) shoreline features observable from remote sensing images. An example of the latter is the boundary between water and non-water pixels [13,14,15,16]. The shoreline boundary between coastal land and water is fuzzy since there is a gradual transition from coastal land to water. Given the nature of the fuzzy shoreline and its changing position, detection of shoreline requires dealing with uncertainty. Fisher [17] mentioned three types of uncertainty: (a) errors: if a shoreline is clearly identified, the uncertainty may arise from error, for example in data processing, spatial generalization, and measurement; (b) vagueness: if it is not possible to define the spatial extent of coastal land, water, and the transition zone (Williamson (1994) after Fisher [17]); and (c) ambiguity: relating to the confusion of land and water definition considering a different classification system or a different perception.
Previous studies have proposed several ways of generating shoreline positions. Shoreline survey and photogrammetry have been primary technology for shoreline mapping, yet both methods are time consuming and expensive [18]. Therefore, image classification is used widely nowadays to detect shoreline positions. Most studies regarding shoreline detection have used hard classification such as thresholding, water indices, iterative self-organizing data analysis (ISODATA), binary slicing, maximum likelihood classification (MLC) and manual digitizing [7,15,19,20,21], whereas only a few applied soft classifications [16,22,23].
Due to the fuzziness of shoreline positions, using hard classification for shoreline mapping could produce errors on the classification results, since hard classification assigns a single label to a pixel, based on its highest membership. To overcome this limitation, this paper explores fuzzy classification to detect shoreline positions from a remote sensing image. In our previous work, we proposed two procedures to derive fuzzy shorelines: (a) we derived shorelines by applying a threshold equal to 0.5 to the membership and depicted shorelines as a single line; and (b) we derived shorelines as a margin determined by the choice of thresholds on the membership function [16]. In the current paper, we proposed a third procedure to distinguish shoreline proxies from digital images. A shoreline is represented as the transition zone between water and land. In this case, pixels at which the membership value (μ) exceed 0.99 are the core of a class, whereas pixels with 0.01 < μ < 0.99 belong to transition zones and pixels with μ < 0.01 do not belong to objects. In this way, we can account for the gradual transition between water and land (vagueness of the boundary). Moreover, in change detection, use of transition zones instead of crisp shorelines allows us to account for the influence of ambiguity resulting from comparing images recorded under different circumstances, such as weather, and have a more detailed description, of not only the magnitude and direction of the changes, but also of the related uncertainty.
Various change detection techniques have been developed. They can be divided into two groups, namely bi-temporal change detection and temporal trajectory analysis [24,25]. The former measures changes based on two separate time periods, for example image differencing and post classification comparison. Image differencing does not provide a detailed change matrix while post classification comparison does not allow the detection of subtle changes within a class. The latter, temporal trajectory analysis, is based on the temporal development curve or trajectory for successive times. It focuses both on what has changed between dates, as well as on the trend of the change over the period [24,25,26,27]. Change detection in this research utilizes the second method. To measure the change of the fuzzy shoreline, change vector analysis (CVA) based upon pixel-wise comparison was used to estimate the changes of successive shorelines. CVA identifies changes of features which were acquired at different times. In previous studies, CVA was applied to the brightness and greenness indices [28,29], normalized difference vegetation index [25,30], near infrared band and vegetation index [31], wetness and bare soil index [32], and spectral bands and textural images [33,34]. In this study, CVA was applied to the water membership values of shoreline images. Furthermore, in the earlier studies CVA has been applied in the multi-spectral space [29,35], and then extended to be applied in multi-temporal observation vectors of an indicator variable measured at different times [25]. CVA provides an overall change magnitude and change direction showing the trend of the fluctuating shoreline.
The objective of this study was to develop a method that is useful for monitoring the changes of a fuzzy shoreline. The method is based on fuzzy classification and CVA. A series of Landsat images is used to detect shoreline positions as a transition zone while taking tides into account. For this study, the uncertainty of shoreline positions was estimated by means of confusion indices. We focus on inherent uncertainty caused by continuous variation of a shoreline over time, and on uncertainty as it propagates from extraction and implementation of the shoreline change detection method. The method is applied to an area in Java, in the northern coastal area of the Central Java Province, Indonesia, where extensive shoreline changes associated with coastal inundation have increased in term of frequency and duration.

2. Methodology

2.1. Study Area

The study area is located at the northern part of the Central Java Province in Indonesia. It is characterized by a low-land landscape (Figure 1) with an elevation less than 5 m above mean sea level (AMSL). The extent of the study area is approximately 7.5 km from east to west, and 6.5 km from north to south. The central point of the area is at UTM coordinates 444,243°E and 9,234,731°N zone 49S or geographic coordinates 6°55′S and 110°29′E. Alluvial and sand sedimentation dominate its soil type [36,37]. As a coastal area, this area has a mixed semi-diurnal tide with two high tides and two low tides each day. These two highs and lows differ in height, whereas the average tidal range is 0.6 m. The highest tidal ranges occur in December and June during the rainy and dry seasons, while the lowest tidal ranges occur in March and September during the transitional seasons. Two types of flood regularly occur: (a) floods caused by a tidal flood occurring daily in line with tidal cycles [38,39]; and (b) floods due to poor drainage systems during rainy seasons.
Extensive fishponds and rice fields are covering the study area. Settlements are found bordering the sea and along the riverbanks, which are threatened if tidal floods become higher. This area has suffered from a changing shoreline position leading to severe coastal inundation and erosion. The coastal inundation has increased recently both in terms of frequency and duration. Some factors such as extreme winds and waves contribute to this increase. Furthermore, in the long run, other causes such as land subsidence, sea level rise, mangrove conversion, beach reclamation and a seaport extension are potential causes for increased coastal inundation [38,40,41,42].

2.2. Satellite Images, Data Pre-Processing and Reference Data Generation

2.2.1. Satellite Images and Data Pre-Processing

Multi-temporal images from the Landsat 8 OLI/TIRS (Operational Land Imager/Thermal Infrared Sensor) with 30 m spatial resolution were used to monitor the shoreline change between 2013 and 2015 (Table 1). We obtained terrain corrected Landsat images (L1T product) from USGS EarthExplorer [43]. Those images were acquired at the low tide. Tidal data relating to the time of acquisition of the images were collected from the Indonesian Geospatial Information Agency.
Pre-processing of Landsat 8 OLI/TIRS comprises two steps: (a) histogram minimum adjustment; it was applied to remove the influence of atmospheric path radiance [44,45]; and (b) geo-referencing; it was implemented using >100 ground control points (GCP) collected from road intersections, rivers and other prominent features. The root mean square error (RMSE) values were less than 0.1 pixels. Geo-registration of Landsat images was conducted using geometrically corrected reference images: (1) a Pleiades image at a 2 m spatial resolution; (2) a SPOT 6 (Satellite Pour l’Observation de la Terre) image at a 6 m spatial resolution; and (3) a Sentinel 2 image at a 10 m spatial resolution. The spectral band information for each reference image including Landsat 8 OLI/TIRS is available in Table 2.

2.2.2. Reference Data Generation

To evaluate the accuracy of a fuzzy classification, it is necessary to use soft reference data [46,47]. We generated soft reference data from available fine resolution datasets [48,49]. These datasets (Pleiades, Spot 6 and Sentinel 2) were rectified using a 2015 orthoimage. To reduce the variance of the Pleiades image, smoothing was performed using the average filter applied to a 3 × 3 window size. Afterwards, we applied fuzzy c-means (FCM) with the number of classes c = 2 and the fuzzy weight m = 1.7 [16]. Further, membership images generated using FCM classification from these high resolution datasets were used as reference images.
For accuracy assessment purpose, the pixel size of Spot 6 image was resampled to 10 m using nearest neighbour resampling, so that the spatial resolution of Pleiades, Spot 6, Sentinel 2 and Landsat images were in the ratio 15:3:3:1. Hence, 225 pixels (15 × 15) of Pleiades, 9 pixels (3× 3) of Spot 6, and 9 pixels (3 × 3) of Sentinel 2 were combined (pixel values averaged) to achieve the pixel dimension of Landsat images. Furthermore, an effective comparison could be made between images of different resolutions.
For the alternative methods, MLC and hardened classification, we visually interpreted Pleiades and Spot 6 images as hard reference data for the year 2013 and 2014 respectively, whereas ground data were used as the 2015 reference data.

2.3. FCM Classification

To discriminate water classes from non-water, we applied a fuzzy c-means (FCM) classification [50]. FCM iteratively separates data clusters with fuzzy means and fuzzy boundaries and the results assign each pixel to a partial membership of land cover classes. The membership values (μ) range from 0 to 1, and add up to 1 for each pixel. In this work, the membership values of the classification follow the trapezoidal membership function. In the literature, there are two possible ways of generating the membership function: Similarity Relation Model (SRM) and Semantic Import Model (SIM) [51,52,53]. The former derives the membership function using classifiers like for example fuzzy k-means, fuzzy c-means, and neural networks [53,54]. The first two are data driven, partitioning the observations based on multivariate attributes. The latter, SIM, generates membership based on expert knowledge [51,55].
The FCM results assign each pixel to membership of the two classes. The value of m determines the level of fuzziness in FCM classification. If m = 1, FCM is hard classifier. FCM was carried out by labeling two membership images resulting from each FCM classification as the water and non-water images. To do so, the combination of near infrared (NIR) and shortwave infrared (SWIR) of Landsat bands were used. The water label was given to the class which has the minimum value of the sum of the cluster means in the infrared bands. Detailed descriptions regarding the FCM algorithm are available in Bezdek et al. [50], whereas detailed explanations regarding membership function, pixel labeling, and parameter estimation for FCM classification can be found in Dewi et al. [16].

2.4. Validation

To quantify the accuracy of the FCM classifier, a conventional error matrix cannot be used. In this study, we used a fuzzy error matrix which has non-negative real numbers [48,56,57], since pixels have a partial membership to two classes.
For accuracy assessment, soft reference images were generated by applying an FCM classification to Pleiades, Spot 6, and Sentinel 2 images which were all captured during low tides. Let the value of μ i k and μ j l represent membership values of the kth pixel for class i in the classified image and lth pixel for class j in the reference images. It was assumed that the rows of the matrix are classes of the classified image and the columns are classes of the reference image. The fuzzy error matrix (FERM) is obtained using minimum operator showing the maximum possible overlap between the classified and reference images and indicating the agreement between classes in both images [48,49,56]:
A i = j =   M I N ( μ i k ,   μ j l )
To calculate the agreement in FERM, a group of 225 Pleiades pixels (15 × 15), 9 Spot 6 pixels (3 × 3) and 9 Sentinel 2 pixels (3 × 3) were averaged to achieve pixel dimension of Landsat images. Using this reference data, membership of 200 pixels randomly selected from both classified and reference images were computed to obtain the overall accuracies (OA) of the FCM classifications:
O A = i = 1 c A i = j k s
where k s represents the number of pixels used to generate the FERM.
Moreover, we compared the quality of the FCM results with respect to alternative pixel-based classification methods. Firstly, we classified the multi-spectral bands of Landsat using the MLC classifier being the most commonly used supervised classification technique for remote sensing images [58]. Secondly, we classified the multi-spectral bands of Landsat images using FCM and then labelled each pixel to the class to which it has the highest membership. It was assumed that hard output is the highest membership value which is actually computed from the soft output [59,60,61,62]. We called this the hardened classification. After classification, post classification comparisons were applied to detect the changes of the shorelines by superimposing the classification results in GIS.

2.5. Deriving Fuzzy Shoreline

FCM classification derives two raster layers, namely: water and non-water membership images. Each layer consists of fuzzy regions with fuzzy boundaries. Estimation of the spatial extent of objects i.e., water, non-water and shoreline, and their representations is related to the interpretation of the fuzziness of objects [52]. To derive shorelines at the locations where water and non-water objects meet, we modified the fuzzy-crisp object model based upon Cheng [52]. The two classes (water and non-water) are spatially disjoint, but their boundary is vaguely defined, whereas their interiors are crisp. Given this concept, we consider the boundary between water and non-water as fuzzy and form a transition zone that we call shoreline. To determine the spatial extent of water, non-water and shoreline, it is necessary to combine class objects from different layers into a single layer. The decision function d w k assigns pixel k with water membership value μ w k to a sub-area of water class based upon the following conditions:
If   ( μ w k > 0.99 )   then   ( d w k = 1 )
which means that the pixels belong to sub-areas water. Threshold 0.99 was set to represent the highest water membership values indicating the core of water.
If   ( 0.01 < μ w k < 0.99 )   then   ( d w k = μ w k )
This equation classifies pixels as shoreline.
If   ( μ w k < 0.01 )   then   ( d w k = 0 )
Pixels not belonging to water or shoreline constitute non-water. A threshold of 0.01 represents the lowest water membership values. This indicates pixels with membership below that threshold not belong to water or shoreline areas. The results after deriving fuzzy shorelines by applying Equations (3)–(5) were called as shoreline images.

2.6. Uncertainty Estimation

The uncertainty in class assignment was estimated by a measure of the confusion index CI for each pixel resulting from FCM classification as follows [63,64,65]:
C I = 1 ( μ i k 1 μ i k 2 )
If C I approaches 1 then the difference in membership values between the first and the second highest membership values are small meaning that both membership values are almost equal. Thus, it is more likely that the pixel defines a fuzzy boundary and the uncertainty of the pixel to belong to the class with the largest membership is high. If C I approaches 0, however, then the difference in membership values between the first and the second highest membership values are high and the uncertainty of the pixel to belong to the class with the largest membership is low.

2.7. Shoreline Change Detection

For establishing the changes over time, shoreline images obtained using Equations (3)–(5) of the same year were stacked and compared with the stack of shoreline images of the next year with corresponding seasons. If membership values to water ( μ w k ) of shoreline images within year T 1 and T 2 are given by G = ( g 1 ,   g 2 , .. , g z ) T 1 and H = ( h 1 ,   h 2 , .. , h z ) T 2 , respectively, and z is the number of shoreline images, a change vector is defined as:
Δ C V = H G = ( h 1   g 1 h 2   g 2 .. h z   g z )
Here, Δ C V includes all the change information between two years for a given pixel. The final result of CVA is an image of vector changes. The shoreline change is defined as the vector difference between successive time periods and is represented by a vector in a multi-dimensional space. The length of the change vector indicates the magnitude of change and its direction indicates the nature of the change [25,66].

2.7.1. Change Magnitude

The change magnitude Δ C V was derived by determining the Euclidean distance between shoreline images as:
Δ C V = ( h 1   g 1 ) 2 + ( h 2   g 2 ) 2 + ... + ( h z   g z ) 2
Δ C V represents the total membership differences between two years and measures the intensity of the shoreline change. Two categories of change were identified, namely change and no-change. A change category was defined when the water membership difference between T 1 and T 2 is larger than zero, whereas a no-change category is related to a magnitude equal to zero. A higher change magnitude corresponds with a large water membership difference between shoreline images in T 1 and T 2 . When the change magnitude is low, the water membership difference between shoreline images in T 1 and T 2 is small.

2.7.2. Change Direction

For all pixels classified as change, we estimated the change directions. Change direction was determined by evaluating the water membership difference between shoreline images in two successive years. It quantifies the variation of water membership in each pixel and shows how frequent the changes have occurred. Change direction estimation started by calculating the number of change combinations (CC) as:
C C = d p
where d refers to the types of change direction which can be distinguished when comparing the stack of shoreline images from both years for corresponding seasons and p refers to the number of shoreline image pairs. We identified three types of change direction to water: positive change direction (or in short positive direction), negative change direction (negative direction) and unclear change direction (unclear direction).
The change vector ( C V ) showing water membership difference between a pair of shoreline images in T 1 and T 2 from corresponding seasons needs to be estimated: (a) if the water membership difference between pair of shoreline images within years T 1 and T 2 is less than zero then C V = 1 showing a decrease of water membership in T 2 ; (b) if the water membership difference is larger than zero then C V = + 1 showing an increase of water membership in T 2 ; (c) if the water membership difference is equal to zero then C V = 0 showing that the water membership in T 1 and T 2 were the same. The total change vector (TCV) values are defined as:
T C V = C V 1 + C V 2 + ... + C V z
CV1 refers to ( h 1   g 1 ) , CV2 refers to ( h 2   g 2 ) , and CVz refers to ( h z   g z ) . Finally, the change direction (Chg.dir) categories showing the degree of change direction to water membership were obtained by grouping the direction values: (a) TCV values from +1 up to + z were grouped as positive direction; (b) TCV values from −1 up to z were grouped as negative direction; (c) TCV values equal to 0 showing unclear change directions were classified as unclear direction; and (d) TCV values equal to 0 having water membership differences equal to 0 at all time periods were classified as no-change. Table 3 shows the procedure to determine the change direction categories by using four pairs of image used in this study.
Based upon these results, the change area of a specific change direction category (positive direction, negative direction, and unclear direction) and the no-change area were defined as:
A r ( C h g ) =   P i x k ( C h g ) ×   A r ( k )
where P i x k ( C h g ) is the number of pixels belonging to the area of change and no-change, and A r ( k ) is area of pixel k (30 × 30 m2).

2.7.3. Change Uncertainty

Based upon the change detection results, the change uncertainty of related areas was estimated by the confusion index CI. If C I of two images for T 1 and T 2 are given by Q = ( q 1 ,   q 2 , .. , q z ) T 1 and R = ( r 1 ,   r 2 , . , r z ) T 2 , respectively, then the change confusion is derived as:
Δ C U = ( q 1   r 1 ) 2 + ( q 2   r 2 ) 2 + .. + ( q z   r z ) 2
A high Δ C U value is related to a large difference of confusion indices between images for T 1 and T 2 , whereas a low change confusion corresponds to a small difference of confusion indices between images for T 1 and T 2 .

3. Results

3.1. FCM Classification and Accuracy Assessment

Table 4 presents the accuracy assessment of classification results using FCM and alternative classification methods. The FCM classifier outperformed MLC and the accuracy values of FCM are generally higher than the hardened classification.
Figure 2 presents an example of FCM outputs, together with MLC and hardened classification. In the image, MLC overestimated the non-water area shown by the larger area of non-water (see Figure 2d–f, e.g., grid cells B1 and B2), whereas the hardened classification underestimated the non-water area (see Figure 2g–i, e.g., grid cells B1 and B2). Both methods failed to distinguish the gradual transition between water and non-water.
The results of FCM classification are presented in Figure 3 with the values ranging from 0 to 1 for both membership images of water (Figure 3a–c) and non-water (Figure 3d–f). Areas with higher water membership values were located for example in marine areas, fishponds, and water-covered agricultural areas (Figure 3a, e.g., grid cell A2). In Figure 3d, higher non-water membership pixels are located near settlements adjacent to the shorelines, and mangrove forests (Figure 3d, e.g., grid cells B2 and C3).

3.2. Fuzzy Shoreline and Uncertainty Estimation

Figure 3g–i and Figure 4a,d show the results of the fuzzy-crisp objects model to derive shorelines. Figure 4b,e presents the C I . Dark pixels with C I close to 0 indicate the areas classified as lower uncertainty (Figure 4e, e.g., grid cell A1), whereas brighter pixels with confusion index close to 1 indicate areas classified as the fuzzy boundary with higher uncertainty (Figure 4e, e.g., grid cell B2). Figure 4c,f shows shorelines images with fuzziness represented by C I values. These ambiguous areas indicate shoreline positions represented by pixels in grey shades (Figure 4f, grid cells e.g., A2 and B2). The width of these shorelines is determined by natural conditions of the coastal areas, for example, a wider shoreline is more likely to be found in a muddy coastal area or at a gently sloping beach, whereas a narrow shoreline is usually found along a steeper slope beach and coastal area with embankment and other man-made structures.

3.3. Shoreline Change Detection

3.3.1. Change Magnitude and Change Uncertainty

The change magnitude and the change and no-change categories of shoreline are displayed in Figure 5. Low change magnitude values correspond to a small water membership difference between shoreline images at T 1 and T 2 . They cover marine areas (Figure 5a, e.g., grid cells A2 and B3) and a relatively undisturbed coastal land (Figure 5c e.g., grid cells D1 and D2). In addition, high change magnitude values correspond to a large water membership difference. Those pixels cover muddy areas (Figure 5a,c, e.g., grid cells C1 and D1) and coastal land which was highly-influenced by tidal floods (Figure 5b,d, e.g., grid cells B2 and B3).
The fuzziness of shoreline changes is presented in Figure 6. Low change confusion correspond to small C I differences between images in T 1 and T 2 . This indicates a low uncertainty that the changes have occurred as can be seen in Figure 6a,c, e.g., grid cells A1, A2 and B1. High values are associated with large C I differences and indicate a high uncertainty that the changes have occurred (Figure 6b,d, e.g., grid cells B1, B2 and B3).

3.3.2. Change Direction

The representation of change direction of shoreline showing the variation of water membership in each pixel can be seen in Figure 7. Positive directions to water membership correspond to the increase of water membership at time T 2 (Figure 7a, e.g., grid cells A2, B1, and C1). On the contrary, negative directions to water membership were associated with the decrease of water membership values at time T 2 (see Figure 7d, e.g., grid cells A2 and B1). The no-change category indicates an undisturbed environment (see Figure 7a, e.g., grid cell B1), whereas the unclear direction category indicates an ambiguous condition since the changes occurred without an obvious trend (see Figure 7b, e.g., grid cells A2 and C1).
The change areas for each category are presented in Table 5. A positive direction to water covers an area of approximately 1828 ha in the period 2013–2014 and 1120 ha in the period 2014–2015. A negative direction has occurred for 920 ha and 1635 ha in the period 2013–2014 and in the period 2014–2015, respectively. Unclear direction category presented as pink pixels can be seen in Figure 7c, e.g., grid cells A2 and C1 covering an area of 616 and 528 ha in 2013–2014 and in 2014–2015, respectively. No-change direction shows a stable area which is mostly located at the sea or inundated fishponds represented by light yellow color (see Figure 7a,d, e.g., grid cells A1 and B1) covering an area of 1319 and 1403 ha in the period 2013–2014 and in the period 2014–2015, respectively.

3.3.3. Change Confusion

The intensity of the change confusion was identified for each change direction category in the period 2013–2014 and in the period 2014–2015, respectively (Figure 8a,c). Three change confusion values were identified including positive direction, negative direction and unclear direction.

3.3.4. Comparison with Alternative Change Detection Methods

Figure 9 and Figure 10 show the change detection of shoreline using post classification comparison of MLC and hardened classification, respectively. Both MLC and hardened classification present shoreline as a single line. The changes of this single shoreline have occurred due to the changes of water and non-water area. Binary images from two dates T 1 and T 2 were superimposed in GIS and four types of change were identified, namely: non-water to water, water to non-water, water to water and non-water to non-water. Figure 9a–c and Figure 10a–c show the changes of shoreline in three consecutive dates in 2013, whereas Figure 9d,e and Figure 10d,e present two examples of shoreline changes from 2013 to 2014 (Figure 9d and Figure 10d) and from 2014 to 2015 (Figure 9e and Figure 10e).
Figure 11 shows the comparison between the proposed method and the alternative method at the selected study area. In this example, both methods agree on the results of change detection as can be seen in Figure 11, e.g., grid cells B3 and C2. From CVA results (Figure 11c,d), the area in yellow polygons shows a negative change to water membership with high change magnitudes as shown in Figure 11e,f. The negative change to water membership means a change towards non-water. Similarly, post classification results also denote that these yellow polygon sites experienced a change from water to non-water without further information on the intensity of the change (Figure 11a,b).

3.3.5. Multi-Year Pattern of Water Membership Changes

Each pixel from the resulting change vectors provides information regarding its change direction and magnitude. Each combination represents specific types of change processes that may occur in the field and shows a multi-year pattern of water membership changes over the observation periods. Four combinations of change and their related processes are interpreted as follows:
(a) High change direction and high change magnitude
The areas with high change direction and high change magnitude values are observed for both positive and negative directions. Both conditions indicate a continuous change of an area to a certain direction with a relatively large intensity. A consistency to positive direction indicates a persistence of enhanced water influence as those pixels show an increase of water membership in multi-temporal images (see Figure 12a–d). This probably corresponds to the land subsidence and coastal inundation. As the land subsides and the water level increases, some mangrove trees located closely to the sea are falling down. The RGB 542 of Landsat images in Figure 12e,f depict these changes indicated by the decrease of vegetation cover between 2013 and 2015.
Figure 13 presents areas characterized by continuously decreasing water membership in multi-temporal images categorized as negative direction (Figure 13a,b) with high change magnitude (Figure 13c,d). This may indicate a success in shoreline protection scheme that caused sediment accretion to occur allowing mangroves to grow as can be seen from the RGB 542 of Landsat images in white-dashed ellipses in Figure 13e–h.
(b) Low change direction and high change magnitude
This category indicates an abrupt change which may be influenced by random events. Figure 14b shows a low positive direction with high change magnitude (Figure 14d) which may result from coastal flooding triggered by spring tides, extreme waves and winds. Since the magnitude of the changes is high and the change is sudden, this type of change may indicate a higher risk. Images made available by Google Earth from 2013 to 2015 in Figure 14e–h show the decrease of mangrove coverage. In fact, mangroves can act as sediment trap and can reduce the energy of the high waves, therefore, when the mangroves disappear, the threat from tidal floods increases. Figure 15a,b shows the embankment which protects settlements from high tide; however, during an extreme event for example when a higher tide combines with an extreme wind, the water level may increase and overflow this embankment.
(c) High change direction and low change magnitude
A gradual, continuous increase of wet conditions was observed by an increase in water membership with low magnitude values. This type of change was categorized as positive direction which could be due to cyclical tides and coastal processes, for example flooded land (Figure 16a,c), and water turbidity (Figure 16b,d). Even though the magnitude of the change is low, the changes occur frequently. Hence, this type of change may give a higher risk. In a longer observation, if the areas persistently become wetter, this location may have a risk of coastal inundation as well.
(d) Low change direction and low change magnitude
This type of change probably indicates an undisturbed environment with a low change magnitude (see black-dashed circle sites in Figure 17c,d). This category mainly occurs in water areas, probably due to the changes in water turbidity (see black-dashed circle sites in Figure 17a,b). In addition, this type of change was observed in small patches of the coastal land probably resulting from changes in soil moisture (see black rectangle sites in Figure 17b,d).

4. Discussion

In this study, the dynamics of fuzzy shorelines have been assessed using fuzzy classification and a raster-based change detection technique. FCM classification was used to discriminate the land and water classes and to estimate their memberships. FCM is a well-known clustering method which is less dependent on the initial state of clustering [62] and capable of describing phenomena such as water and non-water which is changing gradually. Instead of FCM, there are various ways to derive a fuzzy classification for example from fuzzy maximum likelihood classification [59], and artificial neural network fuzzy classification [67]. Membership values obtained by applying FCM are used to deal with the uncertain information on the position of fuzzy shorelines. In FCM classification, we set c = 2 since we were interested in identifying the boundary between water and non-water and because both classes give the largest spectral differences in image [16,68]. Finding the suitable number of clusters in the beginning of the classification could be difficult. A priori knowledge regarding the study area, for example by observing an aerial photo, can be used to define the suitable number of clusters [69,70]. For other situations, by assessing the homogeneity measure using a posteriori indicators, the number of clusters could be determined using entropy and non-fuzziness index [70], exponential cluster validity [71], and spatial fuzzy clustering [72].
Shorelines and their changes were presented as fuzzy areas. The fuzzy-crisp object model in this study was successful in identifying the extent of shoreline positions as the transition zone between water and non-water. Setting the threshold to the highest (0.99) and the lowest (0.01) memberships are intended to find the core of water and non-water, respectively. The uncertainty addressed in this research corresponds to the existential and extensional uncertainty of shoreline objects as has been mentioned by Molenaar and Cheng [73] and Cheng [52]. Existential uncertainty expresses the uncertainty of the existence of shoreline in reality. It refers to the possibility of existence of a shoreline to be detected on an image. Extensional uncertainty implies that the area indicated as shoreline can be determined with limited certainty, for example with boundaries that reflect the transition zone between water and non-water. Moreover, when the values of an adjacent grid are very similar, the zones of confusion divide regions indicating the presence of gradual transitions. The extensional uncertainty in shoreline identification includes differences in applied threshold values when defining the core of water and non-water, the applied shoreline definition, tides condition during image acquisition, time series of observation, and the nature of the beaches (such as flat or steep slope beaches, and muddy or rocky beaches). In addition, the changed areas of the fuzzy shoreline are thus associated with the distribution of changed confusion indices. The change uncertainty represented by changed confusion indices shows the degree of uncertainty of the changes that have occurred. It can be seen from the results that a location having a higher change magnitude, has a higher change confusion value as well. It corresponds to the higher differences of both water memberships and confusion indices between corresponding images in T 1 and T 2 . Explicit handling of uncertainty by addressing the shoreline as a transition zone allows decision makers and planners to include this uncertainty in spatial planning. Moreover, it visualizes not only the changes in shoreline, but also the uncertainty of these changes for every location, thereby providing a better base for a debate on the combined effects of land subsidence and sea level rise in this area.
The change of shoreline was explained in terms of change magnitude and change direction using CVA. Information provided by CVA allows us to see the trend of the fluctuating shoreline over time, whereas the change detection results of the alternative method could provide only “from-to” change information and the detail of subtle within-class changes was lost because it only compared images from two dates [24,25]. In our previous study, we used post classification comparison because we were interested in observing the changes of shoreline over a longer period from 1994 up to 2015 [16]. Given different methods have been implemented in monitoring the change position of shoreline, both studies confirmed that shoreline changes associated with coastal inundation have occurred in this study area.
The analysis of information provided by the change magnitude and direction reveals that each change combination represents one specific change process type. The processes could vary depending on the characteristic of the coastal areas. For example, shorelines could change due to floods triggered by land subsidence, and floods caused by seasonal variation, abrupt shoreline changes due to extreme tides and waves, and the changes of water turbidity and soil moisture triggered by daily weather events. These specific type processes were explained on the basis of the analysis of four images for each observation period. In fact, the number of images could easily be extended to more than four images. The seasonal variation of shorelines and other information regarding whether the changes would lead to a permanent coastal inundation could not fully be assessed in this study because only four observations each year were compared. This shows that the change vector analysis is sensitive to the length of the stacked period and the number of images stacked over that period, as is also confirmed by Lambin and Strahler [25].
To have images captured during similar tides and with corresponding seasons for each pair is important, hence shorelines can be compared equally. Images captured during wet seasons gave more uncertainty to the classification results since the rain affected the wetness of soil surfaces, thus the classified images produced false impression of higher flooding. Therefore, to reduce the uncertainty due to seasonal influence, images captured during dry season should be preferred. Furthermore, the uncertainty of the results was also influenced by tide differences. Images captured during high tides and low tides produced different positions of the boundary between water and non-water. Hence, this increased the uncertainty of shoreline position as well, because shoreline by its definition is an intersection of coastal land and water surface. To have images acquired at exactly the same tide level is hardly possible. Therefore, all images were acquired at the low tides with negligible differences. We only considered astronomical tide level assuming it was more influential than meteorological factors, as confirmed by Pugh [74]. Astronomical and meteorological factors have different influences on different slope conditions and the magnitude of the influences may become larger if the slope is gentler. However, if there would be any remaining small influence of meteorological factors, it would be accounted for by the use of fuzzy classification in deriving fuzzy shoreline.
In the accuracy assessment, soft reference data were generated from various higher resolution images. To obtain the required image resolution, resampling and aggregating have been implemented. Resampling of image and aggregation of pixel values were potential sources of error but were ignored in this work since the error was likely to be very small [49,75]. Furthermore, using a soft classifier to generate soft reference data is more likely to reduce the uncertainty due to the vagueness in class definition and mixed pixel problems. Hence, the finer resolution dataset was not assumed to be pure and no information was lost due to the hardening of the soft classification [60,61,76]. This could be an explanation for the higher accuracies obtained by using the proposed classification in comparison with alternative methods, as confirmed by Zhang and Foody [61] and Chawla [75]. Furthermore, although the differences of the accuracy results were only small, the advantage of a fuzzy approach is not only in the improved accuracy of the shoreline, but also in the fact that it makes clear what the margins of uncertainty around the shoreline are, which provides a better basis for decision.
Analysis of the shoreline changes in the northern coastal area of Central Java shows the changes of shoreline positions from 2013 to 2015. This could be related to the processes that shape shorelines determined by the interaction of several factors, including: (a) the change of sea-level; (b) the amount of sediment supplied to the beach by rivers; (c) the movement of the sediment by marine processes; and (d) the role of waves, currents, tides and winds in moving the sediment [77]. Furthermore, sediment transport is not constant, and it is constantly subject to change. The alteration of sediment transport can come from changes in water flow, water level, weather events and human influence. In addition, previous studies have mentioned that this location has suffered from the changing of shorelines for more than 20 years due to coastal inundation accelerated by, for example, land subsidence, sea level rise, seaport development, and ground water extraction [16,20,38]. Many attempts have been made to combat coastal inundation and erosion along the 1.3 km coast in Demak. Elevated roads, raised floor of the houses, breakwater, and mangrove planting have been applied. In 2013, a Dutch–Indonesian consortium agreed to start a pilot project “Building with Nature” building a permeable dam of natural material called “hybrid engineering” [78,79]. This development could be one reason for the increase of negative changes to water membership in the period of 2014–2015. Hybrid engineering is one type of coastal protection combining technical and ecosystem-based solutions referred to as sediment traps [79]. Netherlands-Water-Partnership [78] reported that after a year, new sediment layers were deposited at the surrounding areas.

5. Conclusions

In this article, we present a method to identify shoreline positions and their changes as a fuzzy area including a measure of change confusion. Shoreline changes could be detected, and the method provided information regarding the change magnitude and the trend of water membership in each pixel. Our results reveal that this information represents specific type of change processes showing multi-year patterns of water membership changes over time. These include: (a) high change direction and high change magnitude with a consistent positive change direction, which probably corresponds to land subsidence and coastal inundation, while a consistent negative direction may indicate a success in shoreline protection scheme; (b) low change direction and high change magnitude indicates an abrupt change which may be influenced by random events, such as flooding triggered by spring tides, extreme waves and winds; (c) high change direction and low change magnitude, which, if in positive direction, could be due to cyclical tides and coastal processes, for example flooded land; and (d) low change direction and low change magnitude probably indicates an undisturbed environment, such as water areas with changes in water turbidity and coastal land with changes in soil moisture. Finally, we conclude that the proposed method can assess changes in a shoreline by taking into account that it is a fuzzy boundary.
The change area estimation, change magnitude and direction of the shorelines may support local government and stakeholders in monitoring the change of fuzzy shorelines. Combining information given by this research with other information such as distribution of population could help to determine priority locations prioritized in disaster preparedness and response. To include digital elevation model in the processing phase is important for further research, because it allows the change analysis to focus more on the area affected by tidal floods. We further realized that astronomical and meteorological factors have different influences on different slope conditions. Observation data regarding the wave run-up and other incident wave conditions in the study area were not available for each observation period. This information may be collected and used in a near-shore wave model. This might be included in future studies as well.

Acknowledgments

This study has been conducted at the Faculty of Geo-information Science and Earth Observation (ITC), University of Twente, supported by the scholarship Program for Research and Innovation in Science and Technologies (RISET-PRO), Ministry of Research, Technology and Higher Education of Indonesia. The authors express their gratitude to Valentyn Tolpekin from Faculty of Geo-information Science and Earth Observation (ITC), University of Twente for providing the script code for FCM.

Author Contributions

Ratna Sari Dewi, as the first author, had responsibility for conducting the research including a writing task. The two other authors, Wietske Bijker and Alfred Stein, supervised the work, and revised the manuscript. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zanuttigh, B.; Martinelli, L.; Lamberti, A.; Moschella, P.; Hawkins, S.; Marzetti, S.; Ceccherelli, V.U. Environmental design of coastal defence in Lido di Dante, Italy. Coast. Eng. 2005, 52, 1089–1125. [Google Scholar] [CrossRef]
  2. Jin, D.; Hoagland, P.; Au, D.K.; Qiu, J. Shoreline change, seawalls, and coastal property values. Ocean Coast. Manag. 2015, 114, 185–193. [Google Scholar] [CrossRef]
  3. Baldassarre, G.D.; Schumann, G.; Bates, P.D. A technique for the calibration of hydraulic models using uncertain satellite observations of flood extent. J. Hydrol. 2009, 367, 276–282. [Google Scholar] [CrossRef]
  4. Kaergaard, K.; Fredsoe, J. A numerical shoreline model for shorelines with large curvature. Coast. Eng. 2013, 74, 19–32. [Google Scholar] [CrossRef]
  5. Snoussi, M.; Ouchani, T.; Khouakhi, A.; Niang-Diop, I. Impacts of sea-level rise on the Moroccan coastal zone: Quantifying coastal erosion and flooding in the Tangier Bay. Geomorphology 2009, 107, 32–40. [Google Scholar] [CrossRef]
  6. O’Connor, M.C.; Cooper, J.A.G.; McKenna, J.; Jackson, D.W.T. Shoreline management in a policy vacuum: A local authority perspective. Ocean Coast. Manag. 2010, 53, 769–778. [Google Scholar] [CrossRef]
  7. Tamassoki, E.; Amiri, H.; Soleymani, Z. Monitoring of shoreline changes using remote sensing (case study: Coastal city of Bandar Abbas). In 7th IGRSM International Remote Sensing & GIS Conference and Exhibition; IOP Publishing: Bristol, UK, 2014; Volume 20. [Google Scholar]
  8. Davidson-Arnott, R. An Introduction to Coastal Processes and Geomorphology; Cambridge University Press: Cambridge, MA, USA, 2010; p. 442. [Google Scholar]
  9. Bird, E.C.F. Coastline Changes: A Global Review; John Wiley & Sons: New York, NY, USA, 1985; p. 219. [Google Scholar]
  10. Boak, E.H.; Turner, I.L. Shoreline definition and detection: A review. J. Coast. Res. 2005, 21, 688–703. [Google Scholar] [CrossRef]
  11. Oertel, G.F. Coasts, coastlines, shores, and shorelines. In Encyclopedia of Coastal Science; Schwartz, M.L., Ed.; Springer: New York, NY, USA, 2005; pp. 323–327. [Google Scholar]
  12. Bird, E.C.F. Coastal Geomorphology: An Introduction; John Wiley & Sons: Chicester, UK, 2000. [Google Scholar]
  13. Liu, X.; Xia, J.C.; Blenkinsopp, C.; Arnold, L.; Wright, G. High Water Mark Determination Based on the Principle of Spatial Continuity of the Swash Probability. J. Coast. Res. 2014, 30, 487–499. [Google Scholar] [CrossRef]
  14. Moore, L.J.; Ruggiero, P.; List, J.H. Comparing Mean High Water and High Water Line Shorelines: Should Proxy-Datum Offsets be Incorporated into Shoreline Change Analysis? J. Coast. Res. 2006, 22, 894–905. [Google Scholar] [CrossRef]
  15. Ghosh, M.K.; Kumar, L.; Roy, C. Monitoring the coastline change of Hatiya Island in Bangladesh using remote sensing techniques. ISPRS J. Photogramm. 2015, 101, 137–144. [Google Scholar] [CrossRef]
  16. Dewi, R.S.; Bijker, W.; Stein, A.; Marfai, M.A. Fuzzy Classification for Shoreline Change Monitoring in a Part of the Northern Coastal Area of Java, Indonesia. Remote Sens. 2016, 8, 1–25. [Google Scholar] [CrossRef]
  17. Fisher, P.F. Models of uncertainty in spatial data. Geogr. Inf. Syst. 1999, 1, 191–205. [Google Scholar]
  18. Li, R.; Ma, R.; Di, K. Digital tide-coordinated shoreline. Mar. Geod. 2002, 25, 27–36. [Google Scholar] [CrossRef]
  19. Senthilnatha, J.; H, V.S.; Omkar, S.N.; Mani, V. Spectral-spatial MODIS image analysis using swarm intelligence algorithms and region based segmentation for flood assessment. In Proceedings of the Seventh International Conference on Bio-Inspired Computing: Theories and Applications; Bansal, J.C., Singh, P.K., Deep, K., Pant, M., Nagar, A.K., Eds.; Springer: New Delhi, India, 2012. [Google Scholar]
  20. Marfai, M.A.; Almohammad, H.; Dey, S.; Susanto, B.; King, L. Coastal dynamic and shoreline mapping: Multi-sources spatial data analysis in Semarang Indonesia. Environ. Monit. Assess. 2008, 142, 297–308. [Google Scholar] [CrossRef] [PubMed]
  21. Shenbagaraj, N.; Mani, N.D.; Muthukumar, M. Isodata classification technique to assess the shoreline changes of Kolachel to Kayalpattanam coast. Int. J. Eng. Res. Technol. 2014, 3. [Google Scholar]
  22. Muslim, A.M.; Foody, G.M.; Atkinson, P.M. Localized soft classification for super-resolution mapping of the shoreline. Int. J. Remote Sens. 2006, 27, 2271–2285. [Google Scholar] [CrossRef]
  23. Taha, L.G. E.-d.; Elbeih, S.F. Investigation of fusion of SAR and Landsat data for shoreline super resolution mapping: The Northeastern Mediterranean sea coast in Egypt. Appl. Geomat. 2010, 2, 177–186. [Google Scholar] [CrossRef]
  24. Jianya, G.; Haigang, S.; Guorui, M.; Qiming, Z. A Review of Multi-Temporal Remote Sensing Data Change Detection Algorithms. Int. Arch. Photogramm., Remote Sens. Spat. Inf. Sci. 2008, 37, 757–762. [Google Scholar]
  25. Lambin, E.F.; Strahler, A.H. Change-Vector Analysis in Multitemporal Space: A Tool to Detect and Categorize Land-Cover Change Processes Using High Temporal-Resolution Satellite Data. Remote Sens. Environ. 1994, 48, 231–244. [Google Scholar] [CrossRef]
  26. Singh, A. Digital Change Detection Techniques using Remotely-Sensed Data. Int. J. Remote Sens. 1989, 10, 989–1003. [Google Scholar] [CrossRef]
  27. Lu, D.; Li, G.; Moran, E. Current situation and needs of change detection techniques. Int. J. Image Data Fusion 2014, 5, 13–38. [Google Scholar] [CrossRef]
  28. Kauth, R.J.; Thomas, G.S. The Tasselled Cap—A Graphic Description of the Spectral-Temporal Development of Agricultural Crops as Seen by Landsat; LARS Symposia; Purdue University: Lafayette, IN, USA, 1976; pp. 41–51. [Google Scholar]
  29. Malila, W.A. Change Vector Analysis: An Approach for Detecting Forest Changes with Landsat; LARS Symposia; Purdue University: Lafayette, IN, USA, 1980. [Google Scholar]
  30. Lambin, E.F.; Ehrlich, D. Land-cover Changes in Sub-Saharan Africa (1982–1991): Application of a Change Index Based on Remotely Sensed Surface Temperature and Vegetation Indices at a Continental Scale. Remote Sens. Environ. 1997, 64, 181–200. [Google Scholar] [CrossRef]
  31. Landmann, T.; Schramm, M.; Huettich, C.; Dech, S. MODIS-based change vector analysis for assessing wetland dynamics in Southern Africa. Remote Sens. Lett. 2013, 42, 104–113. [Google Scholar] [CrossRef]
  32. Vorovencii, I. A change vector analysis technique for monitoring land cover changes in Copsa Mica, Romania, in the period 1985–2011. Environ. Monit. Assess. 2014, 186, 5951–5968. [Google Scholar] [CrossRef] [PubMed]
  33. He, C.; Wei, A.; Shi, P.; Zhang, Q.; Zhao, Y. Detecting land-use/land-cover change in rural–urban fringe areas using extended change-vector analysis. Int. J Appl. Earth Obs. Geoinf. 2011, 13, 572–585. [Google Scholar] [CrossRef]
  34. Singh, S.; Talwar, R. Performance Analysis of Different Threshold Determination Techniques for Change Vector Analysis. J. Geol. Soc. India 2015, 86, 52–58. [Google Scholar] [CrossRef]
  35. Özyavuz, M.; Satir, O.; Bilgili, B.C. A Change Vector Analysis Technique to Monitor Land-Use/Land-Cover in the Yildiz Mountains, Turkey. Fresenius Environ. Bull. 2011, 20, 1190–1199. [Google Scholar]
  36. Hartini, S. Flood Risk Modelling on Agricultural Area in the North Coastal Area of Central Java. Ph.D. Thesis, Gadjah Mada University, Yogyakarta, Indonesia, 2015. [Google Scholar]
  37. Marfai, M.A. Preliminary assessment of coastal erosion and local community adaptation in Sayung coastal area, Central Java–Indonesia. Quaest. Geogr. 2012, 31, 47–55. [Google Scholar]
  38. Harwitasari, D.; van Ast, J.A. Climate change adaptation in practice: People’s responses to tidal flooding in Semarang, Indonesia. J. Flood Risk Manag. 2011, 4, 216–233. [Google Scholar] [CrossRef]
  39. Marfai, M.A.; King, L.; Sartohadi, J.; Sudrajat, S.; Budiani, S.R.; Yulianto, F. The impact of tidal flooding on a coastal community in Semarang, Indonesia. Environmentalist 2008, 28, 237–248. [Google Scholar] [CrossRef]
  40. Sarbidi. Rob influence on Coastal Settlements (Case Study: Semarang City). In Losses in Buildings and Areas Due Sea-level Rise on Coastal Cities in Indonesia, Bandung, 19–20 March 2001; Research Center for Settlements, Ministry of Public Works: Bandung, Indonesia, 2001. [Google Scholar]
  41. Kurniawan, L. Tidal flood assessment in Semarang (case: Dadapsari). Alami 2003, 8, 54–59. [Google Scholar]
  42. Marfai, M.A.; King, L. Tidal inundation mapping under enhanced land subsidence in Semarang, Central Java Indonesia. Nat. Hazards 2007, 44, 93–109. [Google Scholar] [CrossRef]
  43. USGS United States Geological Survey. Available online: https://earthexplorer.usgs.gov/ (accessed on 22 November 2015).
  44. Hadjimitsis, D.G.; Papadavid, G.; Agapiou, A.; Themistocleous, K.; Hadjimitsis, M.G.; Retalis, A.; Michaelides, S.; Chrysoulakis, N.; Toulios, L.; Clayton, C.R.I. Atmospheric correction for satellite remotely sensed data intended for agricultural applications: Impact on vegetation indices. Nat. Hazards Earth. Syst. Sci. 2010, 10, 89–95. [Google Scholar] [CrossRef]
  45. Mather, P.M. Computer Processing of Remotely-Sensed Images: An Introduction, 3rd ed.; John Wiley & Sons: Chicester, UK, 2004; p. 352. [Google Scholar]
  46. Lu, D.; Weng, Q. A survey of image classification methods and techniques for improving classification performance. Int. J. Remote Sens. 2007, 28, 823–870. [Google Scholar] [CrossRef]
  47. Congalton, R.G.; Green, K. Assessing the Accuracy of Remotely Sensed Data : Principles and Practices, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2009. [Google Scholar]
  48. Pontius, R.G.; Cheuk, M.L. A generalized cross-tabulation matrix to compare soft-classified maps at multiple resolutions. Int. J. Geogr. Inf. Sci. 2006, 20, 1–30. [Google Scholar] [CrossRef]
  49. Harikumar, A.; Kumar, A.; Stein, A.; Raju, P.L.N.; Murthy, Y.V.N.K. An effective hybrid approach to remote-sensing image classification. Int. J. Remote Sens. 2015, 36, 2767–2785. [Google Scholar] [CrossRef]
  50. Bezdek, J.C.; Ehrlich, R.; Full, W. FCM: The fuzzy c-means clustering algorithm. Comput. Geosci. 1984, 10, 191–203. [Google Scholar] [CrossRef]
  51. Burrough, P.A. Fuzzy mathematical methods for soil survey and land evaluation. J. Soil Sci. 1989, 40, 477–492. [Google Scholar] [CrossRef]
  52. Cheng, T. Fuzzy Objects: Their Changes and Uncertainties. Photogramm. Eng. Remote Sens. 2002, 68, 41–49. [Google Scholar]
  53. Zhu, A.-X.; Yang, L.; Li, B.; Qin, C.; Pei, T.; Liu, B. Construction of membership functions for predictive soil mapping under fuzzy logic. Geoderma 2010, 155, 164–174. [Google Scholar] [CrossRef]
  54. Chuang, K.-S.; Tzeng, H.-L.; Chen, S.; Wu, J.; Chen, T.-J. Fuzzy c-means clustering with spatial information for image segmentation. Comput. Med. Image Graph. 2006, 30, 9–15. [Google Scholar] [CrossRef] [PubMed]
  55. Medasani, S.; Kim, J.; Krishnapuram, R. An overview of membership function generation techniques for pattern recognition. Int. J. Approx. Reason. 1998, 19, 391–417. [Google Scholar] [CrossRef]
  56. Silván-Cárdenas, J.L.; Wang, L. Sub-pixel confusion–uncertainty matrix for assessing soft classifications. Remote Sens. Environ. 2008, 112, 1081–1095. [Google Scholar] [CrossRef]
  57. Binaghi, E.; Brivio, P.A.; Ghezzi, P.; Rampini, A. A fuzzy set-based accuracy assessment of soft classification. Pattern Recognit. Lett. 1999, 20, 935–948. [Google Scholar] [CrossRef]
  58. Richards, J.A. Remote Sensing Digital Image Analysis, 5th ed.; Springer: New York, NY, USA, 2013. [Google Scholar]
  59. Wang, F. Fuzzy Supervised Classification of Remote Sensing Images. IEEE Trans. Geosci. Remote Sens. 1990, 28, 194–201. [Google Scholar] [CrossRef]
  60. Dutta, A. Fuzzy c-Means Classification of Multispectral Data Incorporating Spatial Contextual Information by Using Markov Random Field; ITC: Enschede, The Netherlands, 2009. [Google Scholar]
  61. Zhang, J.; Foody, G.M. A fuzzy classification of sub-urban land cover from remotely sensed imagery. Int. J. Remote Sens. 1998, 19, 2721–2738. [Google Scholar] [CrossRef]
  62. Tso, B.; Mather, P. Classification Methods for Remotely Sensed Data; CRC Press: Boca Raton, FL, USA, 2009; p. 347. [Google Scholar]
  63. Burrough, P.A.; van Gaans, P.F.M.; Hootsmans, R. Continuous classification in soil survey: Spatial correlation, confusion and boundaries. Geoderma 1997, 77, 115–135. [Google Scholar] [CrossRef]
  64. Cheng, T.; Molenaar, M.; Lin, H. Formalizing fuzzy objects from uncertain classification results. Int. J. Geogr. Inf. Sci. 2001, 15, 27–42. [Google Scholar] [CrossRef]
  65. Zhang, J.; Kirby, R.P. Alternative Criteria for Defining Fuzzy Boundaries Based on Fuzzy Classification of Aerial Photographs and Satellite Images. Photogramm. Eng. Remote Sens. 1999, 65, 1379–1387. [Google Scholar]
  66. Lunetta, R.S.; Elvidge, C.D. Remote Sensing Change Detection: Environmental Monitoring Methods and Applications; Taylor & Francis: London, UK, 1999. [Google Scholar]
  67. Zhang, J.; Foody, G.M. Fully-fuzzy Supervised Classification of sub-urban Land Cover from Remotely Sensed Imagery: Statistical and Artificial Neural Network Approaches. Int. J. Remote Sens. 2001, 22, 615–628. [Google Scholar] [CrossRef]
  68. Bijker, W.; Hamm, N.A.; Ijumulana, J.; Wole, M.K. Monitoring a fuzzy object: The case of Lake Naivasha. In Proceedings of the 6th International Workshop on the Analysis of Multi-Temporal Remote Sensing Images (Multi-Temp), Trento, Italy, 12–14 July 2011.
  69. Wang, Q.; Shi, W. Unsupervised classification based on fuzzy c-means with uncertainty analysis. Remote Sens. Lett. 2013, 4. [Google Scholar] [CrossRef]
  70. Fisher, P.F.; Pathirana, S. The evaluation of fuzzy membership of land cover classes in the suburban zone. Remote Sens. Environ. 1990, 34, 121–132. [Google Scholar] [CrossRef]
  71. Zarandi, M.H.F.; Faraji, M.R.; Karbasian, M. An exponential cluster validity index for fuzzy clustering with crisp and fuzzy data. Trans. E Ind. Eng. 2010, 17, 95–110. [Google Scholar]
  72. Delinom, R.M.; Assegaf, A.; Abidin, H.Z.; Taniguchi, M.; Suherman, D.; Lubis, R.F.; Yulianto, E. The contribution of human activities to subsurface environment degradation in Greater Jakarta Area, Indonesia. Sci. Total Environ. 2008, 407, 3129–3141. [Google Scholar] [CrossRef] [PubMed]
  73. Molenaar, M.; Cheng, T. Fuzzy spatial objects and their dynamics. ISPRS J. Photogramm. 2000, 55, 164–175. [Google Scholar] [CrossRef]
  74. Pugh, D.T. Tides, Surges and Mean Sea-Level; John Wiley & Sons: London, UK, 1987. [Google Scholar]
  75. Chawla, S. Possibilistic C-Means—Spatial Contextual Information Based Sub-Pixel Classification Approach for Multi-Spectral Data. Master’s Thesis, University of Twente, Enschede, The Netherlands, 2010. [Google Scholar]
  76. Foody, G.M. Status of land cover classification accuracy assessment. Remote Sens. Environ. 2002, 80, 185–201. [Google Scholar] [CrossRef]
  77. Dean, R.G.; Dalrymple, R.A. Coastal Processes: With Engineering Applications; Cambridge University Press: Cambridge, MA, USA, 2002; p. 473. [Google Scholar]
  78. Netherlands-Water-Partnership, Indonesia—The Netherlands integrated approach of future water challenges. In Water for Food and Ecosystems: Using the Power of Nature; Braak, B.T. (Ed.) Netherlands Water Partnership (NWP): Den Haag, The Netherlands, 2016; pp. 26–30.
  79. Winterwerp, H.; van Wesenbeeck, B.; van Dalfsen, J.; Tonneijck, F.; Astra, A.; Verschure, S.; van Eijk, P. A Sustainable Solution for Massive Coastal Erosion in Central Java; Discussion Paper; Deltares: Delft, The Netherlands, 2014; pp. 1–45. [Google Scholar]
Figure 1. (a) The study area is located in the north of Java, Indonesia displayed using the RGB 542 of Landsat image. Some examples of coastal inundation impact: (b) daily flooded-houses; (c) an abandoned and flooded-paddy field; (d) coastal land embankment; and (e) permanent-inundated houses.
Figure 1. (a) The study area is located in the north of Java, Indonesia displayed using the RGB 542 of Landsat image. Some examples of coastal inundation impact: (b) daily flooded-houses; (c) an abandoned and flooded-paddy field; (d) coastal land embankment; and (e) permanent-inundated houses.
Remotesensing 09 00147 g001
Figure 2. Example of classification results using: FCM (ac); MLC classifier (df); and hardened classification (gi). (ai) are the detail presentations of yellow rectangle site in the insert map. Hard classification resulted from alternative methods are of limited use in identifying the transition zone between water and non-water, for example, see grid cells, e.g., B1 and B2.
Figure 2. Example of classification results using: FCM (ac); MLC classifier (df); and hardened classification (gi). (ai) are the detail presentations of yellow rectangle site in the insert map. Hard classification resulted from alternative methods are of limited use in identifying the transition zone between water and non-water, for example, see grid cells, e.g., B1 and B2.
Remotesensing 09 00147 g002
Figure 3. FCM results show the membership of: water (ac); and non-water (df). To derive shoreline position, we combined both membership images using fuzzy-crisp object model (gi). Blue pixels indicate core of water, orange pixels represent the core of non-water and shoreline is represented by light green pixels.
Figure 3. FCM results show the membership of: water (ac); and non-water (df). To derive shoreline position, we combined both membership images using fuzzy-crisp object model (gi). Blue pixels indicate core of water, orange pixels represent the core of non-water and shoreline is represented by light green pixels.
Remotesensing 09 00147 g003
Figure 4. The representation of fuzzy-crisp object model: (a) the core of water and non-water objects, and shorelines; (b) confusion index values considered for the quantification of classification uncertainty; and (c) shoreline image with fuzziness represented by confusion index. Detailed presentation of shorelines in red rectangle sites are displayed in (df).
Figure 4. The representation of fuzzy-crisp object model: (a) the core of water and non-water objects, and shorelines; (b) confusion index values considered for the quantification of classification uncertainty; and (c) shoreline image with fuzziness represented by confusion index. Detailed presentation of shorelines in red rectangle sites are displayed in (df).
Remotesensing 09 00147 g004
Figure 5. The magnitude of shoreline change during: 2013–2014 (a,b); and 2014–2015 (c,d). The magnitude values vary from high magnitude represented by dark blue pixels up to low magnitude represented by light blue pixels, whereas light yellow pixels show the no-change areas.
Figure 5. The magnitude of shoreline change during: 2013–2014 (a,b); and 2014–2015 (c,d). The magnitude values vary from high magnitude represented by dark blue pixels up to low magnitude represented by light blue pixels, whereas light yellow pixels show the no-change areas.
Remotesensing 09 00147 g005
Figure 6. The fuzziness of the shoreline is represented by change confusion values in the periods: 2013–2014 (a,b); and 2014–2015 (c,d). The change confusion values vary from high values represented by dark orange pixels up to low values represented by light orange pixels, whereas light yellow pixels show the no-change areas.
Figure 6. The fuzziness of the shoreline is represented by change confusion values in the periods: 2013–2014 (a,b); and 2014–2015 (c,d). The change confusion values vary from high values represented by dark orange pixels up to low values represented by light orange pixels, whereas light yellow pixels show the no-change areas.
Remotesensing 09 00147 g006
Figure 7. The representation of shoreline change direction: in the period 2013–2014 (a,b); and in the period 2014–2015 (c,d). Darker color pixels show a higher frequency of change to a certain direction. Shades of violet pixels depict a positive direction to water membership while shades of green pixels illustrate a negative direction to water membership. Figures in the second row show the magnitude of each change direction category in the period 2013–2014 (e,f); and in the period 2014–2015 (g,h). Darker color pixels represent a higher change magnitude while lighter color pixels show a lower change magnitude.
Figure 7. The representation of shoreline change direction: in the period 2013–2014 (a,b); and in the period 2014–2015 (c,d). Darker color pixels show a higher frequency of change to a certain direction. Shades of violet pixels depict a positive direction to water membership while shades of green pixels illustrate a negative direction to water membership. Figures in the second row show the magnitude of each change direction category in the period 2013–2014 (e,f); and in the period 2014–2015 (g,h). Darker color pixels represent a higher change magnitude while lighter color pixels show a lower change magnitude.
Remotesensing 09 00147 g007
Figure 8. Total intensity of confusion indices for each change direction category: in the period 2013–2014 (a,b); and in the period 2014–2015 (c,d). Shades of orange pixels represent change confusion values for the area with positive direction, and shades of grey pixels show change confusion values for the area with negative direction. The change confusion values for the unclear direction category are represented by shades of green, whereas no-change category is depicted by light yellow color.
Figure 8. Total intensity of confusion indices for each change direction category: in the period 2013–2014 (a,b); and in the period 2014–2015 (c,d). Shades of orange pixels represent change confusion values for the area with positive direction, and shades of grey pixels show change confusion values for the area with negative direction. The change confusion values for the unclear direction category are represented by shades of green, whereas no-change category is depicted by light yellow color.
Remotesensing 09 00147 g008
Figure 9. Change detection of shorelines using post classification comparison of MLC results; (ac) the change of shorelines in three consecutive dates in 2013; (d) shoreline changes from 2013 to 2014; and (e) shoreline changes from 2014 to 2015. Blue polygons show the changes of non-water to water and red polygons display the changes from water to non-water. No-change areas of water and non-water are represented by white and black polygons, respectively.
Figure 9. Change detection of shorelines using post classification comparison of MLC results; (ac) the change of shorelines in three consecutive dates in 2013; (d) shoreline changes from 2013 to 2014; and (e) shoreline changes from 2014 to 2015. Blue polygons show the changes of non-water to water and red polygons display the changes from water to non-water. No-change areas of water and non-water are represented by white and black polygons, respectively.
Remotesensing 09 00147 g009
Figure 10. Post classification comparison of shoreline as the results of hardened classification: (ac) the changes of shoreline in three consecutive dates in 2013; (d) shoreline changes from 2013 to 2014; and (e) from 2014 to 2015.
Figure 10. Post classification comparison of shoreline as the results of hardened classification: (ac) the changes of shoreline in three consecutive dates in 2013; (d) shoreline changes from 2013 to 2014; and (e) from 2014 to 2015.
Remotesensing 09 00147 g010
Figure 11. An example of comparison results between post classification comparison and CVA method. Both methods agree on change results of the area in yellow polygons that show a change from water to non-water (a,b) which equal to negative direction (c,d) with high change magnitude (e,f).
Figure 11. An example of comparison results between post classification comparison and CVA method. Both methods agree on change results of the area in yellow polygons that show a change from water to non-water (a,b) which equal to negative direction (c,d) with high change magnitude (e,f).
Remotesensing 09 00147 g011
Figure 12. Multi-year pattern of water membership changes showing a high change direction to positive direction (see dark violet pixels in black-dashed polygons in (a,b)) and high change magnitude (see dark blue pixels in black-dashed polygons in (c,d)). RGB 542 of Landsat images show a decrease of vegetation coverage from 2013 to 2015 (see white-dashed ellipses in (eh)).
Figure 12. Multi-year pattern of water membership changes showing a high change direction to positive direction (see dark violet pixels in black-dashed polygons in (a,b)) and high change magnitude (see dark blue pixels in black-dashed polygons in (c,d)). RGB 542 of Landsat images show a decrease of vegetation coverage from 2013 to 2015 (see white-dashed ellipses in (eh)).
Remotesensing 09 00147 g012
Figure 13. Water membership changes showing a continuous change to negative direction (see dark green pixels in black-dashed circles in (a,b)) with high change magnitude (see dark red pixels in black-dashed circle in (c,d)). RGB 542 of Landsat images show an increase of sediment and mangrove coverage from 2013 to 2015 (see white-dashed circle in (eh)).
Figure 13. Water membership changes showing a continuous change to negative direction (see dark green pixels in black-dashed circles in (a,b)) with high change magnitude (see dark red pixels in black-dashed circle in (c,d)). RGB 542 of Landsat images show an increase of sediment and mangrove coverage from 2013 to 2015 (see white-dashed circle in (eh)).
Remotesensing 09 00147 g013
Figure 14. The location shows an unclear direction (see pink pixels in black-dashed circle in (a)), while in the period 2014-2015 the location shows a low change to positive direction (light violet pixels in black-dashed circle in (b)). The change magnitude values were low in the period 2013-2014 (c), while in the period of 2014-2015 the values were high (see dark blue pixels in black-dashed circle in (d)). Images made available by Google Earth (eh) show the decrease of mangrove coverage.
Figure 14. The location shows an unclear direction (see pink pixels in black-dashed circle in (a)), while in the period 2014-2015 the location shows a low change to positive direction (light violet pixels in black-dashed circle in (b)). The change magnitude values were low in the period 2013-2014 (c), while in the period of 2014-2015 the values were high (see dark blue pixels in black-dashed circle in (d)). Images made available by Google Earth (eh) show the decrease of mangrove coverage.
Remotesensing 09 00147 g014
Figure 15. The embankment (shown by red arrows) for protecting the settlements (ac). (c) shows the river and settlements built on the river banks prone to high tides from both the sea and the river.
Figure 15. The embankment (shown by red arrows) for protecting the settlements (ac). (c) shows the river and settlements built on the river banks prone to high tides from both the sea and the river.
Remotesensing 09 00147 g015
Figure 16. Water membership changes showing a higher change to positive direction (see dark violet pixels in black rectangle sites in (a,b)) with low change magnitude values (see dark blue pixels in black rectangle sites in (c,d)).
Figure 16. Water membership changes showing a higher change to positive direction (see dark violet pixels in black rectangle sites in (a,b)) with low change magnitude values (see dark blue pixels in black rectangle sites in (c,d)).
Remotesensing 09 00147 g016
Figure 17. Water membership changes show a low change direction to positive direction in the period 2013-2014 (a); and to negative direction in the period 2014-2015 (b). The change magnitude values were low in both periods (c,d). This type of change was also observed as small patches of the coastal land (see black rectangle sites in (b,d)).
Figure 17. Water membership changes show a low change direction to positive direction in the period 2013-2014 (a); and to negative direction in the period 2014-2015 (b). The change magnitude values were low in both periods (c,d). This type of change was also observed as small patches of the coastal land (see black rectangle sites in (b,d)).
Remotesensing 09 00147 g017
Table 1. Landsat 8 OLI/TIRS images captured in the low tides supplemented by tide level and reference images used in the accuracy assessment purpose for each period.
Table 1. Landsat 8 OLI/TIRS images captured in the low tides supplemented by tide level and reference images used in the accuracy assessment purpose for each period.
Acquisition DateAstronomical Tide Level (m)Reference DataAcquisition DateAstronomical Tide Level (m)Reference Data
23 May 2013−0.1Pleiades (27 February 2013)29 May 2015+0.04Sentinel 2 (26 December 2015)
12 September 2013−0.118 September 2015−0.1
14 October 2013−0.320 October 2015−0.3
1 December 2013−0.321 November 2015−0.3
10 May 2014−0.01Spot 6 (5 October 2014)
15 September 2014−0.2
1 October 2014−0.2
18 November 2014−0.3
Table 2. The spectral band information of Landsat 8 OLI/TIRS used in image classifications, and Pleiades, SPOT 6 and Sentinel 2 used as reference images.
Table 2. The spectral band information of Landsat 8 OLI/TIRS used in image classifications, and Pleiades, SPOT 6 and Sentinel 2 used as reference images.
SatelliteBandsWavelength (μm)SatelliteBandsWavelength (μm)
Landsat 8 OLI/TIRSCoastal and Aerosol0.43–0.45SPOT 6Blue0.45–0.52
Blue0.45–0.51Green0.53–0.59
Green0.53–0.59Red0.625–0.695
Red0.64–0.67NIR0.76–0.89
NIR0.85–0.88
SWIR 11.57–1.65
SWIR 22.11–2.29
PleiadesBlue0.43–0.55Sentinel 2Blue0.49
Green0.50–0.62Green0.56
Red0.59–0.71Red0.665
NIR0.74–0.94SWIR0.842
Table 3. The procedure to estimate change directions of shoreline. It quantifies the variation of water membership in each pixel and shows how frequent the changes have occurred.
Table 3. The procedure to estimate change directions of shoreline. It quantifies the variation of water membership in each pixel and shows how frequent the changes have occurred.
CCCVTCVChg.DirCCCVTCVChg.Dir
CV1CV2CV3CV4CV1CV2CV3CV4
100000No-change410+10−10Unclear direction
2+1+1+1+1+4Positive direction42+100−10Unclear direction
3+10+1+1+3Positive direction4300−1+10Unclear direction
40+1+1+1+3Positive direction440−10+10Unclear direction
5+1+10+1+3Positive direction45−1+1+1−10Unclear direction
6+1+1+10+3Positive direction46+1−1+1−10Unclear direction
7+1−1+1+1+2Positive direction47+1−1−1+10Unclear direction
8−1+1+1+1+2Positive direction48−1+1−1+10Unclear direction
9+1+1−1+1+2Positive direction49+1+1−1−10Unclear direction
10+1+1+1−1+2Positive direction50−1−1+1+10Unclear direction
110+1+10+2Positive direction510−100−1Negative direction
12+10+10+2Positive direction52−1000−1Negative direction
13+100+1+2Positive direction5300−10−1Negative direction
140+10+1+2Positive direction54000−1−1Negative direction
15+1+100+2Positive direction550+1−1−1−1Negative direction
1600+1+1+2Positive direction560−1+1−1−1Negative direction
17+1000+1Positive direction57−1+10−1−1Negative direction
180+100+1Positive direction58−10+1−1−1Negative direction
19000+1+1Positive direction59+10−1−1−1Negative direction
2000+10+1Positive direction60+1−10−1−1Negative direction
210+1−1+1+1Positive direction610−1−1+1−1Negative direction
220−1+1+1+1Positive direction62−1−1+10−1Negative direction
23−1+10+1+1Positive direction63+1−1−10−1Negative direction
24−10+1+1+1Positive direction64−1+1−10−1Negative direction
25+10−1+1+1Positive direction65−1−10+1−1Negative direction
26+1−10+1+1Positive direction66−10−1+1−1Negative direction
27+10+1−1+1Positive direction6700−1−1−2Negative direction
28+1−1+10+1Positive direction68−10−10−2Negative direction
29+1+10−1+1Positive direction690−1−10−2Negative direction
30+1+1−10+1Positive direction70−1−100−2Negative direction
310+1+1−1+1Positive direction710−10−1−2Negative direction
32−1+1+10+1Positive direction72−100−1−2Negative direction
330+1−100Unclear direction73+1−1−1−1−2Negative direction
340−1+100Unclear direction74−1+1−1−1−2Negative direction
35−1+1000Unclear direction75−1−1−1+1−2Negative direction
36−10+100Unclear direction76−1−1+1−1−2Negative direction
37+10−100Unclear direction770−1−1−1−3Negative direction
38+1−1000Unclear direction78−10−1−1−3Negative direction
3900+1−10Unclear direction79−1−1−10−3Negative direction
40−100+10Unclear direction80−1−10−1−3Negative direction
81−1−1−1−1−4Negative direction
Notes: C V : Change vector (based on Equation (7); C C : Change combinations number (based on Equation (9)); T C V : Total change vector (based on Equation (10)); and C h g . D i r : Change direction.
Table 4. Summary of the overall classification accuracy using FCM, MLC and hardened classification.
Table 4. Summary of the overall classification accuracy using FCM, MLC and hardened classification.
Classified ImagesOverall Accuracy
FCMMLCHardened Classification
23 May 20130.870.720.86
12 September 20130.850.760.85
14 October 20130.860.730.86
1 December 20130.860.730.84
10 May 20140.890.780.87
15 September 20140.900.760.88
1 October 20140.900.780.88
18 November 20140.910.790.90
29 May 20150.840.750.84
18 September 20150.880.790.87
20 October 20150.890.790.86
21 November 20150.890.800.88
Table 5. Change area (in ha) for each change category in the period of 2013–2014 and 2014–2015.
Table 5. Change area (in ha) for each change category in the period of 2013–2014 and 2014–2015.
Change Category2013–20142014–2015
Positive direction18281120
Negative direction9201635
Unclear direction616528
No-change13191403

Share and Cite

MDPI and ACS Style

Dewi, R.S.; Bijker, W.; Stein, A. Change Vector Analysis to Monitor the Changes in Fuzzy Shorelines. Remote Sens. 2017, 9, 147. https://doi.org/10.3390/rs9020147

AMA Style

Dewi RS, Bijker W, Stein A. Change Vector Analysis to Monitor the Changes in Fuzzy Shorelines. Remote Sensing. 2017; 9(2):147. https://doi.org/10.3390/rs9020147

Chicago/Turabian Style

Dewi, Ratna Sari, Wietske Bijker, and Alfred Stein. 2017. "Change Vector Analysis to Monitor the Changes in Fuzzy Shorelines" Remote Sensing 9, no. 2: 147. https://doi.org/10.3390/rs9020147

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop