forked from kilianv/CoastSat_WRL
You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
47 lines
1.0 KiB
Python
47 lines
1.0 KiB
Python
# -*- coding: utf-8 -*-
|
|
|
|
# Preamble
|
|
import ee
|
|
import matplotlib.pyplot as plt
|
|
import matplotlib.cm as cm
|
|
import numpy as np
|
|
import pandas as pd
|
|
from datetime import datetime
|
|
import pickle
|
|
import pdb
|
|
import pytz
|
|
from pylab import ginput
|
|
import scipy.io as sio
|
|
import scipy.interpolate
|
|
import os
|
|
|
|
# image processing modules
|
|
import skimage.filters as filters
|
|
import skimage.exposure as exposure
|
|
import skimage.transform as transform
|
|
import sklearn.decomposition as decomposition
|
|
import skimage.morphology as morphology
|
|
import skimage.measure as measure
|
|
|
|
# my functions
|
|
import functions.utils as utils
|
|
import functions.sds as sds
|
|
|
|
np.seterr(all='ignore') # raise/ignore divisions by 0 and nans
|
|
ee.Initialize()
|
|
|
|
with open('data\wl_final.pkl', 'rb') as f:
|
|
wl = pickle.load(f)
|
|
|
|
i = 0
|
|
x = wl[i]['quad_data']['x']
|
|
y = wl[i]['quad_data']['y']
|
|
z = wl[i]['quad_data']['z']
|
|
x = x.reshape(x.shape[0] * x.shape[1])
|
|
y = y.reshape(y.shape[0] * y.shape[1])
|
|
z = z.reshape(z.shape[0] * z.shape[1])
|
|
|
|
idx_nan = np.isnan(z)
|
|
x = x[~idx_nan]
|
|
y = y[~idx_nan]
|
|
z = z[~idx_nan] |