Computational modeling of skin: Using stress profiles as predictor for tissue necrosis in reconstructive surgery

Exploring uncertainty in hyper-viscoelastic properties of scalp skin through patient-specific finite element models for reconstructive surgery

Understanding skin responses to external forces is crucial for post-cutaneous flap wound healing. However, the in vivo viscoelastic behavior of scalp skin remains poorly understood. Personalized virtual surgery simulations offer a way to study tissue responses in relevant 3D geometries. Yet, anticipating wound risk remains challenging due to limited data on skin viscoelasticity, which hinders our ability to determine the interplay between wound size and stress levels. To bridge this gap, we reexamine three clinical cases involving scalp reconstruction using patient-specific geometric models and employ uncertainty quantification through a Monte Carlo simulation approach to study the effect of skin viscoelasticity on the final stress levels from reconstructive surgery. Utilizing the generalized Maxwell model via the Prony series, we can parameterize and efficiently sample a realistic range of viscoelastic response and thus shed light on the influence of viscoelastic material uncertainty in surgical scenarios. Our analysis identifies regions at risk of wound complications based on reported threshold stress values from the literature and highlights the significance of focusing on long-term responses rather than short-term ones.

Sensory neuron activation from topical treatments modulates the sensorial perception of human skin

Neural signaling of skin sensory perception from topical treatments is often reported in subjective terms such as a sensation of skin "tightness" after using a cleanser or "softness" after applying a moisturizer. However, the mechanism whereby cutaneous mechanoreceptors and corresponding sensory neurons are activated giving rise to these perceptions has not been established. Here, we provide a quantitative approach that couples in vitro biomechanical testing and detailed computational neural stimulation modeling along with a comprehensive in vivo self-assessment survey to demonstrate how cutaneous biomechanical changes in response to treatments are involved in the sensorial perception of the human skin. Strong correlations are identified between reported perception up to 12 hours post treatment and changes in the computed neural stimulation from mechanoreceptors residing deep under the skin surface. The study reveals a quantitative framework for understanding the biomechanical neural activation mechanism and the subjective perception by individuals.

The multiscale characterization and constitutive modeling of healthy and type 2 diabetes mellitus Sprague Dawley rat skin

In type 2 diabetes mellitus (T2DM), elevated glucose level impairs the biochemistry of the skin which may result in alteration of its mechanical and structural properties. The several aspects of structural and mechanical changes in skin due to T2DM remain poorly understood. To fill these research gaps, we developed a non-obese T2DM rat (Sprague Dawley (SD)) model for investigating the effect of T2DM on the in vivo strain stress state, mechanical and structural properties of skin. In vivo strain and mechanical anisotropy of healthy and T2DM skin were measured using the digital imaging correlation (DIC) technique and DIC coupled bulge experiment, respectively. Fluorescence microscopy and histology were used to assess the collagen and elastin fibers microstructure whereas nanoscale structure was captured through atomic force microscopy (AFM). Based on the microstructural observations, skin was modeled as a multilayer membrane where in and out of plane distribution of collagen fibers and planar distribution of elastin fibers were cast in constitutive model. Further, the state of in vivo stresses of healthy and T2DM were measured using model parameters and in vivo strain in the constitutive model. The results showed that T2DM causes significant loss in in vivo stresses (p < 0.01) and increase in anisotropy (p < 0.001) of skin. These changes were found in good correlation with T2DM associated alteration in skin microstructure. Statistical analysis emphasized that increase in blood glucose concentration (HbA1c) was the main cause of impaired biomechanical properties of skin. The presented data in this study can help to understand the skin pathology and to simulate the skin related clinical procedures. STATEMENT OF SIGNIFICANCE: Our study is significant as it presents findings related to the effect of T2DM on the physiologic stress strain, structural and mechanical response of SD rat skin. In this study, we developed a non-obese T2DM SD rat model which mimics the phenotype of Asian type 2 diabetics (non-obese). Several structural and mechanical characterization techniques were explored for multiscale characterization of healthy and T2DM skin. Further, based on microstructural information, we presented the constitutive models that incorporate the real microstructure of skin. The presented results can be helpful to simulate the realistic mechanical response of skin during various clinical trials.

Geometrical model establishment and preoperative evaluation on A-T flap design: Finite element method-based computer-aided simulation on surgical operation processes

Objective

