【GIS开发】批量地图瓦片转mbtiles文件(Python)
文章目录1、简介2、下载和编译3、运行后记1、简介MBUtil 是用于导入和导出 MBTiles格式的实用程序,通常使用 Mapbox [TileMill](http:// mapbox.com/tilemill/)。2、下载和编译https://github.com/mapbox/mbutilhttps://pypi.org/project/mbutil/将代码解压如下:在命令行运行如下代码:g
文章目录
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_mbtiles
def 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 sqlite3
def writeTofile(data, filename):
# Convert binary data to proper format and write it on Hard Disk
with 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???
谢谢各位童鞋们啦( ´ ▽´ )ノ ( ´ ▽´)っ!!!
开放原子开发者工作坊旨在鼓励更多人参与开源活动,与志同道合的开发者们相互交流开发经验、分享开发心得、获取前沿技术趋势。工作坊有多种形式的开发者活动,如meetup、训练营等,主打技术交流,干货满满,真诚地邀请各位开发者共同参与!
更多推荐
所有评论(0)