Great Circle Earthworks

../../_images/earthworks.png

Obj file was obtained from Miniminuteman (Milo Rossi, creator) on Patreon regarding the The Spectacular Archaeology of the Great Circle Earthworks video.

For the above surface, a section of the ‘terrain’ colormap was created using the cmap_section method, producing the following colormap:

../../_images/cmap_eworks.png 
import numpy as np
import matplotlib.pyplot as plt
import s3dlib.surface as s3d
import s3dlib.cmap_utilities as cmu

# Great Circle Earthworks Obj datafile
title = 'Great Circle Earthworks'
def trim_bottom(xyz) :
    x,y,z = xyz
    return x,y,np.where(z<-0.001,-0.001,z)

# Setup surface ................................................
cmap = cmu.section_cmap('terrain',0.23,0.8,name='eworks')

surface = s3d.get_surfgeom_from_obj("earthworks/Great Circle Earthworks 2.obj", cmap=cmap)
surface.transform( s3d.axisRot(-90,[1,0,0]) )
surface.map_geom_from_op(trim_bottom)
surface.transform(scale=[1,1,5])
surface.map_cmap_from_op()

# Construct figure, add surface, plot ..........................
fig = plt.figure(figsize=(7,7))
fig.text(0.975,0.975,str(surface), ha='right', va='top', fontsize='medium')
fig.text(0.5,0.8, title, ha='center', va='top', fontsize='xx-large')
ax = plt.axes(projection='3d', aspect='equal',proj_type='ortho' )
ax.set_axis_off()

s3d.auto_scale(ax,surface,uscale=0.7)
ax.add_collection3d(surface.shade())

fig.tight_layout(pad=1)
plt.show()