A-T flap has been extensively applied to repair dermal soft tissue defects. The flap design completely depends on the experience of doctors. Herein, we explored the approach of analyzing the reasonability of A-T flap design and performed a simulation of operation processes by computer-aided technology. Afterward, the finite element analysis software (MSC.Marc/Mentat) was used to establish the simulation model, based on which the computer simulation of flap suturing and release state in A-T flap surgery was performed.

Methods

A geometrical model of the A-T flap was established, and the length-width ratio of the flap, maximum suture distance, and suture area that could influence the postoperative suture effects of the flap were analyzed. The reasonable surgical planning for A-T flap design based on the crossing constraint relationship was achieved. The simulation model was established by the finite element analysis software (MSC.Marc/Mentat), based on which computer simulation of flap suture and release state of A-T flap in surgery processes were performed. The flap's stress and deformation distribution results confirmed the applicability of the A-T flap design method proposed in the present study.

Results

When the apex angle of the A-T flap was 60 degrees, the suture area was the smallest, and the flap design had the highest practicability.

Conclusion

Computer-assisted preoperative assessment, which has high clinical value, could provide a theoretical basis for A-T flap design in clinical practice.

Advances in non-invasive biosensing measures to monitor wound healing progression

Impaired wound healing is a significant financial and medical burden. The synthesis and deposition of extracellular matrix (ECM) in a new wound is a dynamic process that is constantly changing and adapting to the biochemical and biomechanical signaling from the extracellular microenvironments of the wound. This drives either a regenerative or fibrotic and scar-forming healing outcome. Disruptions in ECM deposition, structure, and composition lead to impaired healing in diseased states, such as in diabetes. Valid measures of the principal determinants of successful ECM deposition and wound healing include lack of bacterial contamination, good tissue perfusion, and reduced mechanical injury and strain. These measures are used by wound-care providers to intervene upon the healing wound to steer healing toward a more functional phenotype with improved structural integrity and healing outcomes and to prevent adverse wound developments. In this review, we discuss bioengineering advances in 1) non-invasive detection of biologic and physiologic factors of the healing wound, 2) visualizing and modeling the ECM, and 3) computational tools that efficiently evaluate the complex data acquired from the wounds based on basic science, preclinical, translational and clinical studies, that would allow us to prognosticate healing outcomes and intervene effectively. We focus on bioelectronics and biologic interfaces of the sensors and actuators for real time biosensing and actuation of the tissues. We also discuss high-resolution, advanced imaging techniques, which go beyond traditional confocal and fluorescence microscopy to visualize microscopic details of the composition of the wound matrix, linearity of collagen, and live tracking of components within the wound microenvironment. Computational modeling of the wound matrix, including partial differential equation datasets as well as machine learning models that can serve as powerful tools for physicians to guide their decision-making process are discussed.

Multiscale Computational Model Predicts Mouse Skin Kinematics Under Tensile Loading

Skin is a complex tissue whose biomechanical properties are generally understood in terms of an incompressible material whose microstructure undergoes affine deformations. A growing number of experiments, however, have demonstrated that skin has a high Poisson's ratio, substantially decreases in volume during uniaxial tensile loading, and demonstrates collagen fiber kinematics that are not affine with local deformation. In order to better understand the mechanical basis for these properties, we constructed multiscale mechanical models (MSM) of mouse skin based on microstructural multiphoton microscopy imaging of the dermal microstructure acquired during mechanical testing. Three models that spanned the cases of highly aligned, moderately aligned, and nearly random fiber networks were examined and compared to the data acquired from uniaxially stretched skin. Our results demonstrate that MSMs consisting of networks of matched fiber organization can predict the biomechanical behavior of mouse skin, including the large decrease in tissue volume and nonaffine fiber kinematics observed under uniaxial tension.

Improving reconstructive surgery design using Gaussian process surrogates to capture material behavior uncertainty

