Mathematica code for least-squares cone fitting and equal-area stereonet representation

Kieran F. Mulchrone, Daniel Pastor-Galán, Gabriel Gutiérrez-Alonso

Research output: Contribution to journalArticlepeer-review

6 Citations (Scopus)


In structural geology it is often assumed that folds are cylindrical. However, most structures are conical to some degree. Due to the lack of software capable of accurately estimating the best fit cone from a set of oriented data, we developed a Mathematica application capable of (1) plotting oriented data (lines and planes) on an equal area stereonet, (2) calculating the orientation matrix, the distribution shape and intensity parameters, (3) plotting the eigenvectors and (4) estimating and plotting the best fit cone, a small circle. We present both synthetic and natural data demonstrating its robustness and accuracy calculating the best fit cone.

Original languageEnglish
Pages (from-to)203-210
Number of pages8
JournalComputers and Geosciences
Publication statusPublished - 2013 Apr
Externally publishedYes


  • Complex conical shapes
  • Conical best fit
  • Stereogram plotting software
  • Stereographic projection
  • Stereographic projection software

ASJC Scopus subject areas

  • Information Systems
  • Computers in Earth Sciences


Dive into the research topics of 'Mathematica code for least-squares cone fitting and equal-area stereonet representation'. Together they form a unique fingerprint.

Cite this