R/closeMeshes.r

In zarquon42b/mesheR: Meshing Operations on Triangular Meshes

#' close two meshes where one is only an offset of the other
#'
#' close two meshes where one is only an offset of the other
#'
#' @param mesh1 triangular mesh in first state
#' @param mesh2 triangular mesh in second state
#' @param invert logical: if TRUE the orientation of mesh1 will be flipped
#' @importFrom Morpho mergeMeshes invertFaces
#' @importFrom Rvcg vcgGetEdge vcgClean
#' @export
closeMeshes <- function(mesh1,mesh2,invert=TRUE) {
    if (invert)
        mesh1 <- invertFaces(mesh1)
    nvb <- ncol(mesh1$vb)
    merge <- mergeMeshes(mesh1,mesh2)
    gE <- vcgGetEdge(mesh1)
    gE <- gE[gE$border == T,]
    itnew <- rbind(gE$vert1,gE$vert2,gE$vert1+nvb)
    itnew <- cbind(itnew,rbind(gE$vert1+nvb,gE$vert2+nvb,gE$vert2))
    merge$it <- cbind(merge$it,itnew)
    merge <- vcgClean(merge,sel=7)
    merge <- vcgUpdateNormals(merge)
    return(merge)
}