THUMS RBF Morphing Pipeline (Standing ↔ Seating)

Radial Basis Function (RBF) morphing pipeline for posture transfer of the THUMS human body model.

This repository provides a reproducible workflow to transform a standing THUMS mesh into seated or alternate postures while preserving finite element mesh quality and anatomical consistency. The method uses landmark-based control points and PyGeM RBF interpolation to generate full-body LS-DYNA keyword meshes.

This pipeline was developed for computational biomechanics and spinal injury risk simulations under high-G loading conditions.

---

Overview

The workflow:

1. Select anatomical landmarks on source (THUMS) and target geometry
2. Generate RBF control point parameter file (.prm)
3. Compute deformation using PyGeM RBF
4. Apply deformation to full mesh
5. Export:
   • morphed .k (LS-DYNA keyword)
   • .vtu files for visualization

---

Percentile setting: Each script contains a percentile = ... variable near the top.
Update this value in each file you run (e.g., 50 / 80 / 98) so the input/output paths match your target model folder.

Percentile

The percentile is configurable in each script via:

percentile = 98

Scripts:

    Standing posture morphing: 

        File: `creating_RBF_standing.py`
        - Creates RBF parameter file using selected landmarks
        - Outputs `.prm` and visualization `.vtu`

        File: `RBF_Standing.py`
        - Applies RBF to full THUMS mesh
        - Exports:
            Full Body model: `.k`

    Seating posture morph
        File: `creating_RBF_seating.py`
        - Builds control points for seated configuration
        - Writes seating `.prm` and landmark `.vtu`

        File: `RBF_Seating.py`
        - Computes displacement field
        - Applies RBF morph to full mesh
        - Exports:
            - Visualizations:`.vtu`
            - Full Body model: `.k`