地图下载核心，如何使用Python下载高德地图瓦片并拼接成大图？

最新推荐文章于 2024-06-12 16:14:43 发布
最新推荐文章于 2024-06-12 16:14:43 发布
阅读量2.2k 收藏 17
点赞数 6
分类专栏： SworldMap 文章标签： python 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/liyajun54/article/details/130474016
版权

地图下载核心，如何使用Python下载高德地图瓦片并拼接成大图？

WMTS地图瓦片

  WMTS (Web Map Tile Service) 是一种提供基于瓦片的地图数据服务的网络协议。它通过将地图数据分割成相同大小的瓦片，并预先生成这些瓦片以供客户端请求，从而提高了地图数据呈现速度和效率。WMTS使用RESTful Web服务和基于HTTP/HTTPS协议进行数据交换，可通过Web浏览器、GIS软件和其他客户端应用程序访问和使用。WMTS广泛用于Web地图、桌面应用程序、移动应用程序等。 它是Open Geospatial Consortium（OGC）发布的开放标准之一。

Requst请求

  下面我们对高德地图瓦片进行分析，瓦片是通过X,Y,Z形式的行列号来请求一张图片的方式获取瓦片数据的，请求返回的数据为一个进制的图片数据，所以我们要建立一个图片请求，并将请求到的图片返回。

在这里插入图片描述

request请求并返回瓦片代码：

def request_tile(self, tile_url):
    headers = {
        "Connection": "keep-alive",
        "User-Agent":
        "Mozilla/5.0 (Windows NT 10.0; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/73.0.3683.86 Safari/537.36",
        "Cookie":
        "HWWAFSESTIME=1681889500446; HWWAFSESID=53a7d799024cfb3f627"
    }

    blank_img = Image.new(mode="RGB", size=(256, 256))
    try:
        r = requests.get(url=tile_url, headers=headers, timeout=20)
        if r.status_code == 200:
        	#将请求的内容转为图片
            img = Image.open(BytesIO(r.content))
            return img
        else:
            return blank_img
    except:
        print("failed")
        return blank_img

影像瓦片与道路瓦片合并

  在我们使用的示例下载地图数据（高德地图）中，影像数据包括：影像地图、标签地图，两者结合才是我们所看到的影像地图，标签地图与影像地图的请求方式一样，为一张256*256的透明背景PNG。

在这里插入图片描述

下载影像地图的时候，我们要将标签地图叠加在影像地图之上，下面代码中vec_tile_url为影像地图地址，cta_tile_url 为标签地图地址，xyz是瓦片的行列号，我们根据传入的行列号与地图级别对瓦片进行下载并重叠。

def get_tile(self, col, row, zoom):
    vec_tile_url = "https://webst01.is.autonavi.com/appmaptile?style=6&x={}&y={}&z={}".format(
        col, row, zoom)
    cta_tile_url = "https://wprd03.is.autonavi.com/appmaptile?x={}&y={}&z={}&lang=zh_cn&size=1&scl=1&style=8&ltype=11".format(
        col, row, zoom)

    vec_img = self.request_tile(vec_tile_url)
    cta_img = self.request_tile(cta_tile_url)

    new_vec_img = Image.new(mode="RGBA", size=(256, 256))
    new_vec_img.paste(vec_img, box=(0, 0))

    new_cta_img = Image.new(mode="RGBA", size=(256, 256))
    new_cta_img.paste(cta_img, box=(0, 0))

    merge_img = Image.alpha_composite(new_vec_img, new_cta_img)
    return merge_img

行列号的计算

  在下面计算瓦片行列号的时候我们传入的参数为经度（l）纬度（b）级别（zoom）,由于我们要下载的瓦片数据为墨卡托投影坐标，所以在计算之前我们要先将经纬度坐标通过lonLat2Mercator方法转换为默卡托投影坐标，grithHalf：20037508.3427892（地球周长的一半）。

