形状文件：将R中的点（多点）几何体转换为点（点Z）

我有一个多点几何图形中的形状文件（

mult\u point\u example.shp

）：

# Packages
require(sf)

# get AOI
download.file(
  "https://github.com/Leprechault/trash/raw/main/mult_point_example.zip",
  zip_path <- tempfile(fileext = ".zip")
)
unzip(zip_path, exdir = tempdir())

# Open the files
setwd(tempdir())
my_multi_points <- sf::st_read("mult_point_example.shp") 

#Reading layer `mult_point_example' from data source `C:\Users\fores\AppData\Local\Temp\RtmpmQybFP\mult_point_example.shp' using driver `ESRI Shapefile'
#Simple feature collection with 8 features and 10 fields
#Geometry type: MULTIPOINT
#Dimension:     XYZ
#Bounding box:  xmin: -52.73354 ymin: -19.79479 xmax: -52.72586 ymax: -19.79067
#z_range:       zmin: 0 zmax: 0
#Geodetic CRS:  WGS 84

#包
要求（sf）
#获取AOI
下载文件(
"https://github.com/Leprechault/trash/raw/main/mult_point_example.zip",
zip_path使用
st_cast（，“POINT”）
。但是，请注意，当您转换到
POINT
时，功能会重复，即在一个
多点的情况下（一行具有多个功能）包含20个点，当您投射到
点
时，您将获得20行
点
，其特征与原始的
多点
相同：

#包
要求（sf）
#>装载所需包裹：sf
#>链接到GEOS 3.9.0、GDAL 3.2.1、项目7.2.1
#获取AOI
下载文件(
"https://github.com/Leprechault/trash/raw/main/mult_point_example.zip",
zip_路径几何图形类型：多点
#>尺寸：XYZ
#>边界框：xmin:-52.73354 ymin:-19.79479 xmax:-52.72586 ymax:-19.79067
#>z_范围：zmin:0 zmax:0
#>大地测量CRS:WGS 84
#>Talhao Uso特别是数据事件点
#>1142工厂商业EUUR 2012-2017-06-15 Pequeno-52.72675
#>2 142工厂商业EUUR 2012 2017-06-09 Pequeno-52.72843
#>3142工厂商业EUR 2012-2017-06-12 Pequeno-52.73070
#>4142工厂商业区2012-2017-06-09梅奥区-52.72847
#>5142工厂商业区2012-2017-06-12梅奥区-52.73066
#>6142工厂商业区2012-2017-06-09四分之一-52.72868
#>7142工厂商业区2012-2017-06-12 Um_象限-52.73065
#>8142工厂商业区2012-2017-06-12卡雷罗斯-52.73217
#>点Y Fazenda TS_m2几何体
#>1-19.79296 Lobo 1多点Z（-52.72692-1。。。
#>2-19.79197 Lobo 1多点Z（-52.731-19…）。。。。
#>3-19.79308 Lobo 1多点Z（-52.73354-1。。。
#>4-19.79192 Lobo 6多点Z（-52.73086-1。。。
#>5-19.79317 Lobo 6多点Z（-52.73345-1。。。
#>6-19.79181 Lobo 10多点Z（-52.73075-1。。。
#>7-19.79280 Lobo 10多点Z（-52.73343-1。。。
#>8-19.79275 Lobo 1多点Z（-52.73217-1。。。
#解决方案：按点转换
st_cast.sf中的my_points警告（my_multi_points，“POINT”）：重复所有属性
#>它们可能不是常数的子几何体
我的观点
#>具有416个功能和10个字段的简单功能集合
#>几何图形类型：点
#>尺寸：XYZ
#>边界框：xmin:-52.73354 ymin:-19.79479 xmax:-52.72586 ymax:-19.79067
#>大地测量CRS:WGS 84
#>前10个功能：
#>Talhao Uso特别是数据事件点
#>1142工厂商业EUUR 2012-2017-06-15 Pequeno-52.72675
#>1.1 142工厂商业EUUR 2012 2017-06-15 Pequeno-52.72675
#>1.2 142工厂商业EUUR 2012 2017-06-15 Pequeno-52.72675
#>2 142工厂商业EUUR 2012 2017-06-09 Pequeno-52.72843
#>2.1 142工厂商业EUR 2012 2017-06-09 Pequeno-52.72843
#>2.2 142工厂商业EUUR 2012 2017-06-09 Pequeno-52.72843
#>2.3 142工厂商业EUUR 2012 2017-06-09 Pequeno-52.72843
#>2.4 142工厂商业EUR 2012 2017-06-09 Pequeno-52.72843
#>2.5 142工厂商业EUUR 2012 2017-06-09 Pequeno-52.72843
#>2.6 142工厂商业EUUR 2012 2017-06-09 Pequeno-52.72843
#>点Y Fazenda TS_m2几何体
#>1-19.79296 Lobo 1点Z（-52.72692-19.7930。。。
#>1.1-19.79296 Lobo 1点Z（-52.72664-19.7929。。。
#>1.2-19.79296 Lobo 1点Z（-52.72657-19.7928。。。
#>2-19.79197 Lobo 1点Z（-52.731-19.79140）
#>2.1-19.79197 Lobo 1点Z（-52.73098-19.7912。。。
#>2.2-19.79197 Lobo 1点Z（-52.73075-19.7914。。。
#>2.3-19.79197 Lobo 1点Z（-52.73059-19.7913。。。
#>2.4-19.79197 Lobo 1点Z（-52.73052-19.7920。。。
#>2.5-19.79197 Lobo 1点Z（-52.73047-19.7909）
#>2.6-19.79197 Lobo 1点Z（-52.73046-19.7909。。。

由（v2.0.0）于2021年5月26日创建，谢谢@dieghernan！！很好的解决方案和结果，尽管原始的
多点
具有相同的功能。