Table of Contents¶
- Necessary imports
- Connect to your GIS
- Get the data for analysis
- Prepare data for analysis
- Preprocessing of the data
- Band ratio technique
- Postprocessing of results
- Data resources
Due to anthropogenic activities and natural processes i.e. changes in sea level, sedimentation and wave energy coastlines are changing throughout the world. Coastline is contact line between sea and land, it is an important linear feature on earth's surface with dynamic nature.
Traditionally, the coastlines were manually digitized which was time-consuming and labour intensive. Remote sensing is a good alternative to extract coastlines using satellite imagery, this way both temporal and spatial aspects can be covered.
Satellite imagery of visible range can be used for interpretation and can be easily obtained. But the imageries covering infrared wavelength is best to extract boundary between land and water. So, the satellites which covers both visible and infrared spectrum are widely accepted for coastline extraction and mapping.
Landsat-8 multispectral imagery is used in the current study as it covers a wavelength ranging from 0.43 to 12.51 micrometers, and hence suitable for coastal and aerosol studies.
import os import glob from zipfile import * import arcgis import arcpy from arcpy.management import PolygonToLine from datetime import datetime import pandas as pd from arcgis.features import GeoAccessor, GeoSeriesAccessor from arcgis.raster.analytics import convert_feature_to_raster, convert_raster_to_feature from arcgis.geoanalytics.manage_data import clip_layer from arcgis.raster.functions import equal_to, greater_than, clip, apply, extract_band, stretch
from arcgis import GIS gis = GIS('home') gis2 = GIS('https://pythonapi.playground.esri.com/portal', 'arcgis_python', 'amazing_arcgis_123')
Multispectral Landsat includes Landsat GLS and Landsat 8 imagery for use in visualization and analysis. This layer is time enabled and includes a number band combinations and indices rendered on demand. The Landsat 8 imagery includes eight multispectral bands from the Operational Land Imager (OLI) with 30m spatial resolution and two bands from the Thermal Infrared Sensor (TIRS) of 100m spatial resolution. It is updated daily with new imagery directly sourced from the Landsat on AWS collection.
landsat_item = gis.content.search('title:Multispectral Landsat owner:esri_livingatlas tags:Multitemporal, imagery, landsat 8, temporal, MS', 'Imagery Layer', outside_org=True) landsat = landsat_item.layers landsat_item
The buffer feature of 20 km was created from USA boundaries. This feature layer geometry will be used to get the Landsat-8 tiles of coastal areas.
aoi = gis.content.search('usa_coast_buff_f', 'feature layer') aoi