[python]天地地图瓦片下载

之前下载了一个.net GMap来下载地图发现速度太慢。而且经常会闪退。所以写了个脚本去获取。（亲测可用，直接复制吧。）

#!/usr/bin/python
# _*_coding:utf-8_*_

from urllib3 import request

import re

import urllib.request

import os

import random

import math

import sys

count = 0

agents = [

'Mozilla/5.0 (Windows NT 6.1; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/60.0.3112.101 Safari/537.36',
'Mozilla/5.0 (Windows; U; Windows NT 6.1; en-US) AppleWebKit/532.5 (KHTML, like Gecko) Chrome/4.0.249.0 Safari/532.5'
'Mozilla/5.0 (Windows; U; Windows NT 5.2; en-US) AppleWebKit/532.9 (KHTML, like Gecko) Chrome/5.0.310.0 Safari/532.9',
'Mozilla/5.0 (Windows; U; Windows NT 5.1; en-US) AppleWebKit/534.7 (KHTML, like Gecko) Chrome/7.0.514.0 Safari/534.7',
'Mozilla/5.0 (Windows; U; Windows NT 6.0; en-US) AppleWebKit/534.14 (KHTML, like Gecko) Chrome/9.0.601.0 Safari/534.14',
'Mozilla/5.0 (Windows; U; Windows NT 6.1; en-US) AppleWebKit/534.14 (KHTML, like Gecko) Chrome/10.0.601.0 Safari/534.14',
'Mozilla/5.0 (Windows; U; Windows NT 6.1; en-US) AppleWebKit/534.20 (KHTML, like Gecko) Chrome/11.0.672.2 Safari/534.20',
'Mozilla/5.0 (Windows NT 6.1; WOW64) AppleWebKit/534.27 (KHTML, like Gecko) Chrome/12.0.712.0 Safari/534.27',
'Mozilla/5.0 (Windows NT 6.1; WOW64) AppleWebKit/535.1 (KHTML, like Gecko) Chrome/13.0.782.24 Safari/535.1']
resolution={18:5.36441802978515E-06,
17:1.07288360595703E-05,
16:2.1457672119140625E-05,
15:4.29153442382814E-05,
14:8.58306884765629E-05,
13:0.000171661376953125,
12:0.00034332275390625,
11:0.0006866455078125,
10:0.001373291015625,
9:0.00274658203125,
8:0.0054931640625,
7:0.010986328125,
6:0.02197265625,
5:0.0439453125,
4:0.087890625,
3:0.17578125,
2:0.3515625,
1:0.703125}

# 经纬度反算切片行列号 3857坐标系

def deg2num(lat_deg, lon_deg, zoom):
lat_rad = math.radians(lat_deg)
n = 2.0 ** zoom
xtile = int((lon_deg+180.0) / 360.0 * n)
ytile = int((1.0 - math.log(math.tan(lat_rad)+(1 / math.cos(lat_rad))) / math.pi) / 2.0 * n)
# xtile = math.floor((lon_deg+180.0) / (resolution[zoom]*256))
# ytile = math.floor((90.0-lat_deg) / (resolution[zoom]*256))
return (xtile, ytile)
# 下载图片
def getimg(Tpath, Spath, x, y):
global count
count = 0 #清零
try:
f = open(Spath, 'wb')
req = urllib.request.Request(Tpath)
req.add_header('User-Agent', random.choice(agents)) # 换用随机的请求头
pic = urllib.request.urlopen(req, timeout=60)
f.write(pic.read())
f.close()
print(str(x),'_',str(y),'下载成功')
except Exception:
print(str(x),'_',str(y),'下载失败,重试')
getimg(Tpath, Spath, x, y)
count = count + 1
print(count)
if count > 5:
return

#zoom =16 # 下载切片的zoom
#lefttop = deg2num( 28.317231,115.846340, zoom) # 下载切片的左上角角点
#rightbottom = deg2num( 23.567119,120.722095, zoom)

#def GetMapTile(lat1, lon1, lat2, lon2):
# """根据经纬度范围获取背景"""
# zooms = []
# for i in range(1, 19):
# l = deg2num(lat1, lon1, i)
# r = deg2num(lat2, lon2, i)
# if l[0] - r[0] == 0 or l[1] - r[1] == 0:
# continue
# else:
# zooms.append(i)
# 根据有瓦片数据的地方下载瓦片
#zoom = zooms[-1]
# 根据经纬度确定瓦片位置
#lefttop = deg2num(lat1, lon1, zoom)
#rightbottom = deg2num(lat2, lon2, zoom)

def download(k,LTlat,LTlon,RBlat,RBlon):
zoom = k
LT_lat = LTlat
LT_lon = LTlon
RB_lat = RBlat
RB_lon = RBlon
lefttop = deg2num(LT_lat,LT_lon, zoom)
rightbottom = deg2num(RB_lat,RB_lon, zoom)
print(str(zoom))
print(str(lefttop[0]))
print(str(rightbottom[0]))
print(str(lefttop[0]))
print(str(rightbottom[0]))
print("共" + str(lefttop[0] - rightbottom[0]))
print("共" + str(lefttop[1] - rightbottom[1]))
#rootDir = "/root/ditupython/World_Annotation/TiandituSatelliteMap/"
rootDir = "D:\ditu"
for x in range(lefttop[0], rightbottom[0]):
path = rootDir + str(zoom) + "/" + str(x)
if not os.path.exists(path):
os.makedirs(path)
for y in range(lefttop[1], rightbottom[1]):
#天地图-影像
#tilepath = "http://t0.tianditu.gov.cn/DataServer?T=img_w&x=" + str(x) + "&y=" + str(y) + "&l=" + str(zoom) + "&tk=28b495e4df789d971d2ae77b01a55a55"
tilepath = 'http://t0.tianditu.gov.cn/DataServer?T=img_w&x=' + str(x) + '&y=' + str(y) + '&l=' + str(zoom) + '&tk=28b495e4df789d971d2ae77b01a55a55'
#tilepath = "http://t6.tianditu.gov.cn/DataServer?T=cia_w&x=" + str(x) + "&y=" + str(y) + "&l=" + str(zoom) + "&tk=5d22d49fdc586cb5caed68bfb12d1e6b"
filepath = path + "/" + str(y) + ".png"
if not os.path.isfile(filepath):
getimg(tilepath,os.path.join(path,str(y) + ".png"),x,y)

