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'