Georeferenzierung Parzellenplan

geoRefPz.R

# ===
# Funktionen für den Geodatenimport in FastWood
# ===


# lib
# install.packages("dplyr")
# install.packages("tidyr")
# install.packages("shapes")


# Procrustes Analyse für Georeferenzierung von Parzellenplänen --------------------------------

pzGeoRef <- function(origDat, referMat, transMat, procScale) {
  require(dplyr)
  require(tidyr)
  require(shapes)
  
  try(if(!is.matrix(referMat)) stop("Referenz ist keine Matrix"))
  try(if(!is.matrix(transMat)) stop("Zu transformierendes Objekt ist keine Matrix"))
  try(if(ncol(origDat) != 2) stop("Bitte nur Datensatz mit X und Y Koordinaten als Spalten übergeben"))
  
  procAna <- 
    procOPA(
      A = referMat,
      B = transMat,
      scale = procScale,
      reflect = FALSE
    )
  
  origDat <- as.matrix(origDat)
  
  pzGridT <- (((
    origDat * procAna$s -
      matrix(
        colMeans(transMat),
        byrow = TRUE,
        ncol = 2,
        nrow = nrow(origDat)
      )
  )  %*% procAna$R)  +
    matrix(
      colMeans(referMat),
      byrow = TRUE,
      ncol = 2,
      nrow = nrow(origDat)
    ))
  
  pzGeoRef <- vector("list", length = 2)
  
  pzGeoRef[[1]] <- procAna
  pzGeoRef[[2]] <- pzGridT
  
  names(pzGeoRef) <- c("procAna", "pzGridT")
  
    return(pzGeoRef)
}