多进程下载谷歌地图和切片地图合并

1、下载代码

from urllib import request
import re
import urllib.request
import os
import random
import math
import multiprocessing
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']


# 经纬度反算切片行列号 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)
    return (xtile, ytile)


# 下载图片
def getimg(Tpath, Spath, x, y):
    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)

if __name__ == "__main__":
    for zoom in range(14,20):
        # zoom = 19  # 下载切片的zoom
        path = r"E:\地图切片\\"+str(zoom)
        if not os.path.exists(path):
            os.mkdir(path)
        lefttop = deg2num(27.016010, 102.446480, zoom)  # 下载切片的左上角角点
        rightbottom = deg2num(26.822850, 102.675480, zoom)

        print(str(lefttop[0]))
        print(str(rightbottom[0]))
        print(str(lefttop[1]))
        print(str(rightbottom[1]))
        print("共" + str(lefttop[0] - rightbottom[0]))
        print("共" + str(lefttop[1] - rightbottom[1]))

        cpus = multiprocessing.cpu_count()
        # 定义进程数量，设置为cpus核心的20倍
        p = multiprocessing.Pool(cpus *30)
        results = 0
        for x in range(lefttop[0], rightbottom[0]):
            for y in range(lefttop[1], rightbottom[1]):
                # https: // mt3.google.cn / maps / vt?lyrs = y & gl = cn & x = 205719 & y = 110741 & z = 1
                tilepath = "https://mt3.google.cn/maps/vt?lyrs=y&gl=cn&x=" + str(x) + "&y=" + str(y) + "&z=" + str(zoom)
                print(tilepath)
                # 天地图的url 可以换成谷歌地图的
                # getimg(tilepath, os.path.join(path, str(x) + "_" + str(y) + ".png"), x, y)
                p.apply_async(getimg, args=(tilepath,os.path.join(path, str(x) + "_" + str(y) + ".png"), x, y))
                results=results+1
                total=(lefttop[1]-lefttop[0])*(rightbottom[1]-rightbottom[0])
                print('完成'+str(results/total*100)+'%')
        p.close()
        p.join()

2、合并代码

import glob
import re
from PIL import Image
path='E:/地图切片/14/'
files = glob.glob(path+'*.png')
x_offset = 0
y_offset = 0
x1=int(re.findall(r'\d{3,9}',files[0])[0])
x2=int(re.findall(r'\d{3,9}',files[len(files)-1])[0])
y1=int(re.findall(r'\d{3,9}',files[0])[1])
y2=int(re.findall(r'\d{3,9}',files[len(files)-1])[1])
new_image = Image.new('RGB', (256*(x2-x1), 256*(y2-y1)))
for x in range(x1, x2+1):
    for y in range(y1, y2+1):
        subitem=path+'%d_%d.png' % (x,y)
        print(subitem)
        im = Image.open(subitem)
        new_image.paste(im, (x_offset, y_offset))
        x_offset += 0
        y_offset += 256
    x_offset += 256
    y_offset = 0
new_image.save(path+'merge.jpg', quality = 90)