Questions tagged «gdal»

我正在关注Mike Bostock的教程“让我们绘制地图”。他介绍了如何ogr2ogr根据特征的属性来过滤shapefile 。我可以在QGIS中打开shapefile，然后查看属性表以找出那些功能。有没有办法从命令行中找到该信息？

9 shapefile  gdal  ogr2ogr 

我正在使用python和QGIS 2.0。我正在尝试通过一个面要素来裁剪文件夹中的栅格。这是我第一次使用（例如）“ PyQGIS”，以前我习惯使用arcpy。无论如何，我没有让我的简单脚本起作用，任何建议将不胜感激！ import qgis.core, qgis,utils QgsApplication.setPrefixPath("C:/OSGeo4W64/apps/qgis", True) QgsApplication.initQgis() CLIP= "C:/Users/unim/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/study_area_foscagno.shp" INPUT_FOLDER="C:/Users/unimi/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/LE71930282000259EDC00" OUTPUT= "C:/Users/unim/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/foscagno_pyqgis/" for RASTER in INPUT_FOLDER.tif do echo "Processing $RASTER" gdalwarp -q -cutline CLIP -crop_to_cutline -of GTiff RASTER OUTPUT+ "clip_"+ RASTER done QgsApplication.exitQgis() 以下是我从现在开始所做的改进，虽然没有使脚本起作用，但是我想我可能会越来越近... import qgis.core, qgis.utils, os, fnmatch from osgeo import gdal CLIP= "C:/Users/unimi/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/study_area_foscagno.shp" INPUT_FOLDER= "C:/Users/unimi/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/LE71930282000259EDC00/DNs2Reflectance_LE71930282000259EDC00" OUTPUT= …

9 gdal  pyqgis  gdalwarp 

我正在尝试剪切较大的ECW（详细信息如下），但是光栅文件太大而无法完全处理。 以下ECW的一些详细信息 驱动程序：ECW / ERDAS压缩小波（SDK 5.0） 文件大小：50gb大小为450000、565081像素大小：0.15 0.15 COLORSPACE = RGB COMPRESSION_RATE_TARGET = 9 VERSION = 2波段数：4 我要裁剪的区域大约是原始文件的1/5。 这是我尝试不成功的方法： 使用Arcgis将ecw保存为tiff /其他格式...（我很快就放弃了） 使用了Qgis及其裁剪器工具...文件创建停留在40％左右。 除了Qgis以外，还使用gdal_translate从OSGeo4W中退出。（试图以为可以不使用Qgis释放一些内存就可以了） 使用gdal_retile的想法是，我会将图像切成小块，然后抓住想要的图像。命令“ gdal_retile -ps 10000 10000 -of ecw -tileIndex tile.shp -targetDir input.ecw这崩溃得更快” 有人有主意吗？ 有关信息，我在具有16GB内存的i5-3470 3.2Ghz上运行Windows 7 64位。

9 raster  gdal  clip  ecw  big-data 

