default="default",
)
o = parser.parse_args(sys.argv[1:])
slope_value = None
if o.slope != "default":
slope_value = float(o.slope)
if o.site != "default" and o.start != "default" and o.end != "default":
filepath_data = os.path.join(os.getcwd(), "data")
sitename = o.site
metadata = []
kml_polygon = os.path.join(filepath_data, sitename, sitename + ".kml")
polygon = SDS_tools.polygon_from_kml(kml_polygon)
dates = [o.start, o.end]
sat_list = ["L5", "L7", "L8", "S2"]
pts_sl = np.expand_dims(np.array([np.nan, np.nan]), axis=0)
with open(
os.path.join(filepath_data, sitename,
sitename + "_shoreline.csv")) as csv_file:
csv_reader = csv.reader(csv_file, delimiter=",")
for row in csv_reader:
item = np.array([float(row[0]), float(row[1])])
pts_sl = np.vstack((pts_sl, item))
pts_sl = np.delete(pts_sl, 0, axis=0)
pts_world_interp = utils.get_interpolate_points(pts_sl)
# Kilian Vos WRL 2018
#%% 1. Initial settings
# load modules
import os
import numpy as np
import pickle
import warnings
warnings.filterwarnings("ignore")
import matplotlib.pyplot as plt
from coastsat import SDS_islands, SDS_download, SDS_preprocess, SDS_tools, SDS_transects
# region of interest (longitude, latitude in WGS84), can be loaded from a .kml polygon
polygon = SDS_tools.polygon_from_kml(
os.path.join(os.getcwd(), 'example', 'EVA.kml'))
# or enter the coordinates (first and last pair of coordinates are the same)
# polygon = [[114.4249504953477, -21.9295184484435],
# [114.4383556651795, -21.92949300318377],
# [114.4388731500701, -21.91491228133647],
# [114.4250081185656, -21.91495393621703],
# [114.4249504953477, -21.9295184484435]]
# date range
dates = ['2019-01-01', '2019-02-01']
# satellite missions
sat_list = ['S2']
# name of the site
sitename = 'EVA'