# 获取杆塔编号
def getTileNum(self, l, b, zoom):
    num = int(math.pow(2, zoom))
    cellsize = self.grithHalf * 2 / num
    np = self.lonLat2Mercator(l, b)
    x = math.floor((np.x() + self.grithHalf) / cellsize)
    x = max(0, x)
    x = min(x, int(num - 1))
    y = math.floor((self.grithHalf - np.y()) / cellsize)
    y = max(0, y)
    y = min(y, int(num - 1))
    return x, y
#经纬度转墨卡托投影
def lonLat2Mercator(self, l, b):
    if (b > 85.0511287798066):
        b = 85.0511287798066
    if (b < -85.0511287798066):
        b = -85.0511287798066

    x = l * self.grithHalf / 180
    y = math.log(math.tan((90 + b) * math.pi / 360)) / (math.pi / 180)
    y = y * self.grithHalf / 180
    return QPointF(x, y)

根据给定坐标范围下载地图

  在下载地图时，我们要先给定一个坐标范围，这个范围为地图的左下角坐标及右上角坐标，如:[99,20,105,40]；然后通过范围计算出我们瓦片在X方向上的最小瓦片编号与最大瓦片编号，以及在Y方向上的最小瓦片编号与最大瓦片编号，这样我们就可以得到整个范围的瓦片编号，然后对瓦片进行请求合并导出。

在这里插入图片描述

  由于瓦片编号在X方向是由左向右递增，在Y方向上是由上到下递增，所以我们在取瓦片X与Y的最大最小编号时要以左上角和右下角的坐标来进行瓦片编号计算。

def download_img(self,minx,miny,maxx,maxy,zoom):
	#左上角
    nx1, ny1 = self.getTileNum(minx, maxy, zoom)
    #右下角
    nx2, ny2 = self.getTileNum(maxx, miny, zoom)

    print(nx1, ny1, zoom)
    print(nx2, ny2, zoom)

    tile_w = int(qAbs(nx2 - nx1) * 256 + 256)
    tile_h = int(qAbs(ny2 - ny1) * 256 + 256)
    new_img = Image.new(mode="RGBA", size=(tile_w, tile_h))

    #瓦片总数
    total_count = (nx2 - nx1 + 1) * (ny2 - ny1 + 1)
    current_count = 0
    i = 0
    for col in range(nx1, nx2 + 1):
        j = 0
        for row in range(ny1, ny2 + 1):
            t_img = self.get_tile(col, row, zoom)
            new_img.paste(t_img, box=(i * 256, j * 256))
            current_count += 1
            progress = qAbs(current_count / total_count) * 100
            print("\r Downloading File Map [{:.0f}%] ".format(progress),
                    "\r",
                    end='')
            j += 1
        i += 1

    return new_img

导出图片

mapImage = download_img(103.563246,29.632548,103.686545,29.663265,16)
mapImage.save("D:\\crop_result.png", "PNG")

结果图片展示

在这里插入图片描述

完整代码：

import math
from PyQt5.QtGui import *
from PyQt5.QtCore import *
from PyQt5.QtWidgets import *
import requests
from PIL import Image, ImageDraw, ImageFont
from io import BytesIO
import numpy as np

