第６６回 チョコボール計測
今日の計測報告です。
昨日初めてエンゼルが当たった同僚から２箱もらいました。連続なるか！？
計測結果
| date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
|---|---|---|---|---|---|---|---|---|---|
| 2018-06-07 | 2018-12-01 | 34.087 | 4.788 | 16 | 小山工場 | コンビニ（千代田区） | なし | 29.299 | 1.831 |
| 2018-06-07 | 2018-12-01 | 34.22 | 4.784 | 16 | 小山工場 | コンビニ（千代田区） | なし | 29.436 | 1.840 |
連続ならず！
基礎集計
この集計はピーナツ味のチョコボールの集計結果です。
| 項目 | 数 |
|---|---|
| 計測データ数 | 191 |
| 銀のエンゼル出現数 | 5 |
| 金のエンゼル出現数 | 1 |
| 最小 | 中央値 | 最大値 | 平均 | |
|---|---|---|---|---|
| 正味重量 | 28.913 | 29.387 | 32.232 | 29.481 |
| 個数 | 14.000 | 16.000 | 20.000 | 16.445 |
この図は正味の重量のヒストグラムです。
赤い縦線が仕様（28g）を表しています。
青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
| 項目 | 数 |
|---|---|
| 計測データ数 | 115 |
| 銀のエンゼル出現数 | 7 |
| 金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が3.61%、上側が11.44%という結果です。
だいぶ分布が収束してきました。 現在の期待値は6.97%です。
次に金のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が0.04%、上側が2.67%という予測になっています。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 出版社/メーカー: 森永製菓
- 発売日: 2016/03/01
- メディア: 食品&飲料
- この商品を含むブログを見る
SIMERST 改良版 携帯タイプはかり ポケットデジタルスケール(秤) 0.001g-100g精密 業務用(プロ用) デジタルスケール 電子天秤
- 出版社/メーカー: Simerst
- メディア:
- この商品を含むブログを見る
第６５回 チョコボール計測
今日の計測報告です。
同僚から５箱もらいました。
計測結果
| date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
|---|---|---|---|---|---|---|---|---|---|
| 2018-06-06 | 2018-12-01 | 33.999 | 4.748 | 16 | 小山工場 | コンビニ（千代田区） | なし | 29.251 | 1.828 |
| 2018-06-06 | 2018-12-01 | 34.104 | 4.792 | 16 | 小山工場 | コンビニ（千代田区） | 銀 | 29.312 | 1.832 |
| 2018-06-06 | 2018-12-01 | 34.758 | 4.725 | 18 | 小山工場 | コンビニ（千代田区） | なし | 30.033 | 1.669 |
| 2018-06-06 | 2018-12-01 | 34.258 | 4.896 | 18 | 小山工場 | コンビニ（千代田区） | なし | 29.362 | 1.631 |
| 2018-06-06 | 2018-12-01 | 34.152 | 4.885 | 18 | 小山工場 | コンビニ（千代田区） | なし | 29.267 | 1.626 |
銀のエンゼルさん一枚出ました！
今回チョコボールをカンパしてくれた同僚は、ここまで一回もエンゼルが出ていなかったので、これが初出現です。
基礎集計
この集計はピーナツ味のチョコボールの集計結果です。
| 項目 | 数 |
|---|---|
| 計測データ数 | 189 |
| 銀のエンゼル出現数 | 5 |
| 金のエンゼル出現数 | 1 |
| 最小 | 中央値 | 最大値 | 平均 | |
|---|---|---|---|---|
| 正味重量 | 28.913 | 29.387 | 32.232 | 29.482 |
| 個数 | 14.000 | 16.000 | 20.000 | 16.450 |
この図は正味の重量のヒストグラムです。
赤い縦線が仕様（28g）を表しています。
青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
| 項目 | 数 |
|---|---|
| 計測データ数 | 113 |
| 銀のエンゼル出現数 | 7 |
| 金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が3.54%、上側が10.84%という結果です。
次に金のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が0.03%、上側が2.73%という予測になっています。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 出版社/メーカー: 森永製菓
- 発売日: 2016/03/01
- メディア: 食品&飲料
- この商品を含むブログを見る
SIMERST 改良版 携帯タイプはかり ポケットデジタルスケール(秤) 0.001g-100g精密 業務用(プロ用) デジタルスケール 電子天秤
- 出版社/メーカー: Simerst
- メディア:
- この商品を含むブログを見る
第６４回 チョコボール計測
今日の計測報告です。
大人買いチョコボールの開封の続きです。
さいたま市のまちおかで買った２０箱セットは今日で全て開封しました。
計測結果
| date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
|---|---|---|---|---|---|---|---|---|---|
| 2018-05-28 | 2018-11-01 | 35.81 | 4.942 | 17 | 小山工場 | スーパー（さいたま市） | なし | 30.868 | 1.816 |
| 2018-05-28 | 2018-11-01 | 37.118 | 4.886 | 16 | 小山工場 | スーパー（さいたま市） | なし | 32.232 | 2.014 |
| 2018-05-28 | 2018-11-01 | 34.486 | 4.939 | 17 | 小山工場 | スーパー（さいたま市） | なし | 29.547 | 1.738 |
| 2018-05-28 | 2018-11-01 | 34.357 | 4.962 | 16 | 小山工場 | スーパー（さいたま市） | なし | 29.395 | 1.837 |
エンゼルさんは現れませんでした。
今回の２０箱セットの中で、エンゼルは銀のエンゼルが一つだけでした。
基礎集計
この集計はピーナツ味のチョコボールの集計結果です。
| 項目 | 数 |
|---|---|
| 計測データ数 | 184 |
| 銀のエンゼル出現数 | 4 |
| 金のエンゼル出現数 | 1 |
| 最小 | 中央値 | 最大値 | 平均 | |
|---|---|---|---|---|
| 正味重量 | 28.913 | 29.389 | 32.232 | 29.483 |
| 個数 | 14.000 | 16.000 | 20.000 | 16.429 |
この図は正味の重量のヒストグラムです。
赤い縦線が仕様（28g）を表しています。
青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
| 項目 | 数 |
|---|---|
| 計測データ数 | 108 |
| 銀のエンゼル出現数 | 6 |
| 金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が2.98%、上側が10.75%という結果です。
次に金のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が0.05%、上側が2.53%という予測になっています。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 出版社/メーカー: 森永製菓
- 発売日: 2016/03/01
- メディア: 食品&飲料
- この商品を含むブログを見る
SIMERST 改良版 携帯タイプはかり ポケットデジタルスケール(秤) 0.001g-100g精密 業務用(プロ用) デジタルスケール 電子天秤
- 出版社/メーカー: Simerst
- メディア:
- この商品を含むブログを見る
第６３回 チョコボール計測
今日の計測報告です。
先週末から福岡に出張に行っていたので、福岡で購入したチョコボールを開封します。
計測結果
| date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
|---|---|---|---|---|---|---|---|---|---|
| 2018-05-25 | 2018-12-01 | 34.236 | 4.791 | 17 | 小山工場 | コンビニ（福岡） | なし | 29.445 | 1.732 |
福岡のチョコボールではエンゼルさん現れませんでした。
基礎集計
この集計はピーナツ味のチョコボールの集計結果です。
| 項目 | 数 |
|---|---|
| 計測データ数 | 180 |
| 銀のエンゼル出現数 | 4 |
| 金のエンゼル出現数 | 1 |
| 最小 | 中央値 | 最大値 | 平均 | |
|---|---|---|---|---|
| 正味重量 | 28.913 | 29.387 | 31.487 | 29.460 |
| 個数 | 14.000 | 16.000 | 20.000 | 16.425 |
この図は正味の重量のヒストグラムです。
赤い縦線が仕様（28g）を表しています。
青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
| 項目 | 数 |
|---|---|
| 計測データ数 | 104 |
| 銀のエンゼル出現数 | 6 |
| 金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が3.13%、上側が11.22%という結果です。
次に金のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が0.04%、上側が2.92%という予測になっています。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
チョコボール画像からチョコボールの個数を自動計測してみる
【update:2021/03/31】 ソースコードの全面的な更新に伴って、解説記事も書き直しました。本記事は記述に不足がありますので、改訂版を参照ください。
【update:2021/03/25】 本記事で掲載しているソースコードに一部不足がありました。修正して整理したソースコードを以下のGithubリポジトリにupしています。
概要
- 当ブログではこれまで、チョコボールを手動で計測してきました
- しかし、毎回チョコボールの数を数えるのはめんどくさい
- そこで、画像から自動でカウントさせよう（チョコボールに特化した一般物体認識のモデルを使ってカウント）
- Faster R-CNNをfine tuneしたところ、粒が密集した画像でもそれぞれの粒を認識できていそう！すごい！
はじめに
当ブログのメインコンテンツは、ピーナツ味のチョコボールの計測です。 重さを測ったり、個数を数えたりしているのですが、従来、それらは全て手動で行われています。
一方、画像処理分野を発端として、深層学習（Deep Learning）の多方面での有効性が知られています。
（深層学習技術の進展はついていくことすら難しいレベルでどんどん進んでいますね。。。）
特に、一般物体認識（Faster R-CNN, YOLO等）は、専門に調査・研究していない人でもテレビ等で見たことがあるのでは無いでしょうか？
この一般物体認識ですが、ザクっと言うと、特定の物体が画像内のどこにいくつ存在するのかを認識するものです。 この技術を使えばチョコボールの個数をわざわざ数える必要がなくなるのでは無いでしょうか？！ 毎回毎回数える時間コストを考えると、かなりの効率化が図れるような図れないような気がしてきました。
問題設定
ということで今回は、 チョコボール画像からチョコボールの個数を自動で計測する仕組みを作ります。
データ
毎回の計測で、以下のような画像を撮影しています。
スマートフォンのカメラを使って手動で撮影しています。 そのため、影（照明環境）や距離が一定ではありません。
アプローチ
深層学習の適用事例として有名なものの一つに「一般物体認識」というタスクがあります。 これは、画像の中に写っている複数の物体を個別に認識する技術です（下図参照）。
この技術を応用し、画像内に写っているチョコボールを個別に認識することで、個数が計量できるというアプローチを取ります。
ただし、一般物体認識と言っても全ての物体を認識できるわけではなく、 予め学習した物体のみを認識できます。 当然、公開されている学習済みモデルではチョコボールを認識してはくれません。 そのため、独自の学習データセットを作成し、チョコボールが認識できるように追加学習（fine tune）することが必要です。
前処理
ということで、まずは学習データセットを作ります。
アノテーションデータ作成
一般物体認識モデルを学習するためには、 学習データとして画像にチョコボールの位置を示す情報を付与したデータセットを作成する必要があります。 このデータセットのことをアノテーションデータセットと呼びます(下図参照)。
例えば、このような画像ファイルに対して
この図のように、「どこ」に「何」があるのかを示すデータを作ります。
アノテーションデータセットは、画像データ、物体位置座標（長方形の対角位置座標）のリスト、物体カテゴリのリストがセットで必要になります。
このデータを作るために、LabelImgというツールを利用しました。
使い方はgithubのページを見ればわかると思います。
なお今回は、上図のように、認識対象の物体カテゴリとして、 choco-ball(チョコボールの粒)とchoco-package(パッケージ画像)の２つの物体を認識するものとしました。
PascalVOC形式ｘｍｌデータのパース
LabelImgというツールでは、 アノテーションデータをPascalVOC形式のxmlファイルで出力できます。
このままでは学習用データセットにはならないので、xmlファイルを読みこんで必要な情報を抜き出す必要があります。 PascalVOC形式のxmlをパースする便利なライブラリは私がちょっと探したところ見当たらなかったので、 以下のように作りました。
xmltodictというxmlファイルをdict形式でパースするライブラリを使い、 画像データ(imgs)、バウンディングボックスの座標リスト(bboxs)、物体カテゴリIDのリスト(obj_ids)をnumpyオブジェクトで書きだしています。
認識モデル
一般物体認識モデルは複数提案されていますが、今回はFaster R-CNNを使います。 Faster R-CNNについては、日本語での解説記事も多数あるのでここでは言及しません。
認識する物体として今回は、choco-ballとchoco-packageの２つだけを認識するようにしました。 今回のタスクにおいては、一般物体を認識する必要など無いのです。
Faster-RCNNの実装はChianerCVを使いました。 ChainerCVとはChainerベースの画像処理ライブラリです。 後ろにコードを載せていますが、かなり簡潔に書けます。
学習実験
いよいよ学習です。
実験設定＆コード
今回の実験では、学習データとしてチョコボール画像22枚を用意しました。 （画像自体はもっとあるのですが、アノテーション作業が大変だったので。。。） このデータを8:2で分割し、学習データと評価データに分けました
- 学習データ：17ファイル
- 評価データ：5ファイル
学習係数などのパラメータはChainerCV公式exampleの設定をそのまま使いました。
コードは以下の通りです。 ChainerCVのexampleを ほぼほぼコピーさせてもらっています。
４セル目から６セル目でnumpyオブジェクトで保存した画像データ、バウンディングボックス座標リスト、物体カテゴリIDリストを
読み込んでいます。
読んだデータは、ChainerCVで扱いやすいようにTupleDatasetという形式のデータにしています（７セル目）。
９セル目でFasterRCNNのオブジェクトを作っています。 引数で今回認識対象とする物体（choco-ballとchoco-package）の指定をしています。 また、pretrainedmodelには、imagenetを指定します。 独自データで学習する際には、imagenetを指定する必要があるらしいです。
optimizerやiteratorはほぼexampleの通りです。
１７〜21セル目でextentionの指定をしています。 公式exampleではiteration数でモデルの書き出し等設定していましたが、 epoch数で制御したほうが分かりやすかったのでepochで規準を作っています。
学習結果
上記のnotebookの23セル目に学習曲線（lossカーブ）がありますが、 20epoch学習させてみた結果です。 想定通りlossが落ちたように見えます。
もう少し学習させても良さそうですが、一旦これで評価してみます。 notebookの24セル以降が評価用のコードです。
評価用の5データのチョコボール数を算出したところ、
MSE（平均二乗誤差）が1なので、平均1個の誤認識があるようです。
これで十分な精度かどうかというと、、、ちょっと微妙かな。。。
（正確にはcross-validation必要）
試してみる
とにかく、学習したモデルを使ってチョコボールの認識をしてみましょう。 アノテーションをつけていない画像を入力して結果を見てみます。
うまく行った例
1枚目はそれぞれの粒がばらついているので簡単そうです。 2,3枚目はかなり密集していますが、それでもそれぞれの粒を認識できているようです。
検出ミスがあった例
やはり、密集していると見落としがちなようです。 ３枚目の歪な形の物も見落としています。
まとめ
２２枚のデータを使ってFaster R-CNNをfine tuneしてみました。 結果としては、チョコボールの粒が密集した画像でも正しく粒を認識できているものがあるので、 この仕組みは有望だと感じます。 しかし、評価データで定量的に誤差を計算すると、平均１粒の計数ミスが生じるようで、まだ実用的というのは難しいかもしれません。
評価に関しては、評価データが５つしか無いのでcross-validationが必要です。 今回ざっくりと作ってみたところ有効そうな結果が見えてきたので、 データをもっと作り、また、cross-validationによって正確な評価をしてみようと思います。（今後の課題）
そういえば、粒のカウント工数を減らすことが目的だったんですが、、、 なかなか手間が減らなさそう（汗）
参考資料
- ChainerCV — ChainerCV 0.13.1 documentation
- chainercvの公式ドキュメント。とりあえあずこれを見れば動かせるようになります
- GitHub - chainer/chainercv: ChainerCV: a Library for Deep Learning in Computer Vision
- chainercvの公式Git
- GitHub - tzutalin/labelImg: 🖍️ LabelImg is a graphical image annotation tool and label object bounding boxes in images
- 画像処理周りを中心に、理論を含めて比較的わかりやすいと思います。 これを読んで、Tensorflowなどのライブラリを中心にした参考書を読んで実装すれば良いんじゃないかな。
- 作者:岡谷 貴之
- 発売日: 2015/04/08
- メディア: 単行本
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
- 発売日: 2016/03/01
- メディア: 食品&飲料
第６２回 チョコボール計測
今日の計測報告です。
大人買いチョコボールの開封の続きです。
これが
こうなりました。
計測結果
| date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
|---|---|---|---|---|---|---|---|---|---|
| 2018-05-22 | 2018-11-01 | 34.389 | 4.931 | 17 | 小山工場 | スーパー（さいたま市） | なし | 29.458 | 1.733 |
| 2018-05-22 | 2018-11-01 | 34.645 | 4.919 | 16 | 小山工場 | スーパー（さいたま市） | なし | 29.726 | 1.858 |
| 2018-05-22 | 2018-11-01 | 33.897 | 4.893 | 16 | 小山工場 | スーパー（さいたま市） | なし | 29.004 | 1.813 |
| 2018-05-22 | 2018-11-01 | 34.428 | 4.89 | 17 | 小山工場 | スーパー（さいたま市） | なし | 29.538 | 1.738 |
今回はエンゼルさん現れませんでした。
基礎集計
この集計はピーナツ味のチョコボールの集計結果です。
| 項目 | 数 |
|---|---|
| 計測データ数 | 179 |
| 銀のエンゼル出現数 | 4 |
| 金のエンゼル出現数 | 1 |
| 最小 | 中央値 | 最大値 | 平均 | |
|---|---|---|---|---|
| 正味重量 | 28.913 | 29.387 | 31.487 | 29.460 |
| 個数 | 14.000 | 16.000 | 20.000 | 16.425 |
この図は正味の重量のヒストグラムです。
赤い縦線が仕様（28g）を表しています。
青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
| 項目 | 数 |
|---|---|
| 計測データ数 | 103 |
| 銀のエンゼル出現数 | 6 |
| 金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が3.21%、上側が10.84%という結果です。
次に金のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が0.05%、上側が3.07%という予測になっています。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
第６１回 チョコボール計測
今日の計測報告です。
同僚から４箱もらいましたので開封します。
計測結果
| date | best_before | weight | box_weight | number | factory | shop | angel | net_weight | mean_weight |
|---|---|---|---|---|---|---|---|---|---|
| 2018-05-19 | 2018-10-01 | 34.16 | 4.773 | 15 | 小山工場 | ドラックストア（横浜市） | なし | 29.387 | 1.959 |
| 2018-05-19 | 2018-11-01 | 34.691 | 4.852 | 15 | 小山工場 | ドラックストア（横浜市） | なし | 29.839 | 1.989 |
| 2018-05-19 | 2018-11-01 | 34.205 | 4.879 | 15 | 小山工場 | ドラックストア（横浜市） | なし | 29.326 | 1.955 |
| 2018-05-19 | 2018-10-01 | 35.001 | 4.769 | 16 | 小山工場 | ドラックストア（横浜市） | なし | 30.232 | 1.889 |
今回はエンゼルさん現れませんでした。
基礎集計
この集計はピーナツ味のチョコボールの集計結果です。
| 項目 | 数 |
|---|---|
| 計測データ数 | 175 |
| 銀のエンゼル出現数 | 4 |
| 金のエンゼル出現数 | 1 |
| 最小 | 中央値 | 最大値 | 平均 | |
|---|---|---|---|---|
| 正味重量 | 28.913 | 29.385 | 31.487 | 29.461 |
| 個数 | 14.000 | 16.000 | 20.000 | 16.423 |
この図は正味の重量のヒストグラムです。
赤い縦線が仕様（28g）を表しています。
青い太線で正規分布と仮定した最尤推定量をプロットしています。
エンゼル出現確率の予測
通常版のエンゼルの予測を行っていきます。 これまでの通常版パッケージの開封結果は次の通りです。
| 項目 | 数 |
|---|---|
| 計測データ数 | 99 |
| 銀のエンゼル出現数 | 6 |
| 金のエンゼル出現数 | 0 |
この結果を使ってベイズ推定によるエンゼルの出現確率推定を行います。
はじめに銀のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が3.29%、上側が11.61%という結果です。
次に金のエンゼルの出現確率の推定です。
90%信用区間（上下それぞれ５％）は、下側が0.06%、上側が3.10%という予測になっています。
広告
Amazonの欲しいものリスト作ってみました。
チョコボールのカンパ募集中です。
チョコボールをカンパする
