【GIS开发】批量地图瓦片转mbtiles文件(Python)

2023-10-11 23:30

本文主要是介绍【GIS开发】批量地图瓦片转mbtiles文件(Python),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

文章目录

  • 1、简介
  • 2、下载和编译
  • 3、命令行方式运行
    • 3.1 瓦片图片转mbtiles文件
    • 3.2 mbtiles文件转瓦片图片
  • 4、脚本方式运行
    • 4.1 安装nose库
    • 4.2 编写脚本
    • 4.3 运行脚本
  • 5 加载mbtiles文件
    • 5.1 osgEarth加载mbtiles文件
  • 6、导出mbtiles中的图片
  • 后记

1、简介

MBTiles文件,是指由MapBox制定的一种将瓦片地图数据存储到SQLite数据库中并可快速使用,管理和分享的规范。MBTiles文件内容数据是按照OSGeo的TMS规范来生成和组织的地图瓦片,且支持的图片存储类型为PNG和JPG。MBTiles实质上,是由一个SQLite包装起256*256大小的瓦片地图图片。透过数据库索引的方式提高瓦片索引的效率。据说,这种方式比文件夹方式的瓦片组织要高效得多。

MBTiles的标准在Github或者OpenStreetWiki或者Mapbox主页上可以找到,具体格式不再这里赘述。
https://github.com/mapbox/mbtiles-spec
https://wiki.openstreetmap.org/wiki/MBTiles
https://docs.mapbox.com/help/getting-started/