To produce functional, aesthetically natural results, reconstructive surgeries must be planned to minimize stress as excessive loads near wounds have been shown to produce pathological scarring and other complications (Gurtner et al., 2011). Presently, stress cannot easily be measured in the operating room. Consequently, surgeons rely on intuition and experience (Paul et al., 2016; Buchanan et al., 2016). Predictive computational tools are ideal candidates for surgery planning. Finite element (FE) simulations have shown promise in predicting stress fields on large skin patches and in complex cases, helping to identify potential regions of complication. Unfortunately, these simulations are computationally expensive and deterministic (Lee et al., 2018a). However, running a few, well selected FE simulations allows us to create Gaussian process (GP) surrogate models of local cutaneous flaps that are computationally efficient and able to predict stress and strain for arbitrary material parameters. Here, we create GP surrogates for the advancement, rotation, and transposition flaps. We then use the predictive capability of these surrogates to perform a global sensitivity analysis, ultimately showing that fiber direction has the most significant impact on strain field variations. We then perform an optimization to determine the optimal fiber direction for each flap for three different objectives driven by clinical guidelines (Leedy et al., 2005; Rohrer and Bhatia, 2005). While material properties are not controlled by the surgeon and are actually a source of uncertainty, the surgeon can in fact control the orientation of the flap with respect to the skin's relaxed tension lines, which are associated with the underlying fiber orientation (Borges, 1984). Therefore, fiber direction is the only material parameter that can be optimized clinically. The optimization task relies on the efficiency of the GP surrogates to calculate the expected cost of different strategies when the uncertainty of other material parameters is included. We propose optimal flap orientations for the three cost functions and that can help in reducing stress resulting from the surgery and ultimately reduce complications associated with excessive mechanical loading near wounds.

Anticipating local flaps closed-form solution on 3D face models using finite element method

A realistic three-dimensional (3D) computational model of skin flap closures using Asian-like head templates from two different genders, male and female, has been developed. The current study aimed to understand the biomechanics of the local flap designs along with the effect of wound closures on the respective genders. Two Asian head templates from opposite genders were obtained to use as base models. A third-order Yeoh hyperelastic model was adapted to characterize as skin material properties. A single layer composed of combined epidermis and dermis was considered, and the models were thickened according to respective anatomical positions. Each model gender was excised with a fixed defect size which was consequently covered by three different local flap designs, namely advancement, rotation, and rhomboid flaps. Post-operative simulation presented various scenarios of skin flap closures. Rotation and rhomboid flaps demonstrated maximal tension at the apex of the flap for both genders as well as advancement flap in the female face model. However, advancement flap closure in the male face model was presented otherwise. Yet, the deformation patterns and the peak tension of the discussed flaps were consistent with conventional local flap surgery. Moreover, male face models generated higher stresses compared to the female face models with a 70.34% mean difference. Overall, the skin flap operations were executed manually, and the designed surgery model met the objectives successfully while acknowledging the study limitations. NOVELTY FILE: 3D head templates were considered to address the gap as 3D face models were uncommonly employed in understanding the biomechanics of the local flaps realistically. Most of the existing studies focus on the 2D and 3D planar geometry in their models. As gender comparison has yet to be addressed, we intended to fill this gap by exploring the stress contours of the local flap designs in different genders. Create a 3D face model from two opposite genders which is capable of simulating closure of wounds using local flaps with a focus on advancement, rotation, and rhomboid flaps.

Personalized Computational Models of Tissue-Rearrangement in the Scalp Predict the Mechanical Stress Signature of Rotation Flaps

OBJECTIVE

To elucidate the mechanics of scalp rotation flaps through 3D imaging and computational modeling. Excessive tension near a wound or sutured region can delay wound healing or trigger complications. Measuring tension in the operating room is challenging, instead, noninvasive methods to improve surgical planning are needed.

DESIGN

Multi-view stereo allows creation of 3D patient-specific geometries based on a set of photographs. The patient-specific 3D geometry is imported into a finite element (FE) platform to perform a virtual procedure. The simulation is compared with the clinical outcome. Additional simulations quantify the effect of individual flap parameters on the resulting tension distribution.

PARTICIPANTS

Rotation flaps for reconstruction of scalp defects following melanoma resection in 2 cases are presented. Rotation flaps were designed without preoperative FE preparation.

MAIN OUTCOME MEASURE

Tension distribution over the operated region.

RESULTS

The tension from FE shows peaks at the base and distal ends of the scalp rotation flap. The predicted geometry from the simulation aligns with postoperative photographs. Simulations exploring the flap design parameters show variation in the tension. Lower tensions were achieved when rotation was oriented with respect to skin tension lines (horizontal tissue fibers) and smaller rotation angles.

CONCLUSIONS

Tension distribution following rotation of scalp flaps can be predicted through personalized FE simulations. Flaps can be designed to reduce tension using FE, which may greatly improve the reliability of scalp reconstruction in craniofacial surgery, critical in complex cases when scalp reconstruction is essential for coverage of hardware, implants, and/or bone graft.

Preoperative Evaluation of V-Y Flap Design Based on Computer-Aided Analysis

