Questions tagged «gdal-merge»

4
如何使用VSIS3从S3存储桶中通过GDAL有效地访问文件?
因此,GDAL最近添加了一项新功能,该功能允许随机读取S3存储桶文件。我希望从图像的多个图块中裁剪GDAL图像,而不必下载整个文件。我只看过关于如何通过GDAL配置和访问S3存储桶的非常稀疏的文档,并且对如何开始感到困惑。是否有人会提供一个非常简短的示例/教程来说明如何为GDAL设置虚拟文件系统以实现这一目标?如果您的解决方案允许通过Python编写脚本,请加分! 需要说明的是:我们已经用Python完成了它。Python的问题在于,您必须下载整个图像才能对其进行操作。GDAL的最新版本支持安装S3存储桶,因此,如果我们需要对图像的一小部分进行裁剪,则可以直接在该小部分上进行操作。las,该功能仅在一月份在稳定分支中发布,因此我没有找到任何文档。因此,该解决方案应使用GDAL的最新版本中的VSI3系统,否则应巧妙地使用该系统,以防止用户需要将整个映像下载到EBS驱动器上以对其进行操作。 也就是说,悬赏将授予使用最新版本GDAL中的VSI API的答案,从而无需将整个文件读取到内存或磁盘中。另外,我们使用的存储桶并不总是公开的,因此发布的许多HTTP技巧在我们的许多情况下都行不通。

3
gdalwarp和gdal_merge进行镶嵌有什么区别
我想合并约20个dem转换,每个约100MB,而我使用了gdalwarp(我知道虚拟文件选项): gdalwarp $(list_of_tiffs) merged.tiff 但是,此命令花费的时间很长,大约15分钟后,它正在处理集合中的第5张图像。我终止它。 然后,我发现gdal_merge.py脚本也可以用于镶嵌,并尝试了它: gdal_merge.py $(list_of_tiffs) 不到3分钟即可完成。 正如我期望的那样,两个命令将产生相同的结果,我想知道两者之间的区别是什么,如果结果相同,为什么gdalwarp需要这么多时间?

1
调用gdal_merge进入python脚本
我是python的新手,并且对如何在(spyder)python脚本中使用gdal_merge不屑一顾。我使用Windows 10,Python 3.6,并从osgeo4w安装了gdal工具。我意识到还有很多其他文章描述了这个问题,但是没有一个可以帮助我解决这个问题。 当我从cmd调用gdal模块时,它的工作原理很像: python "C:\OSGeo4W64\bin\gdal_merge.py" -o merged.tif input1.tif input2.tif 但是,我无法使其在python脚本(spyder)中正常工作。 第一种方法产生一个输出,但是名称不正确:它产生一个“ out.tif”文件,而不是我要求的“ merged.tif”文件: import sys sys.path.append('C:\\OSGeo4W64\\bin') import gdal_merge as gm gm.main(['-o', 'merged.tif', 'N00W078.tif', 'N00W079.tif']) 第二种方法只是不产生任何输出: import subprocess merge_command = ["python", "C:\OSGeo4W64\bin\gdal_merge.py", "-o", "merged.tif", "input1.tif", "input2.tif"] subprocess.call(merge_command,shell=True) 关于如何解决此问题有什么想法?

3
如何使用GDAL为光栅带添加颜色解释?
我有3个单独的1频段GeoTIFF文件。带的颜色解释为灰色。我想要一个3波段RGB文件。我用过gdal_merge.py gdal_merge.py -separate file1.tif file2.tif file3.tif -o output_file.tif 但是,结果output_file.tif的3个波段的颜色解释为灰色,未定义,未定义。除此之外,所有像素均为0。 如果执行pct2rgb.py,我将得到一个RGB文件: pct2rgb.py output_file.tif output_file_rgb.tif 但当然,像素也为0。因此,我有3个问题: gdal_merge.py是将3个1波段文件合并为一个3波段RGB文件的正确工具吗? 为什么我得到乐队未定义的颜色解释? pct2rgb是将具有这种颜色解释的3波段文件转换为RGB文件的正确工具吗? 更新:栅格没有颜色表。只是颜色解释:灰色。 另一方面,像素值从0到1023(这是有意的) 更多数据:它们是旋转的栅格(不向北),但是它们都具有相同的地理变换。 更新2:我可以扭曲图像以使其朝北,构造一个VRT并为每个波段添加ColorInterp,但是在输出结果中我仍然得到颜色解释为灰色,未定义和未定义。 问题是我至少需要在第一个波段中创建一个颜色表。我知道创建它们的方法,但是我不知道我的表应该有多少个条目。为什么在GDAL Raster FAQ示例中有13个条目?如果有帮助,所有像素的值都在0到1023之间。 更新3:显然,TIFF格式无法真正指定每个频段的颜色解释。读取TIFF文件时,GDAL构建颜色解释的方式是PHOTOMETRIC和EXTRASAMPLES标记的值的组合。 阅读这些标签: PHOTOMETRIC表示图像数据的色彩空间。值2表示像素值的组成部分是RGB,但是它假定为Byte像素,并且我有UInt16像素(我尝试过-co "PHOTOMETRIC=rgb",并遇到错误)。因此,我无法为输出文件指定PHOTOMETRIC标签。 EXTRASAMPLES指定每个像素具有N个额外的分量。我不确定如何使用此标记创建合并的文件。或者,如果我需要它。 因此,在更新2中,我建议创建一个ColorTable,但是如何呢?在我的3个输入文件中,像素值从0到1023。是否必须将它们与颜色匹配?我必须创建一个具有1024个输入的ColorTable吗?怎么样? 在更新3中,似乎在创建合并文件时可以使用一些GeoTIFF标签,但是我不确定是否真的可以使用它们,或者如何使用。
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.