Questions tagged «raster»

网格，TIF和IMG是栅格文件的三种格式。 在ArcGIS Desktop中，它们之间显然没有区别。 这三种格式有什么区别？

17 raster  file-formats 

我有一个带有多边形的shapefile。而且我有一个全局栅格文件。我想将shapefile的多边形覆盖到栅格网格上，并计算每个多边形的平均栅格值。 如何使用GDAL将结果写入shapefile？

17 python  raster  shapefile  gdal  spatial-statistics 

我想提取栅格图像的某些值以创建一个新的栅格，然后该栅格应仅包含提取/选择的值。 如何在QGIS中做到这一点？

17 qgis  raster  extract-by-mask 

我正在尝试处理一些生物气候栅格文件，例如可以从http://www.worldclim.org/current（bioclim set）下载的文件。他们似乎没有-3.4e+38根据QGIS 设置的无数据值（查看gdalinfo的输出是-3.39999999999999996e+38）。 似乎gdal工具无法处理此nodata值，而qgis似乎也无法识别它。在图层样式中，有一个-3.4e + 38项设置为100％透明，但是即使“识别要素”选择器将其显示为值-3.4e + 38，它仍然显示此类值。 我试图创建一个vrt来将nodata值转换为-9999，但这还是行不通的。 如何处理此类文件以具有可用的nodata值？

17 qgis  raster  gdal 

这是一个约13公顷的湖泊，在约10条样线中有81个采样深度点： 以前在ARCgis中有一个50公顷的湖泊和大约100个深度采样点，我使用该TopotoRaster工具获得了不错的输出。但是，QGIS或开源软件似乎不存在这样的精确等效项。 在QGIS中，我使用Raster Interpolation插件中的TIN方法获得了下面的插值栅格。但是，这似乎是一个幸运的结果，当我尝试重复时，尽管没有真正改变任何参数，但我得到了不同的结果（不令人满意）。 这些比城市湖泊更多的城市（构造）水箱，因此，它们的河床轮廓相当规则，因此采样点密度低。一些文章似乎建议使用逆距离加权（IDW）作为最合适的技术，但是我似乎获得了最差的结果（可能是由于低点密度）。 考虑到湖泊的大小，其规则轮廓和点采样密度，可以采用任何启发式方法来得出哪种插值方法合适？（在克里金，IDW，双线性，三次卷积，TIN或样条之间） 还是总是反复试验？

16 qgis  raster  interpolation  bathymetry  lakes 

我有一层coord. ref. : NA。 我该如何分配coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0给coord. ref. : NA？ class : RasterLayer dimensions : 127, 146, 18542 (nrow, ncol, ncell) resolution : 0.01, 0.01 (x, y) extent : -122.39, -120.93, 36.72, 37.99 (xmin, xmax, ymin, ymax) coord. ref. : NA data source : C:\2012Tif\2012ASC5min_ppt_spas1306_0001_20121222_0005_UTC.asc.tif …

16 raster  coordinate-system  r 

我在postgresql / postgis中有一个带有以下列的栅格数据库： （ID，rast，data_of_data）。 “ rast”是具有WKT格式的栅格文件的列。以下示例查询用于查找WGS84系统（30.424，-1.66）中某个点的DN值（对于2002-01-09）： SELECT st_value(rast,(st_GeomFromText('POINT(30.424 -1.66)', 4326))) as val FROM my_table WHERE date_of_data='2002-01-09' 是否有一种方法（例如空间索引）来加快此类查询的速度？

16 postgis  raster  query  performance  spatial-index 

我有两层。包含许多图块的多边形形状层和包含CORINE 2006土地覆盖物的栅格层，并在色图中包含许多类别。我想为shapelayer中的每个多边形获取栅格图层的每个土地覆盖类别的总和。 例如，有一个ID为'2'的多边形，我要为此多边形设置这样的属性（以百分比或平方米为单位）： 耕地：15％ 森林：11％ 街道：2％（...等等） 我试图在草，qgis（无功能），saga（将每个总和求和成一个总值）r（总和）中做到这一点，但我仍然没有找到解决方案。大多数插件（qgis中的区域统计信息）仅支持0-1栅格图层。v.rast.stats也没有帮助。Iam接受任何优秀而明智的解决方案！也许我什至使用了错误的方法或犯了错误。 如果记得正确的话，在Arcgis中，此任务非常容易，但是对于您的日常linux用户，我仍然缺少一个好的解决方案。 我正在运行一个debian linux系统，这就是为什么我只能为此操作系统使用程序的原因。 编辑：因为这个问题仍然有这么多的视图和访问者：我写了一个QGIS插件，它也能够计算栅格图层的土地覆盖。我还没有编写多边形叠加层的代码，但是它确实经过了规划。在此处找到插件，然后首先安装Scipy库。

