leonbohmann/fracture-suite

A small toolkit for analyzing fracture patterns on glass plys.

For citation please use:

@software{fracture-suite,
  author       = {Leon Bohmann},
  title        = {fracsuite v1},
  data         = 2024-03,
  version      = {1.0.0},
  doi          = {10.5281/zenodo.10758208},
  url          = {https://github.com/leonbohmann/fracture-suite/tree/masterrelease}
}

How it works

⚠️ This is deprecated. The documentation will be updated in the future!

Several steps are necessary to analyze a fracture scan:

1. Cropping of input image (optional)
   1. Analyze the image and find the biggest rectangular shape
   2. Perspective transform the image, so that the rectangle is filling the extents
2. Preprocessing
   1. Gaussian Blur + (optional) Resize of the input image
   2. Adaptive Threshold
3. Contour detection
   1. Find all contours on the preprocessed image
   2. Filter Contours, remove all:
      1. Whose perimeter is too small
      2. Whose area is way too large (25000px²)
4. Create stencil with the found contours
   1. This helps to quickly remove all contours that lie within a bigger contour
   2. Draw all contours onto a new image (resulting image will display the cracks)
5. Skeletonization #1
   1. Skeletonize the stencil to minimize the crack width to 1px wide lines
   2. Use Erode/Dilate (closing kernel) to connect gaps in contours (this will widen the 1px wide lines)
6. Skeletonization #2
   1. Skeletonize the image again to retrieve the crack middle lines
7. (enabled by default) Remove splinters, whose are is all black in the original image and fill those spaces
8. Contour detection #2
   1. Now with minimal fuzziness, run the splinter detection again
9. Create splinters from resulting contour list

A couple of outputs are saved to a directory next to the input image. These are:

• Contour plot of recognized splinters (you have to zoom in, contour thickness is 1px)
• Filled plot of splinters
• Figures
  • Splinter size distribution
  • Splinter relative share distribution
  • Voronoi overview
  • Fracture Intensity overlay

Installation

pip install fracsuite

Usage (:warning: 0.1.7 only!)

Use the module directly

py -m fracsuite.splinters "path/to/image"

For details see: API Docs

Create a script

from fracsuite.splinters.analyzer import Analyzer, AnalyzerConfig

image = r"Path/to/some/image.bmp"

config = AnalyzerConfig()
# size of cropped image (if cropping is needed)
config.cropsize = (4000,4000)

analyzer = Analyzer(image, config)