tri = Delaunay(top_granitoid_verticesPD.values[:, 0:2])
top_granitoid_vertices = Mesh([xyz, tri.simplices]).c("darkcyan")
top_granitoid_vertices.name = 'Top of granite surface'
plot += top_granitoid_vertices.flag()
####################
## 3. Point objects
####################
printc("...plotting...", invert=1)
#FORGE Boundary
#Since the boundary area did not have a Z column, I assigned a Z value for where I wanted it to appear
border['zcoord'] = 1650
borderxyz = border[['xcoord', 'ycoord', 'zcoord']]
boundary = Points(borderxyz.values).c('k')
boundary.name = 'FORGE area boundary'
plot += boundary.flag()
#Microseismic
microseismicxyz = microseismic[['xloc', 'yloc', 'zloc']]
scals = microseismic[['mw']]
microseismicPts = Points(microseismicxyz.values, r=3).cellColors(scals,
cmap="jet")
microseismicPts.name = 'Microseismic events'
plot += microseismicPts.flag()
####################
## 4. Line objects
####################
#The path of well 58_32
xyz = well_5832_path[['X', 'Y', 'Z']].values
# Top Granite
xyz = top_granitoid_verticesPD.values
xyz[:, 2] = top_granitoid_verticesPD.values[:, 2] - 20
tri = Delaunay(top_granitoid_verticesPD.values[:, 0:2])
top_granitoid_vertices = Mesh([xyz, tri.simplices]).texture('white1')
top_granitoid_vertices.name = "Top of granite surface"
plt += top_granitoid_vertices.flag()
printc("plotting...", invert=1)
# Microseismic
microseismicxyz = microseismic[["xloc", "yloc", "zloc"]].values
scals = microseismic[["mw"]]
microseismicPts = Points(microseismicxyz, r=5).pointColors(scals, cmap="jet")
microseismicPts.name = "Microseismic events"
plt += microseismicPts.flag()
# FORGE Boundary. Since the boundary area did not have a Z column,
# I assigned a Z value for where I wanted it to appear
border["zcoord"] = 1650
borderxyz = border[["xcoord", "ycoord", "zcoord"]]
boundary = Line(borderxyz.values).extrude(zshift=120,
cap=False).lw(0).texture('wood1')
boundary.name = "FORGE area boundary"
plt += boundary.flag()
# The path of well 58_32
Well1 = Line(well_5832_path[["X", "Y", "Z"]].values, lw=2, c='k')
Well1.name = "Well 58-32"
plt += Well1.flag()