V-Y flap is widely used in plastic surgery as an important technique for reconstructing deformities and improving appearance. In this paper, a geometrical parameter model and finite element analysis were used to study the rationale of the proposed V-Y flap design and the preoperative evaluation of the V-Y flap design. First, a geometric parameter model of the V-Y flap was established to analyze the five key geometric relationships affecting the flap structure and obtain a reasonable plan for the V-Y flap design through the crossing constraint relationship. Second, in order to verify the effectiveness of the V-Y flap design, the suture and release states of the V-Y flap during surgery were evaluated based on a simulation model of the V-Y flap generated by finite element analysis software. The results revealed that the approach proposed in this paper provides a feasible method for clinical V-Y flap design.

On the compressibility and poroelasticity of human and murine skin

A total of 37 human and 33 murine skin samples were subjected to uniaxial monotonic, cyclic, and relaxation experiments. Detailed analysis of the three-dimensional kinematic response showed that skin volume is significantly reduced as a consequence of a tensile elongation. This behavior is most pronounced in monotonic but persists in cyclic tests. The dehydration associated with volume loss depends on the osmolarity of the environment, so that tension relaxation changes as a consequence of modifying the ionic strength of the environmental bath. Similar to ex vivo observations, complementary in vivo stretching experiments on human volar forearms showed strong in-plane lateral contraction. A biphasic homogenized model is proposed which allows representing all relevant features of the observed mechanical response.

Multi-view stereo in the operating room allows prediction of healing complications in a patient-specific model of reconstructive surgery

Excessive mechanical stress leads to wound healing complications following reconstructive surgery. However, this knowledge is not easily applicable in clinical scenarios due to the difficulty in measuring stress contours during complex tissue rearrangement procedures. Computational tools have been proposed as an alternative to address this need, but obtaining patient specific geometries with an affordable and flexible setup has remained a challenge. Here we present a methodology to determine the stress contours from a reconstructive procedure on a patient-specific finite element model based on multi-view stereo (MVS). MVS is a noninvasive technology that allows reconstruction of 3D geometries using a standard digital camera, making it ideal for the operating room. Finite element analysis can then be used on the patient-specific geometry to perform a virtual surgery and predict regions at risk of complications. We applied our approach to the case of a 7-year-old patient who was treated to correct a cranial contour deformity and resect two large areas of scalp scarring. The simulation showed several zones of high stress, particularly near the suture lines at the distal ends of the flaps. The patient did show delayed healing and partial flap tip necrosis at one of such predicted regions at the 30-day follow up visit. Our results further establish the application of computational tools in individualized medical scenarios to advance preoperative planing and anticipate regions of concern immediately after surgery.

Mathematical and computational modelling of skin biophysics: a review

The objective of this paper is to provide a review on some aspects of the mathematical and computational modelling of skin biophysics, with special focus on constitutive theories based on nonlinear continuum mechanics from elasticity, through anelasticity, including growth, to thermoelasticity. Microstructural and phenomenological approaches combining imaging techniques are also discussed. Finally, recent research applications on skin wrinkles will be presented to highlight the potential of physics-based modelling of skin in tackling global challenges such as ageing of the population and the associated skin degradation, diseases and traumas.

An Invariant-Based Damage Model for Human and Animal Skins

Constitutive modelling of skins that account for damage effects is important to provide insight for various clinical applications, such as skin trauma and injury, artificial skin design, skin aging, disease diagnosis, surgery, as well as comparative studies of skin biomechanics between species. In this study, a new damage model for human and animal skins is proposed for the first time. The model is nonlinear, anisotropic, invariant-based, and is based on the Gasser-Ogden-Holzapfel constitutive law initially developed for arteries. Taking account of the mean collagen fibre orientation and its dispersion, the new model can describe a wide range of skins with damage. The model is first tested on the uniaxial test data of human skin and then applied to nine groups of uniaxial test data for the human, swine, rabbit, bovine and rhino skins. The material parameters can be inversely estimated based on uniaxial tests using the optimization method in MATLAB with a root mean square error ranged between 2.15% and 12.18%. A sensitivity study confirms that the fibre orientation dispersion and the mean fibre angle are among the most important factors that influence the behaviour of the damage model. In addition, these two parameters can only be reliably estimated if some histological information is provided. We also found that depending on the location of skins, the tissue damage may be brittle controlled by the fibre breaking limit (i.e., when the fibre stretch is greater than 1.13-1.32, depending on the species), or ductile (due to both the fibre and the matrix damages). The brittle damages seem to occur mostly in the back, and the ductile damages are seen from samples taken from the belly. The proposed constitutive model may be applied to various clinical applications that require knowledge of the mechanical response of human and animal skins.

