[58c332]: / RadETL / R / modality_Occurrence_Table.R

Download this file

38 lines (36 with data), 1.3 kB

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
#' 'modality'
#'
#' modality function indicates modality of each shoot
#'
#'
#' @param DICOMList you can put it like this and then run the function : DICOMList<-DICOMHeaderList(DICOMFolderPath)
#' @import dplyr
#' @importFrom magrittr "%>%"
#'
#'
#' @return A dataframe indicating modality of each shoot
#' @examples
#' DICOMList<-DICOMHeaderList(DICOMFolderPath)
#' modality(DICOMList)
#' @export
#modality
modality<-function(DICOMList){
modality<-lapply(DICOMList, function(x){
modality<-as.character(x[[1]] %>% filter(name=='Modality') %>% select(value))
if(modality=="CR" | modality=="DX"){
modality='XR'
}
else if(modality=="character(0)" | modality=="" | modality=="integer(0)"){
modality='NA'
}
return(modality)})
modality<-as.data.frame(do.call(rbind, modality))
colnames(modality)<-'modality'
modality<-cbind(modality, radiologyOccurrenceId(DICOMList))
modality<-as.data.frame(modality %>% group_by(radiologyOccurrenceId) %>% distinct(modality))
modality<-split(modality, modality$radiologyOccurrenceId)
modality<-sapply(modality, function(x){
paste0(x$modality, collapse=', ')})
modality<-data.frame(radiologyOccurrenceId=names(modality), modality=modality, row.names = NULL)
return(modality)
}