Questions tagged «gdal»

我正在尝试使用ogrinfo获取有关我下载的shapefile的一些详细信息。当前，我知道如何执行此操作的唯一方法是将其加载到QGIS中，然后手动单击以查找有关它的任何信息，例如打开属性表。 我只希望能够看到所有带有功能标记的元数据。如果我做： ogrinfo -al USA_adm0.shp 我可以看到一开始有很多有用的信息，但是随后所有的特征数据就飞逝了。 有人可以帮我吗？ 编辑 这就是我使用-ro和-so标志在Mac上获得的内容，似乎并没有太大帮助。 ->ogrinfo -ro -so USA_adm0.shp INFO: Open of `USA_adm0.shp' using driver `ESRI Shapefile' successful. 1: USA_adm0 (Polygon)

12 shapefile  gdal  ogr 

我有一组包含每日土壤湿度数据的栅格（.vrt）。我想对所有栅格的像素求和，以便按月进行度量。但是，每天都没有覆盖整个世界，这导致在前几天存在度量的地方没有数据值。 我想做的是对每个栅格的值求和。但是，似乎每当无数据像素合计时，结果像素就被直接分类为无数据。我想相反：忽略所有nodata值并求和其余值。 我想到了两种解决问题的方法： 汇总栅格而忽略nodata值 将无数据像素转换为值0，然后将所有栅格相加 不幸的是，我找不到任何工具可以做到这一点。 谁能帮我？

12 raster  gdal  qgis-2 

有没有一种方法可以使用GDAL Python绑定执行与gdalbuildvrt实用程序相同的任务？到目前为止，除了创建单个数据集的vrt和手动编辑xml之外，我没有找到任何其他方法。我想从多个栅格创建vrt（本质上是执行镶嵌）。使用纯Python是否可能？我的另一个选择是使用子进程简单地调用gdalbuildvrt。

12 python  raster  gdal  vrt  gdalbuildvrt 

我正在寻找类似于GRASS r.resample的功能。我有一个装满TIF文件的文件夹，并希望降低其分辨率（“按比例缩小”）以减小文件大小。

12 raster  gdal  ogr 

我正在使用一种工具将栅格图层转换为数组以进行NumPy处理，并且理想情况下，我希望能够处理打包在.gdb中的栅格而不导出所有栅格（特别是如果这需要使用ArcGIS或ArcPy的话） ）。 OpenFileGDB驱动程序可以做到这一点吗？ 据我所知，该驱动程序似乎将栅格图层与矢量图层相同，这使您可以访问有关该图层的一些数据，但不具有ReadAsArray功能。

12 python  gdal  file-geodatabase 

我在具有4个内核的debian机器上运行gdal2tiles： sudo /usr/bin/gdal2tiles.py -r cubic -a 0,0,0 -z 10-15 /home/adm/topo/ekb.vrt /var/www/tiles 在htop我看到： 为什么gdal在工作中只使用一个核心？Gdal版本是2.1.2。

11 gdal 

我想加快与此命令相对应的过程： gdalwarp -t_srs EPSG:4326 -overwrite input.ntf output.tif 有没有办法在GDAL中使用并行处理？半场论题：如果不是，您是否建议使用非GDAL解决方案来加快类似galwarp的过程？ 这是我浏览过的网页： 是否可以在GDAL和QGIS中进行并行处理？ [gdal-dev] GDAL栅格处理：并行计算 编辑：我问了这个问题，因为我认为处理gdalwarp时我的CPU％保持在100％以下。但是，再看一看，它达到了555.5％。

11 gdal  parallel-processing 

