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

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

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

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

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

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

然后做一个渲染的话，可以是这样的：

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

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

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

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

最后，给出以上绘制的全部源代码：

import arcpy

importmatplotlib.pyplot as plt

ras = "e:/dem/rkzNorth2.tif"

hy =arcpy.RasterToNumPyArray(arcpy.Raster(ras))

#渲染灰度图

ax =plt.figure(figsize=(9,9))

#用gist_earth色带渲染高程样式图

ax =plt.figure(figsize=(9,9))

#利用ArcGIS的空间分析扩展制作山体阴影分析图

arcpy.CheckOutExtension("Spatial")

#太阳高度角

azimuth = 315.0

#太阳方位角

altitude = 45.0

hs =arcpy.sa.Hillshade(ras, azimuth, altitude)

ax =plt.figure(figsize=(9,9))

#联合渲染

ax =plt.figure(figsize=(9,9))

打完收工