从曲线卫星数据集中提取区域答案

【问题标题】：Extract region from a Curvilinear satellite Dataset从曲线卫星数据集中提取区域
【发布时间】：2017-06-28 17:52:36
【问题描述】：

我有来自 MODIS 的卫星条带数据，需要提取数据的子集（区域）进行分析（不是绘图）。我正在尝试找到使用可能很慢的循环的最佳方法。过去我使用过 set.intersect，但这不适用于 2D 数据。

我的问题是 Lat 和 Lon 都是 2D 的，我需要找到满足我的条件的索引 (lat>=x1)&(lat

纬度样本

经度样本

气溶胶元数据

到目前为止的代码

通常（对于一维纬度/经度）我会使用 Opt_Depth_Land[:,goodlat,goodlon] 来提取我的数据，但这不适用于此类数据集。

任何帮助将不胜感激。

【问题讨论】：

标签： python-2.7 hdf satellite

【解决方案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))

【讨论】：