python實(shí)現(xiàn)DEM數(shù)據(jù)的陰影生成的方法
更新時(shí)間：2019年07月23日 14:26:28 作者：林木聲
這篇文章主要介紹了python實(shí)現(xiàn)DEM數(shù)據(jù)的陰影生成的方法，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧
相關(guān)的依賴庫(kù)在我的github網(wǎng)站上
首先貼出代碼：
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__':
'''
# 另一種從網(wǎng)上直接下載DEM數(shù)據(jù)
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)
# 定義一個(gè)光線向量
# 第一個(gè)表示和豎直方向的夾角，第二個(gè)表示由東向西照射
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') # 不顯示坐標(biāo)軸
plt.subplot(1,2,2)
plt.imshow(shadow,cmap='gray')
plt.axis('off')
plt.show()
生成的結(jié)果如下：左邊是DEM數(shù)據(jù)，右邊是shadow
以上就是本文的全部?jī)?nèi)容，希望對(duì)大家的學(xué)習(xí)有所幫助，也希望大家多多支持腳本之家。
相關(guān)文章
python中的不可變數(shù)據(jù)類型與可變數(shù)據(jù)類型詳解
探尋python的數(shù)據(jù)類型是否可變，也可以更好的理解python對(duì)內(nèi)存的使用情況，下面這篇文章主要給大家介紹了關(guān)于python中不可變數(shù)據(jù)類型與可變數(shù)據(jù)類型的相關(guān)資料，文中通過(guò)示例代碼介紹的非常詳細(xì)，需要的朋友可以參考下2018-09-09
Python 實(shí)現(xiàn)網(wǎng)課實(shí)時(shí)監(jiān)控自動(dòng)簽到、打卡功能
這篇文章主要介紹了Python實(shí)現(xiàn)網(wǎng)課實(shí)時(shí)監(jiān)控自動(dòng)簽到，打卡功能，本文通過(guò)實(shí)例代碼給大家介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2020-03-03
Pytorch統(tǒng)計(jì)參數(shù)網(wǎng)絡(luò)參數(shù)數(shù)量方式
這篇文章主要介紹了Pytorch統(tǒng)計(jì)參數(shù)網(wǎng)絡(luò)參數(shù)數(shù)量方式，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。如有錯(cuò)誤或未考慮完全的地方，望不吝賜教2023-02-02
在Python中，類表示具有相同屬性和方法的對(duì)象的集合。在使用類時(shí)，需要先定義類，然后再創(chuàng)建類的實(shí)例，通過(guò)類的實(shí)例就可以訪問(wèn)類中的屬性和方法了。本文就來(lái)和大家一起聊聊Python中類的定義和使用，需要的可以參考一下2022-08-08
在本篇文章里小編給各位整理的是關(guān)于用vue.js組件模擬v-model指令實(shí)例方法和相關(guān)代碼，需要的朋友們可以跟著學(xué)習(xí)下。2019-07-07
python實(shí)現(xiàn)簡(jiǎn)單的計(jì)時(shí)器功能函數(shù)
這篇文章主要介紹了python實(shí)現(xiàn)簡(jiǎn)單的計(jì)時(shí)器功能函數(shù),涉及Python操作時(shí)間的技巧,具有一定參考借鑒價(jià)值,需要的朋友可以參考下2015-03-03
對(duì)Python游戲有所了解的朋友都知道，在2D的游戲制作中，經(jīng)常會(huì)用到一個(gè)模塊pygame，他能幫助我們實(shí)現(xiàn)很多方便使用的功能,本文就來(lái)用pygame模塊制作代碼雨效果，需要的可以了解一下2023-05-05
這篇文章主要為大家詳細(xì)介紹了python多線程分塊讀取文件，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2019-08-08
這篇文章主要介紹了python可視化爬蟲界面之天氣查詢的相關(guān)資料,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2019-07-07
使用Python實(shí)現(xiàn)Excel表格轉(zhuǎn)圖片
在數(shù)據(jù)處理與信息分享過(guò)程中,Excel表格作為一種強(qiáng)大的數(shù)據(jù)管理工具被廣泛應(yīng)用,這篇文章主要為大家詳細(xì)介紹了如何使用Python將Excel表格轉(zhuǎn)換為圖片,需要的可以參考下2024-04-04