如何通过与特定国家/地区相交的轨道路径过滤轨道路径？

Question

我有一吨的轨道路径，每个包含63000个纬度和经度坐标，这使得像轨道路径列表文件的这个。

这些文件是 HDF 文件，包含两个列表，一个是 63k 纬度点，另一个是 63k 经度点，按顺序排列。 这是一个示例文件。

要选择纬度和经度坐标，请使用以下代码：

hdf = SD(FILEPATH, SDC.READ)
lat = hdf.select('Latitude')
lon = hdf.select('Longitude')

我需要能够过滤掉这些文件，以便我只有轨道路径与格陵兰国家相交的文件。

目前我这样做的效率非常低，我的解决方案现在遍历每个文件路径的每个点并检查该点是否在格陵兰岛内，如果是，我将其添加到文件列表中。 这可能需要几个小时才能完成，这是目前代码的基本内容。

#poly is a polygon around Greenland. 
coords = [(-55, 59.5), (-55,67.5), (-60,67.5), (-60,75), (-73.25,75), (-73.25,79.2),
        (-67.5,79.2), (-67.5,80.82), (-65.34,80.82), (-65.34,81.23), (-62,81.23),
        (-62,82), (-60.25,82), (-60.25,84), (-10,84), (-10,75), (-17,75), (-17,67.5),
        (-30,67.5), (-30,59.5)]
poly = Polygon(coords)

latlon = zip(latitude, longitude)
for index, lalo in enumerate(latlon):
     if (poly.contains(Point(lalo[1], lalo[0]))):
          greenland_files.append(f)
          break

有任何想法吗？

Answer 1

如果您的轨道点按顺序排序（我假设它们必须是），您可以尝试某种搜索技术。 例如，如果您的所有轨道点都是从东到西（反之亦然，但一致），并且如果您可以找到保证位于格陵兰岛以东且保证位于格陵兰岛以西的点，则您应该能够搜索在这些点之间。

我意识到这些情况会变得很奇怪，因为如果你走得足够远，地球上的任何一点都同时位于格陵兰岛的西部和东部，但我也不知道你的轨道是如何构建的（每次的起点都一样吗？他们总是从东到西进行，反之亦然？等）。

编辑：要快速找到两个多边形之间的交点，您可以使用以下算法： https : //dl.acm.org/doi/abs/10.1145/800141.804662