Fault skeletonization for fault identification in a subterranean environment

A system can receive fault likelihood data about a subterranean environment and apply a binary mask filter using a tuning parameter to convert the fault likelihood data to binary distribution data having a plurality of pixels arranged in a plurality of profiles in at least two directions. The system can perform, for each profile of the plurality of profiles, fault skeletonization on the binary distribution data to form fault skeletonization data with pixels connected that represent part of a fracture. The system can convert the fault skeletonization data to seismic volume data and combine and filter the seismic volume data in the at least two directions to form combined seismic volume data. The system can output the combined seismic volume data as an image for use in detecting objects to plan a wellbore operation.