Extrahiere Region von einem kurvilinearen Satelliten Datensatz

Question

Ich habe Satellitenschwaddaten von MODIS und muss eine Teilmenge (Region) von zu analysierenden Daten extrahieren (NOT PLOT). Ich versuche den besten Weg zu finden, dies mit langsamen Schleifen zu tun. In der Vergangenheit habe ich set.intersect verwendet, aber dies funktioniert nicht für 2D-Daten.

Mein Problem ist sowohl Lat und Lon sind 2D und ich muss die Indizes finden, wo meine Bedingungen erfüllt sind (lat> = x1) & (lat < = x2) und ähnlich für lon. und dann diesen 2D-Indizes verwenden, um meinen Hauptdatensatz (AOD)

Latitude Probe

Longitude Probe

Aerosol MetaDaten

-Code so Far

zu schneiden Normalerweise (für 1D lat/lon) würde ich Opt_Depth_Land [:, goodlat, goodlon] verwenden, um meine Daten zu extrahieren, aber dies funktioniert nicht für diese Art von d Ata-Set.

Jede Hilfe würde sehr geschätzt werden.

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))