Osteopetrosis is a bone disease that is marked by an increase in bone density due to genetic defects in cells called osteoclasts. HrpQCT is a common imaging tool used to assess bone diseases such as osteopetrosis because it provides high-resolution, three-dimensional images and can differentiate between cortical (outer) bone and trabecular (inner spongy bone) bone. However, available image analysis methods cannot account for the phenotypical variability that is seen with osteopetrotic bone, and therefore we developed an image analysis method using iterative erosion values in ImageJ. We used this method to quantify bone volume fraction in the trabecular region of healthy bone and compared the results to the measurements of HrpQCT manufacturer (ScanCo) software, to find the erosion value that minimizes error between the two methods. Using our method, we found that higher erosion radius values were associated with males (mean erosion value=1.64457 mm) and younger samples (mean erosion value=1.768078 mm). The findings from this study can be used for the analysis of osteopetrotic bones.