Identification-of-core-XCT-photographs simplified

  1. Extracte the data in the dicom files, such as pixel, rescale slope and rescale intercept.
  2. Use the extracted data to calculate CT value by the formula:

    Hu = pixel * rescale slope + rescale intercept

  3. Scale the value of pixel data into 0~255, so that we can apply image processing algorithms to this gray image.

Use image processing algorithms to find the target

  1. make the image blur
  2. use Canny algorithms to find the edge of the image
  3. find the inner circle of the rock container by Circle Hough Transform
    • Inner circle of the rock container fliter out the rock container when calculating porosity later
  4. Threshold the image to create a binary image
  5. apply findContours algorithm.

Calculate porosity

  1. get CT value of the point inside in circle and contour of the target.
  2. get the CTG value
  3. apply the formula (CTG - Hu) / CTG
  4. sum them up
  5. divide by the number of pixel in the target

Areas for improvement

  1. The efficiency of the whole alogrithm is bad. It take a few second to calculate the porosity of a photo.
  2. CT values of outter cover that protect the core sample and the container of the core sample should be totally cut out.
Select a repo