扫二维码与项目经理沟通
我们在微信上24小时期待你的声音
解答本文疑问/技术咨询/运营咨询/技术建议/互联网交流
基于python代码的3D地图可视化,供大家参考,具体内容如下
成都创新互联是一家专注于网站建设、成都做网站与策划设计,北关网站建设哪家好?成都创新互联做网站,专注于网站建设十年,网设计领域的专业建站公司;建站业务涵盖:北关等地区。北关做网站价格咨询:18980820575介绍
使用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()
另外有需要云服务器可以了解下创新互联scvps.cn,海内外云服务器15元起步,三天无理由+7*72小时售后在线,公司持有idc许可证,提供“云服务器、裸金属服务器、高防服务器、香港服务器、美国服务器、虚拟主机、免备案服务器”等云主机租用服务以及企业上云的综合解决方案,具有“安全稳定、简单易用、服务可用性高、性价比高”等特点与优势,专为企业上云打造定制,能够满足用户丰富、多元化的应用场景需求。
我们在微信上24小时期待你的声音
解答本文疑问/技术咨询/运营咨询/技术建议/互联网交流