r-raster

    0热度

    1回答

    我有一个具有非负值(每个像素都有一个高度值)的栅格矩阵,并且想要使用5x5移动窗口计算一些自定义高度指标,使用对称填充矩阵的边缘。我使用QGIS中的pktools(空间过滤器)做了一些基本的计算(平均值,标准差,最大值),这使我可以做到这一点(使用零填充,复制或对称填充)。 现在我想使用自定义函数(即高度系数,偏度和峰度)来计算一些其他指标,我不能使用pktools。我一直试图使用raster中的

    1热度

    2回答

    我有从日常气候数据创建的栅格堆栈。可以在这里找到: #!/bin/bash wget -nc -c -nd http://northwestknowledge.net/metdata/data/tmmx_1982.nc 目标是从这些日常记录中获得每月95%的温度值。每当我使用raster包中的calc时,它只返回一层而不是12(例如,12个月)我错过了什么? library(raster)

    0热度

    2回答

    我正在使用R(非RStudio)3.3.3在Mac0S10.10.5中的R栅格软件包中处理地理参考图像。使用GUI将工作区保存到xxx.RData文件后,关闭会话重新启动R(无需在启动时读取历史记录文件),然后使用GUI再次加载工作区。当我输入ls()时,所有对象(包括我加载的原始完整全局栅格以及使用crop()创建的子集)都会出现。但是,当我绘制栅格时,只会显示原始非裁剪图像。裁剪的子集不绘制,

    -1热度

    1回答

    我有一个带有2个波段的栅格,我想计算3x3移动窗口中的协方差,并用协方差线代替焦点值。我正在研究栅格包的焦点函数,并且还提供了一些使用getValuesBlock for循环的提示。如果有人指出我解决这个问题,我会非常感激。

    -1热度

    1回答

    我试图覆盖与科罗拉多州(SpatialPolygonsDataFrame)状态相关的蝙蝠位置(SpatialPointsDataFrame)的数据。两个对象的CRS是不同的: crs(colorado) #CRS arguments: # +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 crs(BatRange_

    0热度

    2回答

    我有一个光栅,有50多个波段。我想要搜索栅格第一个波段内的所有像素== 0。以下我想在其他光栅波段中设置所有这些像素NA。因此,我不必再次搜索50次以上。这是我的例子: cl_input <- brick("sometif") for(i in 1:nlayers(cl_input)){ print(names(cl_input[[i]])) cl_input[[i]][cl

    6热度

    2回答

    基本上,我用ASCII的形式计算了一个全局分布概率模型,例如: gdpm。 gdpm的值全部为0和1之间 然后我输入从形状文件的本地地图: shape <- file.choose() map <- readOGR(shape, basename(file_path_sans_ext(shape))) 下一步,我光栅化gdpm,并裁剪使用本地地图: ldpm <- mask(gdpm, ma

    0热度

    1回答

    我已经收到一个.tif文件,其中包含有关每个值的颜色信息(例如,森林是绿色的,云是黑色的).qml文件。 library(raster) library(rasterVis) sm <- raster("UT_classification_coded.tif") ########## Classify raster/build attribute table sm <- rat

    0热度

    1回答

    好吧,我有女佣在科罗拉多州的地图。我喜欢能够遮蔽精选县。我的数据给了我县的名字,但我不知道如何将它应用到我的ggmap。 这是从我的多边形 [email protected][["NAME_2"]] [1] "Adams" "Alamosa" "Arapahoe" "Archuleta" "Baca" [6] "Bent" "Boulder" "Broomfield" "Chaffe

    2热度

    1回答

    library(raster) r <- raster('glc2000_v1_1') # http://forobs.jrc.ec.europa.eu/products/glc2000/products/glc2000_v1_1_Grid.zip extent(r) # class : Extent # xmin : -180.0045 # xmax : 179.9955 #