我有一个简单的脚本,用于绘制3D中的纬度,经度和深度。 是否有一个python模块或一些简单的解决方案我可以用来获得这一点的横截面?
横截面应该是一个垂直平面,我将为其设置两个坐标(在平面视图中),深度和厚度(垂直于我想要数据的横截面每个方向有多少公里)。
类似于point1(lat1,lon1),point2(lat2,lon2),depthmin,depthmax,thicknes
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
lats=[]
lons=[]
depts=[]
file_open=open("catZMAPVpVs")
for fo in file_open:
element=fo.split("\t")
lons.append(float(element[0]))
lats.append(float(element[1]))
depts.append(-float(element[6]))
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.scatter(lons,lats, depts, c='r', marker='o')
ax.set_xlabel('lon')
ax.set_ylabel('lat')
ax.set_zlabel('depth')
plt.show()
谢谢!