将多边形转换为数据框。

huangapple
go评论91阅读模式
英文:

Convert polygon to dataframe

问题

我通过从大型数据框派生的两个多边形相减创建了一个多边形。更具体地说：

  1. plot(Cyp_occ_4326_df$lon, Cyp_occ_4326_df$lat,
  2.      xlim=c(min(Cyp_csl_4326_df$lon),max(Cyp_csl_4326_df$lon)),
  3.      ylim=c(min(Cyp_csl_4326_df$lat),max(Cyp_csl_4326_df$lat)), 
  4.      type='p', col='red', xlab='x', ylab='y')
  5. points(Cyp_csl_4326_df$lon, Cyp_csl_4326_df$lat, col='blue')

将多边形转换为数据框。

其中Cyp_csl_4326_df是一个数据框（蓝色点），Cyp_occ_4326_df是另一个数据框（红色点），它们在上部与第一个数据框重叠。

我通过以下方式创建了一个多边形：

  1. poly1 <- Polygon(as.matrix(Cyp_csl_4326_df))
  2. poly2 <- Polygon(as.matrix(Cyp_occ_4326_df))
  3. p1 <- SpatialPolygons(list(Polygons(list(poly1), "p1"))
  4. p2 <- SpatialPolygons(list(Polygons(list(poly2), "p2"))
  5. res <- gDifference(p1, p2)

然后我得到了以下图片：

  1. plot(res, col="blue")

将多边形转换为数据框。

问题是，我不知道如何将派生的多边形转换为点以供进一步使用。我得到的只是多边形的范围：

  1. > summary(res)
  2. Object of class SpatialPolygons
  3. Coordinates:
  4.        min      max
  5. x 32.27102 34.58771
  6. y 34.56312 35.69493
  7. Is projected: NA 
  8. proj4string : [NA]

我找不到提取点的方法。

英文:

I have created a polygon by substracting two polygons, derived from large dataframes. More specific:

  1. plot(Cyp_occ_4326_df$lon, Cyp_occ_4326_df$lat,
  2.      xlim=c(min(Cyp_csl_4326_df$lon),max(Cyp_csl_4326_df$lon)),
  3.      ylim=c(min(Cyp_csl_4326_df$lat),max(Cyp_csl_4326_df$lat)), 
  4.      type=&#39;p&#39;, col=&#39;red&#39;, xlab=&#39;x&#39;, ylab=&#39;y&#39;)
  5. points(Cyp_csl_4326_df$lon, Cyp_csl_4326_df$lat, col=&#39;blue&#39;)

将多边形转换为数据框。

where Cyp_csl_4326_df is a dataframe (blue points) and Cyp_occ_4326_df
is another dataframe (red points), which overlays with the first dataframe at the upper part.

I have created a polygon by substracting the two dataframes as follow:

  1. poly1 &lt;- Polygon(as.matrix(Cyp_csl_4326_df))
  2. poly2 &lt;- Polygon(as.matrix(Cyp_occ_4326_df))
  3. p1 &lt;- SpatialPolygons(list(Polygons(list(poly1), &quot;p1&quot;)))
  4. p2 &lt;- SpatialPolygons(list(Polygons(list(poly2), &quot;p2&quot;)))
  5. res &lt;- gDifference(p1, p2)

And I got the following pic:

plot(res, col="blue")

将多边形转换为数据框。

The problem is that I do not know how to convert the derived polygon to points for further use. What I get, is simply the extend of the polygon:

  1. &gt; summary(res)
  2. Object of class SpatialPolygons
  3. Coordinates:
  4.        min      max
  5. x 32.27102 34.58771
  6. y 34.56312 35.69493
  7. Is projected: NA 
  8. proj4string : [NA]

I could not find a way to extract points.

答案1

得分: 1

以下是您提供的代码部分的中文翻译：

  1. # 我的答案基于`sf`包提供的数据以进行复制：
  3. # 读取sf包中的`nc.shp`数据文件
  4. library(sf)
  5. nc <- st_read(system.file("shape/nc.shp", package="sf"))
  7. # 将数据转换为与您的类别相同
  8. library(sp)
  9. res <- as_Spatial(nc$geometry)
  10. summary(res)
  1. # 类别为SpatialPolygons的对象
  2. # 坐标范围：
  3. #        最小       最大
  4. # x -84.32385 -75.45698
  5. # y  33.88199  36.58965
  6. # 投影状态：FALSE
  7. # proj4string : [+proj=longlat +datum=NAD27 +no_defs]

通过这个示例数据，以下代码似乎有效：

  1. # 加载data.table用于rbindlist()函数
  2. library(data.table)
  4. # 从SpatialPolygons中提取点
  5. res_points <- lapply(1:length(res@polygons), function(i) {
  6.   
  7.   lapply(1:length(res@polygons[[i]]@Polygons), function(j){
  8.         
  9.     data.frame(res@polygons[[i]]@Polygons[[j]]@coords)
  10.     
  11.   }) %>% rbindlist()
  12.   
  13. }) %>% rbindlist()
  15. head(res_points)
  1. # 前6行提取的点坐标：
  2. #          X1       X2
  3. # 1: -81.47276 36.23436
  4. # 2: -81.54084 36.27251
  5. # 3: -81.56198 36.27359
  6. # 4: -81.63306 36.34069
  7. # 5: -81.74107 36.39178
  8. # 6: -81.69828 36.47178
英文:

I base my answer on data from the sf package for replication:

  1. library(sf)
  2. nc &lt;- st_read(system.file(&quot;shape/nc.shp&quot;, package=&quot;sf&quot;))
  4. # Convert data to the same class as yours
  5. library(sp)
  6. res &lt;- as_Spatial(nc$geometry)
  7. summary(res)
  1. Object of class SpatialPolygons
  2. Coordinates:
  3.         min       max
  4. x -84.32385 -75.45698
  5. y  33.88199  36.58965
  6. Is projected: FALSE 
  7. proj4string : [+proj=longlat +datum=NAD27 +no_defs]

With this example data the following code seems to work:

  1. # Load data.table for the rbindlist() function
  2. library(data.table)
  4. # Extract points from SpatialPolygons
  5. res_points &lt;- lapply(1:length(res@polygons), function(i) {
  6.   
  7.   lapply(1:length(res@polygons[[i]]@Polygons), function(j){
  8.         
  9.     data.frame(res@polygons[[i]]@Polygons[[j]]@coords)
  10.     
  11.   }) %&gt;% rbindlist()
  12.   
  13. }) %&gt;% rbindlist()
  15. head(res_points)
  1.           X1       X2
  2. 1: -81.47276 36.23436
  3. 2: -81.54084 36.27251
  4. 3: -81.56198 36.27359
  5. 4: -81.63306 36.34069
  6. 5: -81.74107 36.39178
  7. 6: -81.69828 36.47178

huangapple
  • 本文由 发表于 2023年2月24日 03:29:40
  • 转载请务必保留本文链接：https://go.coder-hub.com/75549481.html
  • dataframe
  • polygon
  • r
  • r-sp
  • subtraction
匿名

发表评论

匿名网友

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen:

确定