sp

    0熱度

    1回答

    是否可以計算SpatialLinesDataFrame中每個SpatialLines的面積? 我在尋找一種計算個體面積的方法以及總面積的平均值。

    0熱度

    1回答

    我試圖將鄰居空間多邊形識別爲一組空間多邊形,同時考慮鄰居是否存在或爲NA。我使用rgeos包中的gTouches()函數來確定哪些幾何體具有共同的邊界,但是我不知道如何使其考慮鄰居的值是否爲NA,在這種情況下,我希望它找到下一個最近的幾何體。在下面的示例代碼,我想爲tr2的鄰居,其中包括NA,比tr是不同的: library(rgeos) library(sp) grid <- GridTo

    3熱度

    1回答

    我有一個包含在單個stack中的柵格圖層以及包含點座標的SpatialPoints對象。如果單元格包含一個點,我試圖將柵格圖層的值更改爲NA。 我在下面提供了一個可重現的例子。然而,爲了說明這一點,我的真實數據由10層「棲息地」(海拔,樹冠覆蓋等)組成,這些層都包含在堆棧中。 rasters覆蓋面積約爲34 x 26公里。另外,我有大約10,000個來自動物的GPS位置。 因此,對於可重複的例子。

    3熱度

    2回答

    我有一組數據框並希望使用它們來創建一個柵格網格。 這裏有兩個示例文件: dat1dat2 我用sp包繪製空間點: library(raster) library(sp) coordinates(dat1) <- c("x", "y") proj4string(dat1) <- CRS("+proj=longlat") sp1 <- SpatialPoints(coords = dat1)

    1熱度

    1回答

    我有一個我想擴展的shapefile(基本上添加了一些NA行和列)。有沒有像raster::extend那樣的功能?不幸的是,這隻適用於raster對象。 我可以使用raster::crop對其進行裁剪,該方法有SpatialPolygonsDataFrame,但我無法找到任何用於擴展的方法。

    0熱度

    1回答

    我有來自美國census.gov的人口普查數據的Shapefiles。 我已經使用rdgal庫中的readOGR將文件讀入R中。 SP1:投票區 - SpatialPolygonsDataFrame SP2:對於面積較大,比SP1更高分辨率的塊級 - SpatialPolygonsDataFrame 我曾嘗試使用rgeos庫中的gIntersection,但它不是r琢磨我需要什麼。或者我不知道如何

    0熱度

    1回答

    我遇到了一個小問題,我的for循環運行了3.5小時+。我在SQL Server有一個地理數據集,名爲Shape的字段,它存儲數據類型爲geometry的座標。首先,我通過ODBC [R連接到我的數據庫和檢索我想要的信息(Shape列也轉換爲可讀) library(sp) library(rgeos) library(RODBC) ch<-odbcConnect("SpatialAnalysi

    0熱度

    1回答

    我有一個geocoded數據集,我試圖聚合成多邊形,這樣我就可以將結果繪製成一系列不同級別的choropleth地圖(例如郊區,當地政府區域等)。 要做到這一點,我下面的方法 - shown here - 它使用over功能在sp封裝空間對象加入數據,發現其多邊形我的座標(從單獨的文件)裏鑽然後我使用ggplot2加強了空間對象的繪製。 總的來說,我似乎已經掌握了大部分流程,但是正如您從結果圖表中

    2熱度

    1回答

    比方說,我有一些空間數據,我想用spplot功能從SP繪製它包: library('sp') library('lattice') demo(meuse, ask = FALSE, echo = FALSE) spplot(meuse.grid[,'dist']) 可以很容易地更改主題,例如col.regions,只有一個情節: spplot(meuse.grid[,'dist'], c

    1熱度

    1回答

    我正在做幾個空氣質量站之間的idw插值。對於最終的可視化,我想將我製作的所有圖像合併到一個GIF中。一切正常,唯一的問題是每幅圖像的值的大小都會發生變化,並且規模會隨着值(最小 - 最大)而變化。 如何將顏色的比例設置爲固定值? 下面是一些樣本數據: data(meuse.grid) coords <- SpatialPixels(SpatialPoints(meuse.grid[,c('x'