Rで解析：画像のカラーコードを取得する方法！「EBImage」パッケージの紹介

source("http://bioconductor.org/biocLite.R")
biocLite("EBImage")

library("EBImage")
#Imageファイルの読み込み
selectAImage <- paste(as.character(tkgetOpenFile(title = "Imageを選択",filetypes = '{"Imageファイル" {".*"}}',initialfile = "*.*")), sep = "", collapse =" ")
Image <- readImage(selectAImage, mode = "x11")

#読み込んだImageをブラウザで表示
display(Image)

#Imageからカラーコードを取得
#ピクセル数が多いImageほどメモリが必要です
ColCode <- channel(Image, mode = "x11")

#2000×1260ピクセルの処理で私の環境で約1分ほどかかりました。
CountColor <- sort(table(ColCode), decreasing = TRUE)

#特定の色を削除。今回は白色を削除
CountColor <- CountColor[names(CountColor) != "#FFFFFF"]

#検知数TOP5を表示
SelectCol <- 5
head(CountColor, SelectCol)

ColCode
#FAE7D8 #DBD5D4 #CFC8C8 #333335 #5C6062
3422     954     810     779     712 

pie(rep(1, SelectCol),
    labels = names(head(CountColor, SelectCol)),
    col = names(head(CountColor, SelectCol)))

#おまけ
#明度は「+」または「-」で指定します
#明るく
Image1 <- Image + 0.4
display(Image1)
#暗く
Image2 <- Image - 0.4
display(Image2)

#コントラストは「*」で指定します
#明るく
Image3 <- Image * 1.3
display(Image3)
#暗く
Image4 <- Image * 0.7
display(Image4)

#ガンマ値は「^」で指定します
#高い
Image5 <- Image ^ 1.3
display(Image5)
#低い
Image6 <- Image ^ 0.7
display(Image6)

#クリッピングもできます。範囲はdisplayコマンドで画像を表示して座標を取得すると便利です
#指定は[x軸の範囲（幅）, y軸の範囲（高さ）,]です
display(Image[39:159, 87:122,])

#回転も可能です
#translate(rotate(Image, 「+」または「-」で回転), 回転の原点）です
display(translate(rotate(Image, -45), c(0, 0)))

#カラーチャンネルも変更できます
colorMode(Image) <- Grayscale
display(Image)

#カラーにも簡単に戻せます
colorMode(Image) <- Color