Questions tagged «raster»

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



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

5
有哪些可用的栅格平滑/一般化工具?
我有一个DEM,我想对其进行平滑或概括以消除地形的极端情况(切掉峰和填谷)。理想情况下,我还希望控制“模糊性”的半径或水平。最后,我将需要一组从稍微模糊到真正模糊的栅格。(理论上,最模糊的是所有值的算术平均值的恒定栅格)。 我可以使用任何工具或方法(基于Esri,GDAL,GRASS)吗?我需要在家烘烤自己的高斯模糊例程吗?我可以使用低通滤波器(例如ArcGIS的filter),如果是这样,我是否需要运行多次以获取大半径的效果?


2
如何使用R将lat,lon,value形式的数据转换为栅格文件?
我在美国大陆的km网格上有一组数据值,这些列是“纬度”,“经度”和“观测”,例如: "lat" "lon" "yield" 25.567 -120.347 3.6 25.832 -120.400 2.6 26.097 -120.454 3.4 26.363 -120.508 3.1 26.630 -120.562 4.4 或者,作为R数据帧: mydata <- structure(list(lat = c(25.567, 25.832, 26.097, 26.363, 26.63), lon = c(-120.347, -120.4, -120.454, -120.508, -120.562), yield = c(3.6, 2.6, 3.4, 3.1, 4.4)), .Names = c("lat", "lon", "yield"), class = …
40 raster  convert  r 

3
如何代表一段时间内的趋势?
我有一个栅格时间序列,其中每个栅格的像元在特定时间表示某个值。 我想生成一个图表,突出显示每个单元随时间变化的正值或负值趋势的大小。 我有些天真的方法是将简单的线性回归(X =时间和Y =值)拟合到每个像元,并将斜率数组输出到栅格(如下图所示)。只能通过导出重要结果来进行过滤。 在栅格时间序列中,我还能如何表示一段时间内的趋势? 我对通用技术感兴趣,而不对软件专用说明感兴趣。

8
创建shapefile显示栅格的轮廓图?
我有大约1,000张tiff格式的卫星图像,我想创建一个shapefile,将其作为栅格的索引。这类似于栅格目录,但是我不想建立栅格目录。 我可以预见的一些障碍是图像是经过地理反射的,因此它们的形状不是矩形的(我说的是数据区域)。 为了澄清起见,我要求多边形仅覆盖栅格的非零(或非nodata)像素,而不覆盖整个矩形栅格。到目前为止,大多数答案都给出了一个覆盖数据以及非数据像素的矩形多边形。 我的形象 我检查过的工具给出的结果(例如栅格目录,各种Arcscript,在答案之一中给出的自定义Python脚本): 我想要的结果:


2
在R中裁剪栅格
我正在为美国东北部绘制地图。地图背景必须是海拔图或年平均温度图。我有来自Worldclim.org的两个栅格,这些栅格为我提供了这些变量,但是我需要将它们裁剪到我感兴趣的州的范围。有关如何执行此操作的任何建议。这是我到目前为止所拥有的: #load libraries library (sp) library (rgdal) library (raster) library (maps) library (mapproj) #load data state<- data (stateMapEnv) elevation<-raster("alt.bil") meantemp<-raster ("bio_1.asc") #build the raw map nestates<- c("maine", "vermont", "massachusetts", "new hampshire" ,"connecticut", "rhode island","new york","pennsylvania", "new jersey", "maryland", "delaware", "virginia", "west virginia") map(database="state", regions = nestates, interior=T, lwd=2) map.axes() #add …
33 raster  r  clip  maps 

7
在QGIS作曲家图例中使用连续色带作为栅格图例?
我正在寻找一种为栅格图层插入连续颜色图例的方法。我设法用条形图创建带有值的离散符​​号,但我想在条形图上有一个连续的条形图,类似于旧插件“ One-Band Raster”,它不适用于QGIS版本高于1.7 我尝试设置“属性”->“样式”,在其中选择“渲染类型”“单带伪色”,然后看到“生成新色带”框。我尝试了“连续”或“等间隔”。但是,在Composer中,图例始终以离散符号显示。 有什么方法可以使QGIS Composer为栅格图层显示一个不错的连续波段?

5
将栅格值提取到面属性中
我有一个像素为高程数据的州的栅格数据。我也有Vector图层,其中包含300个州的州。我需要的是,我想使用ArcGIS从栅格数据中获取特定分区的平均高程值。 我尝试在ArcGIS中通过遮罩提取来提取单个多边形以将栅格裁剪为300个栅格,然后从单个栅格中收集了平均高程,但是花了几天的时间才能完成工作。还有没有其他简单的方法来获得它。


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 

5
通过R中的许多多边形提高裁剪,遮罩和提取栅格的速度?
我正在从基于数千个多边形边界的栅格中提取不同土地利用类型的面积和覆​​盖率。我发现,如果我遍历每个单独的多边形并进行裁剪,然后将栅格蒙版减小到特定多边形的大小,则提取功能的运行速度会更快。尽管如此,它的运行速度很慢,我想知道是否有人对提高我的代码效率和速度有任何建议。 我发现与此相关的唯一的事情就是这个响应由Roger Bivand使用谁建议GDAL.open(),并GDAL.close()作为以及getRasterTable()和getRasterData()。我研究了这些内容,但过去在gdal上遇到过麻烦,并且对它的了解不够深,不知道如何实现它。 可重现的示例: library(maptools) ## For wrld_simpl library(raster) ## Example SpatialPolygonsDataFrame data(wrld_simpl) #polygon of world countries bound <- wrld_simpl[1:25,] #name it this to subset to 25 countries and because my loop is set up with that variable ## Example RasterLayer c <- raster(nrow=2e3, ncol=2e3, crs=proj4string(wrld_simpl), xmn=-180, xmx=180, ymn=-90, ymx=90) …

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.