Dragon绘制
import shutil
import tarfile
from mayavi import mlab
from os.path import join
dragon_tar_file = tarfile.open('Data/dragon.tar.gz')
dragon_tar_file.extractall('dragon_data')
dragon_tar_file.close()
dragon_ply_file = join('dragon_data', 'dragon_recon', 'dragon_vrip.ply')
mlab.pipeline.surface(mlab.pipeline.open(dragon_ply_file))
mlab.show()
shutil.rmtree('Data/dragon_data')
Canyon地形可视化
import zipfile
import numpy as np
from mayavi import mlab
hgt = zipfile.ZipFile('Data/N36W113.hgt.zip').read('N36W113.hgt')
data = np.fromstring(hgt, '>i2')
data.shape = (3601, 3601)
data = data.astype(np.float32)
data = data[:1000, 900:1900]
data[data == -32768] = data[data