列表以列表或多維名單列表的元素作爲網格

Question

我試圖創建一個包含發現指數在兩個條件得到滿足的經/緯度組合的陣列緯度/經度網格。這種方法可能太複雜，但使用meshgrid或numpy廣播也失敗了。如果有更好的方法，請隨時分享你的知識。 :-)

回合緯度/經度值gridsize的1°的分辨率，但是保留陣列的全長：

x = np.around(lon, decimals=0) 
y = np.around(lat, decimals=0) 

陣列由經度/緯度值的從-180到180和-82°到82°;可能的多douplets

檢查經/緯度多少度量可用於1°/ 1°網格點的每個組合：

a = arange(-180,181)  
b = arange(-82,83) 
totalgrid = [ [ 0 for i in range(len(b)) ] for j in range(len(a)) ] 
for d1 in range(len(a)): 
    for d2 in range(len(b)): 
     totalgrid[d1][d2]=np.where((x==a[d1])&(y==b[d2]))[0] 

此方法將失敗，並只返回一個空列表清單陣列。我無法弄清楚爲什麼它不能正常工作。

totalgrid[d1][d2]=np.where((x==a[0])&(y==b[0]))[0] 

返回所有發現指數從經度/緯度的是存在於-180°/ -82°： 通過更換最後一行。不幸的是它需要一段時間。我在某處丟失了一個for循環嗎？！

更詳細的問題： @askewchan 不幸的是這一個沒有解決我原來的問題。如預期的那樣，結果很好地代表了地面軌道。除了我需要每個網格點的總點數的事實之外，我還需要緯度/經度數組中的每個單獨的緯度/經度組合索引用於進一步的計算。 假設我有一個數組

lat(100000L,), lon(100000L,) and a third one array(100000L,) 

其對應於在每個點處測量。如果滿足條件，我需要經緯度的所有1°/ 1°組合的每個索引來檢查數組中的這個索引（100000L，）。現在讓我們假設lat/lon的指數[10000,10001,10002，...，10025]在同一個網格點上。對於那些索引，我需要檢查array [10000,10001,10002，..，10025]現在是否滿足條件，即np.where（array == 0）。用cts.nonzero（）我只能得到直方圖中的索引。但是，那些對直方圖值有貢獻的每個點的所有信息都會丟失。希望你能得到我最初的問題。

Answer 1

不確定我是否理解這裏的目標，但是你想要計算每個1°部分有多少個經緯度對？這是一個直方圖的作用：

lon = np.random.random(5000)*2*180 - 180 
lat = np.random.random(5000)*2*82 - 82 

a = np.arange(-180,181) 
b = np.arange(-82,83) 

np.histogram2d(lon, lat, (a,b)) 
#(array([[ 0., 0., 1., ..., 0., 0., 0.], 
#  [ 0., 2., 0., ..., 0., 0., 1.], 
#  [ 0., 0., 0., ..., 0., 1., 0.], 
#  ..., 
#  [ 0., 1., 0., ..., 0., 0., 0.], 
#  [ 0., 0., 0., ..., 0., 0., 0.], 
#  [ 0., 0., 0., ..., 0., 0., 0.]]), 

，你有一個非零計數的指標是在：

cts.nonzero() 
#(array([ 0, 0, 0, ..., 359, 359, 359]), 
# array([ 2, 23, 25, ..., 126, 140, 155])) 

您可以繪製太：

cts, xs, ys = np.histogram2d(lon, lat, (a,b)) 

pyplot.imshow(cts, extent=(-82,82,-180,180))