こしあん
2019-01-12

TPUでアップサンプリングする際にエラーを出さない方法

Pocket
LINEで送る
Delicious にシェア

1.3k{icon} {views}



画像処理をしているとUpsamplingが必要になることがあります。Keras/TensorFlowではUpsampling2Dというレイヤーを使ってアップサンプリングができますが、このレイヤーがTPUだとエラーを出すので解決法を探しました。自分でアップサンプリングレイヤーを定義するとうまく行ったので、それを見ていきます。

環境:TensorFlow v1.12.0

TPU使わないとうまくいく

例えば、今簡単な例として、shape=(8,1,1,1)の0~7の数字の配列を画像と見立てます。これを縦横2倍のアップサンプリングをして、shape=(8,2,2,1)という形に変形します。

KerasのUpsampling2Dレイヤーを使った実装ではこうでしょう。

import tensorflow as tf
from tensorflow.keras import layers
from tensorflow.keras.models import Model
import numpy as np

def upsampling_model():
    input = layers.Input((1,1,1))
    x = layers.UpSampling2D(2)(input)
    return Model(input, x)

def upsampling_test():
    model = upsampling_model()

    X = np.arange(8).reshape(-1,1,1,1)
    y = model.predict(X)
    print(y.shape)

if __name__ == "__main__":
    upsampling_test()

これは実際うまく行って、出力のshapeは正しい値になります。

(8, 2, 2, 1)

UpSampling2DはTPUではうまくは行かない

ただこのUpsampling2Dを使った方法はTPUだとエラーを出します。TPUに変換する処理を含めたコードです。

import tensorflow as tf
from tensorflow.keras import layers
from tensorflow.keras.models import Model
import tensorflow.keras.backend as K
from tensorflow.contrib.tpu.python.tpu import keras_support
import numpy as np
import os

def upsampling_model():
    input = layers.Input((1,1,1))
    x = layers.UpSampling2D(2)(input)
    return Model(input, x)

def upsampling_test():
    model = upsampling_model()
    # TPUモデルに変換するためにコンパイルが必要なのでこの値に意味はない
    model.compile(tf.train.GradientDescentOptimizer(0.1), "mean_squared_error")

    tpu_grpc_url = "grpc://"+os.environ["COLAB_TPU_ADDR"]
    tpu_cluster_resolver = tf.contrib.cluster_resolver.TPUClusterResolver(tpu_grpc_url)
    strategy = keras_support.TPUDistributionStrategy(tpu_cluster_resolver)
    model = tf.contrib.tpu.keras_to_tpu_model(model, strategy=strategy)

    X = np.arange(8).reshape(-1,1,1,1)
    y = model.predict(X)
    print(y.shape)

if __name__ == "__main__":
    upsampling_test()

途中コンパイルで適当なオプティマイザーや損失関数を与えていますが、これはTPUのモデルに変換するのにコンパイルが必要なので、この例では、コンパイルで指定した値は全く意味がありません。

ちなみにこのコードはコンパイル失敗します。

RuntimeError: Compilation failed: Compilation failure: Detected unsupported operations when trying to compile graph cluster_1_5423988494562430942[] on XLA_TPU_JIT: ResizeNearestNeighbor (No registered 'ResizeNearestNeighbor' OpKernel for XLA_TPU_JIT devices compatible with node {{node tpu_140457990329008/up_sampling2d_1/ResizeNearestNeighbor}} = ResizeNearestNeighbor[T=DT_FLOAT, align_corners=false, _device="/device:TPU_REPLICATED_CORE"](infeed-infer, tpu_140457990329008/up_sampling2d_1/mul)
    .  Registered:  device='CPU'; T in [DT_DOUBLE]
  device='CPU'; T in [DT_FLOAT]
  device='CPU'; T in [DT_BFLOAT16]
  device='CPU'; T in [DT_HALF]
  device='CPU'; T in [DT_INT8]
  device='CPU'; T in [DT_UINT8]
  device='CPU'; T in [DT_INT16]
  device='CPU'; T in [DT_UINT16]
  device='CPU'; T in [DT_INT32]
  device='CPU'; T in [DT_INT64]
){{node tpu_140457990329008/up_sampling2d_1/ResizeNearestNeighbor}}

