R-CMD检查 CRAN_状态_标签 R构建状态

空间离散余弦变换

sfdct的目标是提供简单特征的约束三角剖分。

限制

根据顶点和边执行三角剖分性能这意味着每个输出功能将由与所有输入边对齐的三角形组成。这也将对应于任何线或多边形的外壳。在GEOMETRYCOLLECTION中,同一路线仅适用于集合中单独的几何体,因此直线、多边形和点集都会进行三角剖分,就像它们是独立的一样。

未来的版本将对GEOMETRYCOLLECTION进行三角剖分,因为它是一组边和顶点。这也将允许将整个数据集三角化为一个数据集。

我还不清楚如何最好地维护“整个数据集”的原始特性标识,或者这是否重要。如果你有兴趣,请联系!

用于处理分组的简单复杂数据结构的更通用的结构正在开发中,但与此包中的简单功能相结合,为这些细微差别以及我们可以在多大程度上推动标准工具提供了有用的说明。

例子

这是一个基本示例,向您展示了如何分解MULTIPOLYGON平方英尺数据帧对象转换为GEOMETRYCOLLECTION平方英尺由三角形构成的数据帧对象:

图书馆(平方英尺)
#>链接到GEOS 3.10.2、GDAL 3.4.3、PROJ 8.2.0;sf_use_s2()为真
图书馆(sfdct)
数控<- 读取sf(system.file(系统文件)(“形状/nc.shp”,包装=“sf”),安静= 真的)
(nc_三角形<- ct_三角形(数控[1:5,c(c)(“名称”)]))
#>具有5个功能和1个字段的简单功能集合
#>几何图形类型:GEOMETRYCOLLECTION
#>尺寸:XY
#>边界框:xmin:-81.74107 ymin:36.07282 xmax:-75.77316 ymax:36.58965
#>大地测量CRS:NAD27
#>#A台:5×2
#>名称几何图形
#>*<chr><几何颜色[°]>
#>1灰烬几何颜色(多边形((-81.54084 36.27251,-81.47276 36.23…
#>2 Alleghany几何选择(多边形(-81.24069 36.37942,-81.23989 36.36…
#>3 Surry几何配色(多边形((-80.87086 36.32462,-80.87438 36.23…
#>4当前几何颜色(多边形((-76.09106 36.50357,-76.15815 36.41…
#>5北安普敦几何颜色(POLYGON((-77.53808 36.30246,-77.58008 36.32…

(阿苏布<- st_几何(nc_三角形)[[4]] )
#>几何配色(POLYGON((-76.09106 36.50357,-76.15815 36.41269,-76.09509 36.34892,-76.09105 36.50355)),POLYOGN,POLYGON((-76.00161 36.41891,-76.04395 36.35359,-76.01735 36.33773,-76.00151 36.4189一),POLYOGON((-75.95126 36.36547,-76.000161 36.4189,-76.01745 36.3377,-75.95126 36.36547),POLLYGON,POLYGON((-76.09106 36.50357,-76.09509 36.34892,-76.00161 36.41891,-76.0910636.50357]),POLYKON((.76.16829 36.42709,-76.09.106 36.50377,-76.1274 36.55716,-76.16829 3642709),POLLYGON,多边形((-76.1274 36.55716,-76.33025 36.55606,-76.16829 36.42709,-76.1274 36.55716)),多边形((-76.09106 36.50357,-76.04596 36.55695,-76.1274 36.55716,-76.09106 36.50357)),多边形((-76.02717 36.55672,-75.97629 36.51793,-75.99866 36.55665,-76.02717 36.55672)),POLYGON((-76.00161 36.41891,-75.97607 36.43621,-76.09106 36.50357,-76.00141 36.4189一),POLYOGON((-76.3321 36.51437,-76.4596 36.55695,-76.09016 36.50357,-76.03321 36.512437),POLLYGON,(-75.95126 36.36547,-76.01735 36.33773,-76.00897 36.3196,-75.951226 36.36547),POLYGON((-75.95751 36.25945,-75.91376 36.2448,-75.94193 36.29434,-75.95752 36.25944),POLYOGN((-75.94193 36.29.434,-75.95126 36.36547,-76.00897 36.3196,-75.94.193 36.2943),POLLYGON,POLYGON((-75.80006 36.11282,-75.91376 36.2448,-75.85516 36.10567,-75.80007 36.11282),POLYOGON((-75.95126 36.36547,-75.94193 36.29434,-75.92459 36.35095,-75.95126 36.36547),POLLYGON,POLYGON((-75.87817 36.55587,-75.78317 36.22519,-75.77316 36.22926,-75.87816 36.55588),POLYKON((7-5.91376 36.2448,-75.80006 36.11282,-75.92459 36.35095,-75.91376 362448),POLLYGON,POLYGON((-75.97607 36.43621,-76.00161 36.41891,-75.96976 36.41512,-75.9760736.4362一),POLYOGON((-75.78317 36.22519,-75.87817 36.55587,-75.90199 36.5562,-75.7831736.22519),POLLYGON,多边形((-75.95126 36.36547,-75.92459 36.35095,-75.92812 36.42324,-75.9512 36.36547)),多边形((7-5.91192 36.54253,-75.99866 36.55665,-75.97629 36.51793,-7591192 36.54 253))

更密集的三角形,以及可选的D类确保满足Delaunay标准。

st_几何(ct_三角形(数控[4, ],一个= 0.0007,D类= 真的))
#>1个特征的几何图形集
#>几何图形类型:GEOMETRYCOLLECTION
#>尺寸:XY
#>边界框:xmin:-76.33025 ymin:36.07282 xmax:-75.77316 ymax:36.55716
#>大地测量CRS:NAD27
#>几何配色(多边形((-76.02771 36.4037。。。

更多示例请参见小插曲。

如果您对该领域的开发感兴趣,请参阅anglr包将空间数据三角剖分为网格的主要工具。开发项目落叶松科旨在使用CGAL为复杂形状的有限元分解提供更强大的工具。有一个夹耳版本,来自决策包在里面硅酸盐


请注意,此项目发布时带有贡献者行为准则。参与此项目即表示您同意遵守其条款。