sp

    1热度

    1回答

    我有一个数据集,包括许多不同物种的个体的GPS收集地点。具体来说,我的数据集有3列:经度,纬度和物种。我想删除每个物种内重复和非常相似的收集地点,同时为此目的独立处理每个物种。在R组件SP1工作,我能够做到这一点使用下面的代码(而忽略了种): data = read.csv(file="Population_locality.csv", header=T) coordinates(dat

    2热度

    2回答

    非常简单的情况:多边形定义了一个地理区域,我想知道它由gps坐标给定的点是否位于该多边形内。 我经历了许多SO问题,并尝试过各种函数和包如sp,但无法弄清楚它为什么失败。 我这个很简单的功能的尝试: https://www.rdocumentation.org/packages/SDMTools/versions/1.1-221/topics/pnt.in.poly install.package

    0热度

    1回答

    我正在研究基于rasterVis::levelplot的绘图函数,用户可以仅通过栅格对象或栅格对象以及多边形对象。 功能是相当复杂的,但显示该问题的最小子集读作: library(sf) library(raster) library(rasterVis) myplot <- function(in_rast, in_poly = NULL) { rastplot <- rast

    2热度

    1回答

    sf包提供了一种处理地理要素的好方法,但是我找不到与poly.counts功能从GISTools包需要sp对象。 poly.counts从SpatialPointsDataFrame秋天SpatialPolygonsDataFrame的多边形内计算的点的数量,并且可以如下使用: 数据每个多边形 我生成两个 ## Libraries library("GISTools") library("ti

    1热度

    1回答

    我需要在CARTO(aka cartodb)中导出等值线图,所以我试图将这个stat2density图表保存为shapefile或geojson等地理数据文件格式。 我可以用ggsave将它保存在SVG中,但将它转换为spdf或sf oblejct会非常有帮助。 library(ggplot2) library(ggmap) data("crime") crime<- head(crime,

    0热度

    1回答

    我试图从包含124段的列表构建一个SpatialLines对象。 列表中的每个段包含每个节点的x,y坐标(见下文)。我使用下面的代码来创建SpatialLines对象,但它只是检索一个段。有什么建议么? test.func = function(x){ for (i in 1:length(x)) { tt[[i]] <- x[i]; tt[[i]] = Line(tt[[i]]);

    0热度

    1回答

    短版:执行以下命令qtm(World, "amount")时,我得到了以下错误消息: Error in $<-.data.frame (*tmp* , "SHAPE_AREAS", value = c(653989.801201595, : replacement has 177 rows, data has 175 免责声明:这是我过去在this question同样的问题,但如果我没有错,其中

    1热度

    2回答

    我想在我设法绘制的网格数据顶部绘制饼图。 数据: nasafile <- "http://eosweb.larc.nasa.gov/sse/global/text/global_radiation" nasa <- read.table(file=nasafile, skip=13, header=TRUE) 下面这个帖子:R plot grid value on maps,我用spplot

    1热度

    1回答

    我有一组点,其指向: > head(Map3) SpatialPoints: lat long Cincinnati OH 39.10312 -84.51202 Dayton OH 39.75895 -84.19161 Lexington KY 38.04058 -84.50372 Zanesville OH 39.94035 -82.01319 Louisvil

    1热度

    1回答

    我想创建一个空的SpatialLines对象。对于多边形很容易: SpatialPolygons(list()) 对于空间线这不起作用: SpatialLines(LinesList = list()) Error in bb[1, ] : incorrect number of dimensions SpatialLines(LinesList = Lines(list(),ID = "