どうもNearestNeighborのアップサンプリングのコンパイルが失敗するみたいですね。

NumpyでのNearestNeighbor法

さて、もっと簡単にNumpyでの例を振り返りましょう。以前こちらの記事で紹介した方法です。

Numpyだけでサクッと画像を拡大する方法
https://blog.shikoan.com/numpy-upsampling-image/

実はNumpyではForループを一切使わずにNearestNeighbor法の画像拡大をすることができます。repeatを2回かませる方法です。

image.repeat(2, axis=0).repeat(2, axis=1)

実はKerasのバックエンド関数に同様のrepeat関数があるため(K.repeat_elements())、このNumpyの方法をそのままTensorFlowに転用することができます。つまり、TensorFlow用のアップサンプリング用の関数を自分で書いてLambdaレイヤーでラップしてあげればOKです。

TPUで動く方法

こうしてみました。「upsampling2d_tpu」というのが自分で定義したアップサンプリング用の関数です。

import tensorflow as tf
from tensorflow.keras import layers
from tensorflow.keras.models import Model
from tensorflow.contrib.tpu.python.tpu import keras_support
import tensorflow.keras.backend as K
import numpy as np
import os

def upsampling_model():
    input = layers.Input((1,1,1))
    x = layers.Lambda(upsampling2d_tpu, arguments={"scale":2})(input)
    return Model(input, x)

def upsampling2d_tpu(inputs, scale=2):
    x = K.repeat_elements(inputs, scale, axis=1)
    x = K.repeat_elements(x, scale, axis=2)
    return x

def upsampling_test():
    model = upsampling_model()
    # TPUモデルに変換するためにコンパイルが必要なのでこの値に意味はない
    model.compile(tf.train.GradientDescentOptimizer(0.1), "mean_squared_error")

    tpu_grpc_url = "grpc://"+os.environ["COLAB_TPU_ADDR"]
    tpu_cluster_resolver = tf.contrib.cluster_resolver.TPUClusterResolver(tpu_grpc_url)
    strategy = keras_support.TPUDistributionStrategy(tpu_cluster_resolver)
    model = tf.contrib.tpu.keras_to_tpu_model(model, strategy=strategy)

    X = np.arange(8).reshape(-1,1,1,1)
    y = model.predict(X)
    print(y.shape)

if __name__ == "__main__":
    upsampling_test()

これはうまくいきます

