先稍微整理一下
※ 引述《MAGICXX (逢甲阿法)》之銘言：
: 大家好
: 1.
: 我現在有一組程式碼有問題
: upperb=NDVI_LST[col[0]]>0.0
: lowerb=NDVI_LST[col[0]]<0.1
: edge_1=NDVI_LST[(upperb & lowerb)]
col_0 = col[0]
col_1 = col[1]
upperb=NDVI_LST[col_0]>0.0
lowerb=NDVI_LST[col_0]<0.1
edge_1=NDVI_LST[(upperb & lowerb)]
: lst_values=[]
: for i in edge_1[col[0]]:
: LST=NDVI_LST.loc[NDVI_LST[col[0]]==i,col[1]]
: NDVI=i # 這行不知道在幹嘛
: for j in LST.values:
: lst_values.append(j)
: df.loc[0.0,'w1月']=max(lst_values)
# 這段看不懂在幹嘛
max_lst_values = edge_1[col_1].max()
min_lst_values = edge_1[col_1].min()
: lst_values=[]
: for i in edge_1[col[0]]:
: LST=NDVI_LST.loc[NDVI_LST[col[0]]==i,col[1]]
: NDVI=i
: for j in LST.values:
: lst_values.append(j)
: df.loc[0.0,'d1月']=min(lst_values)
: 因為我現在會需要有12個月的重複計算以及
: upperb以及lowerb的重複計算(0.0 ,0.1, 0.2, 0.3.....1.0)
: 原本我是用下面的code跑 可是每次都會剩下最後一個值然後把前面的值蓋掉
: NDVI_bound=np.around(np.linspace(0,1,11).tolist(),decimals=2)
: for k in range(1,13):
: for j in range(23):
: for i in range(10):
: upperb=NDVI_LST[col[int(j)]]>NDVI_bound[int(i)]
: lowerb=NDVI_LST[col[int(j)]]<NDVI_bound[int(i)+1]
: edge_1=NDVI_LST[(upperb & lowerb)]
: lst_values=[]
: for l in edge_1[col[int(j)]]:
: LST=NDVI_LST.loc[NDVI_LST[col[int(j)]]==l,col[int(j)+1]]
: NDVI=i
: for m in LST.values:
: lst_values.append(m)
: if lst_values==[]:
: df.loc[NDVI_bound[int(i)+1],'w'+str(k)+'月']=0
: df.loc[NDVI_bound[int(i)+1],'d'+str(k)+'月']=0
: else:
: df.loc[NDVI_bound[int(i)+1],'w'+str(k)+'月']=max(lst_values)
: df.loc[NDVI_bound[int(i)+1],'d'+str(k)+'月']=min(lst_values)
: 我想問一下大家 這樣要怎麼縮減這個程式碼才會是對的?
: 不然我光是用一個一個貼 就要2000多行程式碼 可是明明是重複做一樣的動作...
這邊完全沒有對K做操作 請檢查一下code跟邏輯
: 2.查表
: 我找的值是從上面那張表出來的
: D1 W1
: 0.1 23.7 22.4
: 0.2 25.8 26.88
: 0.3 29.8 30.64
: .
: .
: .
: 1.0 35.26 38.31
: 我現在給他一個值 他要怎麼找到在這張表的哪個範圍?
: 比如說 我現在給他27.33這個數值 他要怎麼找到是在0.2~0.3之間
: 然後跟我說上限是29.8 下限是25.8 這樣?
: 有兩個問題 再麻煩各位了QQ
binary search