了解地图图块的创建属性


10

我正在尝试创建一个瓷砖底图。我有一个要从中创建图块的图像,并且正在使用GDAL。

我设法创建了磁贴,但是当我查看结果时,磁贴地图似乎是:

  1. 它被翻转
  2. 比例不对

这是地图叠加图像: 地图图片

这是切片的结果:

平铺结果

通过比较正射影像上的图像角,我多次检查坐标是否良好。

另一件事是,在具有相同坐标的Google Earth中,地图图像正好适合于Google地图。

谷歌地球覆盖

已编辑

gdalinfo map_overlay.png

输出如下:

驱动程序:PNG /便携式网络图形文件:map_overlay.png大小为3527,2494坐标系为''图像结构元数据:
INTERLEAVE = PIXEL角坐标:左上(0.0,0.0)左下(0.0,2494.0)右上(3527.0 ,0.0)右下(3527.0,2494.0)中心(1763.5,1247.0)Band 1 Block = 3527x1 Type = Byte,ColorInterp = Red Mask Flags:PER_DATASET ALPHA Band 2 Block = 3527x1 Type = Byte,ColorInterp = Green Mask Flags:PER_DATASET ALPHA Band 3 Block = 3527x1 Type = Byte,ColorInterp = Blue
Mask标志:PER_DATASET ALPHA Band 4 Block = 3527x1 Type = Byte,ColorInterp = Alpha

然后翻译:

gdal_translate -of VRT -a_srs EPSG:4326 -gcp 0 0 31.7431761644 35.1680410195 -gcp 3527 0 31.7493769674 35.1784535489 -gcp 3527 2492 31.7431011291 35.1784951643 -gcp 0 2494 31.7431761644 35.1680410195 map_overlay.png map_overlay.png

输出如下:

输入文件的大小是3527、2494

最后:

gdal2tiles.py -z 14-21 map_overlay.vrt

叠加层是简单的.png图像,没有任何特殊数据。

我将不胜感激。

香妮


2
gdal_translate您在此处发布的参数似乎不一致:您正在将(lat,lon)=(31.7431761644,35.1680410195)与两个分别位于(0,0)和(0,2494)的点相关联。这真的是您发出的命令,还是复制粘贴错误?
ub

Answers:


5

尝试将度量标准协调图像平铺为:

   gdalwarp -of GTiff -t_srs EPSG:3857 input.tif output.tif 

更新

注意:即使指定gcp,gdal_translate也不会指定tiff的角坐标。

除此之外,您不需要创建虚拟数据集(.vrt),因为只有一个文件...

使用以下代码:

  gdal_translate -of  GTiff  -gcp 0 0 31.7431761644 35.1680410195 -gcp 3527 0 
       31.7493769674 35.1784535489 -gcp 3527 2492 31.7431011291 35.1784951643
      -gcp 0 2494 31.7431761644 35.1680410195 map_overlay.png map_overlay.tif

然后

  gdalwarp -s_srs epsg:4326 -t_srs epsg:4326  map_overlay.tif warped_map_overlay.tif

现在尝试使用gdalinfo ...如果您看不到任何问题,请运行以下命令:

  gdal2tiles.py -z 14-21 warped_map_overlay.tif

如果您的瓷砖不适合,请在更新之前运行以上代码,我已经给出了。

您可以在这里获得更多信息

我希望它能对您有所帮助。


谢谢。我什么时候应该这样做?就目前而言,我正在执行1.gdalinfo map.png 2.gdal_translate -of VRT -a_srs EPSG:4326 -gcp .... 3. gdal2tiles.py -z 10-16 map.vrt。翻译后需要做吗?我应该在翻译方法中更改EPSG吗?
shannoga 2012年

如果您的地图未协调,尝试使用gdalinfo获取信息是徒劳的。
阿拉贡

所以我需要改用gdalwarp吗?
shannoga 2012年

您的地图是否与任何投影相协调?如果是,请尝试使用gdalwarp重新投影,否则使其与gdal_transform协调...
Aragon

嘿,谢谢你的耐心。如果理解正确,那是不对的。这是根据地图绘制的简单正方形图像,上面没有特殊数据。我已经用我的命令更新了问题。再次感谢您的尝试。
shannoga 2012年
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.