Gis 将MapInfo选项卡转换为MB坐标

我有一个MapInfo选项卡文件，如下所示：

!table
!version 300
!charset WindowsLatin1

Definition Table
  File "hp.tif"
  Type "RASTER"
  (400000,1300000) (0,0) Label "TL",
  (500000,1300000) (10000,0) Label "TR",
  (500000,1200000) (10000,10000) Label "BR",
  (400000,1200000) (0,10000) Label "BL"
  CoordSys Earth Projection 8, 79, "m", -2, 49, 0.9996012717, 400000, -100000
  Units "m"

我还有一个参考区域的图像文件。我正试图将其转换为一个mbtiles文件，该文件几乎可以正常工作，但有点不对齐

我从这个文件中获得了边界框，它是

400000，1200000，500000，1300000

，然后我将它转换成纬度/longtude，它是

-2.0018997431234，60.67074210511，-0.17172476489587，61.580961113449

然后，我使用以下代码使用zoom 12将其转换为MBTiles栅格：

$x = (int)floor((($longitude + 180) / 360) * pow(2, $zoom))
$y = pow(2, $zoom) - (int)floor((1 - log(tan(deg2rad($latitude)) + 1 / cos(deg2rad($latitude))) / pi()) / 2 * pow(2, $zoom)) - 1;

这将生成2025-2046和2921-2943的网格

到目前为止还不错。接下来，我将关联图像缩放为256px*tiles X 256px*tiles Y，并导出映射到网格tiles/列的256x256 tiles

这将生成一个mbtiles文件，该文件可以工作，但和真实地图稍微不对齐

我尝试使用MapTiler Desktop导出相同的文件，这将生成一个mbtiles文件，其网格为

2025-2046

和

2921-2943

，但它完全对齐。我一辈子都不知道为什么

---

CoordSys地球投影

行中的

-100000

是什么？MapTiler将地理位置显示为

40000013000000-10

。我想知道-10是一个偏移量还是我需要考虑的东西。

你可能会更幸运地发布在上。谢谢，我不知道它的存在。我也在那里发帖了。你可能会更幸运地在上面发帖。谢谢，我不知道那里有。我也在那里发帖了。