Questions tagged «raster»

我正在尝试剪切较大的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 

我对创建一个栅格表面感兴趣，该栅格表面描述了一个小花园地块每单位面积的平均太阳强度。我想模型的输入如下： 地块面积-为简单起见，让我们假设花园地块是圆形的。 周围森林边缘的平均树高 各种太阳角下的太阳强度 随时间推移相对于花园地块的平均太阳角度。 我最终想在树木环绕的花园中显示植物生长的最佳位置。谁能指出我的解决方案？尽管我找不到任何参考资料，但我敢肯定有人已经这样做了。 编辑： 出于示例目的，所附图像显示了一个1英亩的花园地块，周围是针叶树的平均年龄。为了进一步简化此分析，假设我们只对6月至10月园艺季节的太阳天数和强度感兴趣。砍树不是一种选择！

9 arcgis-desktop  python  arcpy  raster  r 

我想知道是否有办法找到具有特定值的栅格的最小边界范围。我从全局图像中裁剪了一个栅格，并将范围设置为具有大量NoData区域的全局范围。我想从此栅格中删除NoData区域，并仅保留包含特定值像素的区域。我怎样才能做到这一点？ 这是我的示例：我想提取value = 1（蓝色区域）并使用蓝色区域的范围而不是整个世界进行进一步处理。

9 arcgis-desktop  raster  extents 

因此，看起来很简单，但我似乎无法使其正常工作！我正在尝试使用linux上的curl通过REST API将栅格数据文件（GeoTIFF）添加到geoserver。REST API 可以正常工作，因为我使用它创建了一个新的工作区： curl -u admin:geoserver -v -XPOST -H 'Content-type: text/xml' \ -d '<workspace><name>restProba</name></workspace>' \ http://localhost:8080/geoserver/rest/workspaces 如Geoserver的admin ui所示，它是正确创建的，但是当我尝试通过执行将coveragestore添加到该工作区时 curl -u admin:geoserver -v -XPOST -H 'Content-type: text/xml' -d '<coverageStore><name>int_dec</name><enabled>true</enabled> <type>GeoTIFF</type> <url>$home/int_dec.tif</url></coverageStore>' "http://localhost:8080/geoserver/rest/workspaces/restProba/coveragestores?configure=all" 即使已正确加载，它也会直接进入默认工作空间，这不是我想要的工作空间。是否有任何方法定义我想要的工作空间？我认为只要在指向restProba的工作区URI的URL中提及它就足够了，但似乎忽略了它。 谢谢你的帮助 ：-） pd：与geoserver的新手，我已经尝试过搜索，但是要么我没有使用正确的搜索条件，我就傻了，或者对于任何人来说我都太简单了，我甚至愚蠢地要求它：p

9 geoserver  raster  rest 

我正在使用以下命令将栅格加载到postgres数据库中： raster2pgsql -I -C -e -Y -F -s 3086 -I -C -M myraster.tif myraster -F -t 30x30 | psql -U postgres -d database -h localhost -p 5432 myraster.tif是155M。当我尝试在QGIS中使用栅格插件时，加载时间很长（15-20分钟）。我不确定是否可以采取任何措施来加快此速度，或者不确定是否可以使用其他方法查看栅格。我要做的就是快速查看栅格，并检查栅格是否与数据库中存储的矢量数据正确对齐。（相比之下，矢量数据显示得非常快。）此刻，QGIS冻结并显示栅格。它只是显示一个黑色正方形，并且由于任何交互都会导致程序冻结，因此更改色阶是不可行的。 我不熟悉与加载栅格相关的细微差别，因此也许应该使用不同的标志和参数。我意识到这是一个很大的栅格，所以我很乐意为显示目的重新采样它或类似的东西（尽管我不知道该怎么做）。 我正在wktrasterQGIS中使用最新版本的插件。我正在连接到远程服务器以显示数据。我将OSX与QGIS 1.8一起使用。我愿意使用任何程序在任何OS上查看栅格。

9 qgis  postgis  raster  postgresql 

我在ArcGIS Desktop中有一个GeoTIFF栅格，可以从以下两个链接下载该栅格： http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tfw http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tif 如果将它们加载到ArcGIS Desktop中，将会看到较大的统一区域。使用“信息”工具单击这些区域将显示像素值为“无数据”。但是，当您选择在图层的“符号”首选项中将“无数据”值显示为透明（或任何颜色）时，“无数据”值将不遵循该规则。 有趣的是，图像主要区域之外的“无数据”值确实遵循该规则（将颜色设置为红色即可看到该规则），但图像内部的值则不然。 我真的很难尝试解决此问题-根本无法解决该问题。 有人有什么想法吗？

9 arcgis-desktop  raster  arcmap  arcgis-10.0  spatial-analyst 

