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()