相關的依賴庫在我的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
以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支持腳本之家。
更多文章、技術交流、商務合作、聯系博主
微信掃碼或搜索:z360901061

微信掃一掃加我為好友
QQ號聯系: 360901061
您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。
【本文對您有幫助就好】元
