R/contour.R

In espadon: Easy Study of Patient DICOM Data in Oncology

.contour.level <- function(contour.L){
  rg<- (0:length(contour.L))[-1]
  level <- rep(0,length(rg))
  if (length(rg)>1) {
    level <- sapply (rg, function(j){
      ptj<- contour.L[[j]]
      roi.index.k <- rg[-j]
      # if (length(roi.index.z)!=0) {
      r <- unique (sapply (roi.index.k, function (k) {
        ptk <- contour.L[[k]]
        keep <- .pt.in.polygon (ptj[1,1], ptj[1,2], ptk[ ,1], ptk[ ,2], 1e-20) > 0.5
        return (ifelse (all(keep), k,NA))}))
      r <- r[!is.na (r)]
      return(length(r))
    }) 
  }
  return(level)
}