python实现DEM数据的阴影生成的方法

 更新时间:2019年07月23日 14:26:28   作者:林木声  
这篇文章主要介绍了python实现DEM数据的阴影生成的方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧

相关的依赖库在我的github网站上

首先贴出代码:

import solar
from gradient import *
from shadows import *
import numpy as np

import matplotlib.pyplot as plt # plt 用于显示图片
import matplotlib.image as mpimg # mpimg 用于读取图片


# dem
import srtm
if __name__ == '__main__':

  
  '''
  # 另一种从网上直接下载DEM数据
  geo_elevation_data = srtm.get_data()
  image = geo_elevation_data.get_image((500, 500), (45, 46), (13, 14), 300)
  # the image s a standard PIL object, you can save or show it:
  image.show()

  image = np.asarray(image)

  print(image.shape)

  '''

  # 读入高程信息
  filename = 'dempyrenees.asc'
  dem = np.loadtxt(filename,skiprows=6,delimiter=' ')

  # 高程信息的维度
  # print(dem)
  print(dem.shape)

  # 定义一个光线向量
  # 第一个表示和竖直方向的夹角,第二个表示由东向西照射
  sv = normal_vector(45, 270)

  # 生成阴影
  shadow = project_shadows(dem=dem,sun_vector=sv,dx=30)
  print(shadow[1,:])
  print(shadow.shape)

  # 显示dem和shadow图像
  plt.figure()
  plt.subplot(1,2,1)
  plt.imshow(dem,cmap='gray') # 显示灰度图像
  plt.axis('off') # 不显示坐标轴

  plt.subplot(1,2,2)
  plt.imshow(shadow,cmap='gray')
  plt.axis('off')

  plt.show()

生成的结果如下:左边是DEM数据,右边是shadow

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持脚本之家。

相关文章

  • 使用pyqt5搭建yolo3目标识别界面的方法

    使用pyqt5搭建yolo3目标识别界面的方法

    这篇文章主要介绍了使用pyqt5搭建yolo3目标识别界面的方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2021-03-03
  • Python mplfinance库绘制金融图表实现数据可视化实例探究

    Python mplfinance库绘制金融图表实现数据可视化实例探究

    mplfinance(Matplotlib Finance),它是基于Matplotlib的库,专门用于创建金融图表和交互式金融数据可视化,本文将深入介绍 mplfinance,包括其基本概念、功能特性以及如何使用示例代码创建各种金融图表
    2024-01-01
  • pycharm全局搜索的具体步骤

    pycharm全局搜索的具体步骤

    在本篇文章里小编给大家整理的是一篇关于pycharm全局搜索的具体步骤,有需要的朋友们可以参考下。
    2020-07-07
  • Django传递数据给前端的3种方式小结

    Django传递数据给前端的3种方式小结

    Django从后台往前台传递数据时有多种方法可以实现,下面这篇文章主要给大家介绍了关于Django传递数据给前端的3种方式,文中通过代码介绍的非常详细,需要的朋友可以参考下
    2024-01-01
  • django在接受post请求时显示403forbidden实例解析

    django在接受post请求时显示403forbidden实例解析

    这篇文章主要介绍了django在接受post请求时显示403forbidden实例解析,小编觉得还是挺不错的,具有一定借鉴价值,需要的朋友可以参考下
    2018-01-01
  • Python中列表与元组的乘法操作示例

    Python中列表与元组的乘法操作示例

    这篇文章主要介绍了Python中列表与元组的乘法操作,结合简单实例形式分析了Python中列表、元组的乘法,并附带分析了字符串的乘法及元组乘法操作的注意事项,需要的朋友可以参考下
    2018-02-02
  • Scrapy 配置动态代理IP的实现

    Scrapy 配置动态代理IP的实现

    这篇文章主要介绍了Scrapy 配置动态代理IP的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2020-09-09
  • Python库Theano深度神经网络的设计训练深入探究

    Python库Theano深度神经网络的设计训练深入探究

    Theano是一个用于深度学习的Python库,它提供了高效的数值计算和自动微分功能,使得深度神经网络的设计和训练变得更加容易,本文将深入探讨Theano的功能和用法,并提供丰富的示例代码,帮助大家入门深度学习
    2024-01-01
  • Python 使用 pip 安装 matplotlib 模块的方法

    Python 使用 pip 安装 matplotlib 模块的方法

    matplotlib是python中强大的画图模块,这篇文章主要介绍了Python 使用 pip 安装 matplotlib 模块(秒解版),本文给大家介绍的非常详细,需要的朋友可以参考下
    2023-02-02
  • python3.6根据m3u8下载mp4视频

    python3.6根据m3u8下载mp4视频

    这篇文章主要为大家详细介绍了python3.6根据m3u8下载mp4视频,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2019-06-06

最新评论