从Curvilinear卫星数据集中提取区域

Question

我有来自MODIS的卫星扫描数据，需要提取要分析的数据子集（区域）（NOT PLOT）。我试图找到最好的方法来做出这可能很慢的循环。在过去，我使用过set.intersect，但这不适用于2D数据。

我的问题是Lat和Lon都是2D，我需要找到满足条件的索引（lat＆gt; = x1）＆amp;（lat＆lt; = x2）和lon类似。然后使用这些2D索引切割我的主要数据集（气溶胶光学厚度）

纬度示例

经度示例

气溶胶元数据

代码到目前为止

通常（对于1D lat / lon）我会使用Opt_Depth_Land [：，goodlat，goodlon]来提取我的数据，但这不适用于这种类型的数据集。

非常感谢任何帮助。

Answer 1

    valid_lat=(lat>=user_lat-radius)&(lat<=user_lat+radius)
    valid_lon=(lon>=user_lon-radius)&(lon<=user_lon+radius)
    Valid_Coord=np.where((valid_lat==True)&(valid_lon==True))