利用Python实现山体阴影效果渲染

0
分享 2018-04-04
前段时间老在和Python的matplotlib过不去,虽然这个东西特别简单好用,但是默认绘制出来的图形简单中带着特有的丑,比如我对地震数据情况进行了一个2D散点图绘制,啥参数语句也不加,默认出现的就是这个效果:






好吧……虽然是默认的效果,这样太单薄了……所以就琢磨着要不加个背景图?比如这样的:






嗯,在下面叠个世界地图做背景,立马就有GIS的赶脚了……






所以发现,matplotlib的imshow这个功能确实灰常好用,用了真的可以扶摇直上九万里啊。



所以,我决定用它对dem生成的山体阴影做一个渲染,比如直接出dem是这样的:






然后做一个渲染的话,可以是这样的:








利用ArcGIS的空间分析扩展生成山体阴影分析图,是这样的:








最后,只需要设定两个图形,渲染+阴影,然后就变成这个样子:






是不是瞬间就觉得高大上了……



matplotlib中imshow的功能灰常灰常之强大,能够把任意numpy的数字直接渲染成图,或者可以直接加载图片为任意分析图形的背景,在分析和可视化的那是灰常有用。



最后,给出以上绘制的全部源代码:
import arcpy
import matplotlib.pyplot as plt


ras = "e:/dem/rkzNorth2.tif"
hy = arcpy.RasterToNumPyArray(arcpy.Raster(ras))


#渲染灰度图
ax = plt.figure(figsize=(9,9))
plt.imshow(hy, cmap= plt.cm.gray)


#用gist_earth色带渲染高程样式图
ax = plt.figure(figsize=(9,9))
plt.imshow(hy, cmap= plt.cm.gist_earth)


#利用ArcGIS的空间分析扩展制作山体阴影分析图
arcpy.CheckOutExtension("Spatial")
#太阳高度角
azimuth = 315.0
#太阳方位角
altitude = 45.0
hs = arcpy.sa.Hillshade(ras, azimuth, altitude)
ax = plt.figure(figsize=(9,9))
plt.imshow(arcpy.RasterToNumPyArray(hs), cmap= plt.cm.gray)


#联合渲染
ax = plt.figure(figsize=(9,9))
plt.imshow(arcpy.RasterToNumPyArray(hs), cmap= plt.cm.gray)
plt.imshow(hy, cmap= plt.cm.gist_earth,alpha=0.8)

打完收工

文章来源:http://blog.csdn.net/allenlu2008/article/details/79734221

5 个评论

Merry GISmas!
Merry GISmas!第一个
逮住小尾巴
每个ID只需发现一个彩蛋即可中奖欧!请勿重复回复!
每个ID只需发现一个彩蛋即可中奖欧!请勿重复回复!

要回复文章请先登录注册