python实现dem输出三维模型_python

61.9K

http://www.worldcat.org/wcpa/oclc/868259840?page=frame&url=http%3A%2F%2Fjmlr.org%2Fpapers%2Fvolume14%2Fdemsar13a%2Fdemsar13a.pdf%26checksum%3D3ad7736b9f5fa33f3e6a433eb1ca4a48&title=&linktype=digitalObject&detail=

worldcat.org

全网免费

62.0K

http://www.weblio.jp/redirect?url=http%3A%2F%2Feprints.fri.uni-lj.si%2F2267%2F1%2F2013-Demsar-Orange-JMLR.pdf&etd=8207739801c76b00

weblio.jp

全网免费

61.9K

http://www.jmlr.org/papers/volume14/demsar13a/demsar13a.pdf

jmlr.org

全网免费

61.9K

http://pdfs.semanticscholar.org/37bc/39697da911ccfbe463bad6b207862e55a007.pdf

pdfs.semanticscholar.org

全网免费

61.9K

http://jmlr.org/papers/volume14/demsar13a/demsar13a.pdf

jmlr.org

全网免费

61.9K

http://jmlr.csail.mit.edu/papers/volume14/demsar13a/demsar13a.p

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
要用 Python 实现 DEM 输出三维模型可以使用一些开源的库和工具,比如: 1. GDAL:可以用于读取 DEM 数据,生成高程网格等。 2. PyVista:可以用于可视化 DEM 数据,并生成三维模型。 3. Blender:可以通过 Python 脚本生成三维模型。 以下是使用 PyVista 库生成三维模型的示例代码: ```python import pyvista as pv import rasterio as rio # 读取 DEM 文件 with rio.open('dem.tif') as src: z = src.read(1) transform = src.transform # 创建 PyVista 网格 grid = pv.UnstructuredGrid(pv.vtk.vtkXMLUnstructuredGridReader().GetOutput()) grid.points = pv.vtk.vtkPoints() grid.cells = pv.vtk.vtkCellArray() # 将 DEM 中每个像素转换为网格节点 for i in range(z.shape[0]): for j in range(z.shape[1]): grid.points.InsertNextPoint(transform * (j, i, z[i, j])) # 将节点连接成三角形网格 for i in range(z.shape[0] - 1): for j in range(z.shape[1] - 1): cell = pv.vtk.vtkTriangle() cell.GetPointIds().SetId(0, i * z.shape[1] + j) cell.GetPointIds().SetId(1, i * z.shape[1] + j + 1) cell.GetPointIds().SetId(2, (i + 1) * z.shape[1] + j) grid.cells.InsertNextCell(cell) # 可视化网格 plotter = pv.Plotter() plotter.add_mesh(grid, cmap='terrain') plotter.show() ``` 要用 Python 实现 DEM 数据的阴影生成,可以使用以下步骤: 1. 计算光照方向。可以根据 DEM 数据的法向量和太阳高度角、方位角等参数计算。 2. 计算每个像素的阴影值。可以使用遮蔽法或者投影法等方法。 3. 可视化 DEM 数据,并根据阴影值着色。可以使用 PyVista 或者 Matplotlib 等库进行可视化。 以下是使用遮蔽法计算阴影值的示例代码: ```python import numpy as np import rasterio as rio # 读取 DEM 文件 with rio.open('dem.tif') as src: z = src.read(1) transform = src.transform # 计算 DEM 数据的法向量 dx, dy = np.gradient(z, transform[0], transform[4]) dz = np.ones_like(z) nx, ny, nz = np.cross(dx, dy, dz, axis=0, normalize=True) # 计算光照方向 theta = np.radians(30) # 太阳高度角 phi = np.radians(-45) # 太阳方位角 lx = np.cos(theta) * np.cos(phi) ly = np.cos(theta) * np.sin(phi) lz = np.sin(theta) # 计算每个像素的阴影值 shadow = np.zeros_like(z, dtype=bool) for i in range(z.shape[0]): for j in range(z.shape[1]): p = (j, i, z[i, j]) q = (p[0] - lx, p[1] - ly, p[2] - lz) r = (p[0] + lx, p[1] + ly, p[2] + lz) if np.dot(nz[i, j], q - p) > 0 and np.dot(nz[i, j], r - p) > 0: shadow[i, j] = True # 可视化 DEM 数据,并根据阴影值着色 import matplotlib.pyplot as plt from matplotlib.colors import LightSource fig, ax = plt.subplots() ls = LightSource(azdeg=315, altdeg=45) rgb = ls.shade(z, cmap='terrain', vert_exag=10, blend_mode='soft', dx=transform[0], dy=-transform[4], fraction=1, shade_alpha=0.8) ax.imshow(rgb) ax.imshow(shadow, cmap='gray', alpha=0.5) plt.show() ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值