2つのデータフレームで繰り返し計算した値を別のデータフレームに入れたい

Question

前提・実現したいこと

df_theta_radに22行1列のデータが、df_theta_rad_Tに1行22列のデータが入っています。
これら2つのデータの差の絶対値を繰り返しでdf_thetaijに入れたいのですが、以下のエラーが発生しました。

発生している問題・エラーメッセージ

raise IndexError("single positional indexer is out-of-bounds")

IndexError: single positional indexer is out-of-bounds

該当のソースコード

python

import numpy as np
pi=np.pi
import pandas as pd

df_read_North=pd.read_csv("Uraga_GateCrossingList_20180101NS.csv",encoding="shift-jis")
df_read_South=pd.read_csv("Uraga_GateCrossingList_20180101NS.csv",encoding="shift-jis")
df_read_West=pd.read_csv("Uraga_GateCrossingList_20180101WE.csv",encoding="shift-jis")
df_read_East=pd.read_csv("Uraga_GateCrossingList_20180101WE.csv",encoding="shift-jis")

df_data_North=df_read_North[(df_read_North["i"]==0) & (df_read_North["j"]==0)]
df_data_South=df_read_South[(df_read_South["i"]==1) & (df_read_South["j"]==0)]
df_data_West=df_read_West[(df_read_West["i"]==0) & (df_read_West["j"]==0)]
df_data_East=df_read_East[(df_read_East["i"]==1) & (df_read_East["j"]==0)]

df_data_NS=pd.merge(df_data_North,df_data_South,how="outer")
df_data_WE=pd.merge(df_data_West,df_data_East,how="outer")
df_data_merge=pd.merge(df_data_NS,df_data_WE,how="outer")

df_theta=df_data_merge.groupby('Group').mean()['θ']
df_theta_rad=df_theta*(pi/180)

df_theta_rad=pd.DataFrame(data=df_theta_rad)
df_theta_rad=df_theta_rad.rename(columns={'Sec':'θ'})

df_theta_rad_T=df_theta_rad.T

a=0
b=0
LastNumber=len(df_theta_rad)

df_thetaij=pd.DataFrame()

for a in range(LastNumber):
for b in range(LastNumber):
df_thetaij.iloc[a,b]=np.abs(df_theta_rad.iloc[a,]-df_theta_rad_T.iloc[1,b])
b=b+1
if b==LastNumber:
a=a+1

with pd.ExcelWriter('Tokyo20181.xlsx')as writer:
df_thetaij.to_excel(writer,'θij')

### 試したこと

for文の中のdf_thetaij.iloc[a,b]=np.abs(df_theta_rad.iloc[a,]-df_theta_rad_T.iloc[1,b])でdf_theta_rad_T.iloc[1,b]の1をなくしてみたり、：にしたりしたのですが、その際はtoo many indexerというエラーが発生し、困っています。for文にも少し自信がないのですが、おかしな点があったら教えて頂きたいです。

### 補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

Answer 1

a=df_theta_rad.values.ravel()
b=np.meshgrid(a,a)
c=np.abs(b[0]-b[1])

Comments

[ryosuke197]

こんな方法があるなんて知らなかったです！解決しました！ありがとうございます！