INFO:tensorflow:New input shapes; (re-)compiling: mode=infer (# of cores 8), [TensorSpec(shape=(1, 1, 1, 1), dtype=tf.float32, name='input_12_10')]
INFO:tensorflow:Overriding default placeholder.
INFO:tensorflow:Remapping placeholder for input_12
INFO:tensorflow:Started compiling
INFO:tensorflow:Finished compiling. Time elapsed: 0.20801544189453125 secs
INFO:tensorflow:Setting weights on TPU model.
(8, 2, 2, 1)

TPUでもちゃんとアップサンプリングを行うことができました。NearestNeighbor法での拡大でいいところに、わざわざConv2DTransposeなどのレイヤーをはさむ必要はなさそうです。

まとめ

TPUではどうもUpSampling2Dのコンパイルに失敗するっぽい(今後のバージョンで改善される可能性あり)。K.repeat_elementsなどの関数を使い、自分でアップサンプリングの処理を書くとエラーにならないよ、ということでした。



Shikoan's ML Blogの中の人が運営しているサークル「じゅ~しぃ~すくりぷと」の本のご案内

技術書コーナー

【新刊】インフィニティNumPy――配列の初期化から、ゲームの戦闘、静止画や動画作成までの221問

「本当の実装力を身につける」ための221本ノック――
機械学習(ML)で避けて通れない数値計算ライブラリ・NumPyを、自在に活用できるようになろう。「できる」ための体系的な理解を目指します。基礎から丁寧に解説し、ディープラーニング(DL)の難しいモデルで遭遇する、NumPyの黒魔術もカバー。初心者から経験者・上級者まで楽しめる一冊です。問題を解き終わったとき、MLやDLなどの発展分野にスムーズに入っていけるでしょう。

本書の大きな特徴として、Pythonの本でありがちな「NumPyとML・DLの結合を外した」点があります。NumPyを理解するのに、MLまで理解するのは負担が大きいです。本書ではあえてこれらの内容を書いていません。行列やテンソルの理解に役立つ「従来の画像処理」をNumPyベースで深く解説・実装していきます。

しかし、問題の多くは、DLの実装で頻出の関数・処理を重点的に取り上げています。経験者なら思わず「あー」となるでしょう。関数丸暗記では自分で実装できません。「覚える関数は最小限、できる内容は無限大」の世界をぜひ体験してみてください。画像編集ソフトの処理をNumPyベースで実装する楽しさがわかるでしょう。※紙の本は電子版の特典つき

モザイク除去から学ぶ 最先端のディープラーニング

「誰もが夢見るモザイク除去」を起点として、機械学習・ディープラーニングの基本をはじめ、GAN(敵対的生成ネットワーク)の基本や発展型、ICCV, CVPR, ECCVといった国際学会の最新論文をカバーしていく本です。
ディープラーニングの研究は発展が目覚ましく、特にGANの発展型は市販の本でほとんどカバーされていない内容です。英語の原著論文を著者がコードに落とし込み、実装を踏まえながら丁寧に解説していきます。
また、本コードは全てTensorFlow2.0(Keras)に対応し、Googleの開発した新しい機械学習向け計算デバイス・TPU(Tensor Processing Unit)をフル活用しています。Google Colaboratoryを用いた環境構築不要の演習問題もあるため、読者自ら手を動かしながら理解を深めていくことができます。

AI、機械学習、ディープラーニングの最新事情、奥深いGANの世界を知りたい方にとってぜひ手にとっていただきたい一冊となっています。持ち運びに便利な電子書籍のDLコードが付属しています。

「おもしろ同人誌バザールオンライン」で紹介されました!(14:03~) https://youtu.be/gaXkTj7T79Y?t=843

まとめURL:https://github.com/koshian2/MosaicDeeplearningBook
A4 全195ページ、カラー12ページ / 2020年3月発行

Shikoan's ML Blog -Vol.1/2-

累計100万PV超の人気ブログが待望の電子化! このブログが電子書籍になって読みやすくなりました!

・1章完結のオムニバス形式
・機械学習の基本からマニアックなネタまで
・どこから読んでもOK
・何巻から読んでもOK

・短いものは2ページ、長いものは20ページ超のものも…
・通勤・通学の短い時間でもすぐ読める!
・読むのに便利な「しおり」機能つき

・全巻はA5サイズでたっぷりの「200ページオーバー」
・1冊にたっぷり30本収録。1本あたり18.3円の圧倒的コストパフォーマンス!
・文庫本感覚でお楽しみください

北海道の駅巡りコーナー

日高本線 車なし全駅巡り

ローカル線や秘境駅、マニアックな駅に興味のある方におすすめ! 2021年に大半区間が廃線になる、北海道の日高本線の全区間・全29駅(苫小牧~様似)を記録した本です。マイカーを使わずに、公共交通機関(バス)と徒歩のみで全駅訪問を行いました。日高本線が延伸する計画のあった、襟裳岬まで様似から足を伸ばしています。代行バスと路線バスの織り成す極限の時刻表ゲームと、絶海の太平洋と馬に囲まれた日高路、日高の隠れたグルメを是非たっぷり堪能してください。A4・フルカラー・192ページのたっぷりのボリュームで、あなたも旅行気分を漫喫できること待ったなし!

見どころ:日高本線被災区間(大狩部、慶能舞川橋梁、清畠~豊郷) / 牧場に囲まれた絵笛駅 / 窓口のあっただるま駅・荻伏駅 / 汐見の戦争遺跡のトーチカ / 新冠温泉、三石温泉 / 襟裳岬

A4 全192ページフルカラー / 2020年11月発行


Pocket
LINEで送る
Delicious にシェア

Add a Comment

メールアドレスが公開されることはありません。 * が付いている欄は必須項目です