從rgeos包的writeWKT中設置精度

我從R中的specialPolygons對象創建一個WRT字符串。從rgeos包的writeWKT中設置精度

但是，我對輸出中的數字數量感到驚訝。有沒有辦法減少它？

x = sp::SpatialPolygons(Srl = list(sp::Polygons(srl = list(sp::Polygon(coords = cbind(c(-19.8, -19.9, -19.9, -19.8, -19.8),c(148, 148, 148.2, 148.2, 148)))), ID = "1")), pO = 1:1) 
rgeos::writeWKT(x)

這給了我：

"POLYGON ((-19.8000000000000007 148.0000000000000000, -19.8999999999999986 148.0000000000000000, -19.8999999999999986 148.1999999999999886, -19.8000000000000007 148.1999999999999886, -19.8000000000000007 148.0000000000000000))"

來源

2016-11-23 RockScience

基於聊天我h廣告與[@Spacedman]（http://chat.stackoverflow.com/transcript/message/34211117#34211117），這似乎要求C級修補。 –

具體來說，這行https://gitlab.com/geos/libgeos/blob/svn-trunk/src/geom/Coordinate.cpp#L43似乎將精度設置爲「17」。在那裏很好的硬編碼。 – Spacedman

有人發帖回答？我想一個hacky的替代方案是用正則表達式後處理'writeWKT'以丟棄不需要的有效數字...... –

功能似乎:)工作更像你希望的那樣：

library(sf) 
st_as_text(st_as_sfc(x)) 
## [1] "POLYGON((-19.8 148, -19.9 148, -19.9 148.2, -19.8 148.2, -19.8 148))"

來源

2016-11-24 01:15:45

工作的，煩人需要3個PKGS而不是1，但在優秀的新sf package

library(geojson) 
library(wellknown) 
library(jsonlite) 
geoj <- as.geojson(x) 
geojson2wkt(fromJSON(geoj, FALSE)$features[[1]]$geometry, fmt = 1) 

#> [1] "POLYGON ((-19.8 148.0, -19.9 148.0, -19.9 148.2, -19.8 148.2, -19.8 148.0))"

來源

2016-11-23 15:12:00 sckott

從rgeos包的writeWKT中設置精度

回答

相關問題