我在英国国家网格投影中有一个shapefile： Geometry: 3D Polygon Feature Count: 5378 Extent: (9247.520209, 14785.170099) - (638149.173223, 1217788.569952) Layer SRS WKT: PROJCS["British_National_Grid", GEOGCS["GCS_airy", DATUM["OSGB_1936", SPHEROID["Airy_1830",6377563.396,299.3249646]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",49], PARAMETER["central_meridian",-2], PARAMETER["scale_factor",0.9996012717], PARAMETER["false_easting",400000], PARAMETER["false_northing",-100000], UNIT["Meter",1]] cat: Integer (9.0) 我可以使用GDAL / OGR来获取shapefile中所有多边形的总面积（以公顷为单位）吗？ 我想知道是否可以使用-sql，例如： ogrinfo -sql "SELECT SUM(ST_Area(geom::geography)) FROM mytable" myshapefile.shp 但是，尝试让我得到ERROR 1: Undefined function 'ST_Area' used.。 我想我可以将Shapefile导入QGIS，为每个多边形添加一个area属性，然后对其求和，但是如果可能的话，我宁愿使用命令行工具。

11 gdal  ogr  area  ogrinfo 

我有大量CSV格式的数据，其中包含UTM坐标X和Y和高程值作为Z信息。我需要将这些数据转换为GeoTIFF的DEM以便进一步分析。在这种情况下，大量意味着16 m。线，具有在一个点X，Y和Z每行。这些点是均匀分布的，因此不需要插值；每个点仅需要转换为栅格像元。 原始数据不带分隔符，列宽固定。我已经弄清楚了如何使用流文本编辑器sed将文件语法转换为使用分隔符而不是固定宽度并消除所有空格字符。从这里开始，通常，我的工作流程是通过根据和数据创建要素类，将数据导入ArcGIS X，Y然后Z第二步，使用“ 指向栅格”工具将点shapefile转换为GeoTIFF 。但是，对于该过程，我当前拥有的文件太大。 除了寻找上述工作流程之外，我一直在寻找一种有效的替代方法，并发现了GDAL。但是，在gdal_translate，我可以在支持的文件类型列表中找到的最接近的受支持格式是ASCII网格，但没有逗号分隔的XYZ。另一个困难是，我有UTM坐标，而大多数示例似乎都使用十进制度坐标。但是，我需要停留在UTM系统内（或者至少，我的输出GeoTIFF必须位于UTM坐标系内）。 因此，我正在寻找一种使用GDAL将CSV XYZ转换为GeoTIFF的方法，但到目前为止，尚无法找到处理此确切问题的示例。对于某些提示甚至代码示例，我将感到非常高兴。

11 gdal  geotiff-tiff  csv  gdal-translate  xyz 

为了在子文件夹上递归迭代并将所有shapefile合并为一个脚本，基本脚本是： #!/bin/bash consolidated_file="./consolidated.shp" for i in $(find . -name '*.shp'); do if [ ! -f "$consolidated_file" ]; then # first file - create the consolidated output file ogr2ogr -f "ESRI Shapefile" $consolidated_file $i else # update the output file with new file content ogr2ogr -f "ESRI Shapefile" -update -append $consolidated_file …

11 gdal  ogr2ogr  merge  bash  mergeshapes 

我正在寻找使用GDAL（OSGeo4W）将16位无符号tiff图像缩放为8位图像。但是，我想限制像素值，以使它们保持在8位图像的范围内。我的意思是，我想确保16位图像中的视觉比例在8位图像中表示（0-255）。使用gdal_translate将16位转换为8位会剪切像素值，并且不会生成RGB图像。我不确定从gdal_translate中调用什么函数和设置，或者是否有更好的选择。 GDALINFO C:\>gdalinfo C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Driver: GTiff/GeoTIFF Files: C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif Size is 17310, 11310 Coordinate System is `' Metadata: TIFFTAG_DATETIME=2015:05:03 12:27:13 TIFFTAG_IMAGEDESCRIPTION=UltraCam-Lvl03 -------------- CAM_ID: UC-SXp-1-50215465 [3] IMG_NO: 4761 RECORD_GUID: 8fc3d7a9-d2e5-40e5-babf-939f1f803dcf IMG_GUID: C850D0A4-39CA-4BC8-ACD7-81A363D78C6E FILE_GUID: 86C4F8E4-A2FD-43FB-B3B1-4C4ED8DF0E42 LICENSE_ID: 820342819 SOFTWARE: UltraCam Aerial Radiometry Core 12.6.1408.2501 PIXEL_SIZE_WIDTH: 6 [micron] PIXEL_SIZE_HEIGHT: 6 [micron] APERTURE: F_8 EXPOSURE_TIME: 0.002000000000 …

11 gdal  imagery  gdal-translate 

我正在尝试学习使用Python GDAL绑定和numpy进行遥感图像处理的绳索。作为首次尝试，我正在读取Landsat8 geotiff文件，进行简单的操作并将结果写入新文件。下面的代码看起来工作正常，除了原始栅格转储在输出文件中，而不是转储在操纵的栅格中。 欢迎提出任何意见或建议，但要特别注意为何未在结果中显示操纵的栅格。 import os import gdal gdal.AllRegister() file = "c:\~\LC81980242015071LGN00.tiff" (fileRoot, fileExt) = os.path.splitext(file) outFileName = fileRoot + "_mod" + fileExt ds = gdal.Open(file) band = ds.GetRasterBand(1) arr = band.ReadAsArray() [cols, rows] = arr.shape arr_min = arr.Min() arr_max = arr.Max() arr_mean = int(arr.mean()) arr_out = numpy.where((arr < arr_mean), …

11 python  gdal  geotiff-tiff  landsat  numpy 

我正在努力进行处理，希望可以在此解决。 我使用应用于林业的遥感技术，尤其是处理LiDAR数据。想法是使用Scikit图像进行树顶检测。由于我是Python的新手，因此我认为完成以下工作非常有个人才能： 导入CHM（使用matplotlib）； 运行高斯过滤器（带有scikit-image包）； 运行一个最大值过滤器（带有scikit-image包）； 运行peak_local_max（带有scikit-image包）； 显示带有局部最大值的CHM（带有matplotlib）； 现在我的问题。当我使用matplot导入时，图像会丢失其地理坐标。因此，我拥有的坐标只是基本图像坐标（即250,312）。我需要获取图像中局部最大值点（图像中的红色点）下的像素值。在论坛的这里，我看到一个人在问同样的事情（在没有NumPy的情况下，在OGR点下获取GDAL栅格的GDAL像素的像素值？），但是他已经在shapefile中包含了这些点。在我的情况下，这些点是使用scikit-image计算的（这是一个具有每个树顶坐标的数组）。所以我没有shapefile。 总之，我最后想要的是一个txt文件，其中每个局部最大值的坐标都在地理坐标中，例如： 525412 62980123 1150 ...

11 python  gdal  remote-sensing 

我想在Raspberry Pi上运行GDAL QGIS，以充当负担得起的GIS处理引擎。该计划是购买大约50 Pis。 是否有人在Raspberry Pi上安装和运行GDAL或仅在RASPBIAN环境中安装GDAL方面有任何链接或信息？ 最好的是，Raspbian是我的首选，它是基于Debian的，但是ARCH，RASPBMC，PIDORA，RISC OS和OPENELEC是开箱即用的Linux风格。 希望使用Python。

11 gdal  linux  raspberry-pi 

我有大量的RapidEye图像（+300 Gb），我想用作在线应用程序的底图。到目前为止，我已经按照名为“类固醇中的Geoserver”的演示文稿中显示的步骤设法在geoserver中将它们用作ImagePyramid。问题是我无法进行适当的色彩平衡，因此在查看整个数据集时，马赛克看起来很差，但是当我放大时，马赛克会变得更好（下面的图片）。 我的过程是： 使用Gdal将所有图像转换为8位，epsg 4326，三次插值 使用所有图像生成Gdal虚拟栅格（gdalbuildvrt） 生成图像金字塔（带有压缩和geotiff耕作的gdal_retile）并发布在geoserver上（使用图层样式中的直方图拉伸） 关于如何改善流程并获得更好的色彩平衡的任何技巧？

11 geoserver  gdal  basemap  image-pyramids