Analysis of Coastline Extraction from Landsat-8 OLI Imagery

: Coastline extraction is a fundamental work for coastal resource management, coastal environmental protection and coastal sustainable development. Due to the free access and long-term record, Landsat series images have the potential to be used for coastline extraction. However, dynamic features of different types of coastlines (e.g., rocky, sandy, artiﬁcial), caused by sea level ﬂuctuation from tidal, storm and reclamation, make it difﬁcult to be accurately extracted with coarse spatial resolution, e.g., 30 m, of Landsat images. To access this problem, we analyze the performance of coastline extraction by integrating downscaling, pansharpening and water index approaches in increasing the accuracy of coastline extraction from the latest Landsat-8 Operational Land Imager (OLI) imagery. In order to prove the availability of the proposed method, we designed three strategies: (1) Strategy 1 uses the traditional water index method to extract coastline directly from original 30 m Landsat-8 OLI multispectral (MS) image; (2) Strategy 2 extracts coastlines from 15 m fused MS images generated by integrating 15 m panchromatic (PAN) band and 30 m MS image with ten pansharpening algorithms; (3) Strategy 3 ﬁrst downscales the PAN band to a ﬁner spatial resolution (e.g., 7.5 m) band, and then extracts coastlines from pansharpened MS images generated by integrating downscaled spatial resolution PAN band and 30 m MS image with ten pansharpening algorithms. Using the coastline extracted from ZiYuan-3 (ZY-3) 5.8 m MS image as reference, accuracies of coastlines extracted from MS images in three strategies were validated visually and quantitatively. The results show that, compared with coastline extracted directly from 30 m Landsat-8 MS image (strategy 1), strategy 3 achieves the best accuracies with optimal mean net shoreline movement (MNSM) of − 2.54 m and optimal mean absolute difference (MAD) of 11.26 m, followed by coastlines extracted in strategy 2 with optimal MNSM of − 4.23 m and optimal MAD of 13.54 m. Further comparisons with single-band thresholding (Band 6), AWEI, and ISODATA also conﬁrmed the superiority of strategy 3. For the various used pansharpening algorithms, ﬁve multiresolution analysis MRA-based pansharpening algorithms are more efﬁcient than the component substitution CS-based pansharpening algorithms for coastline extraction from Landsat-8 OLI imagery. Among the ﬁve MRA-based fusion methods, the coastlines extracted from the fused images generated by Indusion, additive à trous wavelet transform (ATWT) and additive wavelet luminance proportional (AWLP) produced the most accurate and visually realistic representation. Therefore, pansharpening approaches can improve the accuracy of coastline extraction from Landsat-8 OLI imagery, and downscaling the PAN band to ﬁner spatial resolution is able to further improve the coastline extraction accuracy, especially in crenulated coasts.