MBUtil程序库,是用于导入和导出 MBTiles格式的实用程序,通常使用 Mapbox [TileMill](http:// mapbox.com/tilemill/)。

在这里插入图片描述
在这里插入图片描述

2、下载和编译

https://github.com/mapbox/mbutil
https://pypi.org/project/mbutil/
在这里插入图片描述

将代码解压如下:
在这里插入图片描述

在命令行(cmd)运行如下代码,进行mbutil库的下载和安装操作:

git clone git://github.com/mapbox/mbutil.git 
cd mbutil
python setup.py install

在这里插入图片描述
在命令行(cmd)运行如下代码,执行如下语句查询mbutil的使用帮助信息。

python mb-util -h

在这里插入图片描述

3、命令行方式运行

3.1 瓦片图片转mbtiles文件

MBUtil 将元数据导入和导出为 JSON,在 tile 目录的根目录中,保存名为 metadata.json 的文件。

{"name": "World Light","description": "A Test Metadata","version": "3"
}

或者

{"name": "ZYX Tiles","type": "baselayer","description": "","version": "1","format": "png"
}

下载网上的瓦片数据,可以通过本人另外一篇文章:【GIS开发】批量下载和拼接地图瓦片(python)
将下载的瓦片图片数据准备好,如下图所示:
在这里插入图片描述
在这里插入图片描述

在命令行运行如下代码:

# 打开 mbutil 所在文件夹
cd E:\Python\mbutil
# TMS瓦片转换为 .mbutiles文件
python mb-util d:\maps\ d:\test.mbtiles

在这里插入图片描述在这里插入图片描述

3.2 mbtiles文件转瓦片图片

在这里插入图片描述

python mb-util C:\Users\tomcat\Desktop\test\test.mbtiles C:\Users\tomcat\Desktop\test2

第一个参数是python的路径(python.exe)
第二个路径是MBUtil工具的路径(mb-util)
第三个参数是压缩的地图切片路径(.mbtiles)
第四个参数是解压后切片存储的路径(必须是不存在的文件夹)

在这里插入图片描述

4、脚本方式运行

4.1 安装nose库

测试框架nose库(nosetests)。

pip install nose

在这里插入图片描述

4.2 编写脚本

test.py:

import os, shutil
import sys
import json
from nose import with_setup
from mbutil import mbtiles_to_disk, disk_to_mbtilesdef setup_test():print('setup_test 将用于 with_setup')clear_data()def teardown_test():print('teardown_test 也将用于 with_setup')def clear_data():try: shutil.rmtree('D:/test/output')except Exception: pass@with_setup(setup_test, teardown_test)
def test_mbtiles_to_disk():mbtiles_to_disk('D:/test/test1.mbtiles', 'D:/test/output')assert os.path.exists('D:/test/output/0/0/0.png')assert os.path.exists('D:/test/output/metadata.json')@with_setup(setup_test, teardown_test)
def test_mbtiles_to_disk_and_back():# 分解文件mbtiles_to_disk('D:/test/test1.mbtiles', 'D:/test/output')assert os.path.exists('D:/test/output/0/0/0.png')assert os.path.exists('D:/test/output/metadata.json')# 合并文件disk_to_mbtiles('D:/test/output', 'D:/test/output/test2.mbtiles')assert os.path.exists('D:/test/output/test2.mbtiles')# 再分解文件mbtiles_to_disk('D:/test/output/test2.mbtiles', 'D:/test/output/test2')assert os.path.exists('D:/test/output/test2/0/0/0.png')assert os.path.exists('D:/test/output/test2/metadata.json')

4.3 运行脚本

nosetests -v  D:\test\test.py

在这里插入图片描述

5 加载mbtiles文件

5.1 osgEarth加载mbtiles文件

将上文生成的mbtiles文件放在osgEarth的指定数据文件夹里,如下图所示:
在这里插入图片描述
新建一个.earth文件,用来加载mbtiles文件,如下图所示:
在这里插入图片描述
编写批处理命令文件,方便osgEarthViewer.exe运行时参数传入,如下图所示:

在这里插入图片描述
osgEarthViewer.exe运行后,结果如下所示:
在这里插入图片描述
在这里插入图片描述

6、导出mbtiles中的图片

由于mbtiles文件实际就是一个sqlite3数据库文件,所以可以按照数据库的方式导出里面的资源。

import sqlite3def writeTofile(data, filename):# Convert binary data to proper format and write it on Hard Diskwith open(filename, 'wb') as file:file.write(data)print("Stored blob data into: ", filename, "\n")def readBlobData(empId):try:sqliteConnection = sqlite3.connect('d:\\out\\MapWorld.mbtiles')cursor = sqliteConnection.cursor()print("Connected to SQLite")sql_fetch_blob_query = """SELECT * from tiles where tile_row = ?"""cursor.execute(sql_fetch_blob_query, (empId,))record = cursor.fetchall()for row in record:print("Id = ", row[0], "Name = ", row[1])tile_column = row[1]tile_row = row[2]tile_data = row[3]print("Storing employee image and resume on disk \n")photoPath = "d:\\out\\earth" + str(tile_column) + "_"+ str(tile_row) + ".tif"writeTofile(tile_data, photoPath)cursor.close()except sqlite3.Error as error:print("Failed to read blob data from sqlite table", error)finally:if sqliteConnection:sqliteConnection.close()print("sqlite connection is closed")readBlobData(1)
readBlobData(2)

后记

如果你觉得该方法或代码有一点点用处,可以给作者点个赞、赏杯咖啡;╮( ̄▽ ̄)╭
如果你感觉方法或代码不咋地//(ㄒoㄒ)//,就在评论处留言,作者继续改进。o_O???
谢谢各位童鞋们啦( ´ ▽´ )ノ ( ´ ▽´)っ!!!

这篇关于【GIS开发】批量地图瓦片转mbtiles文件(Python)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/191620

相关文章

SpringBoot分段处理List集合多线程批量插入数据方式

《SpringBoot分段处理List集合多线程批量插入数据方式》文章介绍如何处理大数据量List批量插入数据库的优化方案:通过拆分List并分配独立线程处理,结合Spring线程池与异步方法提升效率... 目录项目场景解决方案1.实体类2.Mapper3.spring容器注入线程池bejsan对象4.创建

Python的Darts库实现时间序列预测

《Python的Darts库实现时间序列预测》Darts一个集统计、机器学习与深度学习模型于一体的Python时间序列预测库,本文主要介绍了Python的Darts库实现时间序列预测,感兴趣的可以了解... 目录目录一、什么是 Darts?二、安装与基本配置安装 Darts导入基础模块三、时间序列数据结构与

Python正则表达式匹配和替换的操作指南

《Python正则表达式匹配和替换的操作指南》正则表达式是处理文本的强大工具,Python通过re模块提供了完整的正则表达式功能,本文将通过代码示例详细介绍Python中的正则匹配和替换操作,需要的朋... 目录基础语法导入re模块基本元字符常用匹配方法1. re.match() - 从字符串开头匹配2.

基于 Cursor 开发 Spring Boot 项目详细攻略

《基于Cursor开发SpringBoot项目详细攻略》Cursor是集成GPT4、Claude3.5等LLM的VSCode类AI编程工具,支持SpringBoot项目开发全流程,涵盖环境配... 目录cursor是什么?基于 Cursor 开发 Spring Boot 项目完整指南1. 环境准备2. 创建

Python使用FastAPI实现大文件分片上传与断点续传功能

《Python使用FastAPI实现大文件分片上传与断点续传功能》大文件直传常遇到超时、网络抖动失败、失败后只能重传的问题,分片上传+断点续传可以把大文件拆成若干小块逐个上传,并在中断后从已完成分片继... 目录一、接口设计二、服务端实现(FastAPI)2.1 运行环境2.2 目录结构建议2.3 serv

通过Docker容器部署Python环境的全流程

《通过Docker容器部署Python环境的全流程》在现代化开发流程中,Docker因其轻量化、环境隔离和跨平台一致性的特性,已成为部署Python应用的标准工具,本文将详细演示如何通过Docker容... 目录引言一、docker与python的协同优势二、核心步骤详解三、进阶配置技巧四、生产环境最佳实践

Python一次性将指定版本所有包上传PyPI镜像解决方案

《Python一次性将指定版本所有包上传PyPI镜像解决方案》本文主要介绍了一个安全、完整、可离线部署的解决方案,用于一次性准备指定Python版本的所有包,然后导出到内网环境,感兴趣的小伙伴可以跟随... 目录为什么需要这个方案完整解决方案1. 项目目录结构2. 创建智能下载脚本3. 创建包清单生成脚本4

Python实现Excel批量样式修改器(附完整代码)

《Python实现Excel批量样式修改器(附完整代码)》这篇文章主要为大家详细介绍了如何使用Python实现一个Excel批量样式修改器,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一... 目录前言功能特性核心功能界面特性系统要求安装说明使用指南基本操作流程高级功能技术实现核心技术栈关键函

python获取指定名字的程序的文件路径的两种方法

《python获取指定名字的程序的文件路径的两种方法》本文主要介绍了python获取指定名字的程序的文件路径的两种方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要... 最近在做项目,需要用到给定一个程序名字就可以自动获取到这个程序在Windows系统下的绝对路径,以下

使用Python批量将.ncm格式的音频文件转换为.mp3格式的实战详解

《使用Python批量将.ncm格式的音频文件转换为.mp3格式的实战详解》本文详细介绍了如何使用Python通过ncmdump工具批量将.ncm音频转换为.mp3的步骤,包括安装、配置ffmpeg环... 目录1. 前言2. 安装 ncmdump3. 实现 .ncm 转 .mp34. 执行过程5. 执行结