从R中的范围对象列表中获取Mercator坐标

Question

我有一个带有多个多边形的形状文件，我需要在列表文件中的其他文件或中的每个多边形获得Mercator坐标。在这里，列表与我的ShapeFile的每个多边形的范围

a<-extent(659302.4, 659802.4, 9860325, 9860825)
b<- extent(647663, 648163, 10033987, 10034487)
c<-extent(609015.1, 609515.1, 9871028, 9871528)
ab<-c(a,b,c)

我需要在列表>或中为每个文件中包含坐标的文件生成一个Mercator坐标列表。到目前为止，我唯一的选择是将各个程度栅格化，分配CRS坐标，然后转换为Mercator，如下所示

bb<-raster(b)
crs(bb)<- "+proj=utm +zone=17 +south +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"
mbb<-projectRaster(bb, crs="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0
                   +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs")

MBB类：抛光剂尺寸：12，12，144(nrow，ncol，， NCELL(分辨率：50，50.3(x，y(范围：-8869218， -8868618，34170.81，34774.41(Xmin，Xmax，Ymin，Ymax(坐标。参考。： proj = merc a = 6378137 b = 6378137 lat_ts = 0.0 lon_0 = 0.0 x_0 = 0.0 y_0 = 0 y_0 = 0 k = 1.0 k = 1.0 units = m nadgrids =@null null null null null no_defs

这种方法的唯一问题是非常缓慢，特别是当我需要获得数百个多边形范围的默卡托坐标时。

编辑：我想要的是将投影从UTM更改为Mercator，我正在考虑一个允许为 latne List List 执行此操作的函数而不是如上示例所示。

Answer 1

如果要更改范围的投影(不是栅格(，则可以使用：

library(raster)
library(dplyr)
library(purrr)
a<-extent(659302.4, 659802.4, 9860325, 9860825)
b<- extent(647663, 648163, 10033987, 10034487)
c<-extent(609015.1, 609515.1, 9871028, 9871528)
ab<-c(a,b,c)
ab <- map(ab, function(x) {as(x,'SpatialPolygons')}) %>%
      map(.,function(x){x@proj4string <- CRS("+proj=utm +zone=17 +south +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0");
                        spTransform(x,CRS("+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs"))})
# back to extent
new_extent <-  ab %>% map(extent)