グリッド上の選択点を元に曲面を補完するテスト

blender_interpolate_study.py

import bpy
import numpy as np
from scipy.interpolate import griddata

def func():
    bpy.ops.object.mode_set(mode='OBJECT')
    grid = bpy.data.objects['Grid']
    vs = grid.data.vertices
    
    points = []
    values = []
    for v in [v for v in vs if v.select]:
        points.append(v.co.xy)
        values.append(v.co.z)

    xi = np.array([v.co.xy for v in vs])
    zi = griddata(np.array(points), np.array(values),
        xi, method='cubic')

    for v,z in zip(vs, zi):
        v.co.z = z

func()