"""Render a multi-joint URDF at a pose, auto-framed, from several angles. [RECONSTRUCTED 2026-07-06 —
same interface as the lab original] Usage: MUJOCO_GL=egl render_arm.py <urdf> <out.png> [q]"""
import sys
import xml.etree.ElementTree as ET
import mujoco
import numpy as np
from PIL import Image

urdf, out = sys.argv[1], sys.argv[2]
q = float(sys.argv[3]) if len(sys.argv) > 3 else 0.0
root = ET.parse(urdf).getroot()
for l in root.findall("link"):
    if l.find("inertial") is None:
        i = ET.SubElement(l, "inertial")
        ET.SubElement(i, "mass", {"value": "0.1"})
        ET.SubElement(i, "inertia", dict(ixx="1e-4", ixy="0", ixz="0", iyy="1e-4", iyz="0", izz="1e-4"))
mj = ET.Element("mujoco")
c = ET.SubElement(mj, "compiler")
c.set("discardvisual", "false")
c.set("meshdir", ".")
root.insert(0, mj)
m = mujoco.MjModel.from_xml_string(ET.tostring(root, encoding="unicode"))
d = mujoco.MjData(m)
if m.nq:
    d.qpos[:] = q
mujoco.mj_forward(m, d)
print("nq", m.nq, "nbody", m.nbody, "extent", round(float(m.stat.extent), 3))
ims = []
for az, el in [(45, -20), (135, -18), (90, -89), (0, -5)]:
    cam = mujoco.MjvCamera()
    cam.lookat[:] = m.stat.center
    cam.distance = 1.6 * m.stat.extent
    cam.azimuth, cam.elevation = az, el
    with mujoco.Renderer(m, 470, 630) as r:
        r.update_scene(d, cam)
        ims.append(r.render())
Image.fromarray(np.hstack(ims)).save(out)
print("rendered", out)