Reconstructs stains from cast-off event to reproduce the motion of cast-off.
Scott McCleary
Email: scott.thomas.mccleary@gmail.com | daniel.attinger@gmail.com
Phone: (515) 975-5544
Spatter Stains to Cast-off Reconstruction
Center for Statistics and Applications in Forensic Evidence
Department of Mechanical Engineering - Attinger Lab
Iowa State University
Last Updated 01/06/2024 🎉
- Surface Tangential Vector (t̂) :
- Definition : vector tangential to a surface aligned in the direction of gravity (**-**z-direction) for all non-horizonal surfaces
- for perfectly horizontal surfaces, e.g. upward surface (ceiling) and downward surface (floor), the -z-direction aligns with the Surface Normal Vector. In these instances, the Surface Tangential Vector (t̂) is defined in the **+**x-direction
- Definition : vector tangential to a surface aligned in the direction of gravity (**-**z-direction) for all non-horizonal surfaces
- Surface Normal Vector (n̂) :
- Definition : vector normal to a surface pointing outside the room
- Gamma (γ) : stain directional angle (relative to impacted surface):
- Definition : stain glancing angle between the projected trajectory vector onto the impacted surface and the tangent vector to the surface
- Alpha (α) : stain impact angle
- Definition : stain impact angle between the trajectory vector and the impacted surface
-
Spatter Measurement Data, e.g.
Ink_Trial_INPUT.csv,Swineblood_Trial_INPUT.csv -
Castoff_Reconstruction_DRIVER.m -
DRIVER.csvproducesDRIVER.matrequired forCastoff_Reconstruction_MAIN.m -
Castoff_Reconstruction_MAIN.m -
Castoff_Reconstruction_FUNC.m -
lineSegmentIntersect.m -
meshVolume.m -
point_to_line.m -
gauss_distribution.m -
CircleFitByPratt.m -
Castoff_Reconstruction_POST.m -
combvec2.m -
inpolyhedron.m -
triangulateFaces.m -
linecirc.m -
generate_input.mlinecirc.mplane_line_intersect.mrotate_3D.mtriangulateFaces.m
All licenses for third party scripts are included and must be kept with provided scripts. If third party materials were not cited within the repository Licenses folder, this was not intentional by the author.
- Updated GNU Octave (written and tested with 5.2.0)
- Here is the link for previous GNU Octave version downloads (5.2.0) for windows
- All packages can be updated to the latest version by running:
pkg update
- Updated Java
- Psychtoolbox-3 Requirements install prior to Psychtoolbox download
- Psychtoolbox is a free set of Matlab and GNU Octave functions for vision and neuroscience research.
- Subversion 1.7.x command-line client (e.g. Sliksvn) MUST INSTALL COMPLETE SOFTWARE(https://sliksvn.com/download/)
- SlikSVN is a standalone command-line Subversion client for Windows.
- 64-Bit GStreamer-1.16.0 MSVC runtime or later versions MUST INSTALL COMPLETE SOFTWARE(https://gstreamer.freedesktop.org/data/pkg/windows/1.16.0/gstreamer-1.0-msvc-x86_64-1.16.0.msi) - GStreamer is a library for constructing graphs of media-handling components.
- Download Psychtoolbox Toolbox Version 3 (PTB-3)
- Unzip download in new folder C:\toolbox
- Open Octave, copy and paste the following script to the command window to install Psychtoolbox:
cd C:\toolbox DownloadPsychtoolbox C:\toolbox
- Follow command window prompts to finish installation.
- Prior to running Castoff_Recosntruction_DRIVER.m or Castoff_Reconstruction_MAIN.m the io and signal packages need to be loaded.
- Load packages on an as needed basis by copying and pasting the following script to the command window:
pkg load io pkg load signal
- Required packages for download.
- Automatically load packages at Octave startup with the following:
- Go to C:\Octave\Octave-5.2.0\mingw64\share\octave\site\m\startup
- Open octaverc with a text editor.
- Copy and paste the following lines to the end of the document:
pkg load io pkg load signal
- Load packages on an as needed basis by copying and pasting the following script to the command window:
- Save file.
- If octaverc does not exist, copy, paste, and save the following to a text file and save to the directory (C:\Octave\Octave-5.2.0\mingw64\share\octave\site\m\startup):
pkg load io pkg load signal
- If octaverc does not exist, copy, paste, and save the following to a text file and save to the directory (C:\Octave\Octave-5.2.0\mingw64\share\octave\site\m\startup):
Instructions to Run a first spatter pattern example (Ink_Trial_INPUT.csv), with stains measured by hand or with Hemospat:
- Make sure you have Matlab or Octave installed as per the installation instructions above
- Save all files included in the distribution to the same directory.
Run(F5) the DRIVER (Castoff_Reconstruction_DRIVER.m).- Check that the data name line 64 in the MAIN (
Castoff_Reconstruction_MAIN.m) readsINK_Trial_INPUT_DRIVER.mat Run(F5) the MAIN (Castoff_Reconstruction_MAIN.m); producing figure(3) and figure(4+) which simulate Figures 5 and 7 of McCleary et al FSI 2021, respectively.- Maximize figure 4 to observe the reconstructed swing regions in blue, green and red
- MAIN also outputs the
Total Elapsed Cluster Analysis Time:, in seconds (s), with the total program run time excluding any variable time from user input. See User Inputs: (DRIVER) res for estimated runtimes. - MAIN also save the Resultant Variables as a .mat file denoted by Castoff_Reconstruction.mat
- MAIN and Output the Cast-off Reconstruction Results with warnings as a txt-file denoted by
Castoff_Reconstruction_OUTPUT.txtdisplayed in Figures(4+) (figure number is dependent on clustering method). Figure(1) shows the Pratt fit automated reference point. Figure(3) outputs each clustered cast-off reconstructed arc. - Run again after changing the Stain Width Measurement Uncertainty (mm) or Cast-off Reconstruction Resolution (cm) in the input file
Ink_Trial_INPUT.csv
Instructions to Run a second spatter pattern example (FARO_Trial_10_INPUT.csv), with stains measured with FARO:
- Make sure you have Matlab or Octave installed as per the installation instructions above
- Save all files included in the distribution to the same directory.
- Modify the data name line 64 in the MAIN (
Castoff_Reconstruction_MAIN.m) toFARO_Trial_10_INPUT_DRIVER.mat Run(F5) the DRIVER (Castoff_Reconstruction_DRIVER_FARO.m).Run(F5) the MAIN (Castoff_Reconstruction_MAIN.m); producing figure(3) and figure(4+) which simulate Figures 5 and 6 of McCleary et al FSI 2021, respectively.- Maximize figure 4 to observe the reconstructed swing regions in blue, green and red
- MAIN also outputs the
Total Elapsed Cluster Analysis Time:, in seconds (s), with the total program run time excluding any variable time from user input. See User Inputs: (DRIVER) res for estimated runtimes. - MAIN also save the Resultant Variables as a .mat file denoted by Castoff_Reconstruction.mat
- MAIN and Output the Cast-off Reconstruction Results with warnings as a txt-file denoted by
Castoff_Reconstruction_OUTPUT.txtdisplayed in Figures(4+) (figure number is dependent on clustering method). Figure(1) shows the Pratt fit automated reference point. Figure(3) outputs each clustered cast-off reconstructed arc. - Run again after changing the Stain Width Measurement Uncertainty (mm) or Cast-off Reconstruction Resolution (cm) in the input file
Ink_Trial_INPUT.csv
-
Make sure you have Matlab or Octave installed as per the installation instructions above
-
Enter spatter pattern information in an input
_INPUT.csvfile with same format as either of the examples above - -
Modify the data name line 64 in the MAIN (
Castoff_Reconstruction_MAIN.m) toinput_DRIVER.mat -
Follow same instructions (4-10) as either examples above.
- creates a
_DRIVER.matfile as input forCastoff_Reconstruction_MAIN.m
-
figure(2)- inputs overview : an initial look at the inputted room dimensions, stain impact locations, stain trajectories, automated reference point, and best fit plane
-
figure(3)- reconstruction overview : an active display of clustered stains & trajectories and clustered planes with projected stains reference during spatter reconstruction
-
figure(4)- reconstruction results : inputs overview overlaid with user-defined three-dimensional region of statistical likelihood
- saved simulation workspace and textual overview of
INPUTS(resolution, user preferences, inputed stain details) andCast-off Reconstruction Results:(High, Medium, and Low Percentile reconstructed XYZ-Vertices & Faces)
- Enter known cast-off motion path x, y, and z-coordinates into
_MOTION.csvwith same format asInk_Trial_MOTION.csv - Follow same instructions (4-10) as either examples above.
Run(F5) the Post-processor (Castoff_Reconstruction_POST.m)
- Show and hide legends from resultant figures using:
legend show legend hide
- Rotate, zoom in/out, and pan figure controls are located on the top of the figures to allow for other viewing angles. Also, for all figures besides Figure(1), view(az,el) can be used for precise three-dimensional viewing angles where az is the azimuth angle (in degrees) and el is the polar (or elevation) angle (in degrees).