R/equalProj.R

In move: Visualizing and Analyzing Animal Track Data

setGeneric("equalProj", function(x) {
  standardGeneric("equalProj")
})
setMethod("equalProj",
  signature = c(x = "list"),
  definition = function(x) {
    if (all(is.na(unlist(lapply(x, proj4string))))) {
      return(T)
    } else {
      if (length(x) == 1) {
        return(T)
      } else {
        return(all(unlist(lapply(x[-1], terra::same.crs, x[[1]]))))
      }
    }
    #             if(!isClass(Class="Raster",x)) stop("The list contains objects that are not Rasters.")
    #           tmp <- unlist(strsplit(gsub(pattern="+", replacement="", unlist(strsplit(as.vector(unlist(lapply(x, projection))), " ")), fixed=T), "="))
    #           if(!all(tmp=="NA")){
    #             nom <- as.character(tmp[seq(1,length(tmp),2)])
    #             dat <- data.frame(nom=rep(unique(nom),length(x)), ID=sort(rep(1:length(x), length(unique(nom)))))
    #             vals <- data.frame(val=as.character(tmp[seq(2,length(tmp),2)]), nom=nom)
    #             vals$ID=rep(1:length(x), unlist(lapply(lapply(lapply(as.vector(lapply(x, projection)), strsplit, " "), unlist),length)))
    #             dat <- merge(dat, vals, by=c("nom", "ID"), all=T)
    #
    #             if(any(is.na(dat$val))){
    #               warning(paste("Projections have different number of arguments (differences in: ",paste(as.character(dat$nom[is.na(dat$val)]), collapse=", "),").", sep=""))}
    #             if(!nrow(unique(dat[complete.cases(dat), c("val", "nom")])) == length(unique(nom))){
    #               return(FALSE)} else {return(TRUE)}
    #             } else {return(TRUE)}
  }
)