Plot Stereonet di Python

Saya teringat ketika sewaktu kuliah dahulu harus secara manual plotting stereonet untuk strike and dip.Dan lagi modalnya cukup banyak yaitu kertas kalkir dan printed template nya. Mungkin di era sekarang sudah cukup menyenangkan harusnya, udah zamannya coding-codingan jadi ya harusnya lebih gampang. Jadi artikel ini hanya sekedar sharing-sharing metode saja. Langsung aja kita go to the coding nya.

Python

Pertama kita harus install library terpentingnya

pin install mplstereonet

Kemudian baru kita import library standar untuk plotting lah

import numpy as np
import mplstereonet
import matplotlib.pyplot as plt
%matplotlib inline

Setelah selesai import bahan-bahannya, mari kita masuk ke input datanya. Untuk tulisan kali ini kita hanya memakai data sedikit saja, nanti bisa diatur-atur pakai pandas atau iterasi untuk plot banyak.

#bedding
strike1, dip1 = 100, 45
#fault
strike2, dip2 = 320, 78

Nah kita sudah input data untuk perlapisan dan sesar. Selanjutnya tinggal kita visualisasikan saja dengan coding di bawah ini.

fig = plt.figure(figsize=(10,10))
ax = fig.add_subplot(111, projection='stereonet')
ax.plane(strike1, dip1, c='b', label='Perlapisan Arah %03d/%02d' % (strike1, dip1))
ax.plane(strike2, dip2, c='r', label='Sesar Arah %03d/%02d' % (strike2, dip2))
ax.legend()

Terima kasih

http://geologyandpython.com/structural_geology.html