Build Status Coverage

Desurvey and composite drill hole tables from the mining industry.


Get the latest stable release with Julia's package manager:

] add DrillHoles


Given a collar table, a survey table and at least one interval table (such as assay and lithology), the function desurvey can be used for desurveying and compositing. Examples of these tables are shown bellow:


  • Collar table: stores the coordinates (X, Y, Z) of each drill hole with given ID (HOLEID).
  • Survey table: stores the arc length (AT) and azimuth (AZM) and dip (DIP) angles along the drill hole trajectory. Together with the collar table it fully specifies the trajectory.
  • Interval table: stores the actual measurements taken on cylinders of rock defined by an interval of arc lenghts (FROM and TO). Usually, there are multiple interval tables with different types of measurements.

Assuming that each of these tables was loaded into a Tables.jl table (e.g. CSV.File, DataFrame), we can use the following constructors to automatically detect the columns:

using DrillHoles
using CSV

collar = Collar(CSV.File("collar.csv"))
survey = Survey(CSV.File("survey.csv"))
assay  = Interval(CSV.File("assay.csv"))
litho  = Interval(CSV.File("litho.csv"))

If the columns of the tables follow an exotic naming convention, users can manually specify the names with keyword arguments:

# manually specify column with hole ID
Collar(CSV.File("collar.csv"), holeid = :MYID)

Please check the documentation of Collar, Survey and Interval for more details.

By default, the desurvey function returns a geospatial table compatible with GeoStats.jl. It supports different dip angle conventions from open pit and underground mining as well as different stepping methods:

samples = desurvey(collar, survey, [assay, litho])

The option geom can be used to control the output format, and the option len can be used for compositing. Please check the documentation for more details.