2024-10-23 16:30:58 +08:00

6 lines
163 B
R

goArea <- function(fn, spx = rgdal::readOGR(fn)){
crs.pcs = crs.Albers(spx=spx)
wbd.pcs = spTransform(spx, crs.pcs)
AA = rgeos::gArea(wbd.pcs)
return(AA)
}