Questions tagged «raster»

是否有一个词可用于描述两个或多个具有相同分辨率和范围的栅格网格？我到处都是“巧合”和“共同延伸”之类的词，但是没有什么感觉很正确。而且，我宁愿找到比“相同范围和分辨率的栅格”更紧凑的东西。

9 raster  terminology 

我有几种tiff格式的栅格，这些栅格具有分配给它们的颜色图： 然后，我创建一个规格符合影像的Mosaic数据集，并向其中添加所有tiff（您可以忽略警告，这只是为了说明我的所作所为）： Colormap不会显示，也不允许我从tiff导入我的Colormap： 有什么方法可以添加/导入我的颜色图，以便我们可以根据这些值查看数据？

9 arcgis-desktop  raster  arcmap  mosaic-dataset  colormap 

我正在使用levelplot超赞rasterVis包的功能来创建一个地图，其值在零附近变化。我正在使用红色到蓝色调色板绘制栅格文件，但是我试图在地图中将灰色分配为零值。 具体来说，我想复制此图的颜色： 请注意，这里使用的是红色到蓝色的比例，但零值已用灰色上色。 目前，我的地图如下所示： 这是再现它的代码（文件可在https://www.dropbox.com/s/cypfdu1eaz2fuok/r.annual.tif?dl=0上获得）： # Load required packages library(rasterVis) # open file r.annual <- raster("Downloads/r.annual.tif") # Set color palette myTheme=rasterTheme(region=brewer.pal('RdBu', n=11)) # Plot levelplot(annual.mask, par.settings=myTheme, margin=F) 如何在上面的地图中为所有零值分配灰度？

9 raster  r  legend  visualisation  color-ramp 

我正在将向量转换为R中的栅格。但是，过程太长了。是否可以将脚本放入多线程或GPU处理中以使其更快地执行？ 我的脚本栅格化矢量。 r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do pixel setor.r = rasterize(setor, r.raster, 'dens_imov') 光栅 类：RasterLayer尺寸：9636、11476、110582736（nrow，ncol，ncell）分辨率：10、10（x，y）范围：505755、620515、8555542、8651792（xmin，xmax，ymin，ymax）坐标。参考 ：+ proj = longlat + datum = WGS84 + ellps = WGS84 + towgs84 = 0,0,0 设定者 类：SpatialPolygonsDataFrame功能：5419范围：505755，620515.4，8555429，8651792（xmin，xmax，ymin，ymax）坐标。参考 ：+ proj = utm …

9 raster  r  rasterization  rgdal  parallel 

使用QGIS 2.4.0 Chugiak我有一个关于裁剪光栅图像的问题。我有一个航拍图片，我想使用shapefile的范围进行裁剪。到目前为止，我已经使用光栅/提取/剪切器工具来做到这一点。问题是我的输出栅格的文件大小几乎是原始文件的三倍（即分别为725.283和249.693 kb），而其覆盖面较小。 为什么我的文件更大？有什么办法可以“设置”输出文件的大小？ Clipper的帮助功能参考以下页面：http : //www.gdal.org/gdal_translate.html尽管我不能从中做很多。 我对QGIS比较陌生。

9 qgis  raster  clip  file-size 

我正在尝试投影栅格。在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 

我要将.adf文件加载到R中。数据来自以下页面：http : //www.fao.org/geonetwork/srv/en/metadata.show?id= 14057 我尝试了以下代码，这些代码是在互联网上进行一些研究后发现的。问题是，在RasterLayer类中，我得到了不应该存在的负值。我不知道为什么会这样，所以希望有人可以帮助我！ 码： library(rgdal) library(RColorBrewer) dpath<- path... x <- new("GDALReadOnlyDataset", dpath) getDriver(x) getDriverLongName(getDriver(x)) xx<-asSGDF_GROD(x) r <- raster(xx) “ r”的输出为： r类：RasterLayer尺寸：2160、4320、9331200（nrow，ncol，ncell）分辨率：0.08333333、0.08333333（x，y）范围：-180、180，-90、90（xmin，xmax，ymin，ymax）坐标。参考 ：+ proj = longlat + ellps = WGS84 + towgs84 = 0,0,0,0,0,0,0 + no_defs数据源：在内存名称中：band1值：-997，16（min，max） 值中的“ 16”表示生长期长度的16类。但是我不知道那些“ -997”来自哪里。坐标可能有问题。参考？ 这也是“ xx”的数据摘要： 数据摘要：最小值 第一区 中位数第三区 最高 不适用-997 3 5 -9 …

