使用Python多线程下载卫星地图

import math
import os
import random
import requests
import time
import random
import threading

经纬度反算切片行列号 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)

zoom = 17
x1,y2=deg2num(32.2359,89.3848,zoom)
x2,y1=deg2num(42.0821,106.8311,zoom)
print(x1,x2,y1,y2)

def xiazai(url,path):
head = {‘user-agent’:‘Mozilla/5.0 (Windows NT 10.0; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/107.0.0.0 Safari/537.36 Edg/107.0.1418.42’}
with open(path, “wb”) as f:
img=requests.get(url=url,headers=head)
f.write(img.content)
print(path)

def xia(url_1,path_1,i):
for j in range(len(url_1[i])):
xiazai(url_1[i][j],path_1[i][j])

z_url = []
z_path =[]
for x in range(x1,x2+1):

for y in range(y1,y2+1):
    path = "D:\\GoogleImg\\" + str(zoom) + "\\" + str(x)
    if not os.path.exists(path):
        os.makedirs(path)
    Tpath = 'https://webst02.is.autonavi.com/appmaptile?style=6&x={0}&y={1}&z={2}'.format(x,y,zoom)
    z_url.append(Tpath)
    Spath="D:\\GoogleImg\\{}\\{}\\{}.jpg".format(zoom,x,y)
    z_path.append(Spath)

index = len(z_url)//50
url_1 = []
path_1 = []
for i in range(0,len(z_url),index):
b = z_url[i:i+index]
c = z_path[i:i+index]
url_1.append(b)
path_1.append©

thread = []
for i in range(len(url_1)):
t = threading.Thread(target=xia,args=(url_1,path_1,i))
thread.append(t)
for t in thread:
t.start()
for t in thread:
t.join()