'''
for x in range(lefttop[0], rightbottom[0]):
for y in range(lefttop[1], rightbottom[1]):
path = "/root/ditupython/World_Annotation/" + str(zoom) + "/" + str(x)
if not os.path.exists(path):
os.makedirs(path)
# tilepath = "http://t0.tianditu.gov.cn/DataServer?T=cia_c&x=" + str(x) + "&y=" + str(y) + "&l=" + str(zoom) + "&tk=28b495e4df789d971d2ae77b01a55a55"
tilepath = "http://t0.tianditu.gov.cn/DataServer?T=img_w&x=" + str(x) + "&y=" + str(y) + "&l=" + str(zoom) + "&tk=5d22d49fdc586cb5caed68bfb12d1e6b"
# 天地图的url 可以换成谷歌地图的
getimg(tilepath, path + "\\" + str(y) + ".png", x, y)
'''

print('完成')
arg1 = sys.argv[1]
arg2 = sys.argv[2]
#湖南
#LT_lat = 30.2
#LT_lon = 108.7
#RB_lat = 24.5
#RB_lon = 114.7

#全国
LT_lat = 53.64464
LT_lon = 73.82812
RB_lat = 24.53713
RB_lon = 114.78516

level1 = int(arg1)
level2 = int(arg2)
for i in range(level1,level2):
delta_lat = LT_lat - RB_lat
delta_lon = RB_lon - LT_lon
if i > 15:
print(str(i))
LT_lat = LT_lat -delta_lat * 1/4
LT_lon = LT_lon + delta_lon * 1/4
RB_lat = RB_lat + delta_lat * 1/4
RB_lon = RB_lon - delta_lon * 1/4
print(str(LT_lat))
print(str(LT_lon))
print(str(RB_lat))
print(str(RB_lon))
download(i,LT_lat,LT_lon,RB_lat,RB_lon)

[python]天地地图瓦片下载相关推荐

谷歌地图离线地图瓦片下载_如何下载Google地图以供离线使用
谷歌地图离线地图瓦片下载 If you've ever wanted to be able to download Google Maps data for offline use, you shou ...
离线版百度地图瓦片下载html版
离线版百度地图瓦片下载html版:http://download.csdn.net/detail/wyazyf/9858560
离线地图瓦片下载高德地图谷歌地图天地图 Mapbox地图矢量卫星地图下载
用于地图瓦片下载离线,支持高德地图瓦片谷歌地图 Mapbox bingMap OpenStreeMap 自定义瓦片地址下载 gitee项目地址 https://gitee.com/houlaid ...
1.Electron-封装地图瓦片下载小工具(MapDownload) [配置开发环境]
本软件完成状况未完成~~ 目的之前写了一个地图瓦片下载小工具地图地图瓦片下载器(有兴趣可以去看看)由于是"黑框"版本的.想把它封装成界面版本.接下来我会用一段时间把我学习El ...
基于GMap.Net的地图瓦片下载工具
关键字:GMap.Net.地图瓦片下载.高德瓦片该地图瓦片下载器是根据GMap.Net源码改编,主要以高德地图源作以演示,提供了矩形.圆形.多边形.图标.报警图标.图表.测距等元素的使用与编辑:提供 ...
python下载谷歌地图瓦片_使用 Python 合并地图瓦片
前文提到了合并瓦片图,而瓦片图应用比较多的则是瓦片地图.对地图本就感兴趣的我,也想试试合并互联网地图的某个范围内的地图图层. 随着技术的发展,国内的地图服务商相继将地图瓦片更新为矢量瓦片[1],这下想 ...
OSM地图瓦片下载器1.0版介绍（win64）
简介为方便在工作中随时使用OSM瓦片数据(公开链接的),特编写此下载工具,并会一直更新,欢迎关注.如果需要了解基本的瓦片知识,请移步地图瓦片讲解注意OSM坐标系是WGS84-webMercato ...
百度离线地图开发，node实现地图瓦片下载
最近有个Web项目要用离线地图,项目是在内网环境,找了很多资料,踩了很多坑,好在已经实现了,下面把资料整理一下. 首先是百度离线地图开发包,原文地址是http://www.xiaoguo123.com ...
轻松实现离线地图-离线地图-地图瓦片下载
此demo为前端页面离线地图,不连外网,不用与后端交互,瓦片下载的全就显示的全: 代码示例下载: 链接:https://pan.baidu.com/s/1We86NFDNBCoD9JPz3z_iyw ...

[python]天地地图瓦片下载

[python]天地地图瓦片下载相关推荐

最新文章

热门文章