16 qgis  raster  polygon  r  statistics 

我们刚刚收到了许多工作中的DEM，我想从中生成轮廓。DEM的分辨率为1m，大小为1kmx1km。 gdalinfo的输出： Driver: AAIGrid/Arc/Info ASCII Grid Files: 380000_6888000_1k_1m_DEM_ESRI.asc Size is 1000, 1000 Coordinate System is `' Origin = (380000.000000000000000,6888000.000000000000000) Pixel Size = (1.000000000000000,-1.000000000000000) Corner Coordinates: Upper Left ( 380000.000, 6888000.000) Lower Left ( 380000.000, 6887000.000) Upper Right ( 381000.000, 6888000.000) Lower Right ( 381000.000, 6887000.000) Center ( 380500.000, 6887500.000) Band …

16 raster  gdal  dem  contour 

给定DEM（数字高程模型）和一条线（x1，y1）->（x2，y2）（成对坐标（x，y）），如何计算在该线上投影的高程轮廓或横截面？ 我正在寻找可以在项目中使用的开源代码或可以用作指南的伪代码。

16 raster  dem  open-source-gis  elevation  profile 

有没有一种方法可以计算到QGIS中到最近点的距离？在ArcGIS中，有一个称为Euclidean Distance的函数，该函数可以获取输入矢量或栅格，并生成距离栅格。 QGIS中有类似的功能吗？我在文档中找不到一个。我确实研究了RasterQGIS中的功能以及MMGIS中的Hub Distance计算（对于非光栅解决方案），但是我无法弄清楚如何使用这两种方法。如果向量方法是唯一的解决方案，那也可以。

16 qgis  raster  distance  point 

我是Python的新手，我想知道是否存在一种快速方法来在ArcGIS 10中使用Python逐个像素获取栅格像元值以及坐标（每个像素中心的地图XY坐标）？ 为了进一步描述这一点，我需要获取第一个像素的地图X，地图Y和像元值，并将这三个值分配给三个变量，并对其余的其他像素重复此步骤（遍历整个栅格）。 我想我需要更多地描述我的问题。问题是，我需要获取第一个栅格的像素的XY位置，并获取与该XY位置相对应的其他几个栅格的像元值。此过程应遍历第一个栅格的每个像素，而无需创建任何中间点shapefile，因为这将真的非常耗时，因为我必须处理具有近80亿像素的栅格。另外，我需要在ArcGIS 10中使用Python进行此操作。 @JamesS：非常感谢您的建议。是的，这将适用于一个栅格，但是我还需要收集其他几个栅格的像元值。问题是，在获取第一栅格的第一个像素的X和Y坐标之后，我需要获取第二栅格的像元值，该值对应于第一栅格的X，Y位置，然后是第三栅格，依此类推。因此，我认为在遍历第一个栅格时，应该同时完成像素的X和Y位置以及与该位置相对应的其他栅格的像元值的获取，但我不确定。这可以通过将第一个栅格转换为点shapefile并在ArcGIS 10中执行将多值提取为点功能来完成，但是我 @hmfly：谢谢，是的，如果我可以获取数组的已知行和列值的坐标，则此方法（RastertoNumpyarray）将起作用。 @whuber：我不想执行任何计算，我要做的就是将XY坐标和单元格值写入文本文件，仅此而已

16 arcgis-10.0  python  raster 

如何在QGIS中获取栅格中像元值计数的摘要？我知道QGIS就是这样做的，以便在层的属性中提供摘要统计信息和直方图。我不需要摘要统计信息或直方图图像，而无需计算单个值。 似乎曾经尝试过为此添加特殊功能（功能＃4321，http://hub.qgis.org/issues/4321）。 我是否真的缺少明显的东西？

16 qgis  raster 

您会建议使用哪种CRS来对历史世界地图进行地理参考，例如伦敦地理学院-人民地图集-1920中的地图，例如 该地图显示一个纬度/经度网格，我想用它来选择控制点。我只是使用WGS84 EPSG：4326设法获得了一些合理的结果，但是想知道是否没有更好的选择。

15 raster  georeferencing 

这似乎应该很简单，但是我无法一生解决。 给定一个包含单个波段的现有GeoTiff，如何将另一个波段添加到该GeoTiff？这样做的唯一方法是使用GDALCreate（）创建具有附加波段的新输出文件，然后将现有波段从原始文件复制到新文件吗？这是我到目前为止提出的唯一解决方案，但似乎应该有一种简单的方法可以将波段添加到现有数据集中。

15 raster  gdal  band