for文の処理について

Question

前提・実現したいこと

輪郭データから輪郭の重心を出し，重心を原点として輪郭の座標とのなす角θを出すプログラムを作成しようとしています．
あらかじめ出してある，輪郭contours1[i]の重心座標x_CoG[i]，y_CoG[i]があり，そちらを用いて以下のプログラムでθを求めようとしています．
θを出すために，重心から輪郭座標までのx方向距離x_d，y方向距離y_dを求め，それらをアークタンジェントを求める関数に入れています．x_dのデータはx_d_listに収納され，アークタンジェントのデータはthetaに収納されていて，それらの要素数は一致するはずなのですが，一致せず以下のような結果になってしまいます．
for文処理を間違えていると思うのですが，改善点が分かりません．

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

print(len(theta))の結果
2932

print(len(x_d_list))の結果
1466

該当のソースコード

x_d_list = []
y_d_list = []
theta = []
for i in range(0, len(contours1)):
    buf_np = contours1[i].flatten()# numpyの多重配列になっているため、一旦展開する。
    for j, elem in enumerate(buf_np):#jはcontours[i]内の座標について，θを出すループ
        if j%2==0:#x,yに分ける
            x_d = elem - x_CoG1[i]#x_d=重心から輪郭座標までのx方向距離
            x_d_list.append(x_d)
        else:
            y_d = elem - y_CoG1[i]
            y_d_list.append(y_d)
        theta.append(math.atan2(y_d, x_d))#thetaθ=atan2(y,x)
        
print(len(theta)) 
print(len(x_d_list))

試したこと

ここに問題に対して試したことを記載してください。

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

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

Answer 1

if文の処理が間違っています。
このままだと、x_dの処理と、y_dの処理、それぞれでtheta.appendするロジックになっています。thetaのリスト長が倍になってしまいます。
theta.appendにインデントをつけて、y_dの処理の時だけtheta.appendが実行されるようにすると、うまく動きます。

Comments

[mk_taro5]

回答ありがとうございます．初歩的なミスでした（thetaの要素数がx_dの要素数の倍なのになぜ気づけない）．
ご指摘通りの修正でうまく動きました，ありがとうございます．