如何在循环中加载多个点的多个点文件并将其转换为CSV文件？

发布于 2025-02-11 08:42:31 字数 790 浏览 4 评论 0原文

我想阅读并将多个Shapefile转换为CSV。该链接包含Shapefile的样本。以下是我正在尝试的代码。

  library(raster)
  setwd("D:/share_1/PAPER_regimes/22_06_22/hp/firepoint_12_veg")
  all.files <- list.files(pattern="\\.shp$")
  out.files <- gsub("\\.shp$", ".csv")
  crs <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")
  for(i in 1:length(all.files)) {
   d <- readOGR(all.files[i], stringsAsFactors = FALSE)
   sp <- SpatialPointsDataFrame(coords = d[c("Longitude", "Latitude")], d, proj4string = crs) 
   utm <- spTransform(sp, CRS("+proj=utm +zone=16 +datum=WGS84"))
   shapefile(utm, out.files[i])
    }

原文

I want to read and convert multiple shapefile into csv. The link contain sample of shapefile. Following is the code I was trying.

  library(raster)
  setwd("D:/share_1/PAPER_regimes/22_06_22/hp/firepoint_12_veg")
  all.files <- list.files(pattern="\\.shpquot;)
  out.files <- gsub("\\.shpquot;, ".csv")
  crs <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")
  for(i in 1:length(all.files)) {
   d <- readOGR(all.files[i], stringsAsFactors = FALSE)
   sp <- SpatialPointsDataFrame(coords = d[c("Longitude", "Latitude")], d, proj4string = crs) 
   utm <- spTransform(sp, CRS("+proj=utm +zone=16 +datum=WGS84"))
   shapefile(utm, out.files[i])
    }

link contain few shapefile and dataset which are in point

分享到QQ

分享到微博

如果你对这篇内容有疑问，欢迎到本站社区发帖提问参与讨论，获取更多帮助，或者扫码二维码加入 Web 技术交流群。

发布评论

需要登录才能够评论，你可以免费注册一个本站的账号。

×眷恋的温暖 2025-02-18 08:42:31

您可以这样做，

library(terra)
all.files <- list.files(pattern="\\.shp$", full=TRUE)
out.files <- gsub("\\.shp$", ".csv", all.files)

for(i in 1:length(all.files)) {
   d <- vect(all.files[i])
   utm <- project(d, "+proj=utm +zone=16 +datum=WGS84")
   xy <- crds(utm)
   write.csv(xy, out.files[i], row.names=FALSE)
}

我不知道为什么您去了属性中的LON/LAT，但是如果您需要的话，您可以做：

for(i in 1:length(all.files)) {
   d <- vect(all.files[i])
   geo <- vect(data.frame(d), c("LONGITUDE", "LATITUDE"), crs="+proj=longlat")
   utm <- project(d, "+proj=utm +zone=16 +datum=WGS84")
   xy <- crds(utm)
   write.csv(xy, out.files[i], row.names=FALSE)
}

You can do that like this

library(terra)
all.files <- list.files(pattern="\\.shpquot;, full=TRUE)
out.files <- gsub("\\.shpquot;, ".csv", all.files)

for(i in 1:length(all.files)) {
   d <- vect(all.files[i])
   utm <- project(d, "+proj=utm +zone=16 +datum=WGS84")
   xy <- crds(utm)
   write.csv(xy, out.files[i], row.names=FALSE)
}

I do not know why you went to the lon/lat in the attributes, but if that is what you need you can do:

for(i in 1:length(all.files)) {
   d <- vect(all.files[i])
   geo <- vect(data.frame(d), c("LONGITUDE", "LATITUDE"), crs="+proj=longlat")
   utm <- project(d, "+proj=utm +zone=16 +datum=WGS84")
   xy <- crds(utm)
   write.csv(xy, out.files[i], row.names=FALSE)
}

回复收藏 0 原文

~没有更多了~