欧美bbbwbbbw肥妇,免费乱码人妻系列日韩,一级黄片

python實現(xiàn)3D地圖可視化

 更新時間:2020年03月25日 17:10:42   作者:機關槍小哥哥  
這篇文章主要為大家詳細介紹了python實現(xiàn)3D地圖可視化,文中示例代碼介紹的非常詳細,具有一定的參考價值,感興趣的小伙伴們可以參考一下

基于python代碼的3D地圖可視化,供大家參考,具體內容如下

介紹

使用Python對地圖進行3D可視化。以地圖為地圖,可以在三維空間對軌跡、點進行可視化。

我們使用了多個庫:

1.gdal;
主要是用于讀取地圖信息,這個庫在GIS中很常用,使用C++代碼編寫的,如果安裝不了需要在pypi里面找一下對應的資源。

2.opencv;
很常用的圖像處理庫。

3.matplotlib;
常用的可視化庫

結果

廢話不多說直接上結果:

代碼

直接上代碼,代碼很簡單。

from osgeo import gdal
import cv2
gdal.UseExceptions()

ds = gdal.Open('E:/Pythoncode/讀取地理信息/無標題.tif')
bandg = ds.GetRasterBand(1)
elevationg = bandg.ReadAsArray()

bandr = ds.GetRasterBand(2)
elevationr = bandr.ReadAsArray()

bandb = ds.GetRasterBand(3)
elevationb = bandb.ReadAsArray()

import matplotlib.pyplot as plt
nrows, ncols = elevationr.shape

elevation= cv2.merge([elevationg,elevationr,elevationb])#
# I'm making the assumption that the image isn't rotated/skewed/etc. 
# This is not the correct method in general, but let's ignore that for now
# If dxdy or dydx aren't 0, then this will be incorrect
x0, dx, dxdy, y0, dydx, dy = ds.GetGeoTransform()

x1 = x0 + dx * ncols
y1 = y0 + dy * nrows

plt.imshow(elevation, cmap='gist_earth', extent=[x0, x1, y1, y0])
plt.show()

from PIL import Image
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
fig = plt.figure()
ax = Axes3D(fig)

img = Image.open('E:/Pythoncode/讀取地理信息/無標題.tif')
xx=[]
yy=[]
colall=[]
x = img.size[0]
y = img.size[1]
for i in range(x):
 for j in range(y):
 
 r = hex(img.getpixel((i, j))[0])[2:]
 b = hex(img.getpixel((i, j))[1])[2:]
 g = hex(img.getpixel((i, j))[2])[2:]
 
 if len(r) == 1:
 r = '0' + r
 if len(b) == 1:
 b = '0' + b
 if len(g) == 1:
 g = '0' + g
 col = '#' + r + b + g
 colall.append(col)
 xx.append(x0 + dx * i)
 yy.append(y0 + dy * j)
 # col = '#FF00FF'
ax.scatter(xx, yy, 5, c=colall, alpha=0.5)
plt.show()

以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支持腳本之家。

相關文章

  • Python中查看文件名和文件路徑

    Python中查看文件名和文件路徑

    本篇文章主要介紹了Python中查看文件名和文件路徑的方法。具有很好的參考價值,下面跟著小編一起來看下吧
    2017-03-03
  • 解決python gdal投影坐標系轉換的問題

    解決python gdal投影坐標系轉換的問題

    今天小編就為大家分享一篇解決python gdal投影坐標系轉換的問題,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧
    2020-01-01
  • Django 查詢數(shù)據(jù)庫返回JSON的實現(xiàn)

    Django 查詢數(shù)據(jù)庫返回JSON的實現(xiàn)

    和前端交互全部使用JSON,如何將數(shù)據(jù)庫查詢結果轉換成JSON格式,本文就來介紹一下,感興趣的小伙伴們可以參考一下
    2021-08-08
  • Python學習筆記之常用函數(shù)及說明

    Python學習筆記之常用函數(shù)及說明

    俗話說“好記性不如爛筆頭”,老祖宗們幾千年總結出來的東西還是有些道理的,所以,常用的東西也要記下來,不記不知道,一記嚇一跳,乖乖,函數(shù)咋這么多捏
    2014-05-05
  • 關于sys.stdout和print的區(qū)別詳解

    關于sys.stdout和print的區(qū)別詳解

    今天小編就為大家分享一篇關于sys.stdout和print的區(qū)別詳解,具有很好的參考價值,希望對大家有所幫助。一起跟隨小編過來看看吧
    2019-12-12
  • 在centos7中分布式部署pyspider

    在centos7中分布式部署pyspider

    PySpider:一個國人編寫的強大的網(wǎng)絡爬蟲系統(tǒng)并帶有強大的WebUI。采用Python語言編寫,分布式架構,支持多種數(shù)據(jù)庫后端,強大的WebUI支持腳本編輯器,任務監(jiān)視器,項目管理器以及結果查看器。
    2017-05-05
  • Python matplotlib可視化繪圖詳解

    Python matplotlib可視化繪圖詳解

    這篇文章主要介紹了Python matplotlib繪圖可視化知識點整理(小結),小編覺得挺不錯的,現(xiàn)在分享給大家,也給大家做個參考。一起跟隨小編過來看看吧
    2021-09-09
  • Python Color類與文字繪制零基礎掌握

    Python Color類與文字繪制零基礎掌握

    這篇文章主要介紹了Python Color類與文字繪制,文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友們下面隨著小編來一起學習學習吧
    2022-08-08
  • 在Python中使用MongoEngine操作數(shù)據(jù)庫教程實例

    在Python中使用MongoEngine操作數(shù)據(jù)庫教程實例

    這篇文章主要介紹了在Python中使用MongoEngine操作數(shù)據(jù)庫教程實例,文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友可以參考下
    2019-12-12
  • Python設計模式之門面模式簡單示例

    Python設計模式之門面模式簡單示例

    這篇文章主要介紹了Python設計模式之門面模式,簡單描述了門面模式的概念、原理,并結合實例形式給出了Python定義與使用門面模式的具體操作技巧,需要的朋友可以參考下
    2018-01-01

最新評論