Questions tagged «r»

如何从堆叠栅格中找到最大像元值。 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 

我已经计算出了物种分布的表面积（合并了shapefile中的多边形），但是由于该区域可以由距离相当远的多边形组成，因此我想计算一些分散度。到目前为止，我所做的就是检索每个多边形的质心，计算它们之间的距离，并使用它们来计算变异系数，如下面的虚拟示例中所示； require(sp) require(ggplot2) require(mapdata) require(gridExtra) require(scales) require(rgeos) require(spatstat) # Create the coordinates for 3 squares ls.coords <- list() ls.coords <- list() ls.coords[[1]] <- c(15.7, 42.3, # a list of coordinates 16.7, 42.3, 16.7, 41.6, 15.7, 41.6, 15.7, 42.3) ls.coords[[2]] <- ls.coords[[1]]+0.5 # use simple offset ls.coords[[3]] <- c(13.8, 45.4, # …

9 polygon  r  distance  centroids 

每个人，都为造成的困扰感到抱歉，但是我刚接触R时遇到了一个非常困难的问题：我想创建一张动画的Russin地图，例如失业率随不同年份的变化。在图像上，您可以看到一年的数据 require(sp) require(maptools) require(RColorBrewer) require(rgdal) rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData") print(load(rus)) unempl <- read.delim2(file="C:\\unempl1.txt", header = TRUE, sep = ";",quote = "", dec=",", stringsAsFactors=F) gadm_names <-gadm$NAME_1 total <- length(gadm_names) pb <- txtProgressBar(min = 0, max = total, style = 3) order <- vector() for (i in 1:total){ order[i] <- agrep(gadm_names[i], unempl$region, max.distance = 0.2)[1] …

9 cartography  r  animation 

我正在尝试从Poly2删除到Poly1的相交。以下是2个多边形的坐标。 > coords1 [,1] [,2] [1,] 992.0161 7462.531 [2,] 950.1962 7413.532 [3,] 902.3632 7421.207 [4,] 706.1985 7472.378 [5,] 654.0139 7711.760 [6,] 657.5960 7726.950 [7,] 786.2667 7844.380 [8,] 789.5935 7842.925 [9,] 981.7046 7498.659 [10,] 983.4246 7493.271 [11,] 990.6680 7469.888 [12,] 992.0161 7462.531 > coords2 [,1] [,2] [1,] 930.3464 7607.313 [2,] 979.3528 …

9 r  polygon  intersection  symmetricdifference 

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

9 arcgis-desktop  python  arcpy  raster  r 

随附的图像显示了一个森林间隙，赤松表示为圆圈，白松表示为十字架。我有兴趣确定两种松树之间是否存在正关联或负关联（即它们是否在同一地区生长）。我知道R spatstat套件中的Kcross和Kmulti。但是，由于我需要分析50个空白并且比R更熟悉python编程，因此我想找到一种使用ArcGIS和python的迭代方法。我也愿意接受R解决方案。 如何实现二元Ripley的K函数？

9 arcgis-desktop  python  arcpy  r  spatial-statistics 

在的文档中，sp我看到coords该SpatialPolygons函数的参数仅接受2列。无法创建3维多边形吗？我正在通过使用package gstat通过将时间添加为a的第三个空间维度来进行时空克里金法SpatialPointsDataFrame。我想将其扩展为以多边形为块的克里金法。有什么办法吗？

9 r  interpolation  spatial-statistics  kriging  gstat 

我需要一个空间网格作为各种专题地图的主网格。如何从丢弃所有NA像素的栅格中生成空间网格？

9 r 

我的工作区中有100个SpatialPointsDataFrames。我以这种方式加载它们： filenames <- list.files(path="", pattern="XYhectareTravelTimes_ez+.*shp") for(i in filenames){ filepath <- file.path("/",i) assign(i, readShapePoints(filepath)) 他们这样称呼： XYhectareTravelTimes_ez10.*shp XYhectareTravelTimes_ez11.*shp XYhectareTravelTimes_ez12.*shp 等等 如何将它们转换为在工作区中循环的栅格？ 我是R用户的新手，希望能找到一些帮助。非常感谢。

9 r  rasterization 

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

9 raster  r  georeferencing 

可能是一个非常简单的问题。我列出了大约一千个潜在地理位置（经纬度），其中我需要选择“分布最广泛”的200个位置。我认为那是总平均距离最高的200分。想想城市中的商店。 有定义的方法可以做到吗？也许在R包中？ 多亏了每个人，使它成为学习它的好地方！ /克里斯

9 r  spatial-statistics  business 

这是stackoverflow的一个交叉点。我对GIS软件了解甚少，如果这是一个太基本的问题，我会预先在R. Apologies中进行所有映射。假设我有两个来自不同来源但具有不同属性的shapefile。可以说，一个是德克萨斯州的行政边界（boundaries.shp），另一个是德克萨斯州河流的（rivers.shp）。我还有第三个文件towns.csv显示州内城镇的位置。读完文件后，我可以在maptools包装的相似边界上绘制城镇位置图： plot(boundaries); points(towns$lon, towns$lat) 但是，如何重叠所有三个呢？当然有一个简单的方法可以做到吗？

9 r  overlay 

我有一个CSV格式的数据集，例如： latitude, longitude, value -45, 45, 10 .... 我想使用ArcGIS或R + grass从这些数据计算县级均值。 我具有ArcGIS USA Counties图层，但是无法成功完成空间连接。 我怎样才能做到这一点？

9 arcgis-desktop  data  geoprocessing  r