python實(shí)現(xiàn)3D地圖可視化
基于python代碼的3D地圖可視化,供大家參考,具體內(nèi)容如下
介紹
使用Python對(duì)地圖進(jìn)行3D可視化。以地圖為地圖,可以在三維空間對(duì)軌跡、點(diǎn)進(jìn)行可視化。
庫(kù)
我們使用了多個(gè)庫(kù):
1.gdal;
主要是用于讀取地圖信息,這個(gè)庫(kù)在GIS中很常用,使用C++代碼編寫的,如果安裝不了需要在pypi里面找一下對(duì)應(yīng)的資源。
2.opencv;
很常用的圖像處理庫(kù)。
3.matplotlib;
常用的可視化庫(kù)
結(jié)果
廢話不多說直接上結(jié)果:
代碼
直接上代碼,代碼很簡(jiǎn)單。
from osgeo import gdal import cv2 gdal.UseExceptions() ds = gdal.Open('E:/Pythoncode/讀取地理信息/無標(biāo)題.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/讀取地理信息/無標(biāo)題.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()
以上就是本文的全部?jī)?nèi)容,希望對(duì)大家的學(xué)習(xí)有所幫助,也希望大家多多支持腳本之家。
- 詳解python 利用echarts畫地圖(熱力圖)(世界地圖,省市地圖,區(qū)縣地圖)
- python使用pyecharts庫(kù)畫地圖數(shù)據(jù)可視化的實(shí)現(xiàn)
- Python學(xué)習(xí)之用pygal畫世界地圖實(shí)例
- python實(shí)現(xiàn)Pyecharts實(shí)現(xiàn)動(dòng)態(tài)地圖(Map、Geo)
- 使用Python實(shí)現(xiàn)畫一個(gè)中國(guó)地圖
- Python地圖繪制實(shí)操詳解
- 代碼分析Python地圖坐標(biāo)轉(zhuǎn)換
- 利用python繪制中國(guó)地圖(含省界、河流等)
- 如何利用Python快速繪制海報(bào)級(jí)別地圖詳解
相關(guān)文章
解決python gdal投影坐標(biāo)系轉(zhuǎn)換的問題
今天小編就為大家分享一篇解決python gdal投影坐標(biāo)系轉(zhuǎn)換的問題,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2020-01-01Django 查詢數(shù)據(jù)庫(kù)返回JSON的實(shí)現(xiàn)
和前端交互全部使用JSON,如何將數(shù)據(jù)庫(kù)查詢結(jié)果轉(zhuǎn)換成JSON格式,本文就來介紹一下,感興趣的小伙伴們可以參考一下2021-08-08Python學(xué)習(xí)筆記之常用函數(shù)及說明
俗話說“好記性不如爛筆頭”,老祖宗們幾千年總結(jié)出來的東西還是有些道理的,所以,常用的東西也要記下來,不記不知道,一記嚇一跳,乖乖,函數(shù)咋這么多捏2014-05-05關(guān)于sys.stdout和print的區(qū)別詳解
今天小編就為大家分享一篇關(guān)于sys.stdout和print的區(qū)別詳解,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2019-12-12在Python中使用MongoEngine操作數(shù)據(jù)庫(kù)教程實(shí)例
這篇文章主要介紹了在Python中使用MongoEngine操作數(shù)據(jù)庫(kù)教程實(shí)例,文中通過示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2019-12-12Python設(shè)計(jì)模式之門面模式簡(jiǎn)單示例
這篇文章主要介紹了Python設(shè)計(jì)模式之門面模式,簡(jiǎn)單描述了門面模式的概念、原理,并結(jié)合實(shí)例形式給出了Python定義與使用門面模式的具體操作技巧,需要的朋友可以參考下2018-01-01