9 raster  data  r  esri-grid-format 

我想从指示的文件夹中导入所有多个镶嵌栅格（每个镶嵌由7层组成）到我的R中。然后将它们作为单个多层栅格进行访问。 所以我做了： # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(), pattern="mos.*.img$") # read files as rasters s <- stack(raster_data) # check my imported rasters p.ex. raster n°8 from "s" raster stack s[[8]] 而我的栅格s [[8]]仅包含1层，因此未导入整个镶嵌！ nlayers(s[[8]]) [[1]] 如果我分别阅读每个镶嵌图，它将起作用： # read 1 mosaic (composed by 7 bands) mosaic1<-brick("mosaic1.img") # extract one …

9 raster  r  remote-sensing  landsat 

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

9 python  raster  shapefile  gdal  ogr 

如何从堆叠栅格中找到最大像元值。 Rmax <- maxValue(RAD1998.all[[1]]) 工作正常，但 Rmax <- maxValue(RAD1998.all[[2]]) 给出NA。 当然不是在堆叠的栅格中。 这是我的代码： RAD1998 <- raster(paste(getwd(), "/1998bil/1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil", sep = "")) list.ras <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names = T, pattern = ".asc.bil")) RAD1998.all <- stack(list.ras)

9 raster  r 

ArcMap中是否有任何工具或方法可以测量内核密度输出的每个级别的面积？ 注意，这只是一个内核密度图（不是形状文件或多边形）。内核密度分析（深绿色（10％）到红色（90％）（1） 我已经尝试了一些方法，但是我需要准确。

9 arcgis-desktop  raster  arcmap  kernel-density 

我像往常一样编译并安装了当前的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 

我在同一区域有两个重叠的.ecw层，尽管其中一个被剪切在中间，并带有白色边缘。当我将背景值设置为255,255,255并将其定义为空心时，它只会空心化该特定的RGB值，而实际上，我的图层中还有其他“白色”值，例如： 254,254,254；253,253,253 ..一直到249,249,249。 有没有一种方法可以定义显示背景值，以包含多个唯一的RGB值集？ 我知道我可以解决此问题（剪切，栅格计算器等），但是我只想保留一个原始的.ecw文件而不产生重复，因此我只想使用Arcmap符号系统或图像渲染工具。 我需要的Custom transparency options与QGIS中的完全一样，但适用于ArcMap。

9 raster  arcmap  arcgis-10.2  ecw  display 

我正在尝试为自然灾害模型改进当前极其繁琐的vector / python处理。目前，我们有一个冗长的脚本，该脚本会从给定点生成距离/方位线，以确定： 与其相交的多边形的类型（例如，森林，草丛，沼泽等） 到该多边形的距离 这些线中有多少条与多边形相交，以确定其“环绕”程度。 涉及的更多，但这就是要点。我正在尝试找到一种方法来改进此问题，目前已在第3部分中介绍。该想法是确定某个点是否完全被多边形包围，例如在200m之内 因此，在我的附件图像中，我希望将点A标记为比点B具有更高的风险，因为它完全被多边形包围了。可以重复执行1300万个点，因此这不是一个小任务，我宁愿有一个表面来从中获取值，而不是运行脚本。我认为必须采用各种水文工具或成本路径，但我似乎无法解决。 我该怎么办？

9 arcgis-10.0  arcgis-desktop  raster  arcpy 

我有一个高程的光栅图像。我添加了包含本地政治分歧的矢量shapefile。我要保留标高颜色，但我还要保留代表政治分歧的黑线。 我降低了政治分歧形状文件的透明度，但是我不知道如何使它不影响政治分歧线。是否建议在保持黑色边框颜色充分的同时降低填充透明度？

9 qgis  raster  vector  elevation  transparency