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 = "