由于data.gc.ca现在正在以FGDB格式分发其某些数据（例如：Plant Hardness Zones），因此我希望能够将其与gdal / QGIS一起使用。但是，编译/安装未按计划进行。 这些是我已采取的步骤： 为ubuntugis-unstable启用了源存储库。 下载了Linux 64位File Geodatabase API 1.3，将其安装在中/usr/local/FileGDB_API，添加了必要的ldconfig设置，并成功构建并运行了测试程序。（如此处建议的过程：[Qgis-user]如何使用ESRI FileGDB编译gdal） 安装了必要的源依赖项： sudo apt-get build-dep gdal 下拉来源： apt-get source gdal 将行添加--with-fgdb=/usr/local/FileGDB_API \到debian/rules文件中的configure命令。 从发出的源目录中dpkg-buildpackage -rfakeroot -uc -b（如此处建议：如何获得通过apt-get安装的软件包的源代码？-Ask Ubuntu） 构建进行得很好，但是在软件包编译阶段失败了，抱怨fakeroot文件系统中缺少Ruby 1.9树： … mv /home/scruss/gdal-1.10.0/debian/tmp/usr/local/lib/site_ruby/1.9/* /home/scruss/gdal-1.10.0/debian/tmp/usr/lib/ruby/1.9/. mv: cannot stat ‘/home/scruss/gdal-1.10.0/debian/tmp/usr/local/lib/site_ruby/1.9/*’: No such file or directory make: *** [install] Error 1 dpkg-buildpackage: …

9 qgis  gdal  ubuntu  linux  file-geodatabase 

几次前，我处理了阿拉斯加大小的每日地面温度卫星图像。 我必须处理大约一千张图像。我下载了它们，进行镶嵌，投影，修剪到shapefile边界，从修剪的栅格中提取均值，然后在时间-温度图上绘制均值。时间-温度图是最终输出。我在ArcGIS中完成了所有这些工作，即使我创建了有助于自动完成一部分工作的ArcPy脚本，也仍然很艰辛。 现在，我想构建一个执行相同操作但自动100％完成的程序。我只想将拐角的坐标输入到程序中，以便定义感兴趣的区域并获得生成的最终输出，即时间-温度图。 您认为使用Python与GDAL库和PostGIS作为栅格数据存储库是否可行？

9 python  postgis  arcgis-desktop  gdal  software-recommendations 

我正在尝试使用精彩的PySAL软件包的映射组件，但是由于某些原因，我没有ogr模块。我觉得它与gdal包装在一起。我想念什么吗？ import gdal import ogr --------------------------------------------------------------------------- ImportError Traceback (most recent call last) <ipython-input-29-ef9eb3418665> in <module>() ----> 1 import gdal 2 import ogr ImportError: No module named gdal 更新： 恐怕我急着赶飞机，所以我没有提供一些关键信息。我已经使用这里描述的apt-get方法安装了python-gdal ，因为我希望没有机会。程序包管理器是否可能没有正确安装？ choct155@choct155-Q550LF:~$ sudo apt-get install python-gdal Reading package lists... Done Building dependency tree Reading state information... Done python-gdal is already the …

9 python  gdal  ogr 

我有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标签，但是我不确定是否真的可以使用它们，或者如何使用。

9 gdal  gdal-merge 

对于那些对GIS编程方面不熟悉的人来说，我觉得许多出色的工具和资源是不可能的。我对基于GUI的向导感兴趣，这些向导用于使用gdal和ogr之类的库处理地理处理。有小费吗？

9 gdal  ogr  gui 

我正在尝试使用python gdal库创建栅格，现在已经到了输出数据的地步，但是输出数据在原点的x轴上翻转了。我知道我一定在忽略某些东西，但是我无法弄清楚哪里出了问题。有任何想法吗？ 创建栅格时，我设置了左上角的x / y值，该数组似乎从左上角开始索引，然后继续向下到右下角。在下面的代码中，我用行的值填充数组。 打印数组时，它看起来像这样： [[ 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. …

9 python  raster  gdal 

我和一位同事最近发布了一个名为Mapoteca的开源项目，该项目是与GDAL / FWTools安装集成的一组Python脚本，允许对目录进行完整的遍历，以查找空间数据并提取范围，几何形状，元素数量，投影和用户将其写入csv文件。我们正在研究描述相关DBF文件的部分，并且我们想生成每个shapefile的缩略图，您知道一种使用Python生成jpeg或png缩略图的方法吗？

9 python  gdal  ogr  jpg 

我正在尝试使用gdal2tiles将20 GB LZW压缩的geotiff转换为tilepyramid。不幸的是，gdal2tiles似乎无法处理此问题，在运行10个小时后，它仍未达到第一个“点”。我用一个较小的文件（1GB）尝试了一下，效果很好。我现在的问题是，是否有其他软件允许我执行此操作（商业或开源）？可以从GlobalMapper创建图块，但是我无法将20GB的文件加载到内存中。

9 gdal  tiles  gdal2tiles 

我有一组正射影像，每个影像都完全充满了数据，但是当放在一起时，coverage不会形成矩形，因此没有数据区域。 在MapGuide中，我需要能够将这些无数据区域显示为白色，以便我的用户在打印时不必浪费黑色墨水。 在大规模情况下，我直接显示原始的全分辨率图像，这没有问题。我只是将地图的背景色设置为白色，而没有正射影像的区域显示背景。 为了提高性能，我需要能够将所有这些源图像合并到重新采样的合成概览图像中，以便以较小的比例显示，从而可以一次查看更多的正射影像。 我一直在尝试使用GDAL对总览进行合并和重新采样，但是默认情况下，它似乎创建了重新采样的复合GeoTIFF图块，其在nodata区域中为黑色，并且MapGuide不允许我将黑色设置为在彩色栅格上透明。 我有办法有效地获得想要的东西吗？ 我已经提供了最终使用GDAL的答案，但是很想看到使用其他图像处理实用程序和GIS应用程序（包括开源和专有）的解决方案。

9 gdal  raster