The Incompatibility of Living Systems: Characterizing Growth-Induced Incompatibilities in Expanded Skin

Skin expansion is a common surgical technique to correct large cutaneous defects. Selecting a successful expansion protocol is solely based on the experience and personal preference of the operating surgeon. Skin expansion could be improved by predictive computational simulations. Towards this goal, we model skin expansion using the continuum framework of finite growth. This approach crucially relies on the concept of incompatible configurations. However, aside from the classical opening angle experiment, our current understanding of growth-induced incompatibilities remains rather vague. Here we visualize and characterize incompatibilities in living systems using skin expansion in a porcine model: We implanted and inflated two expanders, crescent, and spherical, and filled them to 225 cc throughout a period of 21 days. To quantify the residual strains developed during this period, we excised the expanded skin patches and subdivided them into smaller pieces. Skin growth averaged 1.17 times the original area for the spherical and 1.10 for the crescent expander, and displayed significant regional variations. When subdivided into smaller pieces, the grown skin patches retracted heterogeneously and confirmed the existence of incompatibilities. Understanding skin growth through mechanical stretch will allow surgeons to improve-and ultimately personalize-preoperative treatment planning in plastic and reconstructive surgery.

Isogeometric Kirchhoff-Love shell formulations for biological membranes

Computational modeling of thin biological membranes can aid the design of better medical devices. Remarkable biological membranes include skin, alveoli, blood vessels, and heart valves. Isogeometric analysis is ideally suited for biological membranes since it inherently satisfies the C¹-requirement for Kirchhoff-Love kinematics. Yet, current isogeometric shell formulations are mainly focused on linear isotropic materials, while biological tissues are characterized by a nonlinear anisotropic stress-strain response. Here we present a thin shell formulation for thin biological membranes. We derive the equilibrium equations using curvilinear convective coordinates on NURBS tensor product surface patches. We linearize the weak form of the generic linear momentum balance without a particular choice of a constitutive law. We then incorporate the constitutive equations that have been designed specifically for collagenous tissues. We explore three common anisotropic material models: Mooney-Rivlin, May Newmann-Yin, and Gasser-Ogden-Holzapfel. Our work will allow scientists in biomechanics and mechanobiology to adopt the constitutive equations that have been developed for solid three-dimensional soft tissues within the framework of isogeometric thin shell analysis.

Isogeometric Kirchhoff-Love shell formulations for biological membranes

Computational modeling of thin biological membranes can aid the design of better medical devices. Remarkable biological membranes include skin, alveoli, blood vessels, and heart valves. Isogeometric analysis is ideally suited for biological membranes since it inherently satisfies the C¹-requirement for Kirchhoff-Love kinematics. Yet, current isogeometric shell formulations are mainly focused on linear isotropic materials, while biological tissues are characterized by a nonlinear anisotropic stress-strain response. Here we present a thin shell formulation for thin biological membranes. We derive the equilibrium equations using curvilinear convective coordinates on NURBS tensor product surface patches. We linearize the weak form of the generic linear momentum balance without a particular choice of a constitutive law. We then incorporate the constitutive equations that have been designed specifically for collagenous tissues. We explore three common anisotropic material models: Mooney-Rivlin, May Newmann-Yin, and Gasser-Ogden-Holzapfel. Our work will allow scientists in biomechanics and mechanobiology to adopt the constitutive equations that have been developed for solid three-dimensional soft tissues within the framework of isogeometric thin shell analysis.

A mechanistic insight into the mechanical role of the stratum corneum during stretching and compression of the skin

