Questions tagged «gdal»

GDAL(地理空间数据抽象库)是一个开放源代码翻译器库和一组命令行实用程序,用于处理栅格空间数据格式。作为一个库,它为所有受支持的格式向调用应用程序提供一个抽象数据模型。它还带有用于数据转换和处理的各种有用的命令行实用程序。

13
在Windows上使用Python安装GDAL?
谁能解释如何在Windows上使用Python安装GDAL / OGR? 我使用Windows Vista,并且尝试按照网站上的信息进行操作,但似乎没有为我绑定exe文件。 有人可以描述该过程,包括指向我需要的文件/文件夹的链接吗? 我现在尝试使用minGW运行gdal设置,但这也失败了:

23
通过pip安装时,Python GDAL软件包缺少头文件
我意识到这个网站存在之前,我曾在Stack Overflow上问过这个问题。由于该站点更合适,因此我也会在这里提出要求。 我正在尝试从pip install gdal虚拟环境(Ubuntu)中的pip安装gdal 。失败,因为找不到cpl_port.h extensions/gdal_wrap.cpp:2853:22: fatal error: cpl_port.h: No such file or directory compilation terminated 但是GDAL已正确安装,并且头文件位于/usr/include/gdal/cpl_port.h。是否需要为GDAL设置一些环境变量,以便pip查找头文件?

8
是否应将GDAL设置为通过压缩生成GeoTIFF文件?应该使用哪种算法?
我有一个主要由GeoTIFF文件组成的GIS数据文件夹。整套大约重约1.2 GB。我注意到,如果将内容包装到tarball中,它将粉碎到大约82 MB。我想将其检入修订控制系统中,以便其他人可以使用它,并且似乎可以挤出一些空间。 GDAL GeoTIFF 驱动程序页面列出了许多可用于创建压缩的GeoTIFF文件的选项。还有很多选项会影响每种算法的工作方式。 帮助页面在描述选项方面做得很好,但是没有详细说明如何选择与变化的压缩级别相关的算法或权衡取舍。这导致以下问题: 使用压缩的优点是可以节省大量空间。缺点是什么?压缩图像时信息会丢失吗? 如何选择一种算法和压缩级别。某些类型的图像是否适合某种算法?

3
Lanczos重采样在空间环境中有什么用?
GDAL除了近邻,双线性,三次和样条曲线的常规混合以外,还包括一种重采样方法:“兰佐斯窗口正弦重采样”。我知道它是一个卷积滤波器,但与结果倾向于主观的图像不同,用于空间数据的重采样还有其他含义。什么是Lanczos,如何使用Lanczos会影响输出?

1
在没有NumPy的情况下获取OGR点下的GDAL栅格的像素值?
我正在研究整个景观中大量野生传粉媒介的计算模型。该模型本身已经完成,现在我正在努力进行后处理步骤。 我的GDAL授粉媒介供应栅格看起来像这样(较浅的颜色表示授粉媒介对像素的访问量较高): 我有一个OGR形状文件,该文件的点代表景观上的样本位置: 我正在尝试对这些点下的像素进行一些分析,但是要这样做,我需要能够提取一个点下的像素的值。 是否可以通过Python仅使用OGR和GDAL提取点下的像素值?我宁愿避免通过读取整个栅格到内存中ReadAsArray(),因为我的输出栅格非常大(太大而无法容纳到内存中)。 我注意到这篇文章,这很相似,但是需要命令行调用。
45 python  gdal  ogr  pixel 


5
如何从Python代码调用gdal_translate?
是否可以使用某些gdal API gdal_translate从Python代码调用?我的意思不是简单地从文件系统执行gdal_translate.exe,而是以某种方式在代码中调用它,所以我不需要知道gdal_translate可执行文件所在的确切目录?
40 python  gdal 

2
如何检查我拥有哪个版本的GDAL?
我在计算机上安装了以下软件: PosgresSQL 9.3.2.3 PostGIS 2.1.1.1 QGIS 2.0.1 Dufour(通过OSGeo4w软件包安装的一个) QGIS 1.7.0弗罗茨瓦夫 据我所知,它们都以某种方式操作,使用或访问GDAL库。 因此,是否根据我使用的上述软件之一使用了GDAL的独立版本?例如,如果我打开OSGeo4W Shell,则在加载批处理后立即显示: GDAL 1.10.1, released 2015/02/10 那么,如何检查其他GDAL版本并进行比较? 有什么办法可以集中所有这些内容,并强制应用程序使用相同的集中文件夹?

