kvos / coastsat.slope Goto Github PK

Dear Kilian Vos,
I'm trying to run the code,but I don't know how to set up the FES2014 global tide model.

Can you please tell me the solution?

Regards,
TaiQQ

Hey @kvos ,

I am using the CoastSat.slope to get the tide level from the FES 2014 global model as described here. All the setups seem good until I got The kernel appears to have died. It will restart automatically after I run the option 1 cell, I noticed someone as raised the issue once (I have no out of memory issue), but there was no suggested solution.
I did some search on the error message, but solutions are not related. I will appreciate a possible solution to this, thanks.

My Data's Date range is = ['1980-01-01', '2000-12-31']
ROI is (Karachi coast) = [[[67.00833391012465,24.78819451085644],
[67.05107759298598,24.78819451085644],
[67.05107759298598,24.81250397390977],
[67.00833391012465,24.81250397390977],
[67.00833391012465,24.78819451085644]]]
sat_list = ['L5' , 'L7']
I have used all this info to gather the .pkl, transects.geojson and the .output files needed for coastsat.slope module
i am having the error in estimating beach-face slopes along the transect as ' SDS_slope.plot_spectrum_all(dates,composite,tsall,settings_slope, title) ' and further errors related to this

Hello Kilian,
I am trying to install the FES2014 model, but when I get to step 6 Open the Anaconda prompt and activate your environment, then run: conda install -c fbriol fes

I get these messages :
Collecting package metadata (current_repodata.json): done
Solving environment: failed with initial frozen solve. Retrying with flexible solve.
Collecting package metadata (repodata.json): done
Solving environment: failed with initial frozen solve. Retrying with flexible solve.

PackagesNotFoundError: The following packages are not available from current channels:

• fes

Current channels:

• https://conda.anaconda.org/fbriol/win-32
• https://conda.anaconda.org/fbriol/noarch
• https://repo.anaconda.com/pkgs/main/win-32
• https://repo.anaconda.com/pkgs/main/noarch
• https://repo.anaconda.com/pkgs/r/win-32
• https://repo.anaconda.com/pkgs/r/noarch
• https://repo.anaconda.com/pkgs/msys2/win-32
• https://repo.anaconda.com/pkgs/msys2/noarch

To search for alternate channels that may provide the conda package you're
looking for, navigate to

https://anaconda.org

and use the search bar at the top of the page.

