GIS开发:tianditu切片计算

简介: GIS开发:tianditu切片计算

如何根据地图的范围计算其中切片的行列号,各类的地图引擎能加载的,是能计算并显示对应的图片,地图的下载工具,也是计算了切片的等级和行列号,再下载对应的图片,按照顺序存储在本地或者数据库里。

将切片拼接在一起,加上坐标的信息,就变成地图影像了,当然最初可能切片也是根据地图影像,来切割出来的。

下边分享一个简单的tianditu切片的行列号计算方法,对应的是tainditu的“经纬度”地图切片,当然球面墨卡托的计算方式也差不多。

两种地图:
image.png
tianditu经纬度的切片地图范围为[-180,-90,180,90],切片地图的原点是-180,90,地图级别从1开始,第一级别的地图切片尺寸为横向512px,纵向为256px,由两张256*256的切片组成,第1级地图的分辨率为360/512,地图每增加一级,分辨率对应缩小到上一级的1/2。

切片序号的顺序是从原点处开始的,越靠近地图原点,编号就越小。

假如给定一个地图的经纬度范围[xmin,ymin,xmax,ymax],地图的第一级分辨率定义initResolution,要获取地图切片的级别为level。

获取地图当前级别的最小行号和列号:

列号

= ((xmin- (-180)) / (initResolution/2^(level-1)))/256

行号

=-1* ((ymax- 90) / (initResolution/2^(level-1))+0.5)/256

获取地图当前级别的最大行号和列号:

列号

= ((xmax- (-180)) / (initResolution/2^(level-1))+0.5)/256

行号

=-1* ((ymin- 90) / (initResolution/2^(level-1)))/256

256分别为切片的宽高像素值,如有不同,需要根据切片实际尺寸计算,坐标计算的时候,需要减去中心点的值,再根据级别计算当前级别的分辨率。

这样根据最大、最小行列号算出一个数字的范围,上边的值都要转换成int类型的,再加上请求的级别,就可以获取当前的切片数据了。

相关文章
|
25天前
|
存储 数据处理 C语言
NumPy 通用函数(ufunc):高性能数组运算的利器
NumPy的通用函数(ufunc)提供高性能的逐元素运算,支持向量化操作和广播机制,能应用于数组的数学、逻辑和比较运算。ufunc可提高计算速度,避免低效的循环,并允许自定义函数以满足特定需求。例如,ufunc实现加法比循环更高效。通过`frompyfunc`可创建自定义ufunc。判断函数是否为ufunc,可检查其类型是否为`numpy.ufunc`。ufunc练习包括数组的平方、平方根、元素积及性能对比。
22 0
|
2月前
|
存储 数据挖掘 数据处理
地理空间数据分析与NumPy的高效融合
【4月更文挑战第17天】本文探讨了地理空间数据与NumPy的高效结合在数据分析中的应用。地理空间数据具有多维性、空间相关性和复杂性,NumPy通过提供高效的数据存储、强大的数学运算及灵活的数据处理,应对这些挑战。实践案例展示了如何使用NumPy进行坐标处理、统计分析和滤波等,揭示数据规律。未来,随着数据增长和领域拓展,NumPy在地理空间数据分析中的作用将更加显著,期待更多创新方法应对新挑战。
|
12月前
|
编解码 算法 定位技术
GIS开发:切片格式说明(翻译)-wgs84
GIS开发:切片格式说明(翻译)-wgs84
110 0
|
12月前
|
存储 编解码 算法
GIS开发:切片格式说明(翻译)-墨卡托
GIS开发:切片格式说明(翻译)-墨卡托
134 0
|
存储 JavaScript 前端开发
GIS开发:简单的数据差值
GIS开发:简单的数据差值
122 0
GIS开发:简单的数据差值
|
前端开发 定位技术 Python
GIS开发:分享osm切片计算
GIS开发:分享osm切片计算
194 0
|
XML 存储 定位技术
GIS开发:Arcgis的切片格式
GIS开发:Arcgis的切片格式
166 0
|
存储 Java 定位技术
GIS开发:mbtiles数据合并
GIS开发:mbtiles数据合并
294 0
|
定位技术 API 数据库
GIS开发:开源gdal切片
GIS开发中开源gdal切片使用
340 0
|
算法 大数据
在电子表格计算架构上应用稀疏数组技术的设计
在电子表格计算架构上应用稀疏数组技术的设计
62 0
在电子表格计算架构上应用稀疏数组技术的设计