The study of skin biophysics has largely been driven by consumer goods, biomedical and cosmetic industries which aim to design products that efficiently interact with the skin and/or modify its biophysical properties for health or cosmetic benefits. The skin is a hierarchical biological structure featuring several layers with their own distinct geometry and mechanical properties. Up to now, no computational models of the skin have simultaneously accounted for these geometrical and material characteristics to study their complex biomechanical interactions under particular macroscopic deformation modes. The goal of this study was, therefore, to develop a robust methodology combining histological sections of human skin, image-processing and finite element techniques to address fundamental questions about skin mechanics and, more particularly, about how macroscopic strains are transmitted and modulated through the epidermis and dermis. The work hypothesis was that, as skin deforms under macroscopic loads, the stratum corneum does not experience significant strains but rather folds/unfolds during skin extension/compression. A sample of fresh human mid-back skin was processed for wax histology. Sections were stained and photographed by optical microscopy. The multiple images were stitched together to produce a larger region of interest and segmented to extract the geometry of the stratum corneum, viable epidermis and dermis. From the segmented structures a 2D finite element mesh of the skin composite model was created and geometrically non-linear plane-strain finite element analyses were conducted to study the sensitivity of the model to variations in mechanical properties. The hybrid experimental-computational methodology has offered valuable insights into the simulated mechanics of the skin, and that of the stratum corneum in particular, by providing qualitative and quantitative information on strain magnitude and distribution. Through a complex non-linear interplay, the geometry and mechanical characteristics of the skin layers (and their relative balance), play a critical role in conditioning the skin mechanical response to macroscopic in-plane compression and extension. Topographical features of the skin surface such as furrows were shown to act as an efficient means to deflect, convert and redistribute strain-and so stress-within the stratum corneum, viable epidermis and dermis. Strain reduction and amplification phenomena were also observed and quantified. Despite the small thickness of the stratum corneum, its Young׳s modulus has a significant effect not only on the strain magnitude and directions within the stratum corneum layer but also on those of the underlying layers. This effect is reflected in the deformed shape of the skin surface in simulated compression and extension and is intrinsically linked to the rather complex geometrical characteristics of each skin layer. Moreover, if the Young׳s modulus of the viable epidermis is assumed to be reduced by a factor 12, the area of skin folding is likely to increase under skin compression. These results should be considered in the light of published computational models of the skin which, up to now, have ignored these characteristics.

Multi-view stereo analysis reveals anisotropy of prestrain, deformation, and growth in living skin

Skin expansion delivers newly grown skin that maintains histological and mechanical features of the original tissue. Although it is the gold standard for cutaneous defect correction today, the underlying mechanisms remain poorly understood. Here we present a novel technique to quantify anisotropic prestrain, deformation, and growth in a porcine skin expansion model. Building on our recently proposed method, we combine two novel technologies, multi-view stereo and isogeometric analysis, to characterize skin kinematics: Upon explantation, a unit square retracts ex vivo to a square of average dimensions of [Formula: see text]. Upon expansion, the unit square deforms in vivo into a rectangle of average dimensions of [Formula: see text]. Deformations are larger parallel than perpendicular to the dorsal midline suggesting that skin responds anisotropically with smaller deformations along the skin tension lines. Upon expansion, the patch grows in vivo by [Formula: see text] with respect to the explanted, unexpanded state. Growth is larger parallel than perpendicular to the midline, suggesting that elevated stretch activates mechanotransduction pathways to stimulate tissue growth. The proposed method provides a powerful tool to characterize the kinematics of living skin. Our results shed light on the mechanobiology of skin and help us to better understand and optimize clinically relevant procedures in plastic and reconstructive surgery.

Characterization of living skin using multi-view stereo and isogeometric analysis

Skin is our interface with the outside world. In its natural environment, it displays unique mechanical characteristics, such as prestretch and growth. While there is a general agreement on the physiological importance of these features, they remain poorly characterized, mainly because they are difficult to access with standard laboratory techniques. Here we present a new, inexpensive technique to characterize living skin using multi-view stereo and isogeometric analysis. Based on easy-to-create hand-held camera images, we quantify prestretch, deformation and growth in a controlled porcine model of chronic skin expansion. Over a period of 5 weeks, we gradually inflate an implanted tissue expander, take weekly photographs of the experimental scene, reconstruct the geometry from a tattooed surface grid and create parametric representations of the skin surface. After 5 weeks of expansion, our method reveals an average area prestretch of 1.44, an average area stretch of 1.87 and an average area growth of 2.25. Area prestretch is maximal in the ventral region with a value of 2.37, whereas area stretch and area growth are maximal above the center of the expander, with values of 4.05 and 4.81, respectively. Our study has immediate impact on understanding living skin to optimize treatment planning and decision making in plastic and reconstructive surgery. Beyond these direct implications, our experimental design has broad applications in clinical research and basic sciences: it serves as a simple, robust, low cost, easy-to-use tool to reconstruct living membranes, which are difficult to characterize in a conventional laboratory setup.