Blender Pythonでオブジェクトを他のオブジェクトの任意の辺と同じ角度にしたい

Question

Blender、python上での3Dオブジェクトの操作についての質問です。

キューブ(オレンジ)のXZ軸を逆T字のフレーム(赤)に合わせたいです

逆T字は、変形するオブジェクトからガイド用に切り出した物のなので
明確に90度ではありません

X軸を合わせる所まではできたのですが
Z軸をあわせる方法がどうしても分かりませんでした
図中水色の点が実現したい事になります

合わせ具合は厳密ではなくても良いです

詳しい方、どうか教えて頂けないでしょうか
よろしくお願いいたします

以下は上記画像のコードです

import bpy
import bmesh
from mathutils import *

#逆T図形
obj_target=bpy.context.scene.collection.objects["target_line.004"]

"""逆T図形の頂点番号は正面(Y方向)から見るとこんな感じ
           0
           │
        1──3──2
"""
    
bm = bmesh.new()
bm.from_mesh(obj_target.data)

vert_local = [
    Vector((v.co[0], v.co[1], v.co[2]))
    for v in bm.verts
]

#ワールド座標への変換
vert_global = [
    obj_target.matrix_world @ v
    for v in vert_local
]

#逆T図形(水平)のベクトル
v1 = (vert_global[1]) - (vert_global[2])
#逆T図形(垂直)のベクトル(未使用)
v2 = (vert_global[0]) - (vert_global[3])

#シーン内から回転するオブジェクトを探す
obj_to=bpy.context.scene.collection.objects["Cube_to.004"]

#回転を適用
obj_to.rotation_mode = 'QUATERNION'
obj_to.rotation_quaternion = v1.to_track_quat('-X',"Z")

Answer 1

自己解決です

コードに下段部分を追加し対応

２つのオブジェクトのベクトルのなす角度を求め　その角度を一方に適用し
角度が一定以下に合うまでループする力技になりました(最大5回ぐらいで収束)

import bpy
import bmesh
from mathutils import *

#逆T図形
obj_target=bpy.context.scene.objects["target_line.004"]

"""逆T図形の頂点番号は正面(Y方向)から見るとこんな感じ
           0
           │
        1──3──2
"""
    
bm = bmesh.new()
bm.from_mesh(obj_target.data)

vert_local = [
    Vector((v.co[0], v.co[1], v.co[2]))
    for v in bm.verts
]

#ワールド座標への変換
vert_global = [
    obj_target.matrix_world @ v
    for v in vert_local
]

#逆T図形(水平)のベクトル
v1 = (vert_global[1]) - (vert_global[2])
#逆T図形(垂直)のベクトル(未使用)
v2 = (vert_global[0]) - (vert_global[3])

#シーン内から回転するオブジェクトを探す
obj_to=bpy.context.scene.objects["Cube_to.004"]

#回転を適用
obj_to.rotation_mode = 'QUATERNION'
obj_to.rotation_quaternion = v1.to_track_quat('-X',"Z")

#回転を適用
obj_to.rotation_mode = 'QUATERNION'
obj_to.rotation_quaternion = v1.to_track_quat('-X',"Z")

#(New)Z軸がT(垂直)に近くなるまで繰り返す
for i in range(100):
    
    print("Z軸合わせ"+str(i)+"回目")
    
    #この時点の回転結果が反映されてな
    
    #これやるBlender上の回転がデータに
    #オブジェクトをすべて非選択
    bpy.ops.object.select_all(action='DESELECT')
    #キューブを選択
    obj_to.select_set(True)
    #エディットモードにして戻す
    #bpy.ops.object.mode_set(mode='EDIT')
    bpy.ops.object.mode_set(mode='OBJECT')

    #////////////////////////////////
    #キューブのZ軸のベクトルを取得
    #Is there a simpler way to get the Vectors of a local axis?
    #https://blender.stackexchange.com/questions/156170/is-there-a-simpler-way-to-get-the-vectors-of-a-local-axis
    (translation, rotation, scale) = obj_to.matrix_world.decompose()
    local_z_axis = Vector((0.0, 0.0, 1.0))
    local_z_axis_global_coords = rotation @ local_z_axis 
    print("z axis:"+str(local_z_axis_global_coords))
    
    #////////////////////////////////
    #T字垂直と　キューブZ軸
    #2つのベクトルのなす角を求める
    vec_a=v2
    vec_c=local_z_axis_global_coords
    
    #【Python】座標上の3点からなる角度を計算 - Qiita
    #https://qiita.com/hacchi_/items/7e6f433d465df9378d7a
    # コサインの計算
    length_vec_a = np.linalg.norm(vec_a)
    length_vec_c = np.linalg.norm(vec_c)
    inner_product = np.inner(vec_a, vec_c)
    cos = inner_product / (length_vec_a * length_vec_c)

    # 角度（ラジアン）の計算
    rad = np.arccos(cos)

    # 弧度法から度数法（rad ➔ 度）への変換
    degree = np.rad2deg(rad)

    print("degree:"+str(degree))
     
    #一定以下なら終了
    if degree<1.0:
        print("Z軸合わせ完了")
        break
    
    #回転実施
    obj_to.rotation_mode = 'XYZ'
    obj_to.rotation_euler.rotate_axis("X", radians(degree))