Nan values in cross distance about coastsat HOT 4 OPEN

are the transects longer than 1km? there is something hard-coded to be within 1km from transect origin, see issue #499 .
Please try first to compute the intersections with the normal function SDS_transects.compute_intersections(), before using the QA version that has more parameters. Let me know if that fixes the issue.

from coastsat.

Actually no, they are shorter! I'll have a look at that issue and let you know.
Thanks!

from coastsat.

did you try the simple compute_intersections()?

from coastsat.

I found that the issue is related to the clipping operation of the output points I added to the script.
While in Qgis I can see that the Coastsat output points and the clipped ones have the same extension, for example when I plot them by:

_fig = plt.figure()
plt.axis('equal')
plt.xlabel('Eastings')
plt.ylabel('Northings')
plt.grid(linestyle=':', color='0.5')

for i in range(len(output['shorelines'])):
sl = output['shorelines'][i]
date = output['dates'][i]
plt.plot(sl[:,0], sl[:,1], '.', label=date.strftime('%d-%m-%Y'))
#plt.legend();_

the extension of the clipped dataset is smaller and equal to the area where the transect values are computed.
I will try to fix this by analysing which parameter of my clipping operation can affect how the dataset extension is read.
Do you have any suggestions?

Thanks

Enrica

from coastsat.