利用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))

打完收工

  • 发表于:
  • 原文链接http://kuaibao.qq.com/s/20180329G07PUA00?refer=cp_1026
  • 腾讯「云+社区」是腾讯内容开放平台帐号(企鹅号)传播渠道之一,根据《腾讯内容开放平台服务协议》转载发布内容。

扫码关注云+社区

领取腾讯云代金券