class MapDownloader(object):
    def __init__(self):
        self.grithHalf = 20037508.3427892

    #经纬度转墨卡托投影
    def lonLat2Mercator(self, l, b):
        if (b > 85.0511287798066):
            b = 85.0511287798066
        if (b < -85.0511287798066):
            b = -85.0511287798066

        x = l * self.grithHalf / 180
        y = math.log(math.tan((90 + b) * math.pi / 360)) / (math.pi / 180)
        y = y * self.grithHalf / 180
        return QPointF(x, y)

    # 获取杆塔编号
    def getTileNum(self, l, b, zoom):
        num = int(math.pow(2, zoom))
        cellsize = self.grithHalf * 2 / num
        np = self.lonLat2Mercator(l, b)
        x = math.floor((np.x() + self.grithHalf) / cellsize)
        x = max(0, x)
        x = min(x, int(num - 1))
        y = math.floor((self.grithHalf - np.y()) / cellsize)
        y = max(0, y)
        y = min(y, int(num - 1))
        return x, y
    
    def download_img(self,minx,miny,maxx,maxy,zoom):
        nx1, ny1 = self.getTileNum(minx, maxy, zoom)
        nx2, ny2 = self.getTileNum(maxx, miny, zoom)

        print(nx1, ny1, zoom)
        print(nx2, ny2, zoom)

        tile_w = int(qAbs(nx2 - nx1) * 256 + 256)
        tile_h = int(qAbs(ny2 - ny1) * 256 + 256)
        new_img = Image.new(mode="RGBA", size=(tile_w, tile_h))

        #瓦片总数
        total_count = (nx2 - nx1 + 1) * (ny2 - ny1 + 1)
        current_count = 0
        i = 0
        for col in range(nx1, nx2 + 1):
            j = 0
            for row in range(ny1, ny2 + 1):
                t_img = self.get_tile(col, row, zoom)
                new_img.paste(t_img, box=(i * 256, j * 256))
                current_count += 1
                progress = qAbs(current_count / total_count) * 100
                print("\r Downloading File Map [{:.0f}%] ".format(progress),
                      "\r",
                      end='')
                j += 1
            i += 1

        return new_img

    def get_tile(self, col, row, zoom):
        vec_tile_url = "https://webst01.is.autonavi.com/appmaptile?style=6&x={}&y={}&z={}".format(
            col, row, zoom)
        cta_tile_url = "https://wprd03.is.autonavi.com/appmaptile?x={}&y={}&z={}&lang=zh_cn&size=1&scl=1&style=8&ltype=11".format(
            col, row, zoom)

        vec_img = self.request_tile(vec_tile_url)
        cta_img = self.request_tile(cta_tile_url)

        new_vec_img = Image.new(mode="RGBA", size=(256, 256))
        new_vec_img.paste(vec_img, box=(0, 0))

        new_cta_img = Image.new(mode="RGBA", size=(256, 256))
        new_cta_img.paste(cta_img, box=(0, 0))

        merge_img = Image.alpha_composite(new_vec_img, new_cta_img)
        return merge_img

    def request_tile(self, tile_url):
        headers = {
            "Connection": "keep-alive",
            "User-Agent":
            "Mozilla/5.0 (Windows NT 10.0; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/73.0.3683.86 Safari/537.36",
            "Cookie":
            "HWWAFSESTIME=1681889500446; HWWAFSESID=53a7d799024cfb3f627"
        }

        blank_img = Image.new(mode="RGB", size=(256, 256))
        try:
            r = requests.get(url=tile_url, headers=headers, timeout=20)
            # print(r.status_code)
            if r.status_code == 200:
                img = Image.open(BytesIO(r.content))
                return img
            else:
                return blank_img
        except:
            print("failed")
            return blank_img
        
mapDownloader = MapDownloader()
mapImage = mapDownloader.download_img(103.563246,29.632548,103.686545,29.663265,16)
mapImage.save("D:\\crop_result.png", "PNG")

那么，这里一个可以按照地图范围与级别下载地图的方法就做好了；如果大家觉得我的方法还有用的话，也请大家能够点赞鼓励一下，您的鼓励是我分享的动力，谢谢各位！