3
创建透明的山体阴影?
通过组合使用gdal,我得到了优雅的透明山体阴影convert。与通常使用的基于灰色的山体阴影相比,这种透明的山体阴影非常酷,因为它们可以放置在地图背景与其他上层(道路,建筑物)之间,以提供3D感觉,无论背景的类型和颜色如何。 怎么运行的 诀窍:从产生的灰度山丘阴影开始gdal hillshade,诀窍是获取灰度通道的值,将每个值求反,然后将结果流到新的不透明度通道中。黑色像素[0,0,0]变为[0,0,0,255](不透明度= 255),灰色像素[120,120,120]降低至[120,120,120,135](不透明度= 135,又称为255-120),白色像素[255,255,255]变为透明[255,255,255,0](不透明度= 0,又称为255-255),依此类推。丘陵的阴影是不透明的和黑色的,平原逐渐变得(白色)透明。从概念上讲,像素的波段方程是这样的: alpha_value([0]) = 255-0 = 255; // black pixels get full opacity alpha_value([Grey]) = 255-Grey; // grey pixels get 255-GreyValue opacity. alpha_value([255]) = 255-255 = 0 // white pixels get null opacity ! 有关此方法的悠闲视频(由Photoshop设计师解释),请参阅在Photoshop中添加阴影浮雕(16分钟)。 题 将ETOPO或SRTM派生的基于灰色的山体阴影(gdaldem hillshade crop_xl.tmp.tif shadedrelief.tmp.tif -s 111120 -z 5 -az …


5
如何通过切割180子午线将栅格从0 360重新投影到-180 180
我有一个Geotiff栅格图像,它的坐标系的纬度从0到360。图像的水平中心是180度。见下图: 我想将其转换为经度范围为-180 180的EPSG:4326 SRS。我希望图像的中心位于格林威治子午线(0)。我想这个SRS用途非常广泛。我希望结果看起来像这样: 因此,我使用gdalwarp命令重新投影: gdalwarp -s_srs '+proj=latlong +datum=WGS84 +pm=180dW' -t_srs EPSG:4326 test_col.tif test_4326.tif 但是我只会看到尺寸更大(像素更多)和EPSG:4326元数据的tiff。图像本身看起来与初始图像相同。但我希望它能交换半球。 问题是-如何将图像的经度严格调整为-180 180 EPSG:4326,且中心在0度的经度? 这是我的初始文件的gdalinfo: Origin = (-0.102272598067084,89.946211604095552) Pixel Size = (0.204545196134167,-0.204423208191126) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( -0.1022726, 89.9462116) ( 0d 6' 8.18"W, 89d56'46.36"N) Lower Left ( -0.1022726, -89.9462116) ( …

5
如何使用Python GDAL绑定获取栅格角坐标?
有没有办法使用gdal的Python绑定从栅格文件中获取角坐标(以纬度/经度为单位)? 网上的一些搜索使我确信没有这种搜索,因此我通过解析gdalinfo输出进行了开发,它有些基本,但是我认为这可能为那些不太熟悉python的人节省一些时间。仅当gdalinfo包含地理坐标和角坐标时,它也才有效,我认为并非总是如此。 这是我的解决方法,请问有人有更好的解决方案吗? 在适当的栅格上的gdalinfo会在输出中途产生如下所示的结果: Corner Coordinates: Upper Left ( -18449.521, -256913.934) (137d 7'21.93"E, 4d20'3.46"S) Lower Left ( -18449.521, -345509.683) (137d 7'19.32"E, 5d49'44.25"S) Upper Right ( 18407.241, -256913.934) (137d44'46.82"E, 4d20'3.46"S) Lower Right ( 18407.241, -345509.683) (137d44'49.42"E, 5d49'44.25"S) Center ( -21.140, -301211.809) (137d26'4.37"E, 5d 4'53.85"S) 此代码将在gdalinfo看起来像这样的文件上工作。我相信有时坐标将以度和小数表示,而不是度,分和秒。为这种情况调整代码应该是微不足道的。 import numpy as np import subprocess …
30 python  gdal 

5
将numpy数组写入栅格文件
我是GIS的新手。 我有一些代码将火星的红外图像转换成热惯性图,然后将其存储为2D numpy数组。我一直将这些地图另存为hdf5文件,但我真的很想将它们另存为栅格图像,以便可以在QGIS中对其进行处理。我经历了多次搜索,找到了如何执行此操作的方法,但是没有运气。我尝试按照http://www.gis.usu.edu/~chrisg/python/上的教程中的说明进行操作,但是当我将其示例代码生成的文件导入QGIS时,它们以纯灰色框打开。我觉得如果有人可以将最简单的过程建议为我想做的事的简化示例,那么我也许可以取得一些进展。我有QGIS和GDAL,很高兴安装任何人都可以推荐的其他框架。我使用Mac OS 10.7。 因此,例如,如果我有一个numpy的热惯性数组,看起来像: TI = ( (0.1, 0.2, 0.3, 0.4), (0.2, 0.3, 0.4, 0.5), (0.3, 0.4, 0.5, 0.6), (0.4, 0.5, 0.6, 0.7) ) 对于每个像素,我都有纬度和经度: lat = ( (10.0, 10.0, 10.0, 10.0), ( 9.5, 9.5, 9.5, 9.5), ( 9.0, 9.0, 9.0, 9.0), ( 8.5, 8.5, 8.5, 8.5) ) …
30 python  raster  gdal  numpy 

4
学习GDAL的资源?
我经常被困于尝试解决特定问题或正确使用语法。 除了GDAL实用工具上的链接中提供的详细信息外,还有谁知道任何第三方资源/手册以更加用户友好的方式对其进行说明?
29 gdal  references 


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.