Clipping trajectories to eurofirs

Hello, I would like to clip so6 trajectories to specific FIRs. For instance, I have collected all the trajectories crossing a specific FIR:

bdx_so6_LIMM = so6_m1.inside_bbox(eurofirs['LIMM'])

and now I would like to clip them. Unfortunately the clip command does not work with the eurofirs objects (I don’t have AIRAC files):

bdx_so6_LIMM['BCS361'].clip(eurofirs['LIMM'].flatten())

I get the following error: NotImplementedError: Multi-part geometries do not provide a coordinate sequence

I’m probably missing something. Is there a better way to do what I’m trying to do?

Thank you, Rob

[xoolive]

@obbe79 I started to have a look but my faulty flight actually did not intersect the sector (so clip would fail, it will return None in the corrected version).

Could you please try bdx_so6_LIMM['BCS361'].to_pickle('bcs361.pkl') and send the pickled file to my personal email so that I can be sure I debug this particular case?

Thank you for your support

[xoolive]

This is indeed the problem I suspected 😃

I will fix it the same way it is done now with Flight classes, i.e. you will get the trajectory between the first point of entry and the last point of exit of the trajectory. (In practice, this is what makes the most sense from an ATC point of view)

What you noticed with inside_bbox is correct according to the specification. The thing is that since this first check is very fast to compute, it makes sense to apply it before intersects (and not try to check for an intersection if it doesn’t cross the bounding box)