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

前段时间老在和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)

打完收工