Questions tagged «raster»

栅格是一种由规则的值网格组成的数据格式，通常以类似图像的格式存储。

4
使用GDAL平滑/重新插值栅格吗？
简短： 我想更改光栅的分辨率并平滑灰色渐变，如下面的图像所示。首选项是使用GDAL，PIL或Numpy。 说明： 我正在使用High Performance Geostatistical Library将该点克里格化为输出分辨率为20米的栅格。我不想更改输出分辨率，因为插值时间呈指数增长。 在这种分辨率下，输出图像很难看（像素化和混叠）。我不知道它在概念上是否正确，但是我希望图像像下面的示例一样更平滑。就像将图像“重新插值”为分辨率更高的图像一样。我正在使用python，因此我的偏好设置是GDAL，Python Imaging Library或Numpy。答案可能是理论上的，例如指出算法名称或这种操作的概念。 资源： 目的地： 使用gdalwarp三次样条的编辑结果：
15 python  raster  gdal 
4
使用ArcGIS Desktop手动编辑栅格以改善分类？
我有一个来自Landsat TM影像的土地覆盖分类-我使用ENVI创建了该分类。我想在分类的某些区域中编辑单个或一组单元格，以便通过与航空影像进行比较来提高输出的准确性。 有没有一种方法可以使用ArcGIS Desktop通过选择大型单元格并将它们分配给新的或现有的类来手动编辑它们？ 我知道基于类统计信息的合并方法，但是由于云层薄，我要分类的区域存在大量光谱混乱-似乎手动编辑是唯一的方法。
3
使用PostGIS将栅格与多边形相交-Artefact Error
我正在使用PostGIS2.0进行一些栅格/多边形相交。我很难理解应该使用哪种操作以及执行此操作最快的方法是什么。我的问题如下： 我有一个多边形和一个栅格 我想找到落在多边形内的所有像素，并获取像素值的总和 并且（已更新的问题）：执行查询时，我正在获取原始栅格中不存在的一些像素的大量值 我难以理解应该使用ST_Intersects()还是ST_Intersection()。我也不知道对像素求和的最佳方法是什么。这是我尝试过的第一种方法（＃1）： SELECT r.rast FROM raster as r, polygon as p WHERE ST_Intersects(r.rast, p.geom) 这将返回一个rast值列表，我不确定该怎么做。我尝试使用计算汇总统计信息，ST_SummaryStats()但不确定这是否是多边形内所有像素的加权和。 SELECT (result).count, (result).sum FROM ( SELECT ST_SummaryStats(r.rast) As result FROM raster As r, polygon As p WHERE ST_Intersects(r.rast, p.geom) ) As tmp 我尝试过的另一种方法（＃2）使用ST_Intersection()： SELECT (gv).geom, (gv).val FROM ( SELECT ST_Intersection(r.rast, p.geom) AS …
1
使用R将多边形转换为栅格
我正在尝试将shapefile转换为R中的栅格。我的方法是按如下方式读取栅格： library(rgdal) # Loads SP package by default demo <- readOGR('F:/data/', 'shapefile') # Creates a SpatialPolygonsDataFrame class (sp) 这个工作正常，我可以绘制它。但是，它是一个很大的shapefile，我想将其转换为栅格。我尝试了以下方法： r <- raster(ncol=180, nrow=180) Demo_ras = rasterize(r, demo, 'pop') # pop is an integer here 错误信息： Error in function (classes, fdef, mtable) : unable to find an inherited method for function …
15 raster  shapefile  r 
5
在ArcMap中使用Hillshade选项？
这是我使用的栅格的漂亮图片-没有应用任何效果-只是调整了色带和直方图以显示所有浮雕-尽管是以2D方式进行的！ 我有一个栅格图层，要使用默认的hillshade显示。我可以在代码中执行此操作，实际上与在ArcMap中选择Hillshade选项相同（请参见下文）： 但是，当您执行此操作时-看起来很令人失望： 尤其是如果放大栅格，您会看到各种网格状的人工制品-看起来一点也不好！ 但是，如果我使用带有Spatial Analyst的栅格来创建正确的Hillshade图层，并以85％的透明度对其进行覆盖，则它看起来会非常好并且非常不同！ 如何在ArcMap中实现此效果并覆盖Hillshade选项的默认功能（我认为这是不可能的，或者我错过了某些设置而导致栅格显示不佳）。如果可以覆盖或选择正确的设置，则可以在ArcObjects中进行。 如果无法使用选项1，我可以通过创建一个Hillshade图层（不使用Spatial Analyst）来创建效果并将其有效地融合到原始栅格中吗-即，我不希望人们在ArcMap的目录中看到两层。我只想要一层，但是具有默认的我的选择/配置的hillshade渲染。我也需要以编程方式执行此操作，即ArcObjects。 我了解到的一件事是，您必须在正确的空间参考和测量单位中使用DEM，否则生成的山影看起来太黑且过于夸张-请观看此视频以了解更多信息-这是一个很好的解释。 因此，我想我现在知道在生成一些高质量的阴影时需要做些什么，但是我认为这将需要一些底层的ArcObjects代码调整。我想我需要提出自己的图形渲染并将其插入ArcMap（我有一个dll，可以很酷的方式渲染网格的位图）。我认为显示此内容的一个好地方是重写IRasterRenderer之类的Draw（）事件，并插入我自己的图像渲染。当然，使用重投影等可能会变得有些复杂。 是否有人在ArcMap中完成过这种低级图形渲染？如果可以，您可以提供任何指导或我应该注意的事情吗？
2
传单以错误的顺序渲染TileMill瓷砖
我很困惑，很多搜索都没有帮助。 我处于最终将是具有三个栅格层和三个标记层的地图的第一步（带有在栅格及其关联数据之间切换的控件，以及用于数据的弹出窗口/工具提示），但我不能获取Tilemill瓷砖以正确渲染。（我在QGIS中创建了栅格，然后通过GDAL对栅格进行着色和重新投影，然后将其放入TileMill。） 智慧： 在Chrome中。这是Safari： 这是代码，对我来说似乎是无害的（此时，我已经注释掉了除了磁贴之外的所有内容，所以这就是我所包含的全部内容）： <html> <head> <title> - Leaflet demo</title> <script type='text/javascript' src="http://cdn.leafletjs.com/leaflet-0.7.3/leaflet.js"></script> <style type='text/css'> html, body{ width:100%; height: 100%, margin:0; padding:0; } #map{ width:700px; height:600px } img {margin:0px;padding:0px}; </style> <script type='text/javascript'>//<![CDATA[ window.onload=function(){ //one base layer on map var map = new L.Map( "map" ).setView([38.173,-94.384],5); L.tileLayer('http://a.tiles.mapbox.com/v3/mtgmediagroup.beef2013/{z}/{x}/{y}.png', { attribution: 'Dani …
1
在QGIS中访问栅格属性表？
我有一些使用GDAL为GeoTiffs制作栅格属性表的Python代码。 在ArcGIS中，我可以很好地看到这些属性表，并且它们表现出我喜欢的方式（作为分类），但是在QGIS中，似乎缺少属性表的整个概念。 有什么方法可以访问QGIS中的RAT？
2
了解“镶嵌到新栅格”的最小值和最大值吗？
我必须从以下来源将多个DEM栅格合并为一个：http : //srtm.csi.cgiar.org/SELECTION/inputCoord.asp，因此我使用了“镶嵌到新栅格”工具。 我设置了工具： 我输入了栅格（全部来自相同的源，相同的大小且没有投影）； 设置像素类型（与原始栅格中一样为32位浮点）； 将像元大小设置为原始栅格； 带数= 1作为原始栅格；马赛克运算符，我已经用BLEND和MEAN完成了（我得到了相同的结果）。 我遇到的问题是，生成的镶嵌图显示的最大值和最小值范围与单个栅格的最大值和最小值不同，例如，栅格1（-5123.8，23.25），栅格2（-5974.6，40.09），栅格3 （-57770.2，38），栅格4（-2534.3、23.55）和最终镶嵌栅格（-5975.8，81.1）。 我猜这种解决方案是不正确的，至少我没想到会得到这种解决方案。任何人都知道是否可以，如果不可以，如何解决并获得具有正确的最大和最小值的合适的镶嵌栅格？ 我正在使用ArcGIS 10.2.2 for Desktop。
3
R：如何从RasterLayer获取经度和纬度？
我绝对是地理数据的初学者，因此，如果问题不适当，请原谅我。 我从NCDC NARR下载了数据，并设法使用该raster软件包将其加载到R中。我想获得一个包含纬度，经度和价值的列表。我知道该功能rasterToPoints()可以完全满足我的要求，但是，我的纬度和经度值看起来很奇怪： r <- raster(myfile) data_matrix <- rasterToPoints(r) head(data_matrix) x y value [1,] -5405401 4347242 70 [2,] -5372938 4347242 88 [3,] -5340475 4347242 76 [4,] -5308012 4347242 85 [5,] -5275549 4347242 87 [6,] -5243086 4347242 88 我想我应该对当前为Lambert Conformal Conic（LCC）的投影做些什么。以下是有关栅格的更多信息。 > r class : RasterLayer dimensions : 277, 349, 96673 …
14 raster  r 
1
PostGIS栅格求和（地图代数）
我有一个表示特定日期旅行时间等时线的多边形表。对于每个原点，有五个等时线几何体（存储在单独的行中）。对于每个原点，我要栅格化五个等时线（二进制NULL或1），然后将它们组合成一个栅格图层。此栅格图层需要一个简单的地图代数：sum / 5，以便最终每个原点都将与一个单独的栅格图层相关联，该栅格图层的值分别为[NULL，0.2、0.4、0.6、0.8、1]，具体取决于多少组成层重叠。这是一个概率面。 我的数据全部存储在Postgres 9.3（带有PostGIS）中。我的问题是，尽管我想学习使用PostGIS栅格，但是学习曲线似乎非常陡峭，我可以找到的所有示例都涉及单个栅格图层。在示例中，该图层用作多边形叠加层的一部分，可能平均了每个多边形的栅格值。我还没有找到组合的可复制示例：a）矢量->栅格b）地图代数；c）按照我的第一段的GROUP BY属性。 如果必须执行此任务，我可以使用GDAL或GRASS，但这似乎是PostGIS应该能够处理的。鉴于我的输入数据已经是PostGIS几何图形，这样做会很方便；我真的很想了解PostGIS栅格。 一些样本数据结构： areaid time date isogeom (polygon) 1000 07:15:00 2014-05-05 xxx 1000 07:15:00 2014-05-06 xxy ... 1006 07:15:00 2014-05-05 zzz 我想栅格化，按areaid分组，然后执行地图代数得出： areaid isorast (raster) 1000 aaa 1006 bbb 我没有成功将其包含到PostGIS中。我的方法是将矢量转换为栅格，将栅格转储为数组，并通过psycopg2将其与numpy数组组合，然后再将其写入GeoTIFF（可能会放回PostGIS中）。不理想，但可行。
6
将线shapefile转换为栅格，值=像元中线的总长度
我有一个代表道路网络的shapefile文件。我希望对该数据进行栅格化，栅格中的结果值显示了落入栅格像元内的线的总长度。 数据在英国国家网格投影中，因此单位将是米。 理想情况下，我想使用来执行此操作R，并且我猜测程序包中的rasterize功能raster将在实现此功能中起作用，我只是无法弄清楚应应用的功能。
14 raster  r  line  rasterization 
4
从DEM中提取堤坝？
我确实有1x1米的DEM和.las组成的.las中的原始LiDAR点云。我需要提取河堤（堤坝的最高点）到矢量特征（点，折线）。 对算法或现有工具有任何想法吗？ 堤防在第一个图像上为浅灰色，下面是堤防示例区域的点云图像。在蓝色的河中心线。 从DEM提取线性特征的任何工具吗？ 我相信我需要像ENVI的空间特征提取模块（第7页）之类的东西，但是没有ENVI的部分，因为我负担不起：)
2
使用R？将点数据转换为网格数据框以进行直方图分析。
我对使用GIS数据非常陌生，并且对R的使用不多。我一直在阅读有关如何使用space-analyst.net PDF书籍分析空间数据的知识，所以我并没有完全迷失，但我想我可以描述一下我的问题，人们可能会提出想法。 我有一个数据集，该数据集在不同的纬度/经度坐标下有约2000个测量值，尽管我可能会细分该数据集，因为数据是在3年内收集的，而且条件随时间而变化。我们将被测量的变量称为“ IP”。 我想在样本数据上使用Kriging或其他插值方法在整个问题区域创建IP映射。然后，我想创建一个直方图来测量各种IP桶中的土地数量。我还需要创建一个直方图，以显示每个存储桶中的样本数量（请注意，样本的实际IP可能比克里金预测的土地高或低）。 我遵循如何将数据加载到SpatialPointsDataFrame并运行kriging分析的方法，我遇到的麻烦是如何将数据转换为网格数据框，以便进行直方图分析。 有任何将点转换为网格的建议吗？
14 raster  r  kriging 
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.