2017-06-28 9 views
0

私はMODISの衛星スワスデータを持っており、解析するデータのサブセット(領域)を抽出する必要があります(PLOTではありません)。私は遅くなる可能性のあるループでこれを行う最善の方法を見つけようとしています。過去に私はset.intersectを使用しましたが、これは2Dデータでは機能しません。曲線衛星からの領域抽出データセット

問題はLatとLonの両方が2Dであるため、私の条件が満たされている位置(lat> = x1)&(lat < = x2)とlonの類似点を見つける必要があります。その後、私のメインのデータセット(エーロゾルの光学的深さ)

緯度サンプルLatitude Sample

経度サンプルLongitude Sample

エアゾールメタデータAerosol Metadata

コードこれまでCode

をスライスするために、これらの2Dインデックスを使用します通常、(1D lat/lonの場合)私はデータを抽出するためにOpt_Depth_Land [:、goodlat、goodlon]を使用しますが、これはこのタイプのd ata set。

ご迷惑をおかけして申し訳ございません。

答えて

0
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)) 
関連する問題