我的背景：经验丰富的程序员，业余历史学家，GIS初学者（所以我正在使用QGIS）。我确实学会了使用QGIS地理参考工具。 输入：假设我有一个类似Wikipedia的图像：File：VoyagesOfRabbanBarSauma.jpg。假定它是一个“裸”的JPG文件，没有任何地理元数据，并且作者无法提供任何此类信息。输入图像已获得GNU FDL的许可，我很高兴在此发布我的衍生作品。 目标：将Wikipedia地图中显示的路径导入QGIS，以便我可以更改投影，叠加高程，为自然地理添加标签等。 警告：一种解决方法可能只是将地图中位置的纬度/经度输入到新的GIS项目中。我不想走这条路，因为我想记录该特定作者对地理信息的估计。 我尝试过的操作：我已将QGIS Georeferencer应用于栅格JPG文件，并在各大洲的边缘有很多点。大陆内部的点对齐一般。 我想知道的是：这张特定地图的扭曲在欧亚大陆的北部边缘非常明显，而且看起来像是标准的投影。在我看来，应该有一种算法方法来识别生成该图的投影和投影参数，并恢复图像像素与经纬度之间的准确非插值转换？ 我对此有两个问题： 在未知时识别Shapefile的坐标系？ 旧图纸上的未知坐标系 这些人似乎在谈论涉及GDAL或QGIS的蛮力方法，而对于GIS初学者来说，似乎有些劳动强度（基本上是地理配准，并且使其适合已知的投影），但我愿意尝试一下，看看是否可以使它们自动化。 这两个GIS.SE上的说明都比我的经验所允许的要高一些，关于如何完全使用这种蛮力方法的任何提示将不胜感激。例如，上面的问题2显示了如何使用gdaltransform将坐标列表转换为经/纬度。输入坐标是xy像素位置吗？

9 qgis  raster  coordinate-system  georeferencing 

我有两个栅格图像，按布尔分类（1或2）。它们在给定区域重叠。 我需要只有重叠部分的栅格。 使用Arcinfo10。Answer只能使用栅格方法，不能转换为多边形，点等。

9 arcgis-10.0  raster 

关于我的上一个问题，如何使用R将生成的栅格图像“地理编码” /转换为生成的栅格图像，使其可以显示为地图上的图层？简而言之，我的目标是重新创建商业SpatialKey（屏幕截图）软件的结果。 我正在使用GeoServer实例来托管地图数据，但是如何将所需的投影和坐标信息合并到图像中，以使其正确地与真实地理位置对齐。密度热图是使用空间数据生成的。

9 raster  r  georeferencing 

我正在尝试使用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 

我有一个有关意大利北部土地使用情况的栅格文件（作物/市区等）。 我在矢量层的某些给定点周围定义了一些缓冲区，我想计算栅格文件的每个值所占面积的百分比，以便估算所选区域的平均使用率（例如70％的农作物，30％的城市等）。 我知道在ArcGIS中有一个称为“ 制表区 ” 的Spatial Analyst函数，但是我在QGIS中找不到这样的功能。 如果有人可以帮助我，我将非常感激。

9 qgis  raster 

我有一个裁剪的栅格数据集，上面覆盖了折线（来自轮廓工具）（请参见附件）。我正在尝试将折线裁剪到栅格数据集的范围，这似乎很简单，但是我尝试了很多事情，但似乎无法使其工作。我最终将需要使用arcpy对此进行编码，但是我认为现在，只是朝着正确的方向轻推一下哪种工具可以完成此操作，这对我现在有帮助。 谢谢。 PS我正在使用ArcInfo 10.0

9 arcgis-10.0  raster  arcpy  vector  clip 

我有一个山区的DEM。我想提取DEM的区域，在该区域之上，占DEM总面积的60％。因此，基本上，我想提取代表山脉上部的区域。 我听说过这条曲线称为“测压线”（H60线），但尚未找到有关如何实现此功能的任何文档。 我的直觉是要以某种方式计算DEM的统计数据，以告诉我60％的水平在哪里，然后在此高程处绘制DEM的轮廓。然后，我可以通过该轮廓裁剪DEM，并在输出DEM子集上执行计算。 我如何确定这个60％的水平？

9 arcgis-desktop  raster  arcgis-10.0  dem 

在ArcMap 9.3中，我使用了“内核密度”来绘制各种事件，但是生成的shapefile不显示任何度量单位。是否有很好的非技术来源可以用非专业术语来解释根据输入像元大小和搜索半径对输出值的解释？

9 arcmap  raster  density 

我对一种识别损坏的栅格（例如ESRI GRID格式）或shapefile的方式感兴趣，这将阻止它们在分析工具中使用。我在结构化文件夹树中有大量的文件，并且希望能够编写脚本，使文件树走动并能够识别栅格，例如，如果您尝试在ArcCatalog中查看它，则将不允许其显示并因此会使任何地理处理工具失败。shapefile或要素类的处理相同。 我的喜好是以某种方式在Python中执行此操作（并且我不会反对使用非ESRI python模块或其他模块），但是我愿意接受任何东西。我想您可以通过尝试验证几何形状或利用文件设计的某些方面来一起破解某些东西，但是最有意义的是什么？还是有另一种专门为此目的设计的方法？ 我的目标是仅能够对一些为期一周的大型预处理文件运行一种简单的质量控制机制，以确保在开始进行实际的批处理之前它们都声音良好。 谢谢汤姆

9 raster  shapefile