python教程 | 最标准的地图调用方式(国家测绘局提供数据)
时间:2022-07-26
本文章向大家介绍python教程 | 最标准的地图调用方式(国家测绘局提供数据),主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
天地图是国家测绘地理信息局建设的地理信息综合服务网站,是国家地理信息公共服务平台的公众版。 与常用的谷歌地图、腾讯地图、百度地图、微软地图、必应地图相比,天地图有什么不同呢?主要体现在数据的权威性和准确性。天地图发布的国界线、九段线等是准确无误的;另外国内只有天地图影像的坐标是无偏移的,其余地图的坐标都进行过加密处理。 Cartopy是一个基于Python的制图模块,其提供了加载在线地图的功能,那么如何添加调用天地图服务功能呢? 其实前期已有相关的工作,但是由于天地图服务升级,原先的方法都不再适用,这里给出的是最新的调用方法。
准备工作
网站:天地图 首先需要注册账号,然后选择开发资源→地图API
进入控制台,点击创建新应用,选择服务器端,这样可以得到一个密钥(key)
天地图地图服务采用OGC WMTS标准,但是该方法目前没有试成功,故采用了瓦片地图(XYZ Tiles)的形式。
添加调用方法
添加以下代码,注意把代码中的'your_key'替换成之前得到的key
import cartopy.io.img_tiles as cimgt
# 天地图矢量
class TDT_vec(cimgt.GoogleWTS):
def _image_url(self, tile):
x, y, z = tile
key = 'your_key'
url = 'http://t0.tianditu.gov.cn/DataServer?T=vec_w&x=%s&y=%s&l=%s&tk=%s' % (x, y, z, key)
return url
# 天地图遥感
class TDT_img(cimgt.GoogleWTS):
def _image_url(self, tile):
x, y, z = tile
key = 'your_key'
url = 'http://t0.tianditu.gov.cn/DataServer?T=img_w&x=%s&y=%s&l=%s&tk=%s' % (x, y, z, key)
return url
# 天地图地形
class TDT_ter(cimgt.GoogleWTS):
def _image_url(self, tile):
x, y, z = tile
key = 'your_key'
url = 'http://t0.tianditu.gov.cn/DataServer?T=ter_w&x=%s&y=%s&l=%s&tk=%s' % (x, y, z, key)
return url
调用演示
1、矢量底图、影像底图、地形底图
# 导入模块
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeat
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
from cartopy.mpl.ticker import LongitudeFormatter,LatitudeFormatter
import cartopy.io.shapereader as shpreader
import numpy as np
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False
fig = plt.figure(figsize=(18, 12))
ax = fig.add_subplot(1, 3, 1, projection=ccrs.PlateCarree())
ax.set_extent([118, 122, 28, 32],crs=ccrs.PlateCarree())
request = TDT_vec()
ax.add_image(request, 9)
ax.set_title('天地图矢量底图',fontsize=15)
gl = ax.gridlines(draw_labels=True, linewidth=1, color='k', alpha=0.5, linestyle='--')
gl.xlabels_top = gl.ylabels_right = False
gl.xformatter = LONGITUDE_FORMATTER
gl.yformatter = LATITUDE_FORMATTER
ax = fig.add_subplot(1, 3, 2, projection=ccrs.PlateCarree())
ax.set_extent([118, 122, 28, 32],crs=ccrs.PlateCarree())
request = TDT_img()
ax.add_image(request, 9)
ax.set_title('天地图影像底图',fontsize=15)
gl = ax.gridlines(draw_labels=True, linewidth=1, color='k', alpha=0.5, linestyle='--')
gl.xlabels_top = gl.ylabels_right = False
gl.xformatter = LONGITUDE_FORMATTER
gl.yformatter = LATITUDE_FORMATTER
ax = fig.add_subplot(1, 3, 3, projection=ccrs.PlateCarree())
ax.set_extent([118, 122, 28, 32],crs=ccrs.PlateCarree())
request = TDT_ter()
ax.add_image(request, 9)
ax.set_title('天地图地形底图',fontsize=15)
gl = ax.gridlines(draw_labels=True, linewidth=1, color='k', alpha=0.5, linestyle='--')
gl.xlabels_top = gl.ylabels_right = False
gl.xformatter = LONGITUDE_FORMATTER
gl.yformatter = LATITUDE_FORMATTER
plt.show()
2、不同层级的底图
fig = plt.figure(figsize=(18, 12))
for i in range(6,9):
ax = fig.add_subplot(1, 3, i-5, projection=ccrs.PlateCarree())
ax.set_extent([120.4, 122.1, 30.6, 32.1],crs=ccrs.PlateCarree())
request = TDT_img()
ax.add_image(request, i)
ax.set_title('Level='+str(i),fontsize=15)
gl = ax.gridlines(xlocs=np.arange(120.5, 122.5, 0.5),
ylocs=np.arange(30, 32.5, 0.5),
draw_labels=True,linewidth = 0.5,color='k',
alpha=0.5,linestyle='--')
gl.xlabels_top = gl.ylabels_right = False
plt.show()
参考:
- https://my.oschina.net/u/4581316/blog/4396806
- http://bbs.06climate.com/forum.php?mod=viewthread&tid=89165
- MYSQL官方与主流开源版本基准压测对比
- 基础知识:多人使用一个Github仓库
- 基于Ubuntu16.04的微信小程序搭建过程
- 搭建微信小程序教程:基于CentOS 7.6的教程
- MySQL根据5.1.34之前版本的一个BUG分析O_DIRECT
- Git面试常见问题
- 剑指offer代码解析——面试题17合并两个排序的链表
- MySQL优化方案(一)优化SQL脚本与索引
- 从0学习MySQL系列(二)安装篇
- 从0学习MySQL系列(三)概念篇
- SQL Server基础SQL脚本之主外键约束
- SQL Server基础SQL脚本之Group By
- SQL Server基础SQL脚本之内外连接、交叉连接;函数、子查询
- MySQL 面试选择题15道(单选)
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 『Flutter-绘制篇』实现炫酷的雨雪特效
- 图书管理系统(Servlet+Jsp+Java+Mysql,附下载演示地址)
- Vuex 映射完全指南
- 我们是如何改进YOLOv3进行红外小目标检测的?
- Unity3D网络通讯(三)-- HttpRestful请求的简单封装
- 详解Elasticsearch 的性能优化
- console.log的那件事
- Codeforces Round #668 (Div. 2)A-D
- nodejs源码解析之udp服务器
- LeetCode | 35.搜索插入位置
- Redis | 源码阅读 —— 字符串
- 一个小白的角度看JavaScript Promise 完整指南
- 教会舍友玩 Git (再也不用担心他的学习)
- Python极简入门:数据类型、条件语句、循环语句、异常处理
- 手把手带你入门Python爬虫Scrapy