I have data containing the latitude, longitude and altitude of a flight path sampled at 1 sec intervals in an [n,3] array.
flightPath = [-28.0206900000000,153.379210000000,12300;
I also have a rectangular grid structure within the 3d space this flight path exists, defined simply as:
lat = linspace(latMin,latMax,latDivisions);
lng = linspace(lngMin,lngMax,lngDivisions);
alt = linspace(altMin,altMax,altDivisions);
I want to find every area within this grid structure that the flight path enters. I know I could simply check which grid area each coordinate lies within, however the sample rate would have to be infinitely small to ensure no areas are missed and this is not feasible as I have nearly a million paths to process.
I was wondering if it would be possible to represent the flight path with an equation and then find which faces of the grid the path intersects with. It is not important if the equation representing the flight path is not exact (so long as it could be converted to discrete [lat,lng,alt] positions later) as the purpose of this work is to create a 'truth data set' which can be used to determine the accuracy of sampling using different structures.
Any other ideas are more than welcome!