Rでお遊び：「47都道府県の代表的なアニメ」と「キャラクターの髪の色」のまとめ

#パッケージのインストール
install.packages(”png", "jpeg", "colorspace", "DT", "ggplot2", "XLConnect", "maptools")

#パッケージの読み込み
library("tcltk")
library("png")
library("jpeg")
library("colorspace")

#読み込みフォルダを選択
FoldPath <- paste(as.character(tkchooseDirectory(title = "読み込みフォルダを選択"), sep = "", collapse =""))

#フォルダ内のアイテムを取得
ItemList <- list.files(path = FoldPath)

#データ格納用の変数を初期化
MasterRGBCol <- NULL

#アイテム数の取得
ItemVol <- length(ItemList)

for(i in seq(ItemVol)){
  
  #Imageファイルの読み込み
  selectAImage <- paste(FoldPath, "/", ItemList[i],
                           sep = "", collapse = "")
  
  #jpeg/png/gifファイルの判断
  if (unlist(strsplit(ItemList[i], "\\."))[2] == "png"){
    
    Image <- readPNG(selectAImage)
    LABCol <- as(RGB(as.vector(Image[,, 1]), as.vector(Image[,, 2]), as.vector(Image[,, 3])), "LAB")
    RGBCol <- cbind(LABCol@coords[,1:3], hex(LABCol))
    
  }else{
    
    Image <- readJPEG(selectAImage, native = FALSE)
    LABCol <- as(RGB(as.vector(Image[,, 1]), as.vector(Image[,, 2]), as.vector(Image[,, 3])), "LAB")
    RGBCol <- cbind(LABCol@coords[,1:3], hex(LABCol))
    
  }
  
  #マスターデータのデータフレーム化
  RGBCol <- as.data.frame(RGBCol)
  
  #カラーコードを文字列化
  RGBCol[, 4] <- as.character(RGBCol[, 4])
  #head(RGBCol)
  
  #画像で使用されているカラーコードを集計
  CountColor <- as.data.frame(table(RGBCol[, 4]))
  
  #集計データをカラーコードで並び替え
  CountColor <- CountColor[order(CountColor[, 1], decreasing = TRUE), ]
  head(CountColor)
  
  #マスターデータの重複を削除
  RGBCol <- RGBCol[!duplicated(RGBCol[, 4]),]
  
  #マスターデータをカラーコードで並び替え
  RGBCol <- RGBCol[order(RGBCol[, 4], decreasing = TRUE),]
  
  #結合
  RGBCol <- cbind(unlist(strsplit(ItemList[i], "_"))[2], unlist(strsplit(unlist(strsplit(ItemList[i], "_"))[3], "\\."))[1],
                     ItemList[i], RGBCol, CountColor[, 2])
  
  #プロットデータの作成
  MasterRGBCol <- rbind(MasterRGBCol, RGBCol)
  
  #必要のないImage,LABCol,RGBColの削除
  rm(Image,LABCol,RGBCol)
  
}

#行名の付与
colnames(MasterRGBCol) <- c("ID_1", "NAME_1", "ファイル名", "明度:l", "補色次元:a", "補色次元:b", "カラーコード", "出現数")

#ファイルの保存、解析画像数によってはファイルサイズが非常に大きい場合があります
#また、エクセルなどで読み込めないことがあるのでコメントアウトしています。
#setwd(paste(as.character(tkchooseDirectory(title = "保存場所を選択"), sep = "", collapse ="")))
#write.csv(MasterRGBCol, "分析結果.csv", fileEncoding = "CP932", eol = "\r\n")

#軸の設定
xRange <- seq(0, 1, by = 1/ItemVol)

#グラフ幅の設定
xWidth <- seq(1/ItemVol, 1, by = 1/ItemVol)

#プロット領域の作成
par(bg = "#2E2E2E")
plot(0:1, 0:1, type = "n", axes = FALSE, xlab = "", ylab = "")

#日本地図プロットデータ格納用変数
MapDrowData <- NULL

for(k in seq(ItemVol)){
  
  #データの抽出
  Plot <- subset(MasterRGBCol, MasterRGBCol[, 3] == ItemList[k])
  
  #出現数で並び替えてTOP1を抽出
  Plot <- Plot[order(as.vector(Plot[, 8]), decreasing = TRUE),]
  Plot <- head(Plot, n = 1)
  
  #日本地図プロット用データ
  MapDrowData <- rbind(MapDrowData, Plot)
  
  #データのプロット
  rasterImage(as.raster(Plot[, 7]),
              xRange[k], 0, xWidth[k], 1, interpolate = FALSE)
  
}

###ライブラリーの読み込み#####
library("ggplot2")
library("maptools")
########

#シェイプファイルの読み込み
#ファイルはJPN_adm1.shpです
ReadData <- paste(as.character(tkgetOpenFile(title = "マップデータを選択",
                                                filetypes = '{"shpファイル" {".shp"}}',initialfile = "*.shp")), sep = "", collapse =" ")
MapData <- readShapeSpatial(ReadData)

#日本地図プロット用データをIDで並び替え
MapDrowData <- MapDrowData[order(type.convert(as.character(MapDrowData[, 1])), decreasing = FALSE),]

#沖縄の地図データを準備
#沖縄を抽出
OkinawaMap <- MapData[MapData$NAME_1 == "Okinawa",]
#データフレーム化
OkinawaMap <- fortify(OkinawaMap) 
#抽出
OkinawaMapData <- subset(OkinawaMap, OkinawaMap[,1] &gt; 127.5 &amp; OkinawaMap[,1] &lt; 128.5 &amp; OkinawaMap[,2] &gt; 26 &amp; OkinawaMap[,2] &lt; 26.9)
#X軸の描写位置を変更
OkinawaMapData[, 1] <- OkinawaMapData[, 1] + 6 
#y軸の描写位置を変更
OkinawaMapData[, 2] <- OkinawaMapData[, 2] + 14 

#沖縄以外の地図データの準備
#沖縄を除く
JapanMapData <- MapData[MapData$NAME_1 != "Okinawa",]
#データフレーム化
JapanMapData <- fortify(JapanMapData)
#データの抽出
JapanMapData <- subset(JapanMapData, JapanMapData[,1] &gt; 128.5 &amp; JapanMapData[,2] &gt; 30)

#地図データの統合
#データの結合
AllJapanDrow <- rbind(JapanMapData, OkinawaMapData)
#おまじない
AllJapanDrow[, 7] <- type.convert(AllJapanDrow[, 7])

#境界線データの準備
BorderLine <- data.frame(DrowXaxis = c(133, 134.5, 135, 135), DrowYaxis = c(39, 39, 41, 43))

#日本地図のプロット
#時間がかかります
JapanDrow <- ggplot() +
  geom_polygon(data = AllJapanDrow, aes(x = AllJapanDrow[, 1], y = AllJapanDrow[, 2], group = AllJapanDrow[, 6], fill = AllJapanDrow[, 7])) +
  geom_path(data = BorderLine, aes(x = BorderLine[, 1], y = BorderLine[, 2])) +
  scale_fill_gradientn(colours = MapDrowData[, 7], guide = FALSE) + labs(x = "", y = "") +
  theme(axis.ticks = element_blank(), axis.text = element_blank(),
        panel.grid = element_blank(), panel.background = element_rect(fill = "#2E2E2E")) 
print(JapanDrow)

#パッケージの読み込み
library("DT")
library("XLConnect")

#HTMLファイルの保存場所を指定
setwd(paste(as.character(tkchooseDirectory(title = "保存場所を選択"), sep = "", collapse ="")))

#エクセルデータの読み込み
LoadData <- loadWorkbook(paste(as.character(tkgetOpenFile(title = "xlsxファイルを選択",
                                                             filetypes = '{"xlsxファイル" {".xlsx"}}',
                                                             initialfile = "*.xlsx")), sep = "", collapse =" "))
#シートの読み込み。読み込んだデータはデータフレームになります。
SheetData <- readWorksheet(LoadData, sheet = 1)

#セル用データをファイル名で並び替え
MapDrowData <- MapDrowData[order(type.convert(as.character(MapDrowData[, 3])), decreasing = FALSE),]

#参考URLをHTMLタグへ変換##
SheetData[, 5] <- paste0('&lt;a href="', SheetData[, 5], '" target="_blank"&gt;', SheetData[, 5], '&lt;/a&gt;')

#テーブルの作成
DataTable <- datatable(SheetData, rownames = FALSE, options = list(pageLength = 5, lengthMenu = c(3, 10, 25, 50, 100)),
                          caption = '47都道府県 アニメキャラ髪の色', escape = FALSE)

#キャラ名のセルの色をパイプ"%>%"でつなげてformatStyleで設定します
DataTable <- DataTable %>% formatStyle('キャラ',
                                             backgroundColor = styleEqual(SheetData[, 4], MapDrowData[, 7]))
#作成テーブルをhtmlで出力
saveWidget(DataTable, 'DataTable.html')