GDAL -多线程gdalbuildvrt和gdalwarp
GDAL (Geospatial Data Abstraction Library) 是一个开源的地理空间数据抽象库，它提供了处理栅格和矢量数据的功能。GDAL 支持多种数据格式，包括常见的栅格格式（如GeoTIFF、JPEG、PNG）和矢量格式（如Shapefile、GeoJSON）。GDAL 提供了一系列命令行工具，其中包括多线程的 gdalbuildvrt 和 gdalwarp。
- gdalbuildvrt：
- 概念：gdalbuildvrt 是 GDAL 提供的一个命令行工具，用于创建虚拟栅格文件（VRT）。VRT 是一种描述栅格数据集合的 XML 文件，它可以将多个栅格文件组合成一个逻辑上的单一数据集。
- 优势：gdalbuildvrt 可以快速创建大型栅格数据集合的虚拟表示，而无需实际合并这些数据集。这样可以节省存储空间，并提高数据处理的效率。
- 应用场景：gdalbuildvrt 可以用于将多个栅格数据集合并为一个虚拟数据集，方便进行后续的数据处理和分析。例如，可以将多个遥感影像拼接成一个大区域的影像，或者将多个时间序列的影像合并为一个时间序列数据集。
- 推荐的腾讯云相关产品：腾讯云提供了云服务器（CVM）和云存储（COS）等产品，可以用于存储和处理栅格数据。您可以使用腾讯云的云服务器搭建 GDAL 环境，并使用云存储存储和管理栅格数据。
- gdalwarp：
- 概念：gdalwarp 是 GDAL 提供的一个命令行工具，用于栅格数据的重投影、裁剪和重采样等操作。它可以将栅格数据从一个投影系统转换到另一个投影系统，也可以裁剪数据集的范围或者改变数据集的分辨率。
- 优势：gdalwarp 支持多线程处理，可以加快数据处理的速度。它还提供了丰富的参数选项，可以灵活地控制数据处理的方式。
- 应用场景：gdalwarp 可以用于将栅格数据集转换到特定的投影系统，以满足不同应用的需求。例如，可以将地理坐标系的栅格数据转换为投影坐标系，或者将不同投影系统的数据集进行拼接和融合。
- 推荐的腾讯云相关产品：腾讯云提供了云服务器（CVM）和云存储（COS）等产品，可以用于存储和处理栅格数据。您可以使用腾讯云的云服务器搭建 GDAL 环境，并使用云存储存储和管理栅格数据。
更多关于 GDAL 的信息和使用方法，您可以参考腾讯云的文档：
- GDAL 官方网站：https://gdal.org/
- 腾讯云云服务器（CVM）产品介绍：https://cloud.tencent.com/product/cvm
- 腾讯云云存储（COS）产品介绍：https://cloud.tencent.com/product/cos
相关·内容
我有70+光栅图像的TIFF格式，我试图合并。
原件可在这里找到：
在预处理后(pct2rgb、gdalwarp单个图表、gdal_translate以切割项圈)，我尝试通过gdalwarp运行它们，使用如下命令将它们拼接起来：
gdalwarp --config GDAL_CACHEMAX 3000 -overwrite -wm 3000 -r bilinear -srcnodata 0 -dstnodata 0 -wo "NUM_THREADS=3" /data/aeronav/sec/c/Albuquerque_c.tif .....70 other file names
浏览 9提问于2014-03-21得票数 9
1回答
我需要将一个.tiff文件从EPSG:4326 -wgs84-地理重新投影到EPSG: 32638 (UTM Zone 38N)。当我运行下面的代码时，输出文件是一个重新投影的.tiff，大小是原始文件的一半。
# Reprojection
from osgeo import gdal
filename = r"/home/DATA/Test.tif"
input_raster = gdal.Open(filename)
output_raster = r"/home/DATA/Test_UTM38N.tif"
gdal.Warp(output_rast
浏览 5提问于2020-06-17得票数 0
使用Processbuilder执行Windows程序有问题。下面的代码生成一个命令数组，如下所示：
C:\\Program Files\\Gdal\\bin\\gdal\\apps\\gdalwarp.exe, -t_srs, EPSG:32633, -et, 1, -r, near, H:\\GeoData\\Raster\\1M_WEF30JAN2020_WGS84LL_somlos_8.tif, H:\\GeoData\\Raster\\result.tif
String [] cmd = {"C:\\Program Files\\Gdal\\bin\\gdal\\apps\\
浏览 3提问于2021-04-08得票数 0
我使用gdalwarp得到一个带有EPSG 4326投影(等距柱面投影)的行星。在GDAL 2之前，我可以通过下面的命令使用gdalwarp： gdalwarp -te lonMin latMin lonMax latMax -t_srs EPSG:4326 sphere.tiff planisphere.tiff 示例 gdalwarp -te -100 -10 -60 10 -t_srs EPSG:4326 sphere.tiff planisphere.tiff 但是现在在GDAL 3.0.4中，它返回一个错误，就像我在坐标中犯了一个错误一样。 Creating output file
浏览 480提问于2020-07-21得票数 1
回答已采纳
我有一个数据集的MODIS大气产品(MOD04)的.hdf格式，我想重新设计和保存作为Geotiff。我知道使用gdalwarp是可以做到的，但问题是，在保存映像之前，我需要做一些操作。基本上，我需要的是从hdf文件中读取不同的子数据集，使用numpy进行一些修改，最后将数组保存为Geotiff文件。
我不是使用gdal的专家，但我通常知道它是如何工作的。我必须指出，当我通过gdal打开MOD04 hdf文件时，gdal似乎无法读取其GeoProjection和GeoTransform数据！
>>> src = 'D:/Test/MYD04_L2.A2007231.0
浏览 2提问于2015-01-15得票数 2
回答已采纳
1回答
基于GDAL的图像变换
、、、
我有一个在UTM的图像，并希望能够应用在东移和北移使用gdal工具。有没有人知道gdal的哪个函数，以及如何仅用移位来重新投影图像？
谢谢。
浏览 0提问于2015-12-07得票数 0
嗨:在Windows7上，我从Python调用GDAL_translate.exe和GDALwarp.exe来重新采样和重新投影一个单波段geotiff。我试着自己运行GDALwarp，但它抱怨它需要输入文件中的地理配准信息。因此，我首先运行GDAL_translate，并通过-a_ullr标志指定文件角坐标。
有人知道我做错了什么吗？
# learn file upper left coordinate and lower right coordinate
ulXgeo = geoTrans[0] + 0 * geoTrans[1] + 0 * geoTrans[2]
ulYgeo = g
浏览 2提问于2014-09-28得票数 1
我使用apt-get install gdal-bin命令安装了gdal。下面是我得到的结果：
# ls /usr/bin | grep gdal
gdal-config
gdal_contour
gdal_grid
gdal_rasterize
gdal_translate
gdaladdo
gdalbuildvrt
gdaldem
gdalenhance
gdalinfo
gdallocationinfo
gdalmanage
gdalserver
gdalsrsinfo
gdaltindex
gdaltransform
gdalwarp
没有gdal_proximity，但在文档中
浏览 15提问于2017-01-13得票数 2
回答已采纳
目前，我在R中有一个处理链，它下载MODIS数据，然后从系统调用gdalwarp，将特定的子数据集(例如NDVI)重新投影到WGS1984中。然后将生成的GeoTiffs收集到一个HDF5文件中，以便进一步处理。
现在，我将处理链移到python，我想知道是否有一种方法可以跳过将GeoTiffs编写到磁盘的步骤，使用gdal模块的功能。
首先，问题是：
我是否可以执行gdalwarp，严格使用来自gdal模块的python绑定，而不写入磁盘？
我一直在研究，对我的问题最贴切的答案是这样的：
第一个方法需要一个模板，所以不是我想要的。
第二种方法看起来更有希望，它使用的是函数
浏览 2提问于2018-02-09得票数 11
回答已采纳
1回答
我在一个特定的投影上有困难。它似乎在proj4 ()中得到了支持，但当我在gdal中使用它时，它似乎不存在：
gdalsrsinfo：
gdalsrsinfo -o proj4 "+proj=bipc +ns"
收益率
failed to load SRS definition
gdalwarp：
gdalwarp -overwrite -s_srs EPSG:4326 -t_srs "+proj=bipc +ns" -of GTiff in.tiff out.tiff
收益率
ERROR 1: Translating source or target SRS
浏览 3提问于2017-10-30得票数 1
回答已采纳
我有一些MODIS hdf文件，尝试使用Python.Here中的GDAL模块将投影从正弦更改为朗伯特共形圆锥是我的python代码，它返回错误信息：
import os,sys
import gdal
import numpy
os.sys('gdalwarp test.hdf outfile.hdf -t_srs "+proj=lcc"')
File "projection_cmd.py", line 5, in ?
os.sys('gdalwarp test.hdf outfile.hdf -t_srs "
浏览 0提问于2013-01-25得票数 0
1回答
我一直试图在Windows8环境中使用Python2.7和Anaconda来完成gdal库的工作。
除了gdal，我还安装了libgdal (坦白地说，我不太明白两者之间的区别)。我现在似乎有gdal 2.1.0和2.0.2以及libgdal 2.1.0。
但是，当我运行我的Py代码时，会出现一个gdal错误：
“‘gdalwarp”不被识别为内部或外部命令、可操作的程序或批处理文件。
我已经将GDAL_DATA环境变量设置为
C:\Anaconda\pkgs\libgdal-2.1.0-vc9_0\Library\share\gdal
我也增加了一条道路，虽然我不完全确定这应
浏览 4提问于2017-08-31得票数 4
当我尝试从命令行使用".py“函数时，我的GDAL表现出了奇怪的行为：
例如，当我运行1.时，我得到了标准的响应：
GDAL 1.11.3, released 2015/09/16
2. --而且当我运行gdalwarp时，我得到了标准的响应，显示了标志和一切。
3.问题:但是，当我使用.py扩展运行gdal函数(例如gdal_polygonize.py )时，会得到以下内容：
Traceback <most recent call last>:
File "C:\OSGeoW64\bin\gdal_polygonize.py", line 36,
浏览 9提问于2016-02-19得票数 17
回答已采纳
我有一个png图像，它显示了由地理边界UpperLat，LeftLon，BottomLat，RightLon (在度和-ve为西经度和南纬度)定义的区域在板块Carree EPSG:4326。我想使用GDAL将其转换为Mercator EPSG:3857。
当我尝试这样做的时候：
Step1) gdal_translate -of Gtiff -a_srs EPSG:4326 pc_map_piece.png pc_map_piece.tiff
Step2) gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 pc_map_piece.tiff mer_map
浏览 0提问于2012-09-26得票数 1
回答已采纳
1回答
无法复制gdal输出
、、、
我有一组被翻转的GRIB文件(经度跨度从0到365)，我使用gdal首先将数据转换为GeoTIFF，然后将网格数据扭曲为标准的WGS84经度(-180到180)。到目前为止，我一直在使用命令行中的gdal_translate和gdalwarp的组合，并使用parallel快速地遍历所有文件。这些是bash脚本中的函数：
gdal_multiband_transform(){
FILEPATH=$1
SAVEPATH=$2
NUM_BANDS=$(gdalinfo $FILEPATH | grep 'Band' | wc -l)
if [[ $
浏览 1提问于2020-07-07得票数 3
我正在尝试使用Java GDAL API来转换数据集。但是，GDAL Java API没有文档记录，也不清楚WarpOptions采用哪些选项。我给-tr 15 15喂食就像喂gdalwarp一样。这将导致Option not found error
String gdalSource = "LC81910242015182LGN00_B4.TIF";
Dataset originalData = gdal.Open(gdalSource, gdalconstConstants.GA_ReadOnly);
Vector<String> options = new V
浏览 71提问于2019-10-22得票数 0
回答已采纳
我尝试用Nvidia GeForce 1060在ubuntu16上编译最新的Gdal (安装了Cuda toolkit10)。我还通过apt-get为opencl头安装了nvidia-opencl-dev包 ./configure --with-python=python3 --with-spatialite --with-pg --with-cryptopp --with-curl --with-opencl=yes --with-proj=/usr/local --prefix=/usr --with-opencl-include="/usr/local/cuda-10.0/t
浏览 138提问于2020-06-18得票数 1
我有一个包含很多类的C++ COM应用程序。一个Utils类有一个名为GdalWarp()的函数，我还有一个类有类似的函数，但有一个不同的方法叫做GdalUtils，它有一个名为GDALWarp()的函数。 在我的idl中，这两个函数我都有。第一个是大写的GDAL，另一个不是。 我在C#应用程序中使用此COM应用程序，当实例化GdalUtils()时，我只有GdalWarp(大写)，而没有GdalWarp()。但它仍然在de COM应用程序中运行正确的函数。当我注释掉idl中的GDALWarp函数时，我得到的是GdalWarp(小写)。 为什么会有这样的差异？ idl在GitHub上：htt
浏览 13提问于2019-06-14得票数 0
我想将一个tiff文件更改为USGS DEM文件。为此，我正在采取以下步骤：
使用Gdal Wrap实用程序将tiff文件转换为WGS84协同文件
gdalwarp -t_srs EPSG:4326 -dstnodata -90909 -of GTiff Source.tif DestWGS84.tif 84
然后使用Gdal Translate WgS84 utilitygdalTranslate -of USGSDEM DEM 84.tif Final.dem将utilitygdalTranslate tiff文件转换为DEM文件。
当我在gdalwarp中通过-dstno
浏览 0提问于2016-05-27得票数 0
回答已采纳
1回答
我正在尝试在macOS Sierra (10.12.6)上安装旧版本的gdal (1.10)。之所以需要这样做，是因为我已经对土地覆盖变化进行了评估，其中对数据集的初始处理是在2014年完成的。我已经更新了gdal，并且升级到了2.1。然而，返回到重新运行我的完整工作流程(以检查可重复性)，我发现gdalwarp在2.1 (kyngchaos版本)和1.11.5 (最新的家庭酿造版本)之间产生的结果略有不同，并且在我现在使用gdalwarp 1.11.5重新处理的数据集和我最旧的分析之间仍然有一些非常小的差异，我认为这些差异必须是最初使用gdalwarp 1.10完成的。
现在，冗长的前言已经
浏览 4提问于2017-08-28得票数 2
我想包装(重新投影)此netcdf文件中的一个变量。
D:\ gdalwarp -t_srs EPSG:4326 NETCDF:"C:\fie.nc":var "C:\Desktop\SM.img"
但是我得到了这个错误：
错误4:无法打开EPSG支持文件gcs.csv。尝试将GDAL_DATA环境变量设置为指向包含EPSG csv文件的目录。
错误1:转换源或目标SRS失败: EPSG:4326
我使用的是GDAL 1.9.0，发布于2011/12/29。如果有任何帮助，我将不胜感激。
浏览 138提问于2013-01-22得票数 28
回答已采纳
1回答
我正在使用gdal python重新设计一个grib文件。我之所以做这个python而不是gdalwarp，是因为我不想受到编写新文件和阅读python脚本来进行进一步处理的影响。运行gdalwarp非常慢，并且会创建非常大的文件。我目前正在使用来自的重新投影代码。
osng = osr.SpatialReference ()
osng.ImportFromEPSG ( epsg_to )
wgs84 = osr.SpatialReference ()
wgs84.ImportFromEPSG ( epsg_from )
tx = osr.CoordinateTransformation (
浏览 3提问于2015-07-14得票数 0
回答已采纳
1回答
栅格时光栅不对齐
、、
我有一个shapefile (下面是绿点层)，其中我使用gdal_grid将其网格化为栅格，然后gdalwarp将源参考系统(EPSG:27700)转换为目标参考系统(EPSG:3857)。
但是，正如您可以看到的，一旦网格化，它就会移动位置并变得不对齐。
我怎么才能解决这个问题？
浏览 11提问于2022-07-08得票数 0
1回答
我使用gdal创建不同类型的层，例如色浮物和山坡，并使用Mapnik将它们组合成单个图像，作为从单个.hgt文件获得的3D dem模型的纹理。
我的前提是，我对gdal是新手，我正面临一些问题的山坡层。
通过使用gdal命令：
gdal_translate N44E007.hgt N44E007.tif
我得到了N44E007.tif文件，在IrfanView中如下所示
使用以下gdal命令：
gdaldem hillshade -of PNG .\tif\N44E007.tif .\hillshade_png\N44E007_hillshade.png
我得到的N44E007_hi
浏览 5提问于2015-02-16得票数 2
回答已采纳
GDALWarpDestName返回NULL，但未设置错误
我正在尝试将几个栅格(在tif中)重采样到粗略分辨率，而gdalwarp运行结果为空
def resam(inDs):
gdal.Open(inDs)
(shotname,extension) = os.path.splitext(os.path.basename(inDs)) # get the filename and extension
outDs = shotname+'_232_ave'+extension #result tif of resample
gdal.Warp(outDs, i
浏览 2提问于2018-12-28得票数 1
好了，经过一些修改之后，我在第二个注释行中调整了站点超链接的脚本。该脚本的目的是使用一个包含多个多边形(每个多边形都有一个“名称”记录)的shapefile来裁剪/遮罩一个GTiff格式的大型栅格(即不适合32位Python2.7.5应用程序的栅格)，并将裁剪的栅格保存到一个"clip“子目录中，其中每个被遮罩的栅格都以每个多边形的”名称“命名。与原始脚本一样，它假设GTiff和shapefile在相同的投影中并正确重叠，并且每秒处理约100个掩码。然而，我已经修改了该脚本: 1)使用浮点值高程网格，2)仅将较大网格的窗口加载到由当前多边形限制的内存中(即减少内存负载)，2)导出具有正
浏览 6提问于2014-05-21得票数 6
回答已采纳
我有两个不同分辨率的重叠光栅文件(A=25m和B=5m)。我想以最快最有效的方式(如果可能的话使用python绑定实现自动化)用A中的值填充B (-999)中的nodata值。
我测试了不同的方法：
用gdalwarp将光栅A重采样到B的分辨率，然后使用gdal_calc表示公式gdal_calc.py -A MNT_5m.tif -B BATI_5m.tif --outfile=out.tif --calc="maximum(A,B)"，但这没有任何作用，而且重采样是一项耗时的操作。
再次将gdal_calc与calc公式--calc="A*(B=-999)+
浏览 1提问于2018-01-31得票数 3
1回答
我是python的新手，我是GDAL的新手，所以这是一个双重问题。
为了执行一些重复的操作，我想开始使用QGIS python控制台。我试着用这个简单的脚本
import gdal
gdalwarp -ot Float32 -s_srs EPSG:32630 -t_srs EPSG:3035 -r near -of GTiff /home/giacomo/Scrivania/provas/cccp.dat /home/giacomo/Scrivania/provas/cccp_qgs2.tif
但我得到的是语法错误：
文件""，第1行
SyntaxError:无效语法
浏览 1提问于2016-05-12得票数 0
回答已采纳
1回答
是否有一种使用GDAL (使用Python )的移位向量来扭曲图像的方法？
所谓移位向量，我的意思是。CSV (或numpy)文件具有以下列：starting_x,starting_y,target_x,target_y。这幅画很像。Tiff (或GeoTiff)文件。
我没有发现像文档中支持多个移位向量(，)这样的东西，但是也许有人知道如何做到这一点。
编辑
我用一些样例坐标生成了VRT文件。这是我使用的命令：
gdal_translate -gcp 100 100 200 200 -gcp 200 200 300 300 -gcp 500 500 600 600 -of VRT origi
浏览 0提问于2019-04-03得票数 1
1回答
到目前为止，我还没有在GDAL库上工作过。如果这个问题没有任何意义，请原谅。
我正在尝试实现图像的变形从EPSG:32611到EPSG: 3857。我使用命令行使用以下命令成功地实现了warp。我想在我的C#项目中实现相同的命令集。我不知道从哪里开始。我找不到任何用于C#的GDAL包装库。在这方面的任何帮助都将不胜感激。
这些是我在命令行中运行的命令。Raw.png是源镜像。
gdal_translate -of Gtiff -co "tfw=yes" -a_ullr 273996.1204 5070668.029 274682.6204 5069981.529 -a_srs
浏览 1提问于2014-11-26得票数 0
我有一个10°×10°的非地学HDF5，我想把它转换成地点。我尝试过2种方法(1) gdalwarp和(2) gdal_translate，但是无法找到如何转换文件。
1) gdalwarp
gdalwarp -t_srs '+proj=longlat +datum=WGS84 +no_defs ' -te 40 30 50 40 -tr 0.0089285714 -0.0089285714 -multi HDF5:"/PathToTheHDF/FileName.HDF5"://SUBDATASET /PathToOutput/FileName_Subdatas
浏览 3提问于2014-04-08得票数 4
回答已采纳
2回答
我正在尝试使用GDAL将netCDF文件转换为EPSG:3857以便与Mapbox一起使用。这将是从.nc到.nc的转换。而不是光栅。我愿意使用GDAL或其他方法来做到这一点。这些数据必须在发送到控制台应用程序之前重新投影-这个过程需要数周时间才能找到解决方案-我认为这很简单。 我正在给卫星数据上色。有3个.nc文件(蓝色、红色和红外线)，它们在组合和处理时会生成彩色图像。下载完这3个文件后(从Amazon AWS下载)，python控制台应用程序进行处理并将.jpg转储到相同的文件夹中。该应用程序的源代码是Located here so you may validate the data。
浏览 99提问于2019-03-02得票数 9
回答已采纳
我正在寻找改变GDAL栅格数据集分辨率的最佳方法。
例如，我有一个像素大小为(30，-30)的光栅，我想将像素大小更改为(5，-5)，将给定像素的所有值插值到输出光栅中。
因此，对于输入光栅的每个像素，我希望输出光栅中有36个像素，所有这些像素都共享相同的值。
如果我运行gdalwarp -tr 5 -5 inputRaster.tif outputRaster.tif，我会得到我想要的结果，所以我假设我应该能够用一些GDAL函数来复制这个功能。
如果可能的话，我倾向于避免使用对python的Subprocess类的调用。
浏览 0提问于2011-10-11得票数 8
回答已采纳
我已经使用了C++ gdal库来读取光栅图像。具体来说，我在这一点上读过hgt和geotiffs。这些文件类型的默认垂直基准面通常是EGM96大地水准面。我需要能够使用C/C++ gdal库来转换为高于WGS84椭球的高度，以便与来自我的GPS的数据进行比较。我知道我可以使用gdalwarp ()在命令行中完成这项工作，但我更希望能够使用任何用户输入，而不需要运行预处理脚本。如果有必要，我会执行系统调用，但是不是应该有一种方法可以通过c++库在我的应用程序中完成这项工作吗？
浏览 109提问于2020-06-04得票数 0
通常，我使用以下终端命令来转换我的Sentinel2 Level 1C数据(它是.jp2格式)，以将其转换为EPSG3411格式的.tif
gdalwarp -t_srs EPSG:3411 input.jp2 output.tif
由于在卫星图像上我想要转换的数量，我想有一个python代码。我阅读了文章并尝试应用它(没有成功工作)。过了很长一段时间，我发现对如何使用gdal.Warp的“说明”真的很小。
现在我有了下面的代码，它不能工作：
import gdal
input_file = 'path/test.jp2'
output_file = 'test.ti
浏览 0提问于2019-04-02得票数 0
1回答
有人能告诉我如何使用gdal羽化图像的边缘吗？我一直认为-cblend开关可以做到这一点，但也许我误解了什么。
这是我尝试过的：
gdalwarp -t_srs ESPG:4326 -multi -cblend 200 -of GTiff input.jp2 test.tiff
我不介意在QGIS中或仅在OSGeo4W命令提示符下执行此操作。
浏览 5提问于2014-08-27得票数 1
1回答
我想合并一些1波段的光栅图像。我想要处理重叠，以便选择像素中的最小值。所有图像都有相同的投影。
我尝试查看gdalwarp和gdal_merge (在命令行中)，但在重叠部分，它们只使用上一张图像中的值。
我也看到过使用PIL blend或paste的建议，但这些需要一个alpha层，它指示重叠部分应该如何混合，当然我没有。
有没有人知道如何使重叠依赖于图像中的实际值？
浏览 8提问于2016-11-24得票数 2
我在试着把罗宾逊投影转换成墨卡托。例如，我使用这个image。 首先，我应用地理标记： gdal_translate -a_ullr -180 90 180 -90 -a_srs ESRI:54030 source.jpg source_tagged.tif 并最终将其传递给墨卡托： gdalwarp -t_srs ESRI:54030 -s_srs EPSG:3857 source_tagged.tif target.tif 结果在垂直方向上略有延伸，但与适当的墨卡托相去甚远。我做错了什么？
浏览 86提问于2021-10-12得票数 1
回答已采纳
1回答
我可以在我的终端中成功地执行下面的命令。
命令：
gdalwarp -s_srs "+datum=WGS84 +no_defs +geoidgrids=egm96-15.gtx" -t_srs "+datum=WGS84 +no_def" input.tif output.tif
现在，我想将这个命令存储到一个变量中，并在一个docker容器中展开这个命令。
我的脚本run.sh如下所示。我首先将目标命令存储到mycommand中，然后以命令作为输入运行容器。
mycommand=$@;
docker run -ti --rm osgeo/gdal:ubuntu
浏览 4提问于2022-01-06得票数 0
回答已采纳
我有平铺的png文件，这些投影是EPSG:4326。我使用以下两个命令将投影转换为EPSG:3857：
gdal_translate -of Gtiff -a_ullr 135.00000000000003 36.59788913307022 140.62500000000003 31.952162238024975 -a_srs EPSG:4326 4326.png 4326.tiff
gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 4326.tiff 3857.png
我能用一个命令来完成吗？
浏览 1提问于2020-08-09得票数 1
澄清:我忽略了关键方面:不使用os.system或子流程--只是python。
我正在尝试将NOAA GTX偏移网格的一部分转换为垂直数据转换，而不是完全遵循如何使用python在GDAL中这样做。我想使用一个网格(在本例中是一个测深属性网格，但它可能是一个geotif)，并使用它作为我想要做的模板。如果我能做到这一点，我有一种感觉，它将极大地帮助人们使用这类数据。
这是我所拥有的绝对不起作用的东西。当我在生成的目标数据集(dst_ds)上运行gdalinfo时，它与源网格包不匹配。
from osgeo import gdal, osr
bag = gdal.Open(bag_filena
浏览 3提问于2012-05-04得票数 17
回答已采纳
我不确定这是不是正确的网站，
我正在使用Python2.7上的gdal.Open()模块打开陆地卫星GeoTIFF图像。它们具有相同的UTM地图投影，但图像大小不同。
如何将图像大小固定为相同？我想直接进行逐个像素的比较。
浏览 3提问于2011-08-24得票数 1
2回答
如何在NetCDF中更改经度范围
、、、、
我的NetCDF文件有4个维度(经度、纬度、时间、水平压力)和多个变量。
我想替换经度数据340 342 344 360 0 2 4. 18 20乘-20 18 -16 . -2 0 2 4.18 20。事实上，经度是从0°到360°，但我需要它从- 180°到180°。
我在网上找到了一些想法，但我不知道如何应用它们。
这里有一个解释：
“一些netcdf文件包含严格的正值0,360，与gdal/QGis对待经度值的方式不兼容--它们不环绕日期线。这也会导致gdalwarp错误，当与其他数据集( - 180 ,180值)组合时，由于180度的经度偏移，它们不能正确对齐。”
简单的修正是测试经
浏览 4提问于2013-11-21得票数 3
回答已采纳
1回答
如何在CMD中打开文件
、、、
我是Java的初学者，也是一名地理测量学的学生。
我想在DOS命令中从QGIS运行gdal单元。我设法在命令中打开了可执行文件。但是该命令显示了几种类型的子文件夹。
C:\Program Files\QGISWien\bin>"C:\Program Files\QGISWien\bin\gdalwarp.exe"
Usage: gdalwarp [--help-general] [--formats]
[-s_srs srs_def] [-t_srs srs_def] [-to "NAME=VALUE"]
[-order n | -tps
浏览 0提问于2015-06-29得票数 0
我有一个应用程序，需要导入美国国家气象局的表面分析，这些分析以文件的形式分发。我想将这些数据放入PostGIS 2.0的栅格中，进行一些计算和建模，并在GeoServer中显示数据和模型结果。
由于grib2是一种支持GDAL的格式，因此提供的raster2pgsql实用程序应该能够将grib2直接插入到与PostGIS兼容的SQL中，并且一旦它存在，GeoServer就应该能够处理它。然而，我遇到了一些问题，这些问题没有明显的解决方案--至少对我来说不明显！Raster2pgsql运行，显然没有错误，生成SQL，然后运行SQL创建看起来非常像栅格的东西。但是GeoServer不能显示它--特
浏览 1提问于2012-05-18得票数 4
回答已采纳
1回答
我在R中为GDAL编写了一个非常简单的包装器，它利用了一个预写的语句，该语句被传递给系统，创建一个输出，然后我想再次读取到R环境中。 它的工作原理是在工作目录中创建一个临时目录，打印出我们感兴趣区域的ESRI形状文件，然后用这个文件剪切一个带有一些预设信息的栅格。 我的问题:成功运行system()调用并创建输出文件后，函数停止。它不会执行下一个调用，也不会将输出读入R环境。 gdalwarpwarp <- function(source_file, source_srs, newfilename, reread=TRUE, clean=TRUE, cutline){
#Crea
浏览 14提问于2021-05-05得票数 0
回答已采纳
1回答
我对堆栈溢出和bash/python脚本非常陌生。
我需要调整一些.tif格式的数据地形模型文件(.tif)的大小，以便能够将它们转换成.hgt，并且我能够使用gdal工具完成所有这些工作，但只能同时使用单个文件。我猜您已经注意到脚本在哪里出现了:需要将300+文件的过程自动化！
因此，我开始研究bash是如何工作的，并得出如下结论：
#!/bin/bash
for filename in "'/home/fenix/1\ Vari\ HDD/MTB/DTM\ Alos/'"*.tif; do
PATH=/usr/bin/ gdalwarp -of Gt
浏览 0提问于2019-02-10得票数 1
回答已采纳
1回答
我正在使用在线查看天气数据，我使用的网站(俄克拉荷马大学)提供了一个链接，指向显示为geoTIF的数据，用于研究目的。
我正在尝试使用Mapbox地图中的图像，但投影似乎有问题。我使用GDAL工具(尽管我是一个新手)，我甚至不知道它是从什么投影开始的。当我使用gdalinfo时，我得到以下结果：
Warning 1: RowsPerStrip not defined ... assuming all one strip.
Raster dataset parameters:
Projection:
RasterCount: 1
RasterSize (7000,3500)
Usin
浏览 0提问于2018-11-16得票数 1
1回答
用GDAL C++裁剪多边形光栅
、、、、
我试图用多边形和GDAL来剪辑光栅。目前，我收到一个错误，即初始化WarpOperation时存在读取访问冲突。我可以访问我的Shapefile和检查的功能的数字，所以访问是好的，我认为。此外，我还可以访问我的栅格数据(GetProjectionRef)。所有文件都在同一个CRS中。有什么方法可以使用GdalWarp和剪刀吗？
const char* inputPath = "input.tif";
const char* outputPath = "output.tif";
//clipper Polygon
浏览 7提问于2022-01-11得票数 1
回答已采纳
1回答
我有一个全局.nc文件的大型数据集，我正在尝试将它们裁剪到较小的区域。我将这个区域存储为.shp文件。 我试过使用Qgis中的gdal，但需要转换每个变量，我必须为所有文件一个接一个地选择每个变量和相同的形状，400个文件通过每个变量似乎不是最好的主意。同样，这会返回分开的.tiff文件，而不是我想要的.nc文件。 我有这个小脚本，但它不能满足我的需求 import glob
import subprocess
import os
ImageList = sorted(glob.glob('*.nc'))
print('
浏览 90提问于2020-10-29得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
