Questions tagged «gdal»

我正在尝试创建一个瓷砖底图。我有一个要从中创建图块的图像，并且正在使用GDAL。 我设法创建了磁贴，但是当我查看结果时，磁贴地图似乎是： 它被翻转 比例不对 这是地图叠加图像： 这是切片的结果： 通过比较正射影像上的图像角，我多次检查坐标是否良好。 另一件事是，在具有相同坐标的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 …

10 gdal  tiles 

我是python地理空间编程的新手。我运行了以下脚本并得到了相应的错误消息 >>> import osgeo >>> import osgeo.ogr >>> shapefile = osgeo.ogr.Open("tl_2009_us_state.shp") >>> numLayers = shapefile.GetLayerCount() Traceback (most recent call last): File "<pyshell#5>", line 1, in <module> numLayers = shapefile.GetLayerCount() AttributeError: 'NoneType' object has no attribute 'GetLayerCount'

10 python  gdal 

我尝试在QGIS打印管理器工具中显示一个小分辨率的光栅文件。我想，以显示与每个单元格的值（位置-一个小标签大约每个单元的中间）。 我该怎么办？有没有可以帮助我完成此任务的插件？ 我的想法是使用gdal2xyz，使用CSV插件加载此类文件，然后将其另存为矢量文件。不幸的是，这很麻烦，因为我有很多栅格文件。 另一个想法是使用gdal2xyz.py，然后使用ogr2​​ogr将其另存为矢量文件。 有什么方法可以分批处理吗？有什么开源程序可以让我更轻松地做到这一点？

10 qgis  raster  gdal  ogr 

我已经成功地将一堆shapefile加载到PostGIS数据库中。我希望能够基于SQL查询以dxf文件的形式检索特定的几何图形和图层。 我正在使用Windows上OSGeo4W套件中的PostGIS，GDAL，PostgreSQL和Python。 我在想应该像这样使用ogr2​​ogr： $ ogr2ogr -sql "SELECT * FROM some_table" -f "dxf" outputfile.dxf -overwrite "PG:dbname=mydb user=dbusername password=dbpassword" 但是当我使用适当的值执行此操作时，会出现一系列如下错误： ERROR 1: DXF layer does not support arbitrary field creation, field 'objectid' not created. 导致该层中每个属性列都出错，some_table但是我确实得到了dxf文件。 如果已经创建了.dxf文件，并且重新运行了上面的命令（包括该-overwrite选项），则会收到此错误： Layer sql_statement not found, and CreateLayer not supported by driver. ERROR 1: Terminating translation prematurely after …

10 postgis  gdal  ogr  dxf 

有没有办法获取使用python的系统上已安装的gdal格式的列表？在命令行上，可以使用以下命令轻松完成此操作： gdalinfo --formats ...但是可以在python中完成吗？

10 python  raster  gdal 

我一直在努力为我们拥有的高分辨率图像生成图块。我们当前拥有的图像是非常大的（+ 20GB）图像，另存为GeoTiff文件。 我想使用gdal2tiles命令行实用程序生成图块，然后使用TMS图像提供程序提供图块，然后在Cesium中打开并查看它。使用gdalinfo，以下是图片的一些详细信息： Driver: GTiff/GeoTIFF Files: image.tif Size is 52250, 56119 Coordinate System is: PROJCS["WGS 84 / UTM zone 35S", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",27], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",10000000], UNIT["meters",1], AUTHORITY["EPSG","32735"]] Origin = (606276.000000000000000,7197873.000000000000000) Pixel Size = (0.500000000000000,-0.500000000000000) Metadata: AREA_OR_POINT=Area TIFFTAG_MAXSAMPLEVALUE=13165 TIFFTAG_MINSAMPLEVALUE=1 TIFFTAG_RESOLUTIONUNIT=2 …

10 gdal  gdalwarp  cesium  gdal2tiles 

我一直在使用GDAL和python来操纵shapefile。打开源文件时，我先将驱动程序设置为“ ESRI Shapefile”，然后使用所需的shapefile调用驱动程序的“ Open”方法。 我对此没有任何问题；但是我很好奇Driver实例方法“ Open”和ogr函数“ Open”之间的区别是什么。例如，两者之间有什么区别？ driver = ogr.GetDriverByName("ESRI Shapefile") source = driver.Open(“O:\myshp.shp”,0) 和 source = ogr.Open(“O:\myshp.shp”,0) 在特定的场景中，您会使用一个场景吗？

10 python  gdal 

我需要创建多个单层栅格图像的多边形轮廓，而不是范围/边界框，而是要创建没有nodata值的区域，如下所示：创建显示 shapepoint的Raster的shapefile？。 在上述问题的答案中，提到了Image Boundary插件，但是，我在运行于Ubuntu的QGIS 1.8.0 Lisboa中找不到它。 该工具仍然可用吗？ 如果没有，是否有办法使用gdal，R，QGIS，GRASS或类似工具而不是ArcMap使用开源工具来做到这一点？

10 gdal  pyqgis  qgis-plugins  footprint 