However, I found the fbriol page within Anaconda.org ( https://anaconda.org/fbriol/fes ) with the exact same instruction

To install this package with conda run:
conda install -c fbriol fes

Still not working.

Do you have any idea of how to override this obstacle?

Thanks in advance.

Hi kvos：
For slope estimation I used shoreline data from 2000-2023 for prediction and satellites used 'L5,L7,L8,L9',but the estimation was not too good. I searched for the reason and found that the revisit period is 16 days instead of 8 days, is it because of the reason that I chose the data after 2020? But another problem arises, as in figure 4: the range of years selected all have two satellites in operation, I hope you can help me where exactly the problem occurs.
Figure1

Figure2

Figure3

Figure4

Dear All,
I have properly installed pyfes using:
conda install pyfes -c fbriol
I checked that pyfes is in the conda list
then I entered "python" and I got:
(py39) PS C:\WINDOWS\system32> python
Python 3.9.18 | packaged by conda-forge | (main, Dec 23 2023, 16:29:04) [MSC v.1929 64 bit (AMD64)] on win32
Type "help", "copyright", "credits" or "license" for more information.

import pyfes

Now I would like to know how to run pyfes with the FES2024 Database

Could you please clarify ?

Thank you.

CoastSat Slope Website is Down

I went to go check out the CoastSat slope website and discovered it was down. Just thought I would give a heads up.

Thank you for your work on shoreline monitoring. I have been following the coastsat notebook. I need to do tidal correction, so I want to use the fes2014 tide model using this CoastSat.slope repository. After downloading and extracting load_tide and ocean_tide files, I can't find the .ini files in the folder. See the content of my directory in the image below. Your guide is highly appreciated. Thanks

Hello

I have installed coastsat and been using it. Then for Tidal correction, I was trying to install Coastsat.slope but was not able to figure out the steps. Can someone please help with the steps for installing Coastat.slope and FES2014.

Hello,
while running the CoastSat code a RunTimeError. So which cause can't be identified. The absolute path in both ini files are correctly set. Down below is the error I get while running the code:

`RuntimeError Traceback (most recent call last)
in
6 config_load = os.path.join(filepath, 'load_tide.ini') # change to load_tide.ini
7 print(config_ocean,config_load)
----> 8 ocean_tide = pyfes.Handler("ocean", "io", config_ocean)
9 load_tide = pyfes.Handler("radial", "io", config_load)
10

RuntimeError: ConfiguRADIAL_2N2__FILE, RADIAL_2N2__LATITUDE, RADIAL_2N2__LONGITUDE, RADIAL_2N2__AMPLITUDE, RADIAL_2N2__PHASE, RADIAL_EPS2__FILE, RADIAL_EPS2__LATITUDE, RADIAL_EPS2__LONGITUDE, RADIAL_EPS2__AMPLITUDE, RADIAL_EPS2__PHASE, RADIAL_J1__FILE, RADIAL_J1__LATITUDE, RADIAL_J1__LONGITUDE, RADIAL_J1__AMPLITUDE, RADIAL_J1__PHASE, RADIAL_K1__FILE, RADIAL_K1__LATITUDE, RADIAL_K1__LONGITUDE, RADIAL_K1__AMPLITUDE, RADIAL_K1__PHASE, RADIAL_K2__FILE, RADIAL_K2__LATITUDE, RADIAL_K2__LONGITUDE, RADIAL_K2__AMPLITUDE, RADIAL_K2__PHASE, RADIAL_L2__FILE, RADIAL_L2__LATITUDE, RADIAL_L2__LONGITUDE, RADIAL_L2__AMPLITUDE, RADIAL_L2__PHASE, RADIAL_LA2__FILE, RADIAL_LA2__LATITUDE, RADIAL_LA2__LONGITUDE, RADIAL_LA2__AMPLITUDE, RADIAL_LA2__PHASE, RADIAL_M2__FILE, RADIAL_M2__LATITUDE, RADIAL_M2__LONGITUDE, RADIAL_M2__AMPLITUDE, RADIAL_M2__PHASE, RADIAL_M3__FILE, RADIAL_M3__LATITUDE, RADIAL_M3__LONGITUDE, RADIAL_M3__AMPLITUDE, RADIAL_M3__PHASE, RADIAL_M4__FILE, RADIAL_M4__LATITUDE, RADIAL_M4__LONGITUDE, RADIAL_M4__AMPLITUDE, RADIAL_M4__PHA`

Hi Kilian,

I'm getting this issue while running Option 1: get tide levels from FES2014
I must be doing something wrong with FES2014 data.

Got this message:

ValueError Traceback (most recent call last)
in
15 dates_fes, tide_fes = SDS_slope.compute_tide(coords,settings_slope['date_range'],time_step,ocean_tide,load_tide)
16 # get tide level at time of image acquisition
---> 17 tide_sat = SDS_slope.compute_tide_dates(coords, dates_sat, ocean_tide, load_tide)
18
19 # plot tide time-series

~\anaconda3\CoastSat\CoastSat.slope\SDS_slope.py in compute_tide_dates(coords, dates, ocean_tide, load_tide)
482 lats = coords[1]*np.ones(len(dates))
483 # compute heights for ocean tide and loadings
--> 484 ocean_short, ocean_long, min_points = ocean_tide.calculate(lons, lats, dates_np)
485 load_short, load_long, min_points = load_tide.calculate(lons, lats, dates_np)
486 # sum up all components and convert from cm to m

ValueError: Incompatible arguments:

1. calculate(self, numpy.ndarray[m, 1], numpy.ndarray[m, 1], numpy.ndarray[datetime.datetime[m, 1]])
2. calculate(self, numpy.ndarray[m, 1], numpy.ndarray[m, 1], numpy.ndarray[datetime64[us][m, 1]])

Thank you!

Juliana

Hi! I've been using CoastSat and now I'm trying to improve my results with tidal correction, so now I'm using the results given by coastsat slope. The thing is, it doesn't seem to save the .csv file obtained by FES 2014, wich I need as an input at coastsat in tidal correction section. So the question is, how could I edit the coastsat slope code so it saves the .csv file of tide data from FES?

Thanks a lot,
Paula.

Hi Kilian,

I'm trying to run CoastSat.slope code, but i get stucked in ### Option 1: get tide levels from FES2014###

I get...

RuntimeError Traceback (most recent call last)
in
5 config_ocean = os.path.join(filepath, 'ocean_tide.ini') # change to ocean_tide.ini
6 config_load = os.path.join(filepath, 'load_tide.ini') # change to load_tide.ini
----> 7 ocean_tide = pyfes.Handler("ocean", "io", config_ocean)
8 load_tide = pyfes.Handler("radial", "io", config_load)
9

RuntimeError: No such file or directory : C:\Users\Juliana\anaconda3\CoastSat-master\CoastSat-master\fes\data\fes2014\ocean_tide\ocean_tide\mm.nc

But the file really exists.

Followed the steps...but i must be doing something wrong.

Thank you.

Juliana

Hi @kvos, it is a very nice job. I have some questions for you. Why is my estimated slope all one value (set maximum value)? Sometimes it calculates the correct value, sometimes it doesn't. This confuses me very much. I only extract one shoreline each year for analysis.

Hellol Everyone, Hi @kvos
i installed the FES 2004, created output data from 2000 to 2020 (L5, L7, L8), 2 transects. I have a polar coastline so I have only more than 40 satellite images. Its depend of the day/night time in Arctic (polar day). What do You think about this results? Is it measurable? The frequency images of the tides is not similar to your, exampled situation. Thanks for Your answer, hope so to get good tips or something like that.
slope_kaff.pdf

i have encountered this error and cant seem to resolve this , have looked it up as weel and it stated to take a min of 20 years of data and that i have but still same error

Hello @kvos! Firstly, I appreciate your work with the CoastSat and CoastSat.slope! I'm using these techniques in my master study about some beaches in the Brazilian south coastal. Here the tidal is microtide, my beaches have short backshore and they can be high slope. How can I know wether my results (with large confidence band) are affected by this metodology error or my imput datas are insufficient? Can I somehow identify the influence of this problem (TR/tanβ<10) on my results?

Here is an example of one of my results:

This line in coastsat.slope keeps returning an error:
settings_slope['freqs_max'] = SDS_slope.find_tide_peak(dates_sat,tide_sat,settings_slope)

  File "C:\Users\RDCHLNRO\AppData\Local\Continuum\anaconda3\envs\coastsat3152021\lib\site-packages\coastsat\SDS_slope.py", line 506, in find_tide_peak
    idx_max = idx_peaks[freqs[idx_peaks] > settings['freqs_cutoff']][0]

IndexError: index 0 is out of bounds for axis 0 with size 0

These are the variables being sent to the function:

([datetime.datetime(2018, 6, 8, 15, 32, 54, tzinfo=<UTC>),
  datetime.datetime(2018, 6, 15, 15, 39, 9, tzinfo=<UTC>),
  datetime.datetime(2018, 6, 24, 15, 33, 4, tzinfo=<UTC>),
  datetime.datetime(2018, 7, 1, 15, 39, 19, tzinfo=<UTC>),
  datetime.datetime(2018, 7, 10, 15, 33, 13, tzinfo=<UTC>),
  datetime.datetime(2018, 7, 17, 15, 39, 27, tzinfo=<UTC>)],
 array([0.367, 0.358, 0.363, 0.358, 0.363, 0.358]),
 {'slope_min': 0.035,
  'slope_max': 0.2,
  'delta_slope': 0.005,
  'date_range': [datetime.datetime(2018, 6, 1, 0, 0, tzinfo=<UTC>),
   datetime.datetime(2018, 8, 1, 0, 0, tzinfo=<UTC>)],
  'n_days': 8,
  'n0': 50,
  'freqs_cutoff': 3.8580246913580245e-07,
  'delta_f': 1e-08})

Maybe there is an issue with the small amount of data being sent? That was not an issue for me previously, but I recently updated all of the coastsat code.

Thank you for your time.

Hello Kilian,
I just downloaded the load_tide and ocean_tide which contain the 34 files, however, in your example there are two files (load_tide.ini and ocean_tide.ini)
My question is: where these two files downloaded from the Aviso FTP ( I could not find them) or did you made them this using the info provided in the image?

best regards

Tona

Hi Kilian,

Thanks for these helpful tools, I tried to use measured tide data to estimate the slope
I got errors on steps number 4, 5

4. Peak tidal frequency

ValueError Traceback (most recent call last)
in
4 fig, ax = plt.subplots(1,1,figsize=(12,3), tight_layout=True)
5 ax.grid(which='major', linestyle=':', color='0.5')
----> 6 bins = np.arange(np.min(delta_t)/seconds_in_day, np.max(delta_t)/seconds_in_day+1,1)-0.5
7 ax.hist(delta_t/seconds_in_day, bins=bins, ec='k', width=1);
8 ax.set(xlabel='timestep [days]', ylabel='counts',

~\anaconda3\envs\coastsat\lib\site-packages\numpy\core\fromnumeric.py in amin(a, axis, out, keepdims, initial)
2616 """
2617 return _wrapreduction(a, np.minimum, 'min', axis, None, out, keepdims=keepdims,
-> 2618 initial=initial)
2619
2620

~\anaconda3\envs\coastsat\lib\site-packages\numpy\core\fromnumeric.py in _wrapreduction(obj, ufunc, method, axis, dtype, out, **kwargs)
84 return reduction(axis=axis, out=out, **passkwargs)
85
---> 86 return ufunc.reduce(obj, axis, dtype, out, **passkwargs)
87
88

ValueError: zero-size array to reduction operation minimum which has no identity

5. Estimate the beach slope

ValueError Traceback (most recent call last)
in
10 tsall = SDS_slope.tide_correct(composite,tide,beach_slopes)
11 title = 'Transect %s'%key
---> 12 SDS_slope.plot_spectrum_all(dates,composite,tsall,settings_slope, title)
13 slope_est[key] = SDS_slope.integrate_power_spectrum(dates,tsall,settings_slope)
14 print('Beach slope at transect %s: %.3f'%(key, slope_est[key]))

~\OneDrive - Universidade do Algarve\New coastsat\CoastSat.slope-master\CoastSat.slope-master\SDS_slope.py in plot_spectrum_all(dates_rand, composite, tsall, settings, title)
580 days_in_year = 365.2425
581 time_step = settings['n_days']*seconds_in_day
--> 582 freqs = frequency_grid(t,time_step,settings['n0'])
583 beach_slopes = range_slopes(settings['slope_min'], settings['slope_max'], settings['delta_slope'])
584

~\OneDrive - Universidade do Algarve\New coastsat\CoastSat.slope-master\CoastSat.slope-master\SDS_slope.py in frequency_grid(time, time_step, n0)
529 def frequency_grid(time,time_step,n0):
530 'define frequency grid for Lomb-Scargle transform'
--> 531 T = np.max(time) - np.min(time)
532 fmin = 1/T
533 fmax = 1/(2*time_step) # Niquist criterium

~\anaconda3\envs\coastsat\lib\site-packages\numpy\core\fromnumeric.py in amax(a, axis, out, keepdims, initial)
2503 """
2504 return _wrapreduction(a, np.maximum, 'max', axis, None, out, keepdims=keepdims,
-> 2505 initial=initial)
2506
2507

~\anaconda3\envs\coastsat\lib\site-packages\numpy\core\fromnumeric.py in _wrapreduction(obj, ufunc, method, axis, dtype, out, **kwargs)
84 return reduction(axis=axis, out=out, **passkwargs)
85
---> 86 return ufunc.reduce(obj, axis, dtype, out, **passkwargs)
87
88

ValueError: zero-size array to reduction operation maximum which has no identity

FYI: here the code I wrote to write the .pkl file
that might be the problem? (I am not a python user, I collected from videos)
import datetime
import pytz
import numpy as np
data_dict = {'dates': [datetime.datetime(1997,1,8,10,32,6,0,pytz.UTC),
datetime.datetime(1997,1,24,10,32,50,0,pytz.UTC),
datetime.datetime(1997,2,9,10,33,32,0,pytz.UTC),
datetime.datetime(1997,2,25,10,34,12,0,pytz.UTC),
datetime.datetime(1997,3,29,10,35,27,0,pytz.UTC),
datetime.datetime(1997,5,16,10,37,7,0,pytz.UTC),
datetime.datetime(1997,6,1,10,37,42,0,pytz.UTC),
datetime.datetime(1997,6,17,10,38,14,0,pytz.UTC),
datetime.datetime(1997,7,3,10,38,48,0,pytz.UTC),
datetime.datetime(1997,8,20,10,40,23,0,pytz.UTC),
datetime.datetime(1997,9,5,10,40,51,0,pytz.UTC),
datetime.datetime(1997,9,21,10,41,18,0,pytz.UTC),
datetime.datetime(1997,10,7,10,41,46,0,pytz.UTC),
datetime.datetime(1997,10,23,10,42,12,0,pytz.UTC),
datetime.datetime(1997,11,8,10,42,37,0,pytz.UTC),
datetime.datetime(1997,12,26,10,43,45,0,pytz.UTC),
datetime.datetime(1998,6,4,10,46,29,0,pytz.UTC),
datetime.datetime(1998,7,6,10,46,53,0,pytz.UTC),
datetime.datetime(1998,7,22,10,47,2,0,pytz.UTC),
datetime.datetime(1998,8,7,10,47,6,0,pytz.UTC),
datetime.datetime(1998,8,23,10,47,17,0,pytz.UTC),
datetime.datetime(1998,9,24,10,47,34,0,pytz.UTC),
datetime.datetime(1998,11,27,10,47,44,0,pytz.UTC),
datetime.datetime(1998,12,13,10,47,40,0,pytz.UTC)],
'tide': np.array([0.1418,-0.8282,-1.6682,-1.4182,-0.9782,0.6418,1.0418,0.7118,0.1118,-1.3282,-1.0282,-1.2682,-0.7882,
0.1418,0.3018,0.6418,0.6918,0.7318,0.3518,0.2618,-0.6582,-1.0482,-0.0382,0.5018])}

with open ('NARRA_test_tide.pkl','wb') as f:
pickle.dump(data_dict , f)

with open ('NARRA_test_tide.pkl','rb') as f :
newdata = pickle.load(f)

Thanks
Ahmed 

Hi @kvos , thanks very much for your help up to this point with navigating issues with running the script and getting FES model to work. I am able to run the slope estimation script end-to-end now, but I am concerned that the final output is not accurate.

Of particular concern to me are 1) that the Power Spectrum and "Inset into the tidal peak frequency bands" plots look very different from the example script outputs, 2) the confidence bands on the slope estimates are very large, 3) the estimate is the same for each transect.