四维云测
关注
  • 6
    点赞
  • 踩
  • 17
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
地图瓦片下载地址汇总（leaflet 百度地图等）
博客简介
11-16 3685
地图瓦片下载地址汇总（leaflet 百度地图等）
Google 谷歌地图（高德地图）瓦片下载/缓存，使用Python实现、arcgis.js（三）
qq_45062700的博客
04-19 936
通过arcgis for JavaScript 技术与flask、requests实现地图在本地服务器进行缓存与显示。 高德地图瓦片缓存、谷歌地图（google）
Google 谷歌地图（高德地图）瓦片下载/缓存，使用Python实现、自定义上传地图瓦片资源（四）
qq_45062700的博客
05-06 868
自定义上传已经下载好的地图瓦片资源压缩包，可以合并类型，也可以创建新的地图类型。使用arcgis for JavaScript，python flask后端，Google 谷歌地图，高德地图
百度地图瓦片下载地址
最新发布
chenchuang0128的博客
06-12 500
请注意，这些网站可能会随时间变化而更新或改变服务内容，因此建议直接访问上述网址以获取最新的信息和操作指导。
高德地图内网加载地图瓦片、如何下载地图瓦片资源
weixin_44640245的博客
10-29 5813
高德地图离线（内网）加载地图瓦片、如何下载地图瓦片资源
高德、百度瓦片地图免费下载
子建莫敌
05-26 770
内嵌高德、百度离线地图，都需要瓦片地图。然后免费、好用的下载工具太少。
高德谷歌腾讯天地图地图瓦片url
热门推荐
吴道简的博客
05-06 2万+
以下是各地图的瓦片地图地址，可以在QGIS中加载，也可以在openlayer或leaflet前端页面中加载。 腾讯、百度、高德，加上谷歌和天地图应该足够用了。 腾讯、高德是GCJ02坐标系，百度是BD09坐标系，谷歌、必应是WGS84坐标系，天地图是CGCS2000坐标系，瓦片地图都是平面墨卡托投影。WGS84和CGCS2000坐标系，近似认为它们相等就可以了。 高德路网： https://w...
全功能瓦片地图下载（百度、高德、谷歌、腾讯、微软等等）.zip
11-14
全功能瓦片地图下载（百度、高德、谷歌、腾讯、微软等等），亲测可用 解压密码123456 解压密码123456
Python：多线程批量下载高德地图最新的省、市、县行政边界地图数据
guoyc439的博客
06-30 1414
用到mapinfo制作地图数据经常因为行政区划图的变更而烦恼，因此，获取实时最新的行政区划图边界地理数据非常有必要。通过调用上述自定义的类，就可以自动下载所有省、市、线边界的json数据，并应用geopandas将其转换成mapinfo、arcgis可应用的shp格式
Google（高德）地图瓦片python代码下载
04-25
Google（高德）地图瓦片使用python工具类代码下载，包含经纬度转换为瓦片x,y轴。使用简单方便，python3.x版本，使用requests进行下载，可以自定义保存路径，自己设置proxy代理端口，实现谷歌地图瓦片的加载。
python下载谷歌瓦片地图
11-02
6. 可能还会涉及处理下载的瓦片，如拼接成大图、裁剪或转换格式等操作。 为了实现这个功能，开发者需要对Web服务API有一定了解，同时掌握Python的网络请求和文件操作。此外，熟悉GIS概念，如投影、瓦片系统和地图...
3.（地图数据篇）高德地图、google、openstreenmap地图瓦片数据爬取--java代码.zip
07-22
下载如有问题，可私信博主。下载前建议先查看博客内容，其地址为：https://blog.csdn.net/QQ98281642/article/details/116652717
高德地图瓦片下载，只针对高德
03-17
高德地图瓦片下载，只针对高德
单线程与多线程python爬虫地图瓦片源码
06-17
标题中的“单线程与多线程Python爬虫”涉及的核心是如何利用线程技术优化地图瓦片的下载过程。单线程爬虫适用于小型项目，其代码简单易懂，但当需要下载大量瓦片时，速度可能较慢。多线程爬虫则可以显著提升下载速度...
python实现高德地图
07-24
本文将深入探讨如何使用Python与高德地图API交互，包括定位、路径规划、地图展示等核心知识点。我们将涵盖以下几个方面： 1. **高德地图API介绍**： 高德地图API提供了一系列接口，用于获取地理位置信息、进行路线...
TileFetchor地图瓦片下载器
10-06
TileFetchor是用Python编写的互联网地图瓦片下载工具，其中支持多种互联网地图服务器，用户可以自行添加下载地址URL。下载代码后，直接运行main.py，具体参数参见main.py中的usage说明。用户可以进github下载：...
tiler，一个超级厉害的 Python 库！
m0_59236602的博客
01-31 895
我们学习Python必然是为了找到高薪的工作，下面这些面试题是来自阿里、腾讯、字节等一线互联网大厂最新的面试资料，并且有阿里大佬给出了权威的解答，刷完这一套面试资料相信大家都能找到满意的工作。Python所有方向的技术点做的整理，形成各个领域的知识点汇总，它的用处就在于，你可以按照上面的知识点去找对应的学习资源，保证自己学得较为全面。观看零基础学习视频，看视频学习是最快捷也是最有效果的方式，跟着视频中老师的思路，从基础到深入，还是很容易入门的。最后，可以将生成的地图瓦片显示在应用程序或网页中。
C++ 如何根据地理坐标范围获取瓦片地图并使用CImage库实现多张图片（瓦片地图）的快速合并
编辑编辑器
05-25 1775
作者：虚坏叔叔 博客：https://xuhss.com 早餐店不会开到晚上，想吃的人早就来了！😄 C++ 如何获取瓦片地图并使用CImage库实现多张图片（瓦片地图）的快速合并 一、C++ 如何获取瓦片数据 C++ 如何获取瓦片数据，这里通过向谷歌地图发起http请求，获取多个地图保存到本地，部分关键源码如下所示 int CGenerateJpgAccessImpl::FetchTiles(int zoomLevel, int x, int y, int width, int height, int
android高德地图加载瓦片,高德地图开发（二）加载瓦片数据
weixin_34338504的博客
05-28 1237
package com.example.maptest;import android.app.Activity;import android.os.Environment;import android.support.v7.app.AppCompatActivity;import android.os.Bundle;import android.view.View;import android.w...
python爬取高德地图xyz瓦片
02-03
要使用Python爬取高德地图的XYZ瓦片，可以使用以下步骤进行操作： 1. 导入必要的库和模块，包括`requests`用于发送HTTP请求，`Pillow`用于处理图像数据。 2. 构建请求URL。高德地图的XYZ瓦片URL格式为：`http://wprd03.is.autonavi.com/appmaptile?style=6&x={x}&y={y}&z={z}`，其中`{x}`表示瓦片的x坐标，`{y}`表示y坐标，`{z}`表示缩放级别。 3. 使用循环遍历需要下载的瓦片的x和y坐标，并发送GET请求获取瓦片数据。可以根据具体情况设定循环范围和缩放级别。 4. 将获取到的瓦片数据保存为图片文件。可以使用`Pillow`库的`Image`模块将数据解码为图片，然后保存到本地文件夹中。 下面是一个简单的Python代码示例： ```python import requests from PIL import Image # 构建请求URL url_template = "http://wprd03.is.autonavi.com/appmaptile?style=6&x={x}&y={y}&z={z}" # 设置需要下载瓦片的范围和缩放级别 min_x = 0 max_x = 10 min_y = 0 max_y = 10 zoom_level = 10 # 循环遍历需要下载的瓦片 for x in range(min_x, max_x+1): for y in range(min_y, max_y+1): # 构建请求URL url = url_template.format(x=x, y=y, z=zoom_level) # 发送GET请求获取瓦片数据 response = requests.get(url) # 将获取到的瓦片数据保存为图片文件 file_name = f"tile_{zoom_level}_{x}_{y}.png" with open(file_name, "wb") as f: f.write(response.content) print("瓦片下载完成。") ``` 上述代码将会下载高德地图的XYZ瓦片，保存为相应的图片文件，图片文件的命名格式为`tile_{缩放级别}_{x坐标}_{y坐标}.png`。根据实际需求，可以调整代码中的循环范围、URL模板、保存文件的路径等。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值