sf - 優文庫

sf

15熱度

3回答

這裏的玩具例子中，我一直在努力與 # Make points point1 <- c(.5, .5) point2 <- c(.6, .6) point3 <- c(3, 3) mpt <- st_multipoint(rbind(point1, point2, point3)) # create multipoint # Make polygons square1 <- rbind

0熱度

2回答

使用geom_sf（）繪製除了latlong以外的任何投影的sf對象

我想繪製一個多邊形geom_sf（）在任何投影其他latlong。我使用的manual pages for geom_sf() 發現輸入數據集的例子： nc <- sf::st_read(system.file("shape/nc.shp", package = "sf"), quiet = TRUE) 從經緯度到EPSG轉化：3857 nc_3857 <- sf::st_transform

0熱度

1回答

地圖3.2.0包（r）：如何在同一個圖中放置數據集的多個區域？

我需要意大利的lombardy region的地圖。當我運行map("italy")但是我可以看到數據集是在省級而不是區域級別上劃分的。倫巴第由12個省組成。我可以通過做map("italy","Milan")得到米蘭（市和省）的地圖，但是當我做map("italy","Lombardy")時，它給了我一個錯誤。我的問題是：我怎麼能在一個地圖多個省（在倫巴第地區）？是否可以刪除內部邊框？

6熱度

1回答

將ggplot2 geoms添加到簡單特徵圖

sf包似乎是一種更方便用戶使用空間數據的方法，比如說sp。 library(sf) devtools::install_github("tidyverse/ggplot2") library(ggplot2) # generate some data set.seed(123) y = rnorm(10, mean=40, sd=20) x = rnorm(10, mean=-10

0熱度

1回答

在柵格中繪製具有覆蓋多邊形的柵格時確定問題的範圍

我正在研究基於rasterVis::levelplot的繪圖函數，用戶可以僅通過柵格對象或柵格對象以及多邊形對象。功能是相當複雜的，但顯示該問題的最小子集讀作： library(sf) library(raster) library(rasterVis) myplot <- function(in_rast, in_poly = NULL) { rastplot <- rast

3熱度

1回答

多邊形

我有一個data.frame與lat S和lng s表示定義的矩形框的界限，像這樣 geohash north_lat south_lat east_lng west_lng 1 gbsuv 48.69141 48.64746 -4.306641 -4.350586 2 gbsuy 48.69141 48.64746 -4.262695 -4.306641 什麼是這個轉換爲持有列的sf對象

2熱度

1回答

`poly.counts`的等價物用來計算落在多邊形內的sf包中的緯度/長度對

sf包提供了一種處理地理要素的好方法，但是我找不到與poly.counts功能從GISTools包需要sp對象。 poly.counts從SpatialPointsDataFrame秋天SpatialPolygonsDataFrame的多邊形內計算的點的數量，並且可以如下使用：數據每個多邊形我生成兩個 ## Libraries library("GISTools") library("ti

0熱度

1回答

將多個`.shp`文件作爲一個對象讀取？

如何將多個.shp文件作爲一個對象讀取？我想簡單地閱讀代碼。 nc <- st_read(dsn = "nc", layer = c("nc1","nc2")) 將多個文件作爲對象讀取的最佳方法是什麼？ library(sf) nc <- st_read(system.file("shape/nc.shp", package="sf")) nc1 <- nc[1:50,

4熱度

1回答

無法獲得alpha與sg linestring對象上的ggplot2 :: geom_sf配合使用

這需要使用R中的庫sf和ggplot2。我有一個帶有4個線串的sf對象。這些3是相同的，一個被擴展： a <- st_linestring(rbind(c(2,2), c(3,3), c(3,2))) b <- st_linestring(rbind(c(2,2), c(3,3))) c <- st_linestring(rbind(c(2,2), c(3,3))) d <- st_line

3熱度

1回答

如何在地圖上標記個別狀態，而其他人在細分水平

我已經設法生成地圖，但是我需要爲包含分區（級別3）的州（級別2）添加一個標籤，，而不是標記每個細分（僅限於此狀態）。在數據「newpak」中，行641-664對應於這個狀態，有沒有辦法在這個狀態之上只放一個名字。 library(dplyr) library(raster) library(sf) library(tidyverse) library(ggrepel) devtools