I have attached the full output of my script here
slope_code_outpu.pdf

As well as the Output.pkl file
Pekalongan_longer_time_output.zip

and the tidal data CSVs
tide_fes.csv
tide_sat.csv
date_fes_ts.csv
date_sat_ts.csv

I am continuing to read through the documentation and publications, but do you have any intuition as to what might be the issue here @kvos to point me in the right direction?

Hi, may I verify where the load_tide.ini and ocean_tide.ini files are located? Visiting the bitbucket URL mentioned in a previous (now closed) issue mentions that the directory doesn't exist.

Hi,
I'm trying to run CoastSat on Colab. I managed to get it running until the point where pyfes is needed. I managed to install pyfes but when I run the following line:

# Option 1. if FES2014 global tide model is setup
# point to the folder where you downloaded the .nc files
filepath = r'/content/drive/MyDrive/aviso-fes-main/data/fes2014'
config_ocean = os.path.join(filepath, 'ocean_tide.ini') # change to ocean_tide.ini
config_load =  os.path.join(filepath, 'load_tide.ini')  # change to load_tide.ini
ocean_tide = pyfes.Handler("ocean", "io", config_ocean)
load_tide = pyfes.Handler("radial", "io", config_load)

I get the following error:

AttributeError                            Traceback (most recent call last)
[<ipython-input-14-f2ca7493019e>](https://localhost:8080/#) in <cell line: 7>()
      5 config_ocean = os.path.join(filepath, 'ocean_tide.ini') # change to ocean_tide.ini
      6 config_load =  os.path.join(filepath, 'load_tide.ini')  # change to load_tide.ini
----> 7 ocean_tide = pyfes.Handler("ocean", "io", config_ocean)
      8 load_tide = pyfes.Handler("radial", "io", config_load)
      9 

AttributeError: module 'pyfes' has no attribute 'Handler'

Any idea?
Thanks

We see the basic workflow is documented here requires the follow files:

1. netcdf files downloaded from AVISO - this part is really straightforward using CyberDuck FTP, but (apparently like a lot of users) I couldn't figure it out via the website
2. 'ocean_tide.ini' - where is this file? It is not in the downloaded tar.xz file from AVISO. what should it include? do you have examples?
3. 'load_tide.ini' - where is this file? It is not in the downloaded tar.xz file from AVISO
4. Why not make 'ocean_tide.ini' and 'load_tide.ini' available through CoastSat or CoastSat.slope? I assume they are small files?

I'm also curious as to why this is only documented here in CoastSat.slope, not CoastSat - isn't this a necessary step for any location with a measurable tide? Therefore why is it only documented and implemented here in CoastSat.slope?

(sorry for the newbie question, I still haven't gotten around to actually using CoastSat.slope, but we're trying to figure out how to wrap the functionality of vanilla CoastSat in CoastSeg while providing the user option to input slope and provide tide csv file)

thanks!

While running the code of CoastSat.slope, I reiceved an error of "boolean index did not match indexed array along dimension 0; dimension is 0 but corresponding boolean dimension is 3".

The dates for the shoreline in Coastsat is from '2020-01-01' to '2020-02-28', and I have used 3 shorelines of date 2020-01-16, 2020-02-01 and 2020-02-17.
These 3 shoreline were delineated from Landsat 8.

In the error they have highlighted "Cross_distance" is it something related to transect.

Does anyone have any solutions why this error is occurring !

@kvos I have also attached the screenshot of the error.

I would like to ask if it is possible to export the simulated tidal data values from the FES2014 simulation process? I would like to use the simulated tidal data as input for the coastsat tidal correction, is it possible?

Hello my friend!
Thank you for your excellent work. And I used your many jobs to help me do something important.
I have a new question about the CoastSat.slope, I want to get the slope, but the beach has many architectures and other things, so can this job get the correct estimated slope of the beach?

Hi Kilian, So sorry.

I have encountered some problems now, the date I set is [2015-01-01', '2023-01-01'], but the tide level correction step has the following error.

Later, I set the period length to one month （2017-12-01,2018-01-01）, but it did not report an error. I tried many reasons but failed. So I would like to ask why this happens. Attached below are my steps to install FES2014.

Thank you Kilian.

Hey,

I have problems with the beach slope. When I enter the first command for the initial settings, I get a ModuleNotFoundError: No module named SDS_slope. Does anyone know what the problem could be and how to solve it? I am unfortunately a total beginner in this area

Thanks a lot

Hi Kvos, I follow your instruction till to step6. after I download anaconda, I open 'Anaconda prompt' and type $ conda create -n fes2014 python=3.9. and then type $ conda activate fes2014. I follow typing $ conda install -c fbriol pyfes. the error msg shows as

I also tried to reset the env with different versions of python, still face the similar error. Kindly seek your help
thanks
jie

Hi,
I have loaded both examples and get a similar false output - my slopes are the same for each transect, and also, it is the upper bound that was specified (0.2). I have not changed any of the code in the example file; therefore, I do not understand where this error could have come in. I believe this might be a similar issue to #18 or #17 however, I could not spot an error (especially since I have not changed the original example code). Do you perhaps have any idea why this might occur @kvos ?

My code & Output:
Example_Error_slope_Narrabeen.pdf

Hello all,
I am currently attempting to download the above mentioned auxiliary files as per FES2014 installation instructions (https://github.com/kvos/CoastSat.slope/blob/master/doc/FES2014_installation.md), however, I am having trouble actually downloading them. They seem to be approximately 1.3GB and my downloads seem to be failing. I am accessing them through the TDS access service on aviso. The FTP seems to be unresponsive and I am downloading on google chrome. I have been trying for days to no avail. Any advice? What am I missing?

Hi Kilian,
I have the same issue as the question#4,
When i try to plot the transects and shorlines, only transects with no shorelines. But i can plot the example correctly.

Hi there!

I'm unable to successfully estimate beach slopes for my region of interest (error code 'ValueError: zero-size array to reduction operation maximum which has no identity'). When looking at the cross_distance dictionary output, many (but not all) of my transect arrays are are full of NaN values, and I'm not sure why this is happening. I tried a few basic fixes, including swapping land/seaside coordinates for my transects, confirming matching coordinate systems, etc. but haven't had much luck.

Any suggestions/clarification on where I am going wrong would be great! Thanks in advance.

** Update
I have been exploring my issue further but haven't made much progress. The empty arrays always occur for some transects, but not for others. I've confirmed that the transect file is properly formatted with correct coordinates, projections, etc. I would assume I would have an empty array if a given transect didn't cross through any shorelines, but this also doesn't appear to be my problem (see last figure below for an isolated example). What are some other potential reasons I could be getting empty arrays here?

Transect and shoreline locations

Example output of cross_distance dict

Example of a single transect that appears to cross multiple shorelines but still produces empty array

Hi @kvos, thank you for all the help with regards to making sense of the script and getting the FES2014 model to run.
I managed to get the diagrams for the energy tidal frequency band, but I'm confused by why I can't find the slope of the section of coast I used.

The energy in the energy tidal frequency band diagrams, for all transects drawn, appears to taper off without a defined minimum, continuously decreasing even when high values of max_slope were used for testing.

Attached is the notebook and the outputs, including the diagrams:
notebook_1999-2020_slopecalculation.pdf
Apologies for the extra comments and cells, I had to make these to successfully run the FES2014 model in a separate notebook and use the outputs from that to complete the rest of the processing in this notebook.

I'm looking for publications or documentation about how to interpret the power spectrum and tidal peak frequency bands to deduce what I've done incorrectly, but do you have any suggestions on what might've caused this?