我想.tiff从.png文件创建地理参考文件。因此，我.tiff使用以下命令行将.png转换为普通文件： convert image.png image.tiff 在这里，一切正常，我有一个有效的.tiff文件，显示正常。然后，我想在其中添加一些地理配准数据。为此，我找到了一个名为的工具tiff2geotiff，我正在这样使用它： tiff2geotiff -4 "+proj=latlong +ellps=WGS84 +datum=WGS84 +no_defs" -c none -n "-122.5575664 38.5818201 -122.5449425 38.5896175" notgeoreferenced.tiff georeferenced.tiff 这给了我黑色图像，该图像在QGis中导入，在适当的地方显示良好，除了在各处都是黑色（证明已正确添加了地理参考数据）。 因此，当我这样做时： gdalinfo georeferenced.tiff 我得到： Warning 1: TIFFReadDirectory:Bogus "StripByteCounts" field, ignoring and calculating from imagelength Driver: GTiff/GeoTIFF Files: araujo_tmp.tiff Size is 1585, 979 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS …

10 qgis  gdal  php  geotiff-tiff 

我有以下格式的纬度，经度和data_val的ASCII文件。 35-13.643782N, 080-57.190157W, 118.6 ... 我有一个GeoTiff图像文件，可以轻松查看它。 我想在ASCII文件中找到的特定纬度/经度位置上，在图像上放置“大头针”（可以是点/标记/星号，或者最简单的东西）。 到目前为止，这是我设法做到的： 我的源图像如下所示： Driver: GTiff/GeoTIFF Files: /tmp/Charlotte SEC 100.tif Size is 16867, 12358 Coordinate System is: PROJCS["Lambert Conformal Conic", GEOGCS["NAD83", DATUM["North_American_Datum_1983", SPHEROID["GRS 1980",6378137,298.2572221010042, AUTHORITY["EPSG","7019"]], AUTHORITY["EPSG","6269"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4269"]], PROJECTION["Lambert_Conformal_Conic_2SP"], PARAMETER["standard_parallel_1",38.66666666666666], PARAMETER["standard_parallel_2",33.33333333333334], PARAMETER["latitude_of_origin",34.11666666666667], PARAMETER["central_meridian",-78.75], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]]] Origin = (-365041.822331817995291,240536.419747152860509) Pixel Size = (42.334586069440391,-42.334898968590878) …

9 python  gdal  latitude-longitude  geotiff-tiff  ascii 

我正在尝试打包一个将GDAL用于PyPI的Python脚本。我首先在我的文章中加入了直接参考setup.py： install_requires=['GDAL==1.11.2'], 这样，程序包无法在我的测试虚拟环境中安装： extensions/gdal_wrap.cpp:2855:22: fatal error: cpl_port.h: No such file or directory #include "cpl_port.h" ^ compilation terminated. error: Setup script exited with error: command 'x86_64-linux-gnu-gcc' failed with exit status 1 然后，我尝试使用对的引用pygdal，因为它被标记为virtualenv友好版本： install_requires=['pygdal'], 这样，安装完成就不会出现错误（但是会加载通常的编译警告）。但是，当我随后调用脚本时，我得到此错误： Traceback (most recent call last): File "/home/desouslu/.virtualenvs/test_p3/bin/hasc2gml", line 5, in <module> from pkg_resources import load_entry_point File "/home/desouslu/.virtualenvs/test_p3/lib/python3.4/site-packages/pkg_resources.py", …

9 python  gdal 

python GDAL模块是否有官方文档？ 有很多食谱，但是关于每个函数或最基本函数的参数的信息很少。

9 python  gdal  documentation 

我正在尝试投影栅格。在R中有对此projectRaster()功能（在一个完全可复制的示例下面）： # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, crs = CRS(newproj), method = 'bilinear') 哪个工作正常。但是，它很慢。 为了提高速度，我还是gdalwarp改用了（使用SSD时，从磁盘/ R读写数据的成本不是很高）。 但是，我无法重现projectRaster()使用的结果gdalwarp： # using gdalwarp to reproject tf <- …

9 raster  gdal  r  coordinate-system  gdalwarp 

已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 2年前关闭。 我正在使用OGR / GDAL在python中构建脚本。 我有一组shapefile和一组GeoTiff栅格文件。 如果脚本不与栅格区域相交，我希望脚本忽略它们。 shapefile不是矩形，因此我不能简单地比较layer.GetExtent（）返回的xmin / xmax，ymin / ymax值。我需要代表其整体形状的实际多边形，然后通过某种方式确定该多边形是否与栅格正方形相交。 我当时想我可以以某种方式将shapefile中的所有多边形合并为一个要素，然后读取该要素上的几何图形，然后将该信息与栅格范围进行比较。但是，我不确定具体如何执行此操作。 如何从shapefile中提取边界多边形信息？ 如何确定该多边形是否与给定的正方形相交？

9 python  raster  shapefile  gdal  ogr 

我像往常一样编译并安装了当前的QGIS母版，并将一些栅格数据导入了具有postgis 2.0.2扩展名的postgres 9.3数据库中。每当我尝试使用数据库管理器将栅格数据导入qgis时，都会收到以下消息： GDAL provider: Cannot get GDAL raster band: Raster layer: Provider is not valid (provider: gdal, URI: PG: dbname=gis host=127.0.0.1 user=xxx password=xxx port=5432 mode=2 schema=public table=elevation 在命令行中，我得到 ERROR 1: Computed PostGIS Raster dimension is invalid. You've probably specified unappropriate resolution. 我不知道在哪里搜索错误。 我刚搬到Postgis 2.1.2，在那里看到了同样的问题。

9 qgis  postgis  raster  gdal  raster2pgsql