Introduction
Coastline, the boundary of land and sea, is vulnerable to the natural processes, such as coastal erosion-accretion and sea level changes, and human activities. Coastline mapping and change detection are, therefore, becoming a fundamental work for coastal resource management, environmental protection of coastal zone and coastal sustainable development [1][2][3]. Traditional ground surveying is known as an important method to complete coastline maps [4]. However, transferring the information from grounding surveys and photographs to coastline maps is costly and time-consuming [4,5]. With the characteristics of macroscopic, real-time, dynamic and cost-effective, and a large range of spatial and temporal scales, remote sensing data is used widely for coastline mapping and change detection [4,6,7].
In general, the coastline extracted from remote sensing imagery with finer spatial resolution is more likely to achieve higher accuracy [8], and an increasing number of high spatial resolution remote sensing images were used for effective coastline extraction. For example, Bruno et al. [9] exploits SAR image to semi-automatic coastline extraction on a very popular beach in Apulia Region (Italy); Liu et al. [10] presents a segmentation-based image processing method to extracted coastline from LiDAR data. Di et al. [11] employed semi-automatic method based on mean shift segmentation algorithm and a local refinement process to extract coastlines from IKONOS image; Puissant et al. [12] extracted coastlines from Quickbird multispectral (MS) images by using mathematical morphology; Lira et al. [3] used digital aerial photographs and digital orthophotomaps to present the first systematic, national-scale and consistent long-term coastline evolution data of Portuguese mainland low-lying sandy coasts; Sekovski et al. [13] applied a WorldView-2 MS image to delineate coastline by integrating supervised and unsupervised classification methods. Although the high spatial resolution remote sensing images play a significant role in obtaining satisfactory result regarding coastline extraction, their applications for regional-scale mapping in certain areas would be cost-prohibitive [8,13]; moreover, the infrequent revisit coverage of high spatial resolution remote sensing images would make it a challenge to perform coastline extraction in large areas. Valentini [14] proposed a tool for automatic shoreline detection based on remote sensing video systems, which have both high temporal and spatial resolutions. With respect to traditional procedures, video systems are a low-coast method for long-term and continuous coastline monitoring [14], while the storage of continuous video streams for long periods is large.
An alternative remote sensing image used for coastline extraction is the popular Landsat series imagery, because of its free access, large revisit coverage and long-term data record [15]. In real applications, Landsat series images acquired from sensors of Landsat 1-3 Multi-spectral Scanner (MSS), Landsat 4-5 Thematic Mapper (TM), Landsat-7 Enhanced Thematic Mapper Plus (ETM+) and Landsat-8 Operational Land Imager (OLI) are useful data for extracting coastline [16]. Lira and Taborda [17] summarized the advances of coastline extraction from available free-access satellite imagery, and showed that the Landsat imagery offers the best compromise between potential and availability to understand coastal dynamic features. Chen et al. [18] used Landsat MSS and TM MS images to monitor coastline changes in Lingding Bay, and the coastline extends on a coastal stretch of approximately 6 km; Ekercin [19] collected Landsat MSS, TM and ETM+ images between 1975 and 2001 to perform coastline change analysis in a part of the Aegean Sea; Wang et al. [20] used Landsat MSS, TM, ETM+ and OLI images for coastline detection along Ningbo City from 1976 to 2015; Pardo-Pascual et al. [21] evaluated how the storms impact the sandy beaches and the beach recovery process by describing the coastline positions obtained from Landsat TM and ETM+. Coastlines with special characteristics, such as deltas, exhibit less change in position, and they could be monitored by low resolution of Landsat images, while the coastline of sandy areas with small beaches is difficult to be monitored by low resolution images.
The biggest challenge for applying Landsat series images on coastline extraction and monitoring is the limited spatial resolution (30 m) of MS image [22]. Pardo-Pascual et al. [23] proposed an automatic method to extract shorelines from Landsat TM and ETM+ with subpixel precision, where the RMSE obtained in shoreline location ranges from 4.69 to 5.47 m. As the nature beaches are always changing over time, Almonacid-Caballer et al. [24] applied the annual mean shoreline position extracted from Landsat images to largely decline the short-term variability of shoreline, and the extracted shorelines were biased from the seaward by around 4 to 5 m. Foody et al. [25] evaluated a soft classification method in mapping coastline from a degraded simulated Landsat imagery, and the result presented coastline obtained at subpixel scale and the RMSE was 2.25 m. While, these methods are always complex and difficult to be achieved in real application. In the processing, only spectral information of Landsat imagery was utilized and the improvement of accuracy is limited. Notably, besides the 30 m MS image, landsat-7 ETM+ and Landsat-8 OLI sensors simultaneously acquire the 15 m spatial resolution panchromatic (PAN) band. Motivated by this, it is of great interest to improve the coastline extraction accuracy from Landsat ETM+ and OLI images by adding spatial information of PAN band into MS bands. In fact, researchers have developed many algorithms to generate a fine spatial resolution MS image by combining PAN and MS bands, and this concept is known as pansharpening [26].
Pansharpening is a data fusion technique adding high spatial information of a PAN image into a coarse spatial resolution MS image, resulting in a sharpened MS image that has the same spatial resolution of the PAN image. In the process, pansharpening aims to increase the spatial resolution while simultaneously preserving the spectral feature of primary MS image [27,28], and it has been widely used for integrating Landsat ETM+ MS and PAN band to improve the land use/cover mapping, precision agriculture, pollution monitoring, forest resource inventory and urban areas mapping [29][30][31]. However, for Landsat ETM+ satellite sensor, an instrument malfunction occurred on 31 May 2003, with the result that all Landsat-7 scenes acquired since 14 July 2003 have been collected in "SLC-off" mode. The newly launched Landsat-8 OLI image provides seven MS bands with 30 m spatial resolution and a PAN band with 15 m spatial resolution [32]. This makes it readily to obtain a 15 m spatial resolution fused MS image by pansharpening the Landsat-8 30 m MS image and 15 m PAN band. It is noteworthy that the spatial resolution of the sharpened MS image, however, may still not be sufficient to provide subtle variability of coastline retreat and advance information. An alternative solution of this problem is to produce spatially sharpened image with finer spatial resolution beyond both of the Landsat-8 OLI MS and PAN bands [33,34]. For this case, the PAN band is first downscaled to a finer spatial resolution band, and then integrated into the original MS bands by using pansharpening algorithm to produce a finer spatial resolution sharpened MS image.
Once sharpened MS image has been generated, coastline extraction methods should be further applied to extract coastline from the MS image [35]. Previous studies [20,[36][37][38][39][40] have proved that water index method has advantages of universality, user friendliness and low computation cost in coastline extraction, and the common used water index is the normalized difference water index (NDWI) proposed by McFeeters [39]. In this paper, downscaling, pansharpening and water index approaches were used to extract coastlines from Landsat-8 OLI imagery. To evaluate how pansharpening algorithms impact on coastline extraction from fused Landsat-8 MS image, ten popular pansharpening algorithms, including additive à trous wavelet transform (ATWT), additive wavelet luminance proportional (AWLP), Brovey transform (Brovey), Gram Schmidt (GS), high pass filtering (HPF), intensity-hue-saturation (IHS), additive injection model (Indusion), principal component analysis (PCA), partial replacement adaptive component substitution (PRACS) and smoothing-filter-based intensity modulation (SFIM), are used to produce sharpened MS images. Moreover, three strategies are devised to extract coastlines from Landsat-8 OLI imagery. Strategy 1 extracts coastline directly from original Landsat-8 OLI 30 m MS image, and it is a traditional widely used method [41]. Strategy 2 extracts coastlines from pansharpened 15 m MS images generated by integrating Landsat-8 OLI 15 m PAN band and 30 m MS image with different pansharpening algorithms. Strategy 3 first downscales the Landsat-8 OLI PAN band to a finer spatial resolution (e.g., 7.5 m) band, and then extracts coastlines from the fused MS images generated by integrating downscaled PAN band and original MS image with different pansharpening algorithms.

Study Area and Datasets
A portion of Ningbo coast was selected as the study area. As shown in Figure 1, Ningbo locates between 29 • 20 -29 • 32 N and 121 • 50 -122 • 14 E, and is next to the East China Sea. In China, Ningbo is an important port of international economic and cultural exchanges [42]. The increasing connection with the world has pushed the urban development of Ningbo, resulting in rapid coastline position changes during the past decades [20]. The types of Ningbo coast mainly contain bedrock coast, artificial coast and flat sandy coast. Subsites in Figure 1a presents two artificial coasts (subsites 1 and 4), two bedrock coasts (subsites 2 and 5) and two flat sandy coasts (subsites 3 and 6). is an important port of international economic and cultural exchanges [42]. The increasing connection with the world has pushed the urban development of Ningbo, resulting in rapid coastline position changes during the past decades [20]. The types of Ningbo coast mainly contain bedrock coast, artificial coast and flat sandy coast. Subsites in Figure 1a presents two artificial coasts (subsites 1 and 4), two bedrock coasts (subsites 2 and 5) and two flat sandy coasts (subsites 3 and 6). In this study, one Landsat-8 OLI MS image ( Figure 1a) and one PAN image, covering an area of 21 km × 21 km, were used for coastline extraction. The original Landsat-8 OLI imagery located in path 118 and row 039 was acquired at 2:24 a.m. (UTC) on 12 March 2015 and downloaded as the land surface reflectance product in the GeoTIFF format from United States Geological Survey (USGS) earth explore website (http://earthexplorer.usgs.gov/). Landsat-8 OLI imagery is with the Word Geodetic System (WGS84) datum and projected using the Universal Transverse Mercator (UTM) System. The detailed parameters of the Landsat-8 OLI imagery are presented in Table 1. Notably, the wavelength of PAN image only contains the spectral range from band 2 to band 4 of MS image, covering visible wavelengths, while the Landsat-7 ETM+ PAN band covers both visible and NIR wavelengths.
Moreover, one ZiYuan-3 surveying satellite (ZY-3) MS image was used to extract the reference coastline covering the study area as shown in Figure 1. ZY-3 was launched on 9 January 2012, and it is China's first civilian high-resolution stereo mapping satellite [43]. The MS image of ZY-3 has a spatial resolution of 5.8 m. The fine spatial resolution ZY-3 MS image was acquired at 2:45 a.m. (UTC)  Table 1. Notably, the wavelength of PAN image only contains the spectral range from band 2 to band 4 of MS image, covering visible wavelengths, while the Landsat-7 ETM+ PAN band covers both visible and NIR wavelengths.
Moreover, one ZiYuan-3 surveying satellite (ZY-3) MS image was used to extract the reference coastline covering the study area as shown in Figure 1. ZY-3 was launched on 9 January 2012, and it is China's first civilian high-resolution stereo mapping satellite [43]. The MS image of ZY-3 has a spatial resolution of 5.8 m. The fine spatial resolution ZY-3 MS image was acquired at 2:45 a.m. (UTC) on 12 March 2015, which is really close to the acquisition time of Landsat-8 OLI imagery shown in Figure 1. In this case, tide difference can, therefore, be ignored for similar tide levels between Landsat-8 OLI and ZY-3 MS images [20]. Further information about the specifications of ZY-3 MS image is listed in Table 1.
For the acquired Landsat-8 OLI and ZY-3 MS images covering the same geographical region, ZY-3 MS image was first geometrically reprojected into the UTM, WGS84 coordinate system of Landsat-8 OLI image. Image-to-image registration method in ENVI 5.1 was then used to make pixels in the two images precisely coincide to the same points on the ground [27,44]. In the process, Landsat-8 PAN image was chosen as a master and ZY-3 was to be registered. The registration error was less than one pixel of ZY-3.

Methodology
In this study, three strategies were devised for the coastline extraction from Landsat-8 OLI MS and PAN bands. For strategy 1, the original Landsat-8 OLI MS bands were used directly to calculate the NDWI image, and then the dynamic segmentation method of Otsu algorithm [45] was used to separate the NDWI image into water and non-water features. The coastline was extracted finally by converting the water/non-water binary image into vector format with ArcGIS software. For strategy 2, instead of directly using Landsat-8 OLI MS image, pansharpening algorithms were first used to produce 15 m sharpened MS image by integrating Landsat-8 OLI MS and PAN bands, and then the fused MS image was used to calculate 15 m NDWI image and final coastline position map. For strategy 3, the 15 m PAN band was first downscaled to 7.5 m, and then pansharpening algorithms were used to produce 7.5 m sharpened MS image, which then was used to calculate 7.5 m NDWI image and final coastline position map. Downscaling method involved in the whole process of different strategies was performed via a polynomial kernel with 23 coefficients [46]. Ten pansharpening methods were collected to obtain the fine spatial resolution MS images in strategies 2 and 3. With respect to the accuracy assessment, the coastline derived from ZY-3 image with 5.8 m was regarded as reference. Specially, ZY-3 MS image was first used to calculate the 5.8 m NDWI image, and then Otsu segmentation method was applied to get the water/non-water binary map from 5.8 m NDWI, and the final reference coastline was extracted from the water/non-water binary map. The detailed methodology of the whole study is presented in Figure 2.

Water Index
NDWI proposed by McFeeters (1996) is a non-linear conversion of the green and near-infrared band of MS image, and it is one of the most widely used water indices for coastline extraction. For water feature, it reflects the radiation in green band and absorbs in NIR band [35]. As for terrestrial vegetation and soil features, they have higher reflectance in NIR band than in green band. According to the different reflectivity between water and non-water, NDWI was defined as where Green  and NIR  are the reflectance value of green and NIR band of MS image [39]. The results of the index can range from −1 to 1. Water surfaces tend to have positive values and non-water surfaces have negative values. In the case, MS image which has one green band and one NIR band can be used to calculate NDWI.

Water/Non-Water Binary Map Produced by Using Otsu Method
Once the NDWI image was produced, a threshold value is necessary to divide the image into two classes, land and water body. The problem is that the NDWI threshold value is not constant, but varies depending on the land cover components [47]. The Otsu method, proposed by Nobuyuki Otsu [45], selects an optimal threshold automatically from a gray image and the threshold maximize the

Water Index
NDWI proposed by McFeeters (1996) is a non-linear conversion of the green and near-infrared band of MS image, and it is one of the most widely used water indices for coastline extraction. For water feature, it reflects the radiation in green band and absorbs in NIR band [35]. As for terrestrial vegetation and soil features, they have higher reflectance in NIR band than in green band. According to the different reflectivity between water and non-water, NDWI was defined as where ρ Green and ρ NIR are the reflectance value of green and NIR band of MS image [39]. The results of the index can range from −1 to 1. Water surfaces tend to have positive values and non-water surfaces have negative values. In the case, MS image which has one green band and one NIR band can be used to calculate NDWI.

Water/Non-Water Binary Map Produced by Using Otsu Method
Once the NDWI image was produced, a threshold value is necessary to divide the image into two classes, land and water body. The problem is that the NDWI threshold value is not constant, but varies depending on the land cover components [47]. The Otsu method, proposed by Nobuyuki Otsu [45], selects an optimal threshold automatically from a gray image and the threshold maximize the inter-class discrepancy. Pervious researches have proved that Otsu method can successfully be applied to map water bodies from Landsat imagery [37,38].
Assuming t represents the threshold value, ranging from a to b (−1 ≤ a ≤ b ≤ 1), the optimal threshold value t * can be obtained by the following algorithm: where σ is the inter-class variance of non-water and water pixels. P nw is the possibility of a pixel belonging to the non-water and P w means the possibility of a pixel that belongs to water. M nw and M w respectively represent the mean values of all non-water pixels and all water pixels, while M is the mean value of the entire image. As a result, the pixels with values equal or higher the optimal threshold value was divided into water and scored 1. On the contrary, the pixel was classified as non-water and scored 0. The boundary of water and non-water is decided as coastline indicator.

Pansharpening Algorithms
Since the 1980s, studies on pansharpening with multiple source images for better understanding land surface feature has been conducted [48,49]. Many methods have been proposed to realize the fusion process and they can be classified into two categories: component substitution CS-based and multiresolution analysis MRA-based techniques [50]. In order to test the capability of pansharpening for coastline extraction from Landsat-8 OLI imagery, ten widely used standard pansharpening algorithms encompassing different concept were chosen. Specially, five CS-based pansharpening algorithms, including IHS, PCA, GS, Brovey, PRACS, and five MRA-based pansharpening algorithms, including HPF, SFIM, Indusion, ATWT, AWLP, are chosen in the study. The brief introduction of these pansharpening algorithms is described in the following subsections.

Intensity-Hue-Saturation (IHS)
The IHS transformation is a classic image fusion method and it is based on the component substitution to enhance the spatial feature of coarse spatial resolution MS image. Intensity refers to the total brightness of a color image, while hue specifies the dominant or average wavelength of light contributing to a color image and saturation indicates the purity of a color relative to gray [51]. Arbitrary bands of MS image assigned to RGB channel can be transformed to IHS space and replace the first band by PAN image with fine spatial resolution. Then, the substitutive IHS was transformed backward the original RGB space [52]. In the process, the IHS transformation method injects the spatial information of PAN image into MS image and preserves the spectral information of original MS image.

Principal Component Analysis (PCA)
The PCA method is widely used in signal processing, statistics and many other areas, and is also a classic image fusion method [27]. PCA transformation generates a set of uncorrelated components from original MS image. The first component contains the largest information of the color image, and the following ones compose the residual information decreasing successively. Theoretically, the largest information specifies the common feature in all bands, i.e., the spatial feature, while the unique spectral information was mapped to other components [27,53]. Similar to IHS, the first component was substituted by the fine resolution PAN image. Finally, an inverse PCA transform takes the fused image back into the original multispectral feature space [54].

Gram Schmidt (GS)
The widespread GS method was invented by Laben and Brover in 1998 and patented by Eastman Kodak [55]. The GS orthogonalizes the bands of MS image to remove redundant information from imputed coarse spatial resolution MS image [27]. By spatially averaging the multispectral MS bands, a coarse spatial resolution synthetic panchromatic band was generated and used as the first band of GS transform. Then the high spatial resolution PAN image replaces the first GS transform component. Finally, an inverse GS transformation is applied to produce the fine spatial resolution MS image [49].

Brovey Transform (Brovey)
The Brovey transform is a generally used method based on the chromaticity transform [56]. This algorithm considers that the image consists of two parts, hue and brightness [57]. First, each MS band that is chosen to compose a RGB space is divided by the sum of the three bands. The aim of the ratio is to preserve the spectral feature of coarse spatial resolution MS image. Second, multiplying the ratio by the high spatial resolution PAN band obtains the overall brightness of high spatial resolution image. Finally, a new image with the spatial resolution of inputted PAN image is produced. In the process, the Brovey transform technique provides excellent visually contrast of the fusion image.

Partial Replacement Adaptive Component Substitution (PRACS)
The PRACS method based on component substitute was proposed by Choi et al. (2011) [58] to improve the spatial and spectral quality of original MS imagery. To obtain the fine spatial resolution MS image, the method is organized into two steps. The first step is the construction of fine and coarse spatial resolution component images using PAN and MS images by partial replacement. Considering the unfixed relationship varying with different object, area, and circumstance, a linear regression algorithm is performed between PAN image and each MS bands. The second step is find an optimal CS fusion model to minimize the local and global spectral distortion between the PAN image and each MS band while preserving the spatial information within original PAN image [58].

High Pass Filtering (HPF)
The HPF method proposed by Schwengerdt (1980) is regarded as a predecessor of image fusion methods that are based on the concept of adding the spatial information extracted from fine spatial resolution image into the original coarse spatial resolution MS image [27,59]. For HPF, the spatial feature consists with the high frequency information, which can be obtained by filtering the PAN image with a high-pass filter or subtracting the low-pass filter image from the original PAN image. There are many types of filters that can be used in the process, such as Box and Gaussian filters. In this paper, the Box filter was used to realize high spatial resolution information extraction from PAN image, and the final sharpened MS image is produced by injecting the extracted high spatial information into all bands of the initial MS image.

Smoothing-Filter-Based Intensity Modulation (SFIM)
The SFIM is a simple spectral preservation fusion method and it is based on a simplified solar radiation and land surface reflection model [60]. The method was proposed by   [60], and simplified as: where I SFIM means the fused image by SFIM method, I coarse is the coarse spatial resolution MS image and I fine is the fine spatial resolution PAN image. I mean is a smoothed image using spatially averaging filter over the fine spatial resolution PAN image. Notably, I coarse and I mean must have the same pixel size with I fine . In Equation (3), the ratio between I fine and I mean removes the spectral and topographical information of the fine spatial resolution PAN image and maintains the fine spatial resolution textures.
Thus, the pansharpened MS image not only has fine spatial resolution textures of PAN image, but also highly preserves the spectral information of the coarse spatial resolution MS image [60][61][62].

Additive Injection Model (Indusion)
Indusion is proposed by Khan et al. [63] for dealing with the spectral distortion caused by inaccurately registration between PAN and downscaled MS image. The method uses Induction scaling method instead of the bicubic interpolation within SFIM fusion method to produce a better co-registered downscaled MS image. Furthermore, a pair of upscaling and downscaling filters is applied to extract high frequency information from PAN image, and then the extracted spatial information was injected into original MS image, finally producing a fine spatial resolution MS image.
3.3.9. Additive À Trous Wavelet Transform (ATWT) The ATWT method was developed by Núñez et al. [64] for image fusion based on multiresolution wavelet decomposition. In the process of wavelet transform, "à trous" algorithm is applied to obtain a shift-invariant discrete wavelet decomposition for input image. The PAN image is decomposed to several wave planes, and then the wavelet planes including the fine spatial information of PAN image is directly added into the original MS image [64,65]. Each of the fused fine spatial resolution MS bands is finally obtained by the inverse wavelet transform.

Additive Wavelet Luminance Proportional (AWLP)
Similarly to ATWT, AWLP is also a multiresolution wavelet-based image fusion method. To obtain a shift-invariant discrete wavelet decomposition for input image, the same discrete algorithm of the wavelet transform is used. Instead of directly injecting high spatial information into original coarse spatial resolution MS image, AWLP first transforms the RGB components of MS image to IHS space, and then adds the spatial information of PAN image into the first component. Finally, the new IHS space is inversed back into RGB components to produce the fused fine spatial resolution MS image [64,66].

Downscaling Panchromatic (PAN) Band with Interpolation
Downscaling is a potential method in remote sensing that can produce finer spatial resolution imagery than input [67]. In this research, downscaling approach is used to estimate the finer spatial resolution PAN band P 15/s from the observed PAN band P 15 with spatial scale s, and the process is formulated as where ↓ s indicates a downscaling operation. Generally, it can be performed readily by a common interpolation method, such as bilinear and spline interpolations [67]. Bilinear, as an easy-to-operate method, was used in the paper.

Validation of Coastline Extraction
In order to have a quantitative assessment of the capacity of coastline extraction from Landsat-8 OLI imagery with different strategies, six indices, including mean absolute difference (MAD), maximum AD, minimum AD, mean net shoreline movement (MNSM), maximum positive NSM and maximum negative NSM were applied by using the coastline extracted from ZY-3 MS image as reference.
The Digital Shoreline Analysis System (DSAS) [68], a free extent for ArcGIS software, was applied in this study to determine the degree of position differences between extracted coastlines [69], such as NSM and AD. NSM represents a distance between the reference and the extracted coastline at each transect. The positive value of NSM means that the extracted coastline is on the sea side of the reference coastline; otherwise, the extracted coastline is on the land side. AD means the absolute distance between the extracted coastline and the reference coastline at each transect. In the test, AD value equals to the statistic of end point rate provided by DSAS at each transect. In the study, there are total of 866 transects that were created at a horizontal interval of 50 m to calculate the NSM and AD values.

Pansharpened Multispectral (MS) Images
Ten pansharpening algorithms, including five CS-based methods of PCA, GS, IHS, Brovey, PRACS, and five MRA-based method of HPF, SFIM, Indusion, ATWT and AWLP, were used to produce sharpened MS images in strategies 2 and 3.  Figure 3a-c. An artificial coastal subsite along the coastal is magnified 20 times and presents in lower-right corner of each fused MS images. In order to have a quantitative validation of the fused MS images generated by different pansharpening algorithms for strategies 2 and 3, the quality with no reference (QNR) index combined with the D s and D λ were applied in this subsection [70]. QNR, D s and D λ are widely used to examine the spectral and spatial information of the sharpened MS with no reference. Specially, D s is used to validate the spectral information inheriting from the original coarse spatial resolution MS image, and lower D s value means better maintaining of the original MS spectral information. D λ is used to validate the spatial information exploited from the fine spatial resolution PAN band, and lower D λ value indicates better restoring of the spatial information of PAN band. QNR is a combination of D s and D λ , and higher value of QNR means better quality of the sharpened MS image.
In the section, a brief visual and quantitative analysis, including spectral preserving and spatial improvement, is presented. In Figure 3, compared with the original 30 m MS image shown in Figure 3a, the spatial detailed information in all fused images are improved obviously, practically in the results of strategy 3, which shows the promising advance of pansharpening algorithms. Compared with the results of most CS-based pansharpening algorithms, the results of MRA-based pansharpening algorithms of HPF, SFIM, Indusion, ATWT, AWLP can better preserve the spectral information within original 30 m MS image, and they achieve obvious lower D s values and higher QNR values for both strategies 2 and 3 as shown in Table 2. From Figure 3 and Table 2, the CS-based method of PRACS achieves the lowest D s and D λ values, and highest QNR values, and PCA and GS shows obviously spectral distortion. The colors of Sea are moorea blue in the original MS image, while that of PCA and GS fused MS images are cheery cola and medium azul, respectively. The colors of vegetation features in the fused PCA and GS images change from bright red to purple and black red, respectively. This is because PCA and GS algorithms have poor ability to preserve the spectral information of original MS image in the resultant sharpened MS image. Compared with strategy 2, most of the sharpened MS images in strategy 3 achieve an improvement of the QNR, D s and D λ values, and this indicates the superiority of strategy 3 in preserving spectral information of original MS image and exploiting spatial information of PAN band.

Coastline Extraction
Once the fused MS images in strategies 2 and 3 have been produced, they can thus be used to generate corresponding NDWI images, which are shown in Figure 4c-v. Moreover, by applying the Otsu segmentation method, the water and non-water binary maps shown in Figure 5 were generated from the NDWI images. It can be found that more spatial detail information is exploited within the NDWI images of strategies 2 and 3 by comparing with that of the original 30 m NDWI image (Figure 4a), and the corresponding water body maps of strategies 2 and 3 also have more spatially smooth boundaries than that of original 30 m water body map shown in Figure 5a. Compared with the results of strategy 2, the water body maps of strategy 3 contain more spatial details and have spatially smooth boundaries, which shows the efficiency of strategy 3 for coastline extraction. For both strategies 2 and 3, most of the water body maps of MRA-based pansharpening algorithms contain more spatial detail information than that of CS-based pansharpening algorithms. Meanwhile, in Figures 4 and 5, the worst results are associated with the PCA and GS algorithms, and they are almost impossible to distinguish water and non-water features. This is because the fused MS images generated by PCA and GS have a serious spectral distortion of band 3 and band 5, making the calculated NDWI images have a confusion of water and non-water features, and this phenomenon can also be observed in Section 4.1. Although the exploiting of fine spatial information within PAN image is really important, the spectral preserving may be more important for the fused MS image to calculate NDWI image. In the following section, PCA and GS algorithms were, therefore, not applied to produce coastlines. Based on the water body maps shown in Figure 5, twenty vector coastlines were delineated with ArcGIS software. These coastlines were then used to calculate six accuracy indices mentioned in Section 3.5 by using the coastline extracted from the ZY-3 water body map as reference.

Comparising Extracted Coastlines of Different Strategies
To understand the capacity of coastline extraction from Landsat-8 OLI imagery with different strategies, three subsites containing artificial coast, bedrock coast and flat sandy coast, were chosen to analysis how the three strategies impact on the extraction accuracy of coastlines. From the results shown in Figure 6, it can be seen that the results produced in strategies 2 and 3 have better performance than that in strategy 1, and strategy 3 produces the best results. For subsite 1 (artificial coast), the coastlines derived from eight fusion methods in strategy 3 are closest to the reference coastline, followed by the coastline extracted in strategy 2. The coastline extracted from original 30 m MS image presents the largest difference with the reference coastline, and the maximum distance between actual coastline and reference almost reaches more than one pixel in horizon direction. For subsite 2 (bedrock coast), the coastline has twists and turns. It can be seen that the coastline in strategy 1 cannot fully reflect the sinuosity features of bedrock coastline. On the contrary, the coastlines produced in strategies 2 and 3 are able to represent the sinuosity features, and the coastlines in strategy 3 provide much more details about the twists and are the closest to the reference coastline. Coastlines as shown in Figure 6q-x are the results of sandy coast. The coastline of strategy 1 is overestimated by comparing with the reference coastline, and it is on the sea side and almost one pixel (30 m) away from the accurate coastline. It is observed that the coastlines derived from fused MS images of all pansharpening algorithms in strategy 3 present the best results that are approximate to the reference coastline.

Comparising Extracted Coastlines of Different Pansharpening Algorithms
As presented in Section 4.2, eight fusion methods were finally applied for coastline extraction in strategies 2 and 3. In this section, we discussed the capacity of coastline extraction with the above eight pansharpening algorithms. Three subsites, referring to three different coast types, were shown in Figure 7, and the coastlines derived from different fused images are displayed with different colors. For the artificial coast shown in Figure 7a,d, the actual coastline refers to the boundary of construction and it tends to be a straight line. The results of different fusion methods agree well with the fluctuation trend of reference coastline. However, the coastlines of Brovey and PRACS algorithms are farther away from the reference coastline by comparing with the others, and they overlap with each other in both strategies 2 and 3. The coastlines of Indusion and AWLP in both strategies 2 and 3 are closer to the reference coastline than the others. For bedrock coast, as shown in Figure 7b,e, the coastlines extracted from Indusion images both in strategies 2 and 3 produced the most accurate and visually realistic presentation. For the sandy coast shown in Figure 7c,f, the coastline for ATWT in strategy 2 (see Figure 7c) is closest to the reference line, and the following is Indusion. In strategy 3 (see Figure 7f), the coastline extracted from the fused MS image of Indusion produced the most realistic and accurate representation, and the following are AWLP and ATWT.
Water 2017, 9,816 16 of 26 eight pansharpening algorithms. Three subsites, referring to three different coast types, were shown in Figure 7, and the coastlines derived from different fused images are displayed with different colors. For the artificial coast shown in Figure 7a,d, the actual coastline refers to the boundary of construction and it tends to be a straight line. The results of different fusion methods agree well with the fluctuation trend of reference coastline. However, the coastlines of Brovey and PRACS algorithms are farther away from the reference coastline by comparing with the others, and they overlap with each other in both strategies 2 and 3. The coastlines of Indusion and AWLP in both strategies 2 and 3 are closer to the reference coastline than the others. For bedrock coast, as shown in Figure 7b,e, the coastlines extracted from Indusion images both in strategies 2 and 3 produced the most accurate and visually realistic presentation. For the sandy coast shown in Figure 7c,f, the coastline for ATWT in strategy 2 (see Figure 7c) is closest to the reference line, and the following is Indusion. In strategy 3 (see Figure 7f), the coastline extracted from the fused MS image of Indusion produced the most realistic and accurate representation, and the following are AWLP and ATWT.  Table 3 summarizes the coastline extraction accuracy for three strategies. As shown in Table 3 and Figure 8, coastline derived from the original MS image (strategy 1) has the worst accuracy index values by comparing with the extracted coastlines in strategies 2 and 3, which demonstrates the efficiency of pansharpening algorithms for coastline extraction from Landsat-8 OLI imagery. Comparing the MAD values in Table 3 and Figure 8a, the coastline extracted from original MS image has larger MAD value than that of strategies 2 and 3; moreover, it is observed that most of the coastlines extracted in strategy 3 have lower MAD values and maximum AD values than that of strategy 2. The coastline derived from HPF fusion image in strategy 3 has the lowest MAD value of 11.26 m, and followed by ATWT and PRACS. AWLP in strategy 3 achieves the lowest maximum AD value of 124.17 m, while its MAD value is 13.24 m. This demonstrates that strategy 3 has obvious advantage to extract coastlines with lower MAD values by comparing with that of strategies 2 and 1. For the minimum AD values, the coastline can keep very close to the reference coastline. There is no obvious difference between different pansharpening algorithms in the three strategies. Through the coastline extracted from the fused MS image produced by IHS has low minimum AD values, the maximum AD can reach 438.82 m in strategy 2 and 301.81 m in strategy 3. Therefore, it is not a stable algorithm used to extract coastline from Landat-8 OLI imagery.  Table 3 summarizes the coastline extraction accuracy for three strategies. As shown in Table 3 and Figure 8, coastline derived from the original MS image (strategy 1) has the worst accuracy index values by comparing with the extracted coastlines in strategies 2 and 3, which demonstrates the efficiency of pansharpening algorithms for coastline extraction from Landsat-8 OLI imagery. Comparing the MAD values in Table 3 and Figure 8a, the coastline extracted from original MS image has larger MAD value than that of strategies 2 and 3; moreover, it is observed that most of the coastlines extracted in strategy 3 have lower MAD values and maximum AD values than that of strategy 2. The coastline derived from HPF fusion image in strategy 3 has the lowest MAD value of 11.26 m, and followed by ATWT and PRACS. AWLP in strategy 3 achieves the lowest maximum AD value of 124.17 m, while its MAD value is 13.24 m. This demonstrates that strategy 3 has obvious advantage to extract coastlines with lower MAD values by comparing with that of strategies 2 and 1. For the minimum AD values, the coastline can keep very close to the reference coastline. There is no obvious difference between different pansharpening algorithms in the three strategies. Through the coastline extracted from the fused MS image produced by IHS has low minimum AD values, the maximum AD can reach 438.82 m in strategy 2 and 301.81 m in strategy 3. Therefore, it is not a stable algorithm used to extract coastline from Landat-8 OLI imagery. Similar trend as that of MAD can also be observed for NSM index values shown in Table 3 and Figure   Similar trend as that of MAD can also be observed for NSM index values shown in Table 3 and    Figure 9 andFigure 10a,b, it can be seen that the position differences of coastlines in three strategies are mainly within a pixel (30 m) of the original MS image. The coastline in strategy 1 has a more significant position fluctuation than that in strategies 2 and 3. Moreover, both the position change fluctuations of AD and NSM in strategy 3 are slighter than that in strategy 2, which shows the superiority of strategy 3 by comparing with strategy 2. On the other hand, the extracted coastlines   Figures 9 and 10a,b, it can be seen that the position differences of coastlines in three strategies are mainly within a pixel (30 m) of the original MS image. The coastline in strategy 1 has a more significant position fluctuation than that in strategies 2 and 3. Moreover, both the position change fluctuations of AD and NSM in strategy 3 are slighter than that in strategy 2, which shows the superiority of strategy 3 by comparing with strategy 2. On the other hand, the extracted coastlines of different pansharpening algorithms in strategies 2 and 3 change on both sides of the reference coastline, while the MEAN lines change mainly on the ocean side of reference coastline.   To determine the potential use of coastline extracted with the integrated method, the histograms of NSM values (range from −30 m to 30 m) for the extracted coastlines with different methods are presented in Figure 11. For strategies 2 and 3, more than 94% of the NSM values are between −30 m   To determine the potential use of coastline extracted with the integrated method, the histograms of NSM values (range from −30 m to 30 m) for the extracted coastlines with different methods are presented in Figure 11. For strategies 2 and 3, more than 94% of the NSM values are between −30 m To determine the potential use of coastline extracted with the integrated method, the histograms of NSM values (range from −30 m to 30 m) for the extracted coastlines with different methods are presented in Figure 11. For strategies 2 and 3, more than 94% of the NSM values are between −30 m and 30 m. Instead of using histogram, the NSM values of coastline extracted in strategy 1 was drawn as a polyline for better comparison. It can be found that the NSM values of coastline in strategy 1 disperses along the coastline but with an obvious seaward bias. For strategies 2 and 3, the central NSM values of the extracted coastlines gathered closer to the ideal value of 0. For example, the percentage of the NSM values between −12 m and 12 m in strategy 1 is 55.05%. However, for strategy 2, the maximum percentage is 78.13%, and the least is 61.58%. In strategy 3, the percentage of the NSM values between −12 m and 12 m can reach to 88.90%, and the least is 88.13%, which demonstrates the efficiency of strategy 3.

Coastline Extraction Efectiveness
Based on quantitative and visual inspection of the resultant coastlines, it can be seen that the coastlines extracted in strategy 3 provide the most accurate and visually realistic representation, particularly for coasts with complex boundaries, such as bedrock coast. The following is the coastlines extracted in strategies 1 and 2. Increasing the spatial resolution of original MS image by directly using pansharpening algorithms or pansharpening with downscaled PAN image contributes to the accuracy increase of coastline extraction. With the increment of spatial resolution, a mixed pixel may be separated into two classes (water and non-water), and more details hidden in mixed pixels would be exploited in sharpened MS image. This is why coastlines extracted in strategy 3 are better than that in strategy 2. However, it cannot be concluded that with the continuing increase of spatial resolution of fused MS image, the resultant coastline would have continuous accuracy increase. For downscaling PAN band, the uncertainty would increase with the increment of spatial resolution of the downscaled PAN band by using spatial interpolation method.
In the study, ten pansharpening algorithms were used to produce fused image. Among the ten pansharpening algorithms, PCA and GS results have serious spectral distortions, which make them difficult to distinguish water and non-water from fused images. In general, water index is calculated based on spectral ratio and is sensitive to spectral distortion produced within the pansharpening process. CS-based pansharpening algorithms yield sharpened MS images by directly inject PAN band into original MS image, but often results in spectral distortion, whereas MRA approaches present less spectral distortion through a multiscale decomposition of PAN band into MS bands.
Comparing with all of the coastlines extracted from different pansharpening results in strategies 2 and 3, Indusion, ATWT and AWLP are the top three approaches. Instead of using spatial interpolation methods for downscaling original MS image, Induction scaling technique is used in Indusion pansharpening process for producing spectrally consistent downscaled MS image. For ATWT and AWLP, they use the "à trous" algorithm to produce a set of shift-invariant discrete wavelet planes with different degrees of spatial resolution, and additive method is applied to build a fused MS images. In wavelet-based image fusion process, the additive method preserves not only the high spatial information of PAN image but also all coarse spatial information of original MS image which is different from substitution method. In the condition, the ATWT and AWLP algorithms can sharpen the water and non-water boundaries, resulting in higher accuracy coastlines.

Threshold Value Effect on Coastline Extraction
The water index approach used in this study is based on the NDWI, and a key point is to generate water and non-water binary map from the NDWI image. Generally, the NDWI threshold value used to separate water and non-water features is not a constant value, and it varies depending on the land-cover components [47]. This is why the dynamic threshold selection method is needed for water and non-water features separation with NDWI image. Otsu method is a widely used dynamic threshold method to obtain optimal threshold from NDWI image by maximizing between the class various of water and non-water features [45]. For example, Du et al. [71] and Li et al. [38] have proved that the Otsu method can successfully be applied to map water bodies from Landsat imagery. In order to demonstrate the superiority of Otsu method against the constant threshold value of 0 in coastline extraction from Landsat-8 OLI imagery, the following example is presented. Specially, the subarea coastline extracted by using threshold of 0 is presented in Figure 12a, while Figure 12b shows the subarea coastline extracted with optimal threshold by applying Otsu method. Under the same condition, it can be seen that coastlines extracted with 0 threshold retreat toward land obviously and keep away from reference coastline, while that of Otsu is closer to the reference coastline. Moreover, for the whole study area shown in Figure 1, the MNSM of coastline extracted with threshold of 0 is −13.70 m and the MAD is 30.51 m, while that of Otsu result is 13.53 m (MNSM) and 18.62 m (MAD). Therefore, it is suggested to apply the Otsu method to extract coastline from Landsat-8 OLI imagery.

Discussions
Coastline extraction is important for the coastal management and coastline environment protection. The main limitations of coastline extraction from remote sensing image are the spatial resolution and cost. Landsat-8 OLI imagery is free access, but the spatial resolution limits the coastline extraction accuracy. Though several methods have been used and present in many previous research, automatic extraction method at subpixel level are difficult to understand and complex to implement. This study integrated pansharpening and water index approaches on Landsat-8 OLI imagery to discuss and analyze the ability of extracting coastline with higher accuracy. Pansharpening can take advantages of the fine spatial information within Landsat-8 OLI PAN band to increase the spatial resolution of original MS image, and water index was used to extract coastline from sharpened MS image.

Advance Accuracy Evaluation of the Proposed Method
To further evaluate the efficiency of various methods in coastline extraction, three original benchmark methods, including single-band thresholding (Band 6 for Landsat-8 OLI image) [72], classification (ISODATA) [35], water index (AWEI) [73], were used as the comparative methods against the method proposed in the paper. Table 4 shows the coastline extraction accuracy of three original methods. It can be found that the MNSM and MAD values of the three methods are stay at the high level, and much larger than that of strategies 2 and 3, especially the result of coastline extracted by AWEI method. This is because they are all based on the per-pixel scale (at spatial resolution of 30 m), while strategies 2 and 3 can take advantage of the PAN band to increase the accuracy of coastline extraction. In practice, there are many shortcomings of these methods, for example, AWEI combines all possible MS bands, which enhances the spectral characteristics of water and water-like features, but decrease the separability boundaries between them at the same time. On the contrary, single-band thresholding method cannot take full advantage of MS bands, and the extracted coastline position has a large fluctuation. For ISODATA method, the result is unstable, which is highly depends on the section of parameters during classification.

Discussions
Coastline extraction is important for the coastal management and coastline environment protection. The main limitations of coastline extraction from remote sensing image are the spatial resolution and cost. Landsat-8 OLI imagery is free access, but the spatial resolution limits the coastline extraction accuracy. Though several methods have been used and present in many previous research, automatic extraction method at subpixel level are difficult to understand and complex to implement. This study integrated pansharpening and water index approaches on Landsat-8 OLI imagery to discuss and analyze the ability of extracting coastline with higher accuracy. Pansharpening can take advantages of the fine spatial information within Landsat-8 OLI PAN band to increase the spatial resolution of original MS image, and water index was used to extract coastline from sharpened MS image.

Advance Accuracy Evaluation of the Proposed Method
To further evaluate the efficiency of various methods in coastline extraction, three original benchmark methods, including single-band thresholding (Band 6 for Landsat-8 OLI image) [72], classification (ISODATA) [35], water index (AWEI) [73], were used as the comparative methods against the method proposed in the paper. Table 4 shows the coastline extraction accuracy of three original methods. It can be found that the MNSM and MAD values of the three methods are stay at the high level, and much larger than that of strategies 2 and 3, especially the result of coastline extracted by AWEI method. This is because they are all based on the per-pixel scale (at spatial resolution of 30 m), while strategies 2 and 3 can take advantage of the PAN band to increase the accuracy of coastline extraction. In practice, there are many shortcomings of these methods, for example, AWEI combines all possible MS bands, which enhances the spectral characteristics of water and water-like features, but decrease the separability boundaries between them at the same time. On the contrary, single-band thresholding method cannot take full advantage of MS bands, and the extracted coastline position has a large fluctuation. For ISODATA method, the result is unstable, which is highly depends on the section of parameters during classification.

Limitaions and Future Scopes
It is true that uncertainty would be inevitably introduced in the downscaling process, and this is a common but challenging issue in the field of image downscaling [67]. In fact, image downscaling is indeed an ill-posed problem, and the most challenging issue is to reduce the uncertainty [74]. During the past decades, the study of reducing the image downscaling uncertainty is mainly focused on two aspects. Firstly, using auxiliary data, such as multiple shifted images and digital elevation model, to increase the image downscaling accuracy. However, auxiliary data is not always available in real applications, and the collection of auxiliary data is always costly. Secondly, developing some powerful downscaling algorithms, such as machine learning based [75] and deep learning-based [76] approaches, to increase the image downscaling accuracy. In general, downscaling algorithms based on machine learning approaches are always mathematic complicated and time-consuming by comparing with normal image interpolation method. However, there is still of great interest for us to develop the machine learning-based downscaling algorithms, as it has achieved good performance in the field of image super-resolution. In future works, we would focus on the reducing of uncertainty of downscaling-based coastline extraction by proposing novel machine learning-based downscaling algorithms.
Although the coastline extracted from fused MS image have better accuracy, it is expected to add another Landsat OLI image acquired at different dates to prove that the proposed method can be replicated with success. Under this circumstance, we need to ensure that the reference coastline is obtained at almost the same time as that of Landsat OLI imagery, in which to eliminate the tide level difference between Landsat-8 OLI image and reference image. Although Landsat-8 OLI image is readily to be collected, the corresponding high spatial resolution remotely sensed image that used to extract reference coastline is always difficult to be obtained, and the high spatial resolution remotely sensed image, such as ZY-3, IKONOS and Quickbird, is often costly. We also considered to add more Landsat-8 OLI images to validate the performance of the proposed method, but for the study area, it is almost impossible for us to collect other high spatial resolution remotely sensed image that acquired at the same time of Landsat-8 OLI image. Notably, the aim of the paper is to confirm the availability of strategy 2 and strategy 3 and provide a reference for which algorithm is better in coastline extraction. It is, therefore, of great interest for us to apply more datasets and study areas to validate the performance and applicability of the methods in the future work.
Considering the results obtained, the coastline of coast with irregular shape, such as bedrock coast, is more sensitive to the spatial resolution of MS image, and further study should be done to analysis the relationship between the coast type and pansharpening algorithms. For further demonstrating the performance of pansharpening algorithms and downscaling method in coastline extraction, methods based on the integration of pansharpening and other water indices, such as MNDWI, are needed to discuss in the future.

Conclusions
The main purpose of this study was to discuss and analyze the capacity of coastline extraction from Landsat-8 OLI imagery by integrating pansharpening and water index (NDWI) approaches. Based on the spatial resolution of MS images used for coastline extraction, three strategies were proposed; moreover, ten popular pansharpening algorithms, including five CS-based and five MRA-based pansharpening algorithms, were used. From the above results analysis, it is reasonable to conclude that the proposed strategies 2 and 3 are efficient and readily methodologies that can be used for coastline extraction from Landsat-8 OLI imagery.
Specially, the conclusions for this study are summarized into five aspects: (1) For the fusion of Landsat-8 OLI original MS image and PAN band, PCA and GS pansharpening algorithms result in fused MS images with serious spectral distortions, and almost cannot be used to distinguish water and non-water features with NDWI and are also not suggested to be used for coastline extraction. (2) The coastlines produced in strategies 2 and 3 have better performance than that in strategy 1.
With the increase of spatial resolution, subtle coastline changes which are indistinguishable in coarse spatial resolution MS images can be exploited well in the fused MS images in strategies 2 and 3. It indicates that pansharpening approach can improve the coastline extraction from Landsat-8 OLI imagery. (3) Strategy 3 produces coastlines with the best performance, indicating that further downscaling the PAN band is an alternative way to increase the coastline extraction accuracy. (4) For the ten widely used pansharpening algorithms, most of the MRA-based methods presents better results than that of CS-based methods in both strategies 2 and 3. Among the MRA-based methods, Indusion, ATWT and AWLP algorithms are the most efficient for increasing coastline extraction accuracy of Landsat-8 OLI imagery. (5) For the artificial coast, the coastlines extracted from most of the fusion methods agree well with the fluctuation trend of reference coastline. For both bedrock coast and sandy coast, the coastlines extracted from the fused images of Indusion algorithm in strategies 2 and 3 present the most accurate and visually realistic presentation.