JP2013037663A - Ｒｆｉｄタグ移動識別方法及びｒｆｉｄタグ移動識別プログラム - Google Patents

Abstract

【課題】一般的なＲＦＩＤリーダと一般的なアンテナから得られる読取ＲＦＩＤタグの電波受信強度、位相やドップラー周波数の時系列データによって、アンテナ前を移動したＲＦＩＤタグと静止したままのＲＦＩＤタグを識別するＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムを提供する。
【解決手段】ＲＦＩＤリーダによって移動ＲＦＩＤタグ及び静止ＲＦＩＤタグを含む複数のＲＦＩＤタグの時系列読取データを取得し、一定のパラメータに基づく静止ＲＦＩＤタグフィルタによって、読取った複数のＲＦＩＤタグの情報から静止ＲＦＩＤタグを特定する。静止ＲＦＩＤタグフィルタは、チラ読みフィルタ、ＲＳＳＩ絶対値フィルタ、ＲＳＳＩ非連続上昇フィルタ、ＲＳＳＩ谷型フィルタ、過読フィルタ、位相振れ少フィルタ、垂直方向移動フィルタのうち少なくとも一つ以上の個別静止ＲＦＩＤタグフィルタを用いて静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別し、移動ＲＦＩＤタグを識別するＲＦＩＤタグ移動識別方法である。
【選択図】図１

Description

本発明は、ＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムに係り、特に一般的ＲＦＩＤリーダやアンテナを使用したＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムに関する。

最近、ＵＨＦ帯のＲＦＩＤ（Ｒａｄｉｏ Ｆｒｅｑｕｅｎｃｙ Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ）技術を用いて、アンテナから１０ｍ以上といった遠距離に位置する多数のＲＦＩＤタグを一括して瞬時に読取ることが可能となっている。ＲＦＩＤタグチップの性能向上等により読取可能距離は今後更に伸びることが予想されている。

ＲＦＩＤタグは、倉庫・店舗等での高速入出荷検品や棚卸、盗難検知、商品のレジ精算といった多くの用途に使用が広がりつつある。しかしながら、アンテナ周囲に陳列された商品に貼付されたＲＦＩＤタグ等、本来意図した読取対象以外のＲＦＩＤタグも読みとってしまう問題があり、実際にアンテナ前を移動したＲＦＩＤタグのみを識別する技術が必要とされている。

そこで、ＲＦＩＤ タグからの受信情報を経時情報に関連付けるとともに、この受信情報を受信した際のスキャンアンテナのスキャン角と経時情報の組を複数個生成するデータ生成手段と、上記データ生成手段の生成したデータを用いて、上記スキャン角と上記経時情報の間の関係を示す線形近似直線を求め、この線形近似直線の傾きから上記ＲＦＩＤタグが付された移動体の移動方向を検知したり、上記移動体の移動軌跡を求めることにより上記ＲＦＩＤタグが付された移動体の移動方向を検知し、上記ＲＦＩＤタグが付された移動体の移動方向を検知する技術が知られている（特許文献１）。

また、反射体から、互いに異なる複数の搬送周波数によって送信された反射信号を受信する処理を行う距離測定装置において、上記反射体から送信された上記反射信号を解析することによって、該反射体と当該距離測定装置との距離を算出する距離算出部を備え、反射体と距離測定装置との距離をより精度よく算出する技術が知られている（特許文献２）。
さらに、アンテナをＲＦＩＤタグの移動方向に水平ではなく斜めに配置するか出力電波の方向が歪んだ特殊アンテナを使用した技術が知られている（特許文献３）。

特開２００７−３０３９３５号公報 特開２００９−８０１３３号公報 ＵＳ７８１７０１４Ｂ２（Ｆｉｇ．１７ａ，Ｆｉｇ．１８ａ〜１８ｂ）

しかしながら、特許文献１の技術は、スキャンアンテナという特殊なアンテナが必要とされる。そしてその構成が複雑であるため、一般的なＲＦＩＤリーダと一般的なアンテナを用いた汎用的な装置構成には適用できないという不都合があった。
また特許文献２の技術は、ある１つのタグまでの距離を測定するために２つの異なる周波数にて電波送信を行う特殊なリーダが必要とされる。
さらに特許文献３の技術は、アンテナをＲＦＩＤタグの移動方向に水平ではなく斜めに配置するか出力電波の方向が歪んだ特殊アンテナを使用する必要があり、店舗内のゲートといった設置の自由度が少ない場所に適さない場合がある。またＲＦＩＤタグからの電波受信強度は人や物がアンテナ前を移動しただけでも変化してしまう場合があるため、その上昇や下降のみにより移動と判定する手法では誤検知が発生することがある。

電波出力強度やアンテナのハードウェア調整により読取範囲を狭めて所望のＲＦＩＤタグのみを読取る方法が取られることがあるが、遠距離・複数タグ読取のメリットが損なわれてしまう。また、電波のマルチパスフェージングによって予期せず遠くに位置するＲＦＩＤタグまで読まれてしまう場合がある。
またアンテナ２枚を離れた場所に設置したり角度をつけて設置したりすることでそれぞれのアンテナでの読取時間差が一定値以上発生した場合に移動と判定する方法が取られることがあるが、アンテナ間に距離や角度を持たせるための設置スペースが必要とされてしまう。また人や物の移動などによる電波の乱反射によって２つのアンテナ間に位置するタグ等が誤って移動と判定されてしまうことがある。

このように、遠距離・複数ＲＦＩＤタグ読取に優れたＵＨＦ帯ＲＦＩＤには、希望していないＲＦＩＤタグも読んでしまう「読み過ぎ」の問題がある。出力調整等により読取範囲を狭めてＵＨＦ帯のメリットを損ねるのではなく、リーダから取得した膨大な各種無線情報を高速に処理し希望タグ（移動タグ）を識別できるような技術が望まれている。

本発明の目的は、機種等に依存せずに、一般的なＲＦＩＤリーダと一般的なアンテナから得られる読取ＲＦＩＤタグの電波受信強度、位相やドップラー周波数の時系列データによって、アンテナ前を移動したＲＦＩＤタグと静止したままのＲＦＩＤタグを使用アンテナが１枚の場合でも識別することが可能なＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムを提供することにある。
本発明の他の目的は、コンピュータ装置上またはＲＦＩＤリーダ上で動作することが可能なＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムを提供することにある。
本発明の更に他の目的は、他のアプリケーションに識別情報を送信するミドルウエアとしてコンピュータ装置上またはＲＦＩＤリーダ上で動作することが可能なＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムを提供することにある。

前記課題は、本発明のＲＦＩＤタグ移動識別方法によれば、複数のＲＦＩＤタグと、該複数のＲＦＩＤタグからの電波を受信するＲＦＩＤリーダを用いて、複数のＲＦＩＤタグのうち移動ＲＦＩＤタグを識別するＲＦＩＤタグ移動識別方法において、前記ＲＦＩＤリーダによって移動ＲＦＩＤタグ及び静止ＲＦＩＤタグを含む複数のＲＦＩＤタグの時系列読取データを取得し、一定のパラメータに基づく静止ＲＦＩＤタグフィルタによって、前記読取った複数のＲＦＩＤタグの情報から前記静止ＲＦＩＤタグを特定し、前記複数のＲＦＩＤタグから前記移動ＲＦＩＤタグを識別可能とするもので、前記静止ＲＦＩＤタグフィルタは、以下の（１）乃至（６）の工程の何れかからなる個別静止ＲＦＩＤタグフィルタを備え、前記静止ＲＦＩＤタグフィルタは、複数の前記個別静止ＲＦＩＤタグフィルタのうち少なくとも一つ以上の個別静止ＲＦＩＤタグフィルタを用いて静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別し、移動ＲＦＩＤタグを識別することを特徴とするＲＦＩＤタグ移動識別方法；（１）．読取中のタグ枚数を変数とした関数を用いて、秒間読取可能回数を計算し、前記ＲＦＩＤリーダによって読み取った秒間読取回数との比較に基づいて識別する工程；（２）．電波受信強度の時系列分析によって、静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する工程；（３）．前記ＲＦＩＤリーダで読みとられたＲＦＩＤタグの電波受信強度時系列データが、谷型となる場合、すなわち電波受信強度の高い点発生時に所定の時系列データまで遡り、別の電波受信強度が高い点があるかを判定し、高い点があった場合に静止ＲＦＩＤタグと判別する工程；（４）．前記ＲＦＩＤタグの電波受信強度時系列データが示す山型全体、山型の登り部分を線形・非線形曲線として適合する度合い、山型の強度最大値と最小値の差、山型の時系列で前の点から上昇した点の割合のうち少なくとも一つを用いて静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する工程；（５）．前記ＲＦＩＤリーダの送信搬送波とＲＦＩＤタグからの反射で返却された搬送波の位相の時系列分析によって、静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する工程；（６）．ＲＦＩＤタグ応答の位相変曲点またはドップラー周波数０時点と電波受信強度山頂発生時点を比較することによって移動方向を識別する工程、とすることにより解決される。

本発明のＲＦＩＤタグ移動識別方法では、ＲＦＩＤタグ、ＲＦＩＤアンテナとＲＦＩＤリーダは、機種等に依存せずに、時系列データを読取り、識別することが可能となる。
また、静止ＲＦＩＤタグフィルタは、複数の前記個別静止ＲＦＩＤタグフィルタのうち少なくとも一つ以上の個別静止ＲＦＩＤタグフィルタを用いて静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別し、移動ＲＦＩＤタグを識別することにより、複数の観点から静止ＲＦＩＤタグを識別することが可能となり、高精度な移動ＲＦＩＤタグの識別が可能となる。

特に、一つ以上の個別静止ＲＦＩＤタグフィルタを用いて静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別するので、人や物、特に金属製のカゴ車や台車等がアンテナ前を移動すること等よる電波の乱反射やマルチパスフェージング、特に水分が含まれた物や人がアンテナ前を移動することによって発生する電波遮蔽によって電波受信強度が上下するような静止したＲＦＩＤタグにおいても、静止ＲＦＩＤタグと識別することができる。

このとき、前記（１）における読取中タグ枚数を変数とした関数は、読取中タグ枚数を独立変数とした累乗曲線関数または多項式関数であり、この関数を用い計算した計算値に、所定の割合を乗算した閾値を判断基準として、前記判断基準より秒間読取回数が少ない場合に静止タグと識別可能となる。

また、前記各種ＲＦＩＤタグフィルタを使用しアンテナ近くの水平方向移動のみを検知し遠くの水平方向移動や垂直方向移動は検知しないとすることによって、より高精度な移動識別や移動方向判定が行えるように構成すると好適である。

さらに、前記（５）における時系列分析において、位相の標準偏差の値に閾値を設けること、二項検定、カイ二乗検定、Ｇ検定の少なくとも一つを用いて静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する工程を備えるように構成すると好適である。

また、前記一定のパラメータは、予め試験的ＲＦＩＤタグを移動した後に、移動時の開始点及び終了点を選択することで前記パラメータを自動で計算するように構成すると好適である。このように構成すると、ＲＦＩＤリーダやアンテナ、タグ、移動速度といった環境の違いを吸収し最適なパラメータが自動計算されるので、パラメータ設定が容易となる。

前記課題は、本発明のＲＦＩＤタグ移動識別プログラムによれば、コンピュータ装置に、移動ＲＦＩＤタグ及び静止ＲＦＩＤタグを含む複数のＲＦＩＤタグの時系列読取データを、一定のパラメータに基づく静止ＲＦＩＤタグフィルタよって、複数の前記ＲＦＩＤタグの情報から前記静止ＲＦＩＤタグを特定させ、前記複数のＲＦＩＤタグから前記移動ＲＦＩＤタグを識別させるものであって、静止ＲＦＩＤタグフィルタとして、（１）．読取中のタグ枚数を変数とした関数を用いて求められた秒間読取可能回数を計算し、前記ＲＦＩＤリーダによって読み取った秒間読取回数との比較に基づいて識別する手順、（２）．電波受信強度の時系列分析によって、静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する手順、（３）．前記ＲＦＩＤリーダで読みとられたＲＦＩＤタグの電波受信強度時系列データが谷型となる場合、すなわち電波受信強度の高い点発生時に所定の時系列データまで遡り、別の電波受信強度が高い点があるかを判定し、高い点があった場合に静止ＲＦＩＤタグと判別する手順、（４）．前記ＲＦＩＤタグの電波受信強度時系列データが示す山型全体、山型の登り部分を線形・非線形曲線として適合する度合い、山型の強度最大値と最小値の差、山型の時系列で前の点から上昇した点の割合のうち少なくとも一つを用いて静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する手順、（５）．前記ＲＦＩＤリーダの送信搬送波とＲＦＩＤタグからの反射で返却された搬送波の位相の時系列分析によって、静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する手順、（６）．ＲＦＩＤタグ応答の位相変曲点またはドップラー周波数０時点と電波受信強度山頂発生時点を比較することによって移動方向を識別する手順、のうち少なくとも一つ以上の手順からなる個別静止ＲＦＩＤタグフィルタを任意の順でコンピュータ装置に実行させる、ことにより解決される。
本発明のＲＦＩＤタグ移動識別プログラムによれば、プログラムが稼働する装置であれば、ＲＦＩＤタグ移動識別させることが可能となる。

本発明によれば、一般的なＲＦＩＤリーダと一般的なアンテナから得られる読取ＲＦＩＤタグの電波受信強度、位相やドップラー周波数の時系列データによって、アンテナ前を移動したＲＦＩＤタグと静止したままのＲＦＩＤタグを識別することが可能となる。
人や物、特に金属が含まれたカゴ車や台車等がアンテナ前を移動することによる電波の乱反射やマルチパスフェージング、特に水分が含まれた物や人がアンテナ前を移動することによって発生する電波遮蔽によって電波受信強度が一時的に下落するような静止したＲＦＩＤタグにおいても、静止ＲＦＩＤタグと識別することができる。また複数の静止ＲＦＩＤタグフィルタを用いることにより、ＲＦＩＤタグの移動方向を識別でき、高精度な移動識別が可能となる。このように、移動判定の精度を高めることができる。

本発明に係るＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムに使用される構成の概略説明図である。 ＲＦＩＤタグ移動識別方法の説明図である。 移動タグ及び静止タグに関する電波受信強度と時間との関係を示すグラフ図である。 秒間読取可能回数と読取中タグ枚数との関係を示す実測データ及び回帰曲線で表したグラフ図である。 アンテナの電波放射特性とＲＦＩＤタグ位置関係を示す説明図である。 移動ＲＦＩＤタグ及び静止ＲＦＩＤタグに関する電波受信強度と時間との関係を説明するグラフ図である。 静止ＲＦＩＤタグに関する電波受信強度と時間との関係を説明するグラフ図である。 静止ＲＦＩＤタグと移動ＲＦＩＤタグの位相と時間との関係を説明するグラフ図である。 静止ＲＦＩＤタグの位相と時間との関係を説明するグラフ図である。 静止ＲＦＩＤタグの位相と時間との関係を説明するグラフ図である。 ＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムに使用される構成をより具体的にした説明図である。 タグ情報取得スレッド（プログラム）と移動タグ識別スレッド（プログラム）のフローチャートである。 カイ二乗検定のフローチャートである。 判定基データ計算のフローチャートである。 移動判定のフローチャートである。 水平方向移動タグのＲＳＳＩ（電波受信強度）、位相、ドップラー周波数のグラフを計測した時系列データを示すグラフ図である。 垂直方向移動タグのＲＳＳＩ、位相、ドップラー周波数のグラフを計測した時系列データを示すグラフ図である。 移動タグの電波受信強度（上グラフ）と位相（下グラフ）を計測した時系列データを示すグラフ図である。 静止タグの電波受信強度（上グラフ）と位相（下グラフ）を計測した時系列データを示すグラフ図である。 位相の振れの数値化例を示す図である。 移動ＲＦＩＤタグの水平方向と垂直方向とを説明する説明図である。 一般的ＲＦＩＤアンテナの電波放射パターンの説明図である。 水平方向移動と垂直方向移動に関する説明図である。 利用データ期間や電波受信強度最大、最小、期間前最大値を説明するグラフ図である。 過読フィルタの説明図である。 仲上−ライス分布の各種Ｋ及びレイリー分布での累積確率分布に関するもので、累積確率と相対的受信電力との関係を示すグラフ図である。 レイリーフェージングシュミレーションであり、相対的受信電力と試行回との関係を示すグラフ図である。 ＲＦＩＤタグの移動種別を4種類に分けた概念図である。 移動をフィルタするための判断基情報やフィルタ名を示す図表である。 アンテナからの距離及び角度とＲＳＳＩの高低関係を示す説明図である。 垂直移動に関する図１３Ｂと同様な電波放射パターンの説明図である。 アンテナに対するＲＦＩＤタグの移動と角度の関係の説明図である。 ＲＦＩＤタグの移動に伴う動径成分速度の変化を示すグラフ図である。 垂直方向移動時、通常のＲＦＩＤアンテナの放射パターンでのＲＳＳＩ山頂点発生位置に関する説明図である。 パラメータ設定画面である。 移動ＲＦＩＤタグ識別パラメータ設定の説明図である。 移動ＲＦＩＤタグ識別パラメータ設定の説明図である。 移動ＲＦＩＤタグ識別パラメータ設定の説明図である。 移動ＲＦＩＤタグ識別パラメータ設定の説明図である。 パラメータ変更画面である。 移動ＲＦＩＤタグ識別パラメータ変更後の画面の説明図である。 移動ＲＦＩＤタグ識別パラメータ設定画面の説明図である。 パラメータ変更画面である。 アンテナが２つの場合の移動ＲＦＩＤタグ識別の設定画面の説明図である。 ＲＦＩＤタグのＥＰＣにより使用パラメータを分ける場合のフローチャートである。 パラメータ変更・設定のフローチャートである。 ＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムに使用される構成の他の例を示す説明図である。 倉庫内盗難品検知への応用例を示す説明図である。 倉庫内棚卸への応用例を示す説明図である。

以下、本発明の一実施形態について、図を参照して説明する。なお、以下に説明する部材、配置等は、本発明を限定するものではなく、本発明の趣旨に沿って各種改変することができることは勿論である。また、本明細書及び特許請求の範囲において、ＲＦＩＤリーダには、ＲＦＩＤリーダライタを含むものである。

図は本発明に係る実施形態であり、図１はＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムに使用される構成の概略説明図、図２はＲＦＩＤタグ移動識別方法の説明図、図３は移動タグ及び静止タグに関する電波受信強度と時間との関係を示すグラフ図、図４は秒間読取可能回数と読取中タグ枚数との関係を示す実測データ及び回帰曲線で表したグラフ図、図５はアンテナの電波放射特性とＲＦＩＤタグ位置関係を示す説明図、図６は移動ＲＦＩＤタグ及び静止ＲＦＩＤタグに関する電波受信強度と時間との関係を説明するグラフ図、図７は静止ＲＦＩＤタグに関する電波受信強度と時間との関係を説明するグラフ図、図８Ａは静止ＲＦＩＤタグと移動ＲＦＩＤタグの位相と時間との関係を説明するグラフ図、図８Ｂは静止ＲＦＩＤタグの位相と時間との関係を説明するグラフ図、図８Ｃは静止ＲＦＩＤタグの位相と時間との関係を説明するグラフ図、図９ＡはＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムに使用される構成をより具体的にした説明図、図９Ｂはタグ情報取得スレッド（プログラム）と移動タグ識別スレッド（プログラム）のフローチャート、図９Ｃはカイ二乗検定のフローチャート、図１０Ａは判定基データ計算のフローチャート、図１０Ｂは移動判定のフローチャート、図１１ＡはＲＳＳＩ、位相、ドップラー周波数のグラフを計測した時系列データを示すグラフ図、図１１ＢはＲＳＳＩ、位相、ドップラー周波数のグラフを計測した時系列データを示すグラフ図、図１２Ａは 位相の変動についての電波受信強度（上グラフ）と位相（下グラフ）を示すグラフ図、図１２Ｂは位相の変動についての電波受信強度（上グラフ）と位相（下グラフ）を示すグラフ図、図１２Ｃは位相の振れの数値化例を示す図、図１３Ａは移動ＲＦＩＤタグの水平方向と垂直方向とを説明する説明図、図１３Ｂは一般的ＲＦＩＤアンテナの電波放射パターンの説明図、図１３Ｃは水平方向移動と垂直方向移動に関する説明図、図１４Ａは利用データ期間や電波受信強度最大、最小、期間前最大値を説明するグラフ図、図１４Ｂは過読フィルタの説明図、図１５Ａは仲上−ライス分布の各種Ｋ及びレイリー分布での累積確率分布に関するもので、累積確率と相対的受信電力との関係を示すグラフ図、図１５Ｂレイリーフェージングシュミレーションであり、相対的受信電力と試行回との関係を示すグラフ図、図１６ＡはＲＦＩＤタグの移動種別を4種類に分けた概念図、図１６Ｂは移動をフィルタするための判断基情報やフィルタ名を示す図表、図１６Ｃはアンテナからの距離及び角度とＲＳＳＩの高低関係を示す説明図、図１６Ｄは垂直移動に関する図１３Ｂと同様な電波分布の説明図、図１７Ａはアンテナに対するＲＦＩＤタグの移動との関係の説明図、図１７Ｂは動径成分速度の変化を示すグラフ図、図１７Ｃは通常のＲＦＩＤアンテナの放射パターンのRSSI山頂点に関する説明図、図１８はパラメータ設定画面、図１９は移動ＲＦＩＤタグの設定の説明図、図２０は移動ＲＦＩＤタグの設定の説明図、図２１Ａは移動ＲＦＩＤタグの設定の説明図、図２１Ｂは移動ＲＦＩＤタグの設定の説明図、図２２はパラメータ変更画面、図２３は移動ＲＦＩＤタグ識別の設定画面の説明図、図２４は移動ＲＦＩＤタグ識別の設定画面の説明図、図２５はパラメータ変更画面、図２６はアンテナが２つの場合の移動ＲＦＩＤタグ識別の設定画面の説明図、図２７はパラメータ変更のフローチャート、図２８はパラメータ変更・設定のフローチャート、図２９はＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムに使用される構成の他の例を示す説明図、図３０は倉庫内盗難品検知への応用例を示す説明図、図３１は倉庫内棚卸への応用例を示す説明図である。

本発明の本質は、ＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムであり、機種に依存されない一般的に用いられているＲＦＩＤリーダ１０と、機種に依存されない一般的に用いられているアンテナ２０によって読取られたＲＦＩＤタグの電波受信強度や位相及びドップラー周波数の時系列データに基づき、アンテナ２０前を移動したＲＦＩＤタグと静止したままのＲＦＩＤタグを識別するＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムにある。ＲＦＩＤリーダ１０は、その周囲に存在するＲＦＩＤタグとの間において無線電波を送受信することによってＲＦＩＤタグを探索してＲＦＩＤタグの情報（データ）を読み取るものである。

本発明に係るＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムは、ＲＦＩＤリーダ１０を制御するコンピュータ装置Ｃ、或いはＲＦＩＤリーダ１０の動作を制御することによって行われる。これを図１及び図２を用いて概念的に説明すると、ＲＦＩＤリーダ１０で読取られた全てのＲＦＩＤタグの情報から静止ＲＦＩＤタグと判断できるものを、静止ＲＦＩＤタグフィルタＦとしての個別静止ＲＦＩＤタグフィルタで除いていくことにより移動ＲＦＩＤタグを識別するものである。表１はＲＦＩＤタグフィルタ一覧である。

公知のコンピュータ装置Ｃとしては、不図示のＣＰＵ及びメモリ（ＲＯＭ，ＲＡＭ，ＨＤＤ等）、Ｗｉ−Ｆｉ（登録商標）、Ｗｉ−Ｆｉ（登録商標）と接続された出入力アンテナと、電源アダプタ、Ｉ／Ｏ（出入力）手段、Ｉ／Ｏ（出入力）手段と接続されたセンサ（光センサ、）、スイッチ等を備えており、メモリに格納されたＲＦＩＤタグ移動識別プログラムに基づき、ＣＰＵによって、送信出力、もしくはアンテナ２０を制御するものである。本実施形態のＲＦＩＤリーダ１０は、コンピュータ装置Ｃと接続されている。

本実施形態のＲＦＩＤリーダ１０或いはコンピュータ装置Ｃは、アンテナ２０に対して、人や物品（商品）に付されたＲＦＩＤタグが移動した場合に、電波受信強度、位相、ドップラー周波数が時系列で刻々と変化することを検知して、人や物品（商品）に付されたＲＦＩＤタグの移動、静止を認識するように制御する。

人や物品（商品）には、その個々を識別可能とするための商品コードや個品ＩＤなどが記憶されているＲＦＩＤタグが取り付けられている。本発明に係るＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムによって、移動ＲＦＩＤタグと識別されたものを、上記人や物品（商品）の移動、購買、管理等に利用できる。なお、商品の棚などの場合には、棚の各段に、その棚位置およびその位置範囲を特定するために、棚の所定位置（角部など）にはＲＦＩＤタグが取り付けられている。ここで、１つの棚（３次元直方体領域）を構成する所定の部分（例えば角部）には、それぞれＲＦＩＤタグが配置されており、このＲＦＩＤタグによって当該棚の位置範囲を特定するようにしている。

一般に、ＲＦＩＤリーダ１０に備えられたアンテナ２０は、静止したＲＦＩＤタグ３２については、図１で示すように、読取りが常時可能であるが、アンテナ２０前を、ＲＦＩＤタグ３０を持った人や、ＲＦＩＤタグ３０が付着した物が通過することにより、図１で示す静止したＲＦＩＤタグ３１の位置にあるような場合は、通過物の陰に隠れ、ある期間（時間）の間、読取れなくなる場合がある。また静止したＲＦＩＤタグ３１で示すような場所の場合は、人や物に反射した電波により短い期間だけ読みとれる場合がある(図３を参照)。これらを、後述する静止したＲＦＩＤタグと判断する静止ＲＦＩＤタグフィルタＦによって、識別して、静止したＲＦＩＤタグを除き、移動しているＲＦＩＤタグを特定するものである。

このとき、静止ＲＦＩＤタグフィルタＦとして、一定期間中の読取回数に閾値を設ける技術がある。ＲＦＩＤタグを読み続けることができる場合のＲＦＩＤタグ１枚あたり秒間読取可能回数は、ＲＦＩＤリーダ１０の無線通信パラメータにより、例えば読取中ＲＦＩＤタグ枚数を独立変数とした多項式関数や累乗曲線関数で計算することができる（図４を参照）。この計算値に一定割合をかけた閾値より秒間読取可能回数が少ない場合に、静止したＲＦＩＤタグと判定することができる。例えばある時点のＲＦＩＤタグ読取枚数が２０枚である場合の閾値を図４で示すようにＡとして設定できる。

また、静止ＲＦＩＤタグフィルタＦとして、ＲＦＩＤタグからの電波受信強度が距離の４乗に反比例する性質、アンテナ２０の電波放射特性、などから移動したＲＦＩＤタグの電波受信強度時系列グラフは通常山型となるので、このことを利用する技術がある。例えば、図５に示すようなアンテナ２０前を移動したＲＦＩＤタグ及び静止したままのＲＦＩＤタグでは、グラフに表すと、図６のように、移動タグのみが山型になるものである。山型全体やその登り部分を線形・非線形曲線との適合度判定、強度最大値と最小値の差、時系列で前の点から上昇した点の割合といった各種算術的・統計的手法により判定することで静止したＲＦＩＤタグを特定することが可能となる。

さらに、静止ＲＦＩＤタグフィルタＦとして、人や物、特に金属が含まれたカゴ車や台車等がアンテナ２０前を移動することによる電波の乱反射やマルチパスフェージング、前記した図１で示したような静止したＲＦＩＤタグ３１のような場所の場合のように、物陰に隠れ電波が遮蔽され受信強度が一時的に下落するときには、静止したＲＦＩＤタグにおいても、図７で示すように、例えば時間軸点６以降のような山型のグラフとなることがある。この場合、前の時間軸データまで遡り電波受信強度が高い点があるかを判定すること等により図７のような谷型を検出する技術により、静止したＲＦＩＤタグと識別することが可能となる。

さらに、ＲＦＩＤリーダ１０の送信搬送波とＲＦＩＤタグからの反射で返却された搬送波をＲＦＩＤリーダ１０内で合成した強度とＲＦＩＤタグの電波受信強度を比較するなどの技術により、静止したＲＦＩＤタグを特定することが可能である。つまり、２つの搬送波の位相差を０〜１８０度あるいは０〜３６０度の範囲で検知することが可能である。この機能は、ＲＦＩＤリーダ１０の機能として備えられている場合があり、これを利用することが可能である。一般に、アンテナ２０とＲＦＩＤタグとの間の距離が１／４波長または１／２波長変化する毎に位相が１周する。例えば、周波数９５３ＭＨｚの場合で１／４波長は７．８６ｃｍ、１／２波長は１５．７２ｃｍである。

静止したままのＲＦＩＤタグでは位相が変化せず、移動したＲＦＩＤタグは変化を続けることになる。前記周波数帯のように波長が人や物の移動距離に比較し小さい場合、人が手を振る等の少しの移動によって位相は大きく変動し、移動時の位相は、図８Ａの概念グラフのように連続的に変化し一様に分布することになる。
図８Ａで示す２つの時系列データの差異を標準偏差の値に閾値を設けることや二項検定、カイ二乗検定、Ｇ検定等の手法により静止したＲＦＩＤタグを識別することが可能となる。

このとき、標準偏差はデータ数が少ない場合に、変動が大きくなり正しい識別を行えないことがある。二項検定やカイ二乗検定、Ｇ検定は、データ数が少ない場合に発生する識別精度の曖昧さを確率として表現できるため、移動時の読取回数が環境により数回〜数１００回程度と大きく変動するＲＦＩＤタグ読取の場合により有効な方法である。また二項検定、カイ二乗検定、Ｇ検定等を用いた技術により、ＲＦＩＤタグを識別することができる。

図９ＡはＲＦＩＤタグ移動識別方法及びＲＦＩＤタグ移動識別プログラムに使用される構成をより具体的にした説明図であり、コンピュータ装置（ＰＣ）ＣとＲＦＩＤリーダ１０とから構成され、これらはＬＬＲＰ（低レベルのリーダプロトコル）でインターネットやイントラネットで使用されるプロトコル（ＴＣＰ／ＩＰ）を利用して接続されている。これらのコンピュータ装置（ＰＣ）ＣとＲＦＩＤリーダ１０は汎用のものであり、本発明のために特別形成されたものではない。またＲＦＩＤリーダ１０は、汎用入出力（ＧＰＩＯ）によってセンサやスイッチ、パトライト（登録商標）等と接続することができる。

コンピュータ装置Ｃには、識別パラメータ設定ツール、識別パラメータ保存ファイル、上位に位置する上位アプリケーション（上位アプリ：ＲＦＩＤタグデータ表示や保存を行う）、ミドルウエア（ナイスミドル）、ＲＦＩＤタグ識別エンジン等が、ＯＳ上で稼働するように構成されている。図９Ａの例では、ＯＳとしてＷｉｎｄｏｗｓＸＰ（登録商標）、Ｗｉｎｄｏｗｓ７（登録商標）を用いた例を示している。

図９Ｂは、タグ情報取得スレッド（プログラム）と移動タグ識別スレッド（プログラム）のフローチャートである。
本実施形態において、タグ情報取得はステップＳ１０〜Ｓ１５で行われるものであり、コンピュータ装置ＣとＲＦＩＤリーダ１０との接続処理（Ｓ１０）を行い、次に、アンテナ２０から電波を送信し、ＲＦＩＤタグから情報の読取を開始する（ステップＳ１１）。ステップＳ１２〜ステップＳ１５はループとなっており、ステップＳ１２で、読取ったＲＦＩＤタグの情報を取得し、読取ったＲＦＩＤタグの情報をメモリ（記憶媒体）であるハードディスク等に保存する。またステップＳ１４ではメモリ（記憶媒体）であるＲＡＭ等に保持する。このループによりＲＦＩＤタグの枚数とそれぞれの読取回数を掛け合わせた数の読取情報が取得される。ここで、ＲＦＩＤタグ情報は、読取時刻履歴、読取アンテナ履歴、電波受信強度履歴、位相履歴、ドップラー周波数履歴等の情報を取得するものである。

これらのＲＦＩＤタグ情報は、次の移動タグ識別におけるステップＳ２５の判定基データ計算で利用される。
本実施形態において、移動タグ識別はステップＳ２１〜Ｓ３２でおこなわれるもので、ステップＳ２２〜Ｓ３２、ステップＳ２３〜Ｓ３１、ステップＳ２４〜Ｓ２９は、それぞれループとなっており、ステップＳ２１で識別パラメータ設定ファイルを読取り各種パラメータをＲＡＭ等のメモリに格納する処理を行う。

ステップＳ２３では個別ＲＦＩＤタグの移動判定処理を行う。次に、ステップＳ２４で利用データ期間中の移動判定処理を行う。次にステップＳ２５で判定基データの計算処理を行う。このとき、秒間読取回数、最大電波受信強度、最大−最小電波受信強度、電波受信強度上昇割合、期間中最大 - 期間前最大電波受信強度、位相変動カイ二乗片側確率、ドップラー周波数平均値などが計算される。このときのステップＳ２５の処理は後述する。

このようにしてステップＳ２５の処理の後で、ステップＳ２６で移動判定処理を行う。この時のステップＳ２６の処理は後述する。次に、移動ＲＦＩＤタグが移動かどうかを判断する（ステップＳ２７）。移動ＲＦＩＤタグの場合（ステップＳ２７：Ｙｅｓ）、移動タグとして処理を行う（ステップＳ２８）。これにより画面表示、上位アプリケーション・サーバへの送信処理等を行う。ステップＳ２８の後で、判定対象タグがまだある場合にはステップＳ２３〜Ｓ３１のループ処理を継続する。移動ＲＦＩＤタグでない場合（ステップＳ２７：Ｎｏ）、利用データ期間を０．１秒増やして移動判定処理を行う（ステップＳ２４〜Ｓ２９）。移動と判定できない場合には利用データ期間上限に到達するまでステップＳ２４〜Ｓ２９のループを継続する。静止していると判定できる場合には、ステップＳ３０で静止ＲＦＩＤタグとして処理する。判定対象タグがまだある場合にはステップＳ２３〜Ｓ３１のループ処理を継続する。以上のようにして、ステップＳ２２からステップＳ３２までを繰り返す。

ステップＳ２５の判定基データ計算の詳細フローチャートが図１０Ａである。計算の際に用いられる各種パラメータ（山頂期間割合、ドップラー周波数取得期間割合、期間前データ取得期間）はステップＳ２１（図９Ｂ）において設定ファイルからメモリ内に読みこまれたものである。

ステップＳ７０において計算の過程で使用される各種変数を初期化する。数値のカウントや合計を格納する変数は初期値を０とする。ＲＳＳＩの取り得る値の範囲は−１００〜３０程度であるため最大値を格納する変数には―９９９を、最小値を格納する変数には９９９を初期値として格納しておくことにより、後の比較演算の初回は必ずＹｅｓとなり処理を単純化できる。

ステップＳ７１でステップＳ１４でメモリ内に保持された履歴データのうち、対象タグ／対象期間のデータＮ個を計算対象とする。対象期間の全データを走査するためステップＳ７２で変数ｉにＮ−１を初期値として格納し、ステップＳ７３〜Ｓ８４のループ処理でＮ−１から０までｉを１ずつ減算しながら計算を行い、ｉがマイナスになった時点でループを抜ける。

ステップＳ７４において、山頂期間割合パラメータを基にＲｅａｄ＿ｔｉｍｅ［ｉ］が山頂期間内のデータであるか判定を行う。山頂期間内のデータであれば（ステップＳ７４：Ｙｅｓ）、ＲＳＳＩ［ｉ］を用い期間中最大電波受信強度更新処理（ステップＳ７５、Ｓ７６）を行い、山頂期間内のデータでなければ（ステップＳ７４：Ｎｏ）、ステップＳ７７、Ｓ７８で期間中最小電波受信強度更新処理を行う。

ステップＳ７９でＲＳＳＩ［ｉ］とＲＳＳＩ［ｉ−１］を比較し、過去の電波受信強度の方が低いか等しい場合（ステップＳ７９：Ｙｅｓ）、ステップＳ８０で電波受信強度上昇カウンタを１増やす。なお、この比較はｉ＝０の場合には行わない。ステップＳ８１において、ドップラー周波数取得期間割合パラメータを基にＲｅａｄ＿ｔｉｍｅ［ｉ］がドップラー周波数取得期間内のデータか判定する。判定Ｙｅｓの場合には、後に平均値を計算するために値を足し算すると共に、カウンタを１増やす（ステップＳ８２）。ステップＳ８３において、ループ終了判定のためにｉを１減算するとともにデータカウンタを１増やす。ステップＳ８５において、期間前データ取得期間パラメータ分データを遡りその間の電波受信強度最大値を取得する。ステップＳ８６において各種の移動判定基データが計算される。このうち、位相変動カイ二乗片側確率の計算方法は後述する。

ステップＳ２６（図９Ｂ）の移動判定の詳細フローチャートが図１０Ｂである。判定基データＤ１〜Ｄ７はステップＳ２５で計算された値であり、識別パラメータＰ１〜Ｐ９はステップＳ２１において設定ファイルからメモリ内に読みこまれたものである。

ステップＳ６０はチラ読みフィルタの処理であり、読取中ＲＦＩＤタグ枚数を独立変数とした関数で計算された秒間読取可能回数とＰ１との積とＤ１とを比較する処理である。比較判定がＮｏの場合には、静止タグと判定され以後の処理は行わない。Ｙｅｓの場合には次の判定に進む。この分岐は以後の処理で共通である。

ステップＳ６１はＲＳＳＩ絶対値フィルタの処理である。ステップＳ６２及びステップＳ６３はＲＳＳＩ非連続上昇フィルタの処理である。ステップＳ６４はＲＳＳＩ谷型フィルタの処理である。ステップＳ６５は位相振れ少フィルタの処理である。ステップＳ６６で垂直方向移動フィルタ使用の有無を判定し、判定Ｙｅｓの場合には垂直方向移動フィルタの処理であるステップＳ６７を実行する。ステップＳ６７まで全てＹＥＳの場合移動タグと判定される。またステップＳ６６の判定がＮｏの場合にはステップＳ６７を実行せずに移動タグと判定される。
なお、垂直方向移動とは、図１３Ａに図示されるように、アンテナの電波放射面に垂直な方向、水平方向移動とは電波放射面に水平な方向への移動を意味する。

上記した二項検定は一様分布である確率をかなり正確に求めることができるが、データ数が多い場合に計算量が膨大となってしまう欠点がある。カイ二乗検定やＧ検定は計算量が少なく二項検定を高い精度で近似できるため有用な手法である。しかしいずれの手法とも図８Ｂのように人等の移動を契機に位相の値が前半・後半で２つのグループに分かれた場合、一様分布と判定され移動タグと誤識別されてしまう場合がある。

上記状況のような場合においても、期間中のデータを前半・後半に分けそれぞれで観測頻度を計算することにより、静止タグと識別できるように構成することが可能である。図９Ｃは、このように、期間中のデータを前半・後半に分けて観測頻度の計算を行う場合を説明するカイ二乗検定のフローチャートであり、図９Ｃのフローチャートに基づいて説明する。このフローチャートの処理は、それぞれのＲＦＩＤタグの利用データ期間毎に行うものである。

図９Ｃのフローチャートは、図９Ｂにおける移動タグ識別スレッド内の位相変動カイ二乗片側確率を求めるための「判定基データ計算」（ステップＳ２５）の処理に相当する。
この「判定基データ計算」の処理は次のように行われる。
ステップＳ４０において対象タグの対象期間中Ｍ個の位相履歴データをメモリから取得する処理を行う。このデータは図９Ｂにおけるタグ情報取得スレッド（プログラム）処理のステップＳ１４で示す「読取タグ情報メモリ保持」処理によって既に得られたものである。以下、例えば、位相が０〜１８０°の範囲のデータが、ＲＦＩＤリーダ１０から返却されるケースに基づいて説明をする。位相履歴データが配列Phase_angle[0]〜Phase_angle[M-1]，対応する読取時刻履歴データが配列Read_time[0]〜Read_time[M-1]に格納されているものとする。なお以降で角度はラジアンではなく［°］と統一して記述する。Ｃｏｓ，Ｓｉｎの引数やＡｒｃｔａｎの戻り値も［°］で計算されるものとする。

ステップＳ４１〜Ｓ４５は、後に位相平均値を求めるための前処理として位相ベクトル和を計算するためのループである。位相は１８０°の周期データあるためその平均値を算術平均で求めると好ましくない場合がある。例えば静止ＲＦＩＤタグが風等により微小に振れた場合やＲＦＩＤリーダ１０内部のノイズ等によって図８Ｃのように０°（＝１８０°）を挟んだ時系列データとなった場合、その平均値を算術平均で得られる９０°としてしまうと変動の大きいデータと見なされ移動ＲＦＩＤタグと誤認識してしまう場合がある。この場合の平均は０°とした方が望ましい。まず、期間中の前半か後半かで別の平均値とするためにステップＳ４２で期間中の前半かどうか判断する。前半である場合にはステップＳ４３の処理を行い、後半である場合にはステップＳ４４の処理を行う。これらは０〜１８０°の元データを２倍し０〜３６０°に変換してその角度をｘｙ平面上に持つ単位ベクトル和を計算する処理である。前半と後半で別の変数を設け、それぞれ以下式で計算する。ステップＳ４３は
Sum_unit_vector_x_anterior
= Sum_unit_vector_x_anterior + Cos(Phase_angle[x]*2)
Sum_unit_vector_y_anterior
= Sum_unit_vector_y_anterior + Sin(Phase_angle[x]*2)
ステップＳ４４は
Sum_unit_vector_x_posterior
= Sum_unit_vector_x_posterior + Cos(Phase_angle[x]*2)
Sum_unit_vector_y_posterior
= Sum_unit_vector_y_posterior + Sin(Phase_angle[x]*2)
である。

ステップＳ４６において、先に計算された単位ベクトル和のｘ，ｙ座標からその角度を以下式で計算することにより期間前半・後半それぞれの位相平均値を求める。
Average_phase_angle_anterior=Atan2(Sum_unit_vector_y_anterior,
Sum_unit_vector_x_anterior)
Average_phase_angle_posterior=Atan2(Sum_unit_vector_y_posterior,
Sum_unit_vector_x_posterior)
ここでAtan2(y,x)は以下を表す関数である。
x>=0の場合 Arctan(y/x)
x<=0の場合 Sign(y)180+Arctan(y/x)
(Sign(y)はyの符号)

ステップＳ４７〜Ｓ５５は後にカイ二乗値を計算するために各位相を、平均値からの差が９０°未満であるデータとそれ以外（９０〜１８０°）に分けてそれぞれの観測頻度をカウントするためのループである。ステップＳ４８では、期間中の前半か後半かで平均値について別の値を用いるために、期間中の前半かどうか判断する。前半である場合にはステップＳ４９で変数Average_paに前半位相平均値を格納し、後半である場合にはステップＳ５０で変数Average_paに後半位相平均値を格納する。そして、例えば１０°と１７０°の角度差は１６０°ではなく最短距離の２０°として計算するため、ステップＳ５１で以下式を使用する。
Diff_angle =|Average_pa -
Atan2(Sin(Average_pa)+Sin(Phase_angle[x]*2),
Cos(Average_pa)+Cos(Phase_angle[x]*2))|

ステップＳ５１では、角度差が９０°より小さいかどうかを判断する。９０°より小さい場合（ステップＳ５２：Ｙｅｓ）、ステップＳ５３の処理を行い、９０°より大きい場合（ステップＳ５２：Ｎｏ）、ステップＳ５４の処理を行う。これにより、Ｍ個の点を２つのグループに分けそれぞれの観測頻度をカウントする（ステップＳ５３、Ｓ５４）。

ステップＳ４７〜Ｓ５５のループにおいて、もしデータが一様分布であれば、それぞれのグループの発生頻度期待値はM／２となり、ステップＳ５６でその値を計算する処理を行う。
次に、ステップＳ５７で期待値と各観測頻度との差の二乗和を期待値で割ったカイ二乗値（χ２）を以下式で求める。
χ2= ( (Observed[0] − Expected)^2 + (Observed[1] −
Expected)^2 )/Expected
さらに、ステップＳ５８でプログラム起動時にメモリ内配列Chi_probability[]に読みとっていた「自由度１のカイ二乗分布表」からχ２に紐付けられた値を求める。配列Chi_probability[]の添え字は0以上の整数であるためχ2の値も整数に変換する。これが判定基データの「位相変動カイ二乗片側確率」となり、位相時系列データが期間前半にも後半にも偏りのない一様分布である確率を表すことになり、移動ＲＦＩＤタグの場合は、この値が高くなる。

つまり、対象期間中のＭ個のデータの位相平均値を求める処理を行い、平均からの差により全データをＮ個のグループに分ける処理を行う。上記の例では角度差により２個に分けているが、それ以外のＮ個のグループに分けることができる。一様分布のデータであれば全てのグループの頻度の期待値はＭ／Ｎである。さらに、各グループの出現頻度からカイ二乗統計量（χ2）を計算する処理を行う。そして、自由度Ｎ−１でカイ二乗統計量（χ2）の値が発生する片側確率（Ｐ）を求める処理を行う。このとき、プログラム動作時の処理を軽減するために予め計算された代表的な値をメモリに格納しておくことが可能である。そして、この値は、位相データが一様分布である確率とみなすことができる。

次に、予め設定された閾値と、求めた片側確率（Ｐ）とを比較して、片側確率が閾値より大きいかどうか判定する。片側確率（Ｐ）が閾値より小さい場合、静止ＲＦＩＤタグと判断し、大きい場合、移動ＲＦＩＤタグと判断する。

このようにして、移動か静止かの識別を行うことができる。また、実際には移動していない場合に、移動と判定する誤検知をなくすことが必要な防犯ゲートのような応用では、閾値の値を高めに設定することで誤検知を減らすことができる。

ところで、前述した位相を利用した技術だけでは、本来検知したくないアンテナ２０から遠く離れた場所や、アンテナ２０の横や裏での移動についても検知してしまう可能性が残ってしまう。また電波受信強度を利用した技術だけでは、人や物が長時間電波を遮蔽した後に移動することにより移動誤検知が発生する可能性が残ってしまう。そこで、誤検知の可能性を少なくするために、秒間読取可能回数を利用した技術を含めた、複数の技術によって、個別静止ＲＦＩＤタグフィルタを構成させて、これらを組み合わせ、或いは全ての技術から構成される個別静止ＲＦＩＤタグフィルタを用いることにより、ＲＦＩＤタグの高精度な移動識別が可能となる。

図１２Ｂは、図１で示す静止した場所にあるＲＦＩＤタグ３１のような場合を想定し、実際にＲＦＩＤタグ３１を配置し、ＲＦＩＤタグを有しない人がアンテナ２０前に長時間立ち電波遮蔽を行った後に離れた際におけるタグ３１の電波受信強度（上グラフ）及び位相（下グラフ：０〜１８０度。）を計測した時系列データを示すグラフ図である。上グラフの直線で囲まれた期間Ｔ１において、一定以上の秒間読取回数での電波受信強度山型（山の登り）を検知したが、その期間Ｔ１中や期間前後の位相の値は一様に分布していないため、このように位相の値が一様に分布していない場合には、静止したＲＦＩＤタグであると識別することが可能である。

上記位相の変化の傾きは、アンテナ２０に接近している場合と、離れている場合で逆向きとなる。この変曲点発生時と電波受信強度の山頂発生時点を比較することで、ＲＦＩＤタグがアンテナ２０と水平方向に移動したか、垂直方向に移動したかを識別することが可能となる。このため、セキュリティゲート等の応用において重要となるゲート前を通過したＲＦＩＤタグだけを検知することが可能となる。

図１１Ａは人がタグの貼付されたシャツを手に持ちアンテナ近くを水平方向に歩行した実際の時系列データであり、上からＲＳＳＩ、位相、ドップラー周波数のグラフである。x軸の単位は秒であり、左に行くほど過去のデータを示す。ＲＳＳＩグラフ上の2本の垂直線で囲まれた期間(−２．２〜−１．３秒)のデータで移動と識別した。図の−１．５秒付近においてＲＳＳＩ山頂点が発生しているがこの近辺で位相の変曲点とドップラー周波数0点が発生している。

位相は式（５）で示されるようにタグとリーダアンテナ動径方向距離の関数でありアンテナに近づいている時と離れている時ではその変化の方向が反転し、速度0[m/s]地点でグラフに変曲点が発生する。位相変化の情報を得るためには２点の読取データが必要であるが、その読取点間に7.86cm以上の移動(周波数９５３ＭＨｚで位相0〜180°の場合)が発生した場合にはその変化の方向は２通りあり識別できない。同時読取タグ枚数やタグ応答データ量が多い場合のように読取点間隔が空いてしまうと正しい変化の方向を検知できない場合がある。

一方式（６）で表されるドップラー周波数は、１点の読み取りデータで情報を得ることができるため高速に移動した場合でもその符号によって変化の方向を検知できる。式（６）を変形した式（７）によってタグの動径方向移動速度が得られる。

図１１Ａ及び図１１Ｂのドップラー周波数グラフ縦軸の単位は［Ｈｚ］であるが、値に１．１をかけることでほぼ時速［ｋｍ／ｈ］に換算できる。

図１１Ｂは図１１Ａと同環境でアンテナ裏面→表面方向に歩行を行ったデータである。位相のグラフの−１．０秒付近において180°近くから0°近くの点に遷移し下向きの線が引かれているが、位相は周期データであるためこの期間の変化の傾きは実際には正である。その後の−０．９秒付近に固まった２点の傾きは負であることから位相変曲点は−０．９５秒近辺で発生しておりその付近のドップラー周波数は０に近くなっている。−０．３秒付近においてRSSI山頂点が発生しているがこの近辺で位相の傾きは負、ドップラー周波数はマイナスの値になっている。ドップラー周波数０点からＲＳＳＩ山頂点発生までには０．６５秒程度の時間差が発生したことになる。

以上から、ＲＳＳＩ山頂点発生と位相変曲点またはドップラー周波数０点発生の時間差が生じた場合に垂直方向移動と識別可能であることがわかる。なお、図１１Ａ、図１１Ｂのドップラー周波数グラフは揺れが発生しているが、これは人の手の振れや上下方向の揺れ等が影響していると考えられる。この影響を考慮し、図９Ｂで示す判定基データ計算では１点のみの値ではなくある程度の期間中でのドップラー周波数平均値を計算している。

図１２Ａは移動タグの電波受信強度（上グラフ）と位相（下グラフ）を示すグラフ図であり、移動中のＲＦＩＤタグの位相は、通常この図のように変動し続ける。また図１２Ｂのように静止した状態のＲＦＩＤタグは、この図の下グラフの○部分のように一定の値を取り続ける。アンテナ前に人が長時間とどまったあと、移動したために、人の後ろに隠れていた静止タグに電波受信強度変化が発生（直線で囲まれた部分）し移動タグと誤認識された例であり、この期間の位相の振れは上の本当の移動ＲＦＩＤタグに比べて低いことがわかる。

なお、図１２Ｃは、位相の振れの数値化例を示す図で、この図１２Ｃのように、０〜１００の間で数値化される。位相振れ少フィルタの閾値に低い値を設定すると、より多くのＲＦＩＤタグが移動ＲＦＩＤタグと識別される。０と設定することで全てのＲＦＩＤタグが移動候補ＲＦＩＤタグになり（フィルタ無）、１００とすることで全てが静止ＲＦＩＤタグと識別される。デフォルト値は５であるが、識別パラメータ設定ツールにて適宜調整可能である。

以上のように、各種静止ＲＦＩＤタグフィルタＦは、その使用の有無、閾値等が、ＲＦＩＤリーダ１０の電波出力や性能、使用アンテナ２０、ＲＦＩＤタグの性能、周囲の静止ＲＦＩＤタグの想定位置、ＲＦＩＤタグの枚数等の環境により最適なパラメータが異なることになる。このような最適なパラメータが異なるために、実際の読取環境において、予め、確認のために、ＲＦＩＤタグを移動し、その挙動やデータの収集等を行い、これらを基に環境の差が吸収された最適なパラメータを求めることが好ましい。

本発明の応用用途によってはＲＦＩＤタグのアンテナからの距離や移動方向によって移動タグも更にフィルタすることが必要な場合がある。例えば入出荷ゲートが複数存在する倉庫の入出荷検品で本発明が用いられる場合、アンテナが設置されたゲートでのタグ移動のみを検知し隣接する離れたゲートでの移動は検知したくないことがある。また、盗難防止ゲートに本発明が使用される場合にはアンテナ前を水平方向に移動した盗難可能性品のみを検知し店内での買い物客の移動に相当する垂直方向移動は検知したくない場合が多い。図１６Ａが移動種別を４種類に分けた概念図で、図１６Ｂがそれぞれの移動をフィルタするための判断基情報やフィルタ名である。

アンテナ遠くの水平方向移動の場合には、図１６Ｃ上部の破線矢印のようにアンテナ放射パターン上はＲＳＳＩが高い所のみでの移動となり、またアンテナ動径方向の距離減少が小さいことから「ＲＳＳＩ非連続上昇フィルタ」により排除できることが多い。また、フリスの公式から距離が離れた場所でのＲＳＳＩは低くなり、レイリー分布や仲上−ライス分布で示されるようにＲＳＳＩ上振れはほとんど発生しないことから「ＲＳＳＩ絶対値フィルタ」も有効となる。さらに、遠くのタグには十分な電力が供給されにくく読みにくいことから「チラ読みフィルタ」も有用となる場合がある。逆に、アンテナ遠くの移動も検知したい場合にはこれらフィルタの閾値を甘く設定する。

図１６ＤのＣ地点からＤ地点のようなアンテナ表面→裏面方向の垂直方向移動の場合、距離が１／２になることによりＲＳＳＩが＋１２ｄＢとなるがアンテナ放射パターンから−２０ｄＢとなって、合計−８ｄＢとなり「ＲＳＳＩ非連続上昇フィルタ」でフィルタ可能である。
一方、図１６ＤのＥ地点からＦ地点のようなアンテナ裏面→表面方向の垂直方向移動の場合では、アンテナ放射パターンよりＲＳＳＩ上昇は＋３０ｄＢ以上にもなるため「ＲＳＳＩ非連続上昇フィルタ」等でフィルタすることは難しい。そこで垂直方向移動の場合にはＲＳＳＩ山頂点発生時に動径成分速度が０とならないことを利用しフィルタを行う。以下その方法について説明する。

図１７ＡのようにタグがアンテナからＡ［ｍ］地点を速度ｖ［ｍ／ｓ］で水平方向に移動する場合、アンテナ中心に向かう動径方向（図の点線）成分の速度はｖ・ｃｏｓθで与えられ、θ＝９０°となるｘ＝０［ｍ］地点では０［ｍ／ｓ］になる。図１７Ｂは速度１．１１［ｍ／ｓ］ （時速４ｋｍ）で、図１７Ａのｘを−１．５［ｍ］〜１．５［ｍ］の範囲で移動した場合の動径成分速度の変化である。アンテナに向かう方向の速度をプラス、離れる方向の速度をマイナスとした。図１７Ａ中の×点（ｘ＝０）でアンテナ放射パターン及び直線距離上からＲＳＳＩは最大となるがこの点での動径成分速度は０となる。

一方、図１７Ｃで示される垂直方向移動の場合、水平方向移動時と同様にｘ軸上の点で動径成分速度が０［ｍ／ｓ］となるが、アンテナ裏面や側面への電波放射は小さいという通常のＲＦＩＤアンテナの放射パターンによってＲＳＳＩ山頂点はｘ軸上ではなく通常×点のような場所となる。その時点での動径成分速度は表面→裏面移動の場合はプラス、裏面→表面移動の場合はマイナスとなる。

次に、より具体的に説明する。
コンピュータ装置Ｃは、公知の構成からなるもので、ＣＰＵ，記憶装置（外部記憶装置を含むもので、ＨＤＤ，ＲＡＭ，ＲＯＭ、ＣＤ，ＤＶＤ，ＢＬ，フラッシュメモリ等を含むその他のメモリ等）、ディスプレイ、入力部（キーボード、マウス、タッチパネル）、出力部（電波、音波、光その他の通信手段との間の出力関係を含む）等を備えている。

本実施形態では、図１で示すように、ＲＦＩＤリーダ１０と、コンピュータ装置Ｃは接続されている。そして、アンテナ２０で信号を受信し、コンピュータ装置Ｃは、ＲＦＩＤリーダ１０からの情報を受け取り、ＲＦＩＤリーダ１０を制御するように構成されている。なお、この実施形態では、コンピュータ装置ＣによってＲＦＩＤリーダ１０を制御するように構成した例を示しているが、ＲＦＩＤリーダ１０にコンピュータ装置Ｃと同様な構成を搭載することにより、ＲＦＩＤリーダ１０のみで、制御するように構成することができるのは勿論である。

本実施形態におけるＲＦＩＤタグ移動識別プログラムは、以下の３つのモジュール（ソフトウェア：プログラム）から構成されている。
１．識別パラメータ設定ツール
２．移動タグ識別エンジン
３．ＬＬＲＰミドルウエア「ナイスミドル」（商標：以下同じ）
これらのモジュール（ソフトウエア：プログラム）は、ＯＳ上で動作するものであり、本実施形態ではＷｉｎｄｏｗｓＸＰ（登録商標）或いはＷｉｎｄｏｗｓ７（登録商標）で動作するようになっている。

そして、識別パラメータ設定ツールで環境を設定し、次に移動タグ識別エンジンで識別し、ＬＬＲＰミドルウエア「ナイスミドル」で、移動タグ識別エンジンのフロントエンドとして上位アプリとのインターフェース（Ｉ／Ｆ）機能を提供する。これにより、それぞれの目的ごとの機能のみを搭載したアプリケーションを開発するようになっている。アプリケーションとしては、例えば、ＲＦＩＤタグを用いた入出荷検品、在庫管理、入退出管理、盗難品検知、観客管理、商品管理、商品やサービスの販売管理、その他の多くのＲＦＩＤタグの移動によって生じる各種管理等を行うアプリケーションを構築することが可能となる。

まず読取環境（アンテナ２０の構成・使用するＲＦＩＤタグ・移動速度・電波出力等）により最適な識別方法やパラメータが異なるため、実環境において最適となるようにパラメータを設定するが、上記の「識別パラメータ設定ツール」を起動する。
この「識別パラメータ設定ツール」により、所定のＲＦＩＤタグを移動させた後、その環境における判定パラメータを計算し設定することができる。設定内容は「識別パラメータ保存ファイル」に保存され、「移動タグ識別エンジン」において使用されるようになる。なお、パラメータによっては「識別パラメータ設定ツール」でのみ使用され、「移動タグ識別エンジン」で使用されないものも存在する（リーダＩＰアドレスや送信出力等、設定画面で青字の項目）。これらは上位アプリからＬＬＲＰミドルウエア（ナイスミドル）に設定するように構成される。

次に、識別パラメータ設定ツールの使用方法について、画面遷移に基づいて説明する。
まず、PickMovedTagsParamsAutoConfigTool.exeをダブルクリックして識別パラメータ設定ツールを起動する。図１８で示す設定画面において、「移動識別パラメータ」タブで接続するＲＦＩＤリーダ１０の設定を行う。
この設定画面では、上記ＲＦＩＤリーダ１０設定の他に、識別オプション、ログ出力、移動識別パラメータ、識別対象ＥＰＣ、ＥＰＣ別の移動識別パラメータ等を設定する。ここで、ＥＰＣ（Ｅｌｅｃｔｒｉｃ Ｐｒｏｄｕｃｔ Ｃｏｄｅ）とは、ＩＣタグを利用して製品識別を行う場合に利用されるコードであり、商品識別コード、輸送用のコンテナなどの識別コード、事業所コード、返却可能な資産の識別コード、企業内の資産管理コードなどが組み込まれているものである。

ＲＦＩＤリーダ設定では、リーダＩＰアドレス、送信出力（電波出力）、無線チャンネル、アンテナ番号等を特定することができる。
識別オプション（移動識別）では、垂直方向移動フィルタ使用の有無、移動方向表示の有無、アンテナグループを複数設定可能となっており、各グループのアンテナ（本実施形態ではグループＡ，Ｂともにアンテナ１〜４）を指定可能となっている。
ログ出力では、移動識別サマリログ、リプレイ用ログ、使用パラメータバックアップが選択可能となっている。
「移動識別パラメータ」では、設定により移動識別対象となるＥＰＣを絞ることが可能である。対象外のＥＰＣではＯｎＭｏｖｅｄＴａｇＲｅｐｏｒｔイベントが発生しなくなる。

移動識別パラメータでは、利用データ期間、秒間読取回数必要割合（％）、必要ＲＳＳＩ（ｄＢｍ）、必要ＲＳＳＩ差（ｄＢ）、ＲＳＳＩ増加発生必要割合（％）、期間前最大ＲＳＳＩからの必要増分（ｄＢ）、必要位相振れ等の設定が可能となっている。

識別対象ＥＰＣとＥＰＣ別パラメータ設定について説明すると、識別対象ＥＰＣは、記入欄が設けられ、この記入欄が空欄の場合には、全ＥＰＣが対象となる。記入の場合には、ＥＰＣ先頭を１６進数で記述し、複数指定するときは改行して記入する。また、ＥＰＣ毎に別の移動識別パラメータを、設定で４つまで設定可能となっている。ＥＰＣ別の移動識別パラメータは、ＥＰＣ先頭、必要ＲＳＳＩ（ｄＢｍ）、必要ＲＳＳＩ差（ｄＢ）、期間前最大ＲＳＳＩからの必要増分（ｄＢ）の記入欄が設けられている。

識別対象ＥＰＣ及び使用されるパラメータは、図２７で示すフローチャートのようにして決定される。すなわち、ステップＳ１００でタグ読取（ＥＰＣ＝Ｘ）を行う。次に、ステップＳ１０１で、識別対象ＥＰＣが空欄かどうか判断する。空欄である場合（ステップＳ１０１：Ｙｅｓ）、ステップＳ１０２で、ＥＰＣ別の識別パラメータで、ＥＰＣ先頭の何れかを判断する。先頭の何れかである場合（ステップＳ１０２：Ｙｅｓ）、ステップＳ１０３でＥＰＣ別の移動識別パラメータの該当パラメータを使用し移動識別実施処理を行う。そして、ステップＳ１００へ戻り、次の読取タグの処理を同様に行う。空欄でない場合（ステップＳ１０１：Ｎｏ）、ステップＳ１０４で識別対象ＥＰＣの何れかかどうか判断する。識別対象ＥＰＣの何れかでない場合（ステップＳ１０４：Ｎｏ）、ステップＳ１００へ戻り、次の読取タグの処理を同様に行う。

識別対象ＥＰＣの何れかである場合（ステップＳ１０４：Ｙｅｓ）、ステップＳ１０２の判断を行う。ステップＳ１０２においてＥＰＣ別の識別パラメータで、ＥＰＣ先頭の何れかを判断するが、先頭の何れかでない場合（ステップＳ１０２：Ｎｏ）ステップＳ１０５で、移動タグ識別パラメータを使用し移動識別実施処理を行う。そして、ステップＳ１００へ戻り、次の読取タグの処理を同様に行う。

前記した図１８で示す設定画面で、パラメータを設定した場合には、保存ボタンを押下し、識別パラメータ保存ファイルに保存する。
次に、図１９で示す移動タグ識別をする画面では、上下及び二段に分割された画面で、上段側には、移動タグの欄と、静止タグの欄、下段には、ＲＦＩＤタグからのデータの受信記録の時系列グラフが表示される。また、この画面には、「オープン／クローズ」ボタン、「スタート／ストップ」ボタン、「クリア」ボタン、「リプレイ」ボタンが設けられている。「リプレイ」ボタンは、リプレイ機能を行うもので、パラメータ変更後の動作確認（設定内容確認）を行うことができるものである。

そして、「オープン」ボタンを押下後、「スタート」ボタンを押下することで読取が開始され、現在のパラメータで移動と識別されたＲＦＩＤタグが「移動タグ」欄に、静止と識別されたものが「静止タグ」欄に表示される。移動タグは利用データ期間の開始・終了点に垂直な棒がグラフ表示される。この実施形態では、画面に表示されるＲＦＩＤタグの総数は２０４８個までとなっているが、これに限るものではない。

ＲＦＩＤタグが移動したにも拘わらず、静止したＲＦＩＤタグであると識別された場合、図２０で示すように、このＲＦＩＤタグ（図２０の場合は赤シャツ１）を一覧から選択し、マウスの場合、右クリックで「このタグを”移動”と識別させる」を選択する。このようにすると、図２１Ａで示すように、該当するＲＦＩＤタグのグラフだけが、下部に表示される。そして、移動と識別すべき期間の開始点と終了点を選択する。選択はマウスの場合、左クリックで行う。このとき、図２１Ａ、図２１Ｂで示すように、「開始(終了)点を選んで下さい」というメッセージの下部に数値が表示される点を選択することができる。このようにして、安定的に読取が行われている範囲の電波受信強度増加開始点から終了点を選択する。
なお、選択のキャンセルは、グラフの点がないエリアをダブルクリックすることにより行う。

上記のように、選択した範囲のデータを「移動」と識別するときには、図２２で示すように、パラメータ変更の画面が、表示される。このパラメータ変更画面は、「このタグを“移動”と識別するためには以下青色欄のパラメータを変更下さい」とのメッセージと共に、このＲＦＩＤタグの値、“移動”とするための値、青色欄を変更についての欄が表示されており、秒間読取回数必要割合（％）、最大ＲＳＳＩ（ｄＢｍ）、ＲＳＳＩ差（ｄＢ）、ＲＳＳＩ増加発生必要割合（％）、期間前最大ＲＳＳＩからの必要増分（ｄＢ）、利用データ期間、位相振れに対する変更が必要なパラメータが青字で表示される。本実施形態では、「“移動”とするための値」の箇所に変更するための値（本実施形態では期間前最大ＲＳＳＩからの増分を「１２．５以下」）が表示される。そして、青色欄の数値を変更し、「保存」ボタンを押下する。

また、パラメータ変更後の動作確認を行うために、図２３で示すように、「リプレイ」ボタンを押下し、該当リプレイ用ログを選択することにより、パラメータ変更後の動作確認を行うことが可能となる。本実施形態では、「リプレイ」ボタンの下部には、ログファイル名が表示される。また本実施形態では、長期間（１分以上）のログが保存されているデータのリプレイには対応しない設定になっているが、これに限るものではない。読取時のアンテナ番号とリプレイ時の設定が異なる場合には、一覧で選択したＲＦＩＤタグのグラフが選択されないので、設定を合わせて使用することになる。

一方、移動と識別されたＲＦＩＤタグを“静止”とするには、図２４で示すように、該当ＲＦＩＤタグを右クリックして「このタグを“静止”と識別させる」を選択する。これにより、図２５で示すように、「このタグを“静止”と識別するには以下のいずれかのパラメータを変更下さい」とのメッセージと共に、このＲＦＩＤタグの値、“静止”とするための値、いずれかを変更の欄で、秒間読取回数必要割合（％）、最大ＲＳＳＩ（ｄＢｍ）、ＲＳＳＩ差（ｄＢ）、ＲＳＳＩ増加発生必要割合（％）、期間前最大ＲＳＳＩからの必要増分（ｄＢ）、利用データ期間、位相振れに対する変更可能なパラメータが青字で表示される。本実施形態では、これらの項目のうち変更する項目である「いずれかを変更」の欄が青色として表示されると共に、“静止”とするための値についてのガイダンスが表示される。
そこで、いずれか１つ以上の数値を「“静止”とするための値」に変更し、「保存」ボタンを押下する。

アンテナ２０が複数枚設定されている場合は、それぞれのアンテナ２０で移動識別が行われる。図２６で示す例は、同じＲＦＩＤタグがアンテナ２０（二つある場合の両方）の前を通過したケースである。なお、移動識別パラメータ画面でアンテナ２０設定を変更後は、「クローズ」、次に「オープン」ボタンを押下し、ＲＦＩＤリーダ１０へ再接続を行う。

次にログファイルについて説明する。
識別パラメータ設定ツールでは、移動識別サマリログと、リプレイ用ログと、使用パラメータバックアップと、エラーログの４つのログファイルが出力される。移動識別サマリログ、リプレイ用ログと使用パラメータバックアップの出力有無は「移動識別パラメータ」の画面で設定可能となっている。
（１）移動識別サマリログ
「読取状況グラフ表示」の画面で静止や移動と表示されたタイミングでタグ情報が保存される。移動したタグは初めに「静止」と出力された後に「移動」が出力される。
.\LOG ディレクトリ配下にスタートボタンを押下した日付・時刻のファイル名で作成される。
ファイル名例:
SUMMARY_2010_1217_1115_45.csv
出力内容は 「タグ情報表示時PC時刻, 静止・移動種別, EPC, PCBits,
ReaderIp, AntennaId, RSSI,（移動タグの場合は表２の６以降の情報）
出力例:
2010/12/24 19:07:40.108,STATIC,00000000000000000000B302,3000,192.168.55.31,1,-59
2010/12/24 19:07:43.061,MOVED,00000000000000000000B302,3000,192.168.55.31,1,-32,0.8,1293217662.753,1293217663.538,30,1.293,18.0,0.448,16.5,13.80,04
（２）リプレイ用ログ
本実施形態では、リプレイ機能用に全読取データを保存している。長時間読取実施時には出力をOFFとしておくことが好ましい。
ファイル名例:
2010_1217_1115_45.csv
（３）使用パラメータバックアップ
使用した識別パラメータのバックアップ。ツール起動後最初及びパラメータ変更発生後のスタートボタン押下時に作成される。
ファイル名例:
2011_0616_2118_27.xml
（４）エラーログ
リーダとの接続エラーやプログラムエラー等が出力される。
ファイル名: PACT_debug.log

次に、移動タグ識別エンジンについて説明する。
前記したように、識別パラメータ設定ツールで設定された方法で移動タグ識別を行う。
本実施形態で用いるミドルウエア（ナイスミドル）は、移動タグ識別エンジンのフロントエンドとして上位アプリとのインターフェース（Ｉ／Ｆ）機能を提供するものである。上位アプリ開発者が移動タグ識別ソフトのドキュメントとナイスミドルのユーザーマニュアル・ヘルプファイル・サンプルソースを参照して、開発を行う。

本実施形態で使用するＲＦＩＤリーダライタとしては、例えば、Impinj Speedway : F/W version 3.2.4 (LBT無版)、Impinj Speedway Revolution : F/W version 4.6.1以降 (LBT無)などを用いることができる。

本実施形態で動作するＯＳ，上位アプリケーションの開発環境例は、ミドルウエア「ナイスミドル」と同様にＷｉｎｄｏｗｓＸＰ（登録商標）やＷｉｎｄｏｗｓ７（登録商標）ＯＳでＭｉｃｒｏｓｏｆｔ ＶｉｓｕａｌＳｔｕｄｉｏ（登録商標）開発環境である。より具体的な構成ファイルとしては、
PickMovedTagsAgent.dll（移動タグ識別エンジン）
PickMovedTagsAgent.xml（識別パラメータ保存ファイル：編集不可）
PickMovedTagsParamsAutoConfigTool.exe（識別パラメータ設定ツール）
MightyControl.dll（識別パラメータ設定ツール構成ファイル）
EPCtoNAME.csv（識別パラメータ自動設定ツールでこのファイル内のEPCが読まれた場合、「名称 / シリアル」欄に名称を表示。名称,EPCのCSVフォーマットで記述。変更後はツールの再起動が必要。このファイルに含まれないEPCは下位3バイト(EPC欄の右側6文字)を10進数で「名称 / シリアル」欄に表示。）
log4net.dll（識別パラメータ設定ツールでエラーログ出力時に使用されるApache log4netライブラリ）
LICENSE.txt,NOTICE.txt（Apache log4netのLICENSE及びNOTICE）
.\WAVE\notify.wav（識別パラメータ設定ツールで移動タグを識別した場合に鳴る音ファイル。別のWAVEファイルをこのファイルに上書きすることで音を変更することが可能。）
.\Ja-JP\PickMovedTagsParamsAutoConfigTool.resources.dll
MightyControl.resources.dll（識別パラメータ設定ツール用”日本語(日本)”リソースファイル）
.\LOG\（ログファイル格納ディレクトリ）
その他(LLRP.dll, LLRP.Impinj.dll, LLRPNiceMW.dll, SET_READER_CONFIG.xml)はナイスミドルの構成ファイル、
などである。

本実施形態の移動識別方法は、基本的に、全読取タグから静止タグを排除し残ったものを移動タグと認識するものである。つまり、図２で示すように、全てのＲＦＩＤタグを読み取った全読取タグデータ（移動しているＲＦＩＤタグと静止しているＲＦＩＤタグの両方）から複数の個別静止ＲＦＩＤタグフィルタを用いて、静止したＲＦＩＤタグを特定し、取り除くものである。例えば、第１の個別静止タグフィルタＦ（１）によって、移動候補のＲＦＩＤタグを絞りこみ、この絞り込んだ移動候補のＲＦＩＤタグから、第２の個別静止ＲＦＩＤタグフィルタＦ（２）によって、さらに絞りこみ移動候補のＲＦＩＤタグを特定する。そして、次なる個別静止ＲＦＩＤタグフィルタによって、絞りこみ、さらに精度を向上させるためには、第Ｎの個別静止ＲＦＩＤタグフィルタＦ（Ｎ）・・・というように複数の個別静止ＲＦＩＤタグフィルタを用いるもので、最終的に移動しているＲＦＩＤタグを移動ＲＦＩＤタグとして特定するものである。

次に識別パラメータ設定ツールについて説明する。
利用データ期間中の最大値は、図１４Ａで示すように、利用データ期間後半の山頂部から、最小値は前半の山麓部から取得される。
また、期間前最大値取得期間を設定することができ、この期間の値を、大きくすることで、移動判定を厳しくすることが可能となる。つまり、図１４Ａで示すように、最大値は、期間前最大値から所定の増分となるのが通常である。

本実施形態における個別静止タグフィルタの例としては、以下のようなものがある。すなわち、
（１）チラ読みフィルタ（利用データ期間中に秒間読取回数が閾値未満）、
（２）ＲＳＳＩ絶対値フィルタ（最大ＲＳＳＩ（最大電波受信強度）が閾値未満）、
（３）ＲＳＳＩ非連続上昇フィルタ（「最大−最小ＲＳＳＩ（電波受信強度）」が閾値未満）、（ＲＳＳＩ（電波受信強度）（ｔ）＞電波受信強度（ｔ−１）である点或いはＲＳＳＩ（電波受信強度）（ｔ）＜電波受信強度（ｔ+１）である点の発生割合が閾値未満）、
（４）ＲＳＳＩ谷型フィルタ（ＲＳＳＩ（谷）：期間前最大電波受信強度からの増分が閾値未満）、
（５）過読フィルタ（移動候補判定からＸ秒後も読取が存在（移動判定まで＋Ｘ秒必要））、
（６）位相振れ少フィルタ（期間中の位相の変動から計算された検定量が閾値未満）、
である。

上記（１）のチラ読みフィルタは、利用データ期間中に秒間読取回数が閾値未満であるが、ここで、秒間読取回数について説明する。
秒間読取可能回数は、例えばあるＲＦＩＤタグの無線パラメータ（ＩＳＯ／ＩＥＣ１８０００−６Ｃ規格のミラーサブキャリア方式でＭ＝４， Ｔａｒｉ＝２０［ｕｓｅｃ］， ＢＤＲ＝３２０ｋ）において、ＲＦＩＤタグの枚数を変化させて、実際に読取を行い、１枚あたりの秒間読取数を計測する。このときのグラフが、図４で示す実線である。これはタグ枚数が１０枚以下の場合に点線で示された多項式曲線（ｙ＝０．２３６ｘ^２−５．３３３９ｘ＋４６．９７４）、１０枚より上の場合に累乗曲線（ｙ＝１３３．３４ｘ^{−０．８６０５}）によって決定係数０．９９８７、０．９９５４と高い精度で近似することができる。
他の無線パラメータを使用した場合でも同様に近似関数ｆ（ｘ）を求めることができ、使用中の無線パラメータと読取中タグ枚数により１枚あたり秒間読取可能回数を計算することができる。アンテナを１枚以上のN枚使用する場合にはｆ（ｘ／N）／Nで計算される数値を秒間読取可能回数として使用できる。この秒間読取可能回数に設定された割合を掛けることにより秒間読取回数閾値とできる。

上記（２）のＲＳＳＩ絶対値フィルタは、最大ＲＳＳＩ（最大電波受信強度）が閾値未満であるが、ここで、電波受信強度について説明する。
電波受信強度は、ＲＦＩＤリーダ１０が受信するＲＦＩＤタグからの電波受信強度理論値は次の式（１）で表わされる。

上式１をｄＢｍで表記すると、

となる。

ここで、アンテナ２０からの電波放射パターンは通常、図１３Ｂのようになり、図１３ＢのＡ地点とアンテナ２０前のＢ地点を比較すると、直線距離が１／２であるため電波受信強度は＋１２ｄＢとなる。この図１３Ｂのように、Ａ地点からＢ地点までタグが移動すると電波受信強度はｌｏｇ Ｇ_{ｒｅａｄｅｒ}の係数が２０であることから２０ｄＢ高くなる。上述の距離による値を加えると３２ｄＢの大きな差が連続的に傾向を持って発生することになる。

一方、アンテナ２０−タグ間には直接波以外にも床や天井、什器等での反射波・回折波が存在しそれらが全て干渉しあってマルチパスフェージングが発生する。マルチパスフェージングは仲上−ライス分布（直接波の強度が反射波や回折波に比べ強い場合）やレイリー分布（直接波の強度が弱い場合）でモデル化できることが知られており、それぞれの確率密度関数は以下の式（３）（４）で与えられる。

仲上−ライス分布：

ここでＫ＝ａ^２／２σ^２は直接波平均電力（ａ^２／２）と散乱波平均電力（σ^２）の比を与える。ｒは信号振幅である。

レイリー分布：

仲上−ライス分布の各種K及びレイリー分布での累積確率分布は、前記した図１５Ａのようになる。

アンテナ２０前を移動するＲＦＩＤタグは、直接波が支配的となるため仲上−ライス分布(K=10やK=40)となり、この場合の電波受信強度の変動は先述の距離やアンテナ放射パターンにより発生するものに比べかなり小さいものとなることがわかる。また、アンテナ前を人や物が通過することにより静止タグにフェージングが発生した場合の受信強度変化をレイリー分布でシミュレーションしたグラフが、図１５Ｂであり、受信強度の大きな差は発生することがあるがその変化が連続的に傾向を持って発生することは非常に稀であることがわかる。

これらのことからマルチパスフェージングを考慮しても前述手法（山型全体やその登り部分を線形・非線形曲線との適合度判定や強度最大値と最小値の差、時系列で前の点から上昇した点の割合といった各種算術的・統計的手法）が移動タグ識別手法として有効であることがわかる。

上記（３）のＲＳＳＩ非連続上昇フィルタは、「最大−最小ＲＳＳＩ（電波受信強度）」が閾値未満、ＲＳＳＩ（電波受信強度）（ｔ）＞電波受信強度（ｔ−１）である点或いはＲＳＳＩ（電波受信強度）（ｔ）＜電波受信強度（ｔ+１）である点の発生割合が閾値未満である。

上記（４）のＲＳＳＩ谷型フィルタは、期間前最大ＲＳＳＩ（電波受信強度）からの増分が閾値未満である。

上記（５）の過読フィルタは、移動候補判定からＸ秒後も読取が存在（移動判定まで＋Ｘ秒必要）である。図１４Ｂはこれを説明するグラフ図である。

上記（６）の位相振れ少フィルタは、期間中の位相の変動から計算された検定量が閾値未満であるが、ここで、位相について説明する。
位相は次の式（５）で表すことができる。

このうち次の３つの項はある移動判定期間中一定とみなすことができる。
θ_Ｔ： リーダからの電波送信時に変化する位相
θ_Ｒ： リーダでタグ応答波受信時に変化する位相
θ_ＴＡＧ：タグ応答時内部で変化する位相
波長（λ）は使用周波数により定まるので、位相θはリーダ−タグ間距離（ｒ）のみを変数とした関数となり、２π（３６０度）で剰余を取ることで周期関数となる。

ミドルウエアを用いた例について、より具体的に説明すると、移動するＲＦＩＤタグの情報を受け取るには、IsPickMovedTag及びSyncLocalTimeをtrueとしMovedTagReportHandlerを設定する。また、SessionOnプロパティの設定を行わないようにする。
静止・移動の区別無く読まれたＲＦＩＤタグの全てを受け取るには、ミドルウエア（ナイスミドル）のTagReportHandlerデリゲートも設定する。
本実施形態のデリケートとしては、MovedTagReportHandler(MovedTagData)であり、LLRPAgentのOnMovedTagReportイベントを処理するメソッドを表す。
ＬＬＲＰＡｇｅｎｔクラスとしては、イベントがOnMovedTagReportという名前の場合、デリケートはMovedTagReportHandlerであり、移動したタグ情報がMovedTagData型で渡される。
パブリックプロパティとしては、IsPickMovedTagとSyncLocalTimeがあり、ミドルウエア（ナイスミドル）から移動タグ情報を取得する場合trueに設定する。

パブリックプロパティについて、移動するＲＦＩＤタグ情報（ＭｏｖｅｄＴａｇＤａｔａ）とミドルウエア（ナイスミドル）でのＲＦＩＤタグ情報（ＴａｇＤａｔａ）との相違を踏まえて記述する。

次に、サンプルソースについて記述すると、ミドルウエアのソースコードに以下太字の行を追加する。
using mightycard.rfid.llrpmw;
try{
.
.
llrpAgent.OnErrorReport
+=new ErrorReportHandler(llrpAgent_OnErrorReport);
//→移動タグ識別エンジンを使用するための設定
llrpAgent.IsPickMovedTag = true;
llrpAgent.SyncLocalTime = true;
llrpAgent.OnMovedTagReport += new MovedTagReportHandler(llrpAgent_OnMovedTagReport);
//←移動タグ識別エンジンを使用するための設定
.
.
}catch(Exception
e)
MessageBox.Show(e.Message)
}

//OnMovedTagReportイベントの処理関数
void llrpAgent_OnMovedTagReport(MovedTagData tagData)
{
MessageBox.Show(“このタグは移動しました : ” + tagData.Epc)
………………..
}

この他に図２９のようにＲＦＩＤリーダ内のＬｉｎｕｘ ＯＳ（商標）上で「移動タグ識別エンジン」を動作させ、一方、「識別パラメータ設定ツール」については、コンピュータ装置Ｃで動作するような形態もある。この場合ＬＬＲＰミドルウエアを使用しないことになり、移動タグ識別エンジンが直接ＲＦＩＤリーダのSDKからＲＦＩＤタグの情報を取得することになる。
このようにＲＦＩＤリーダ内のＬｉｎｕｘ ＯＳ（商標）上で「移動タグ識別エンジン」を動作させることができる。
次に、識別パラメータ設定ツールでのパラメータ調整手順について、図２８のフローチャートに基づいて説明する。

まずステップＳ１０９で静止ＲＦＩＤタグをＲＦＩＤタグリーダの周囲に配置する。このとき実際の環境で想定される位置にＲＦＩＤタグを配置する。例えば売り場の盗難防止ゲートとしてアンテナが設置される場合、アンテナの一番近くの棚にＲＦＩＤタグ付き商品が置かれた環境で設定ツールを使用する。
次にステップＳ１１０で、不図示の操作手段により、識別パラメータ設定ツールを起動する処理を行う。次に、ステップＳ１２０で、コンピュータ装置Ｃ等の記憶装置に格納されている識別パラメータ保存ファイルを読取り、移動タグ識別エンジン（及びミドルウエア「ナイスミドル」）を使用してＲＦＩＤタグの読取りを開始する処理を行う。
ステップＳ１３０では、周囲に静止タグを配置した環境で移動と認識させたい方法で、ＲＦＩＤタグを移動させる。

ここで、「周囲に静止タグを配置した環境」とは、 実際の環境で想定される位置にＲＦＩＤタグを配置することを意味する。例えば、上述のように、売り場の盗難ゲートとしてアンテナ２０が設置される場合には、アンテナ２０近くの棚にＲＦＩＤタグ付き商品が置かれた環境で設定ツールを使用する。
また、「移動と識別させたい場合」とは、ＲＦＩＤタグが貼付された読取対象物を実際の移動方法（フォークリフト運転や人の歩行）でアンテナ２０前を移動させる。その際に読取条件を通常より厳しい方法（アンテナ２０に近い場所を高速移動、アンテナ２０から遠い場所を低速移動、ＲＦＩＤタグにフォークの爪があたり隠れる、鞄の中に多数のＲＦＩＤタグを隠す等）で移動させ、このような状態においてパラメータを合わせることで通常の移動の場合に余裕をもって移動と識別できるパラメータになる。

このとき、ステップＳ１４０で、静止したＲＦＩＤタグを移動しているＲＦＩＤタグと誤検知あるかどうか判別する。誤検知があると判断された場合（ステップＳ１４０：Ｙｅｓ）、ステップＳ１５０において、画面上で誤認識ＲＦＩＤタグを選択し、パラメータの変更を行う。そして、ステップＳ１６０で、リプレイ機能を起動し、誤検知したものが解消されているかどうかを確認する。そして、ステップＳ１３０の処理へ戻る。

一方、ステップＳ１４０で誤検知がないと判断された場合（ステップＳ１４０：Ｎｏ）、ステップＳ１７０で、移動ＲＦＩＤタグを静止ＲＦＩＤタグと誤検知あるかどうか判断する。ステップＳ１７０で誤検知があると判断された場合（ステップＳ１７０：Ｙｅｓ）、ステップＳ１８０で、画面上で誤認識したＲＦＩＤタグを選択し、グラフ上から移動と識別されるべき期間を選択する。そしてステップＳ１９０でツールで計算された範囲の値にパラメータを変更する処理を行う。この処理の後で、ステップＳ２１０で、リプレイ機能を起動し、誤検知したものが解消されているかどうかを確認する。そして、ステップＳ１３０の処理へ戻る。

また、ステップＳ１７０で、誤検知がないと判断された場合（ステップＳ１７０：Ｎｏ）、ステップＳ２００で識別パラメータ決定し、設定ファイルに保存する処理を行う。

次に、本発明の応用例として倉庫や店舗での入出荷検品、店舗での盗難防止ゲート、倉庫内盗難品検知、倉庫内棚卸について説明する。それぞれの用途の目的や読取環境によって使用するフィルタやそのパラメータ設定が異なってくる。他の応用用途としては空港や工場等でコンベア上を移動するタグ付き荷物や物品の識別が上げられる。電波受信強度変化に加え位相やドップラー情報を使用することによって１枚のアンテナでも複数コンベアの読み分けや移動方向、速度の検知を行うことができる。

まず、倉庫や店舗での入出荷検品への応用例について説明する。図１のようにアンテナ前をタグが貼付された読取対象物（洋服やバックといったアパレル製品やテレビ・携帯電話等の電気製品や部品等）を台車やフォークリフト等で通過させることで高速に大量の検品を行うことができる。この用途の特徴としてはタグが重なる等によって読みにくい（ＲＳＳＩが低いままの）状況が発生し得ること、アンテナのすぐ近くに対象外タグを置かないといった環境コントロールが比較的容易であることがあげられる。

こういった環境では例えばＲＳＳＩ絶対値フィルタを使用せず、またＲＳＳＩ非連続上昇フィルタのパラメータを甘くすることで読みにくいタグも移動と検知させながら、チラ読みフィルタや位相振れ少フィルタによって遠くの静止タグを排除することが可能である。

次に、店舗での盗難防止ゲートへの応用例について説明する。図１３Ｃのようにレジで会計が済んでいないタグ付き商品がＲＦＩＤアンテナ及びリーダが組込まれたゲート前を移動したことを検知しブザーを鳴らすことが可能となる。従来の磁界方式のゲートではゲート間距離が１ｍを超えると検知が難しくなるが、電磁界方式のＵＨＦ帯ＲＦＩＤを使用することにより２ｍ以上といったゲート幅でも検知可能となる。また、長距離通信の特性を活かしアンテナを天井や床下に配置することも可能である。
この応用では、ゲート周りのタグを読んでしまい、盗難を行っていない人が通過したときにブザーを鳴らしてしまうとお客の感情を害し、お店の評価を落としてしまう危険性がある。そのため全ての移動タグを検知できなくても可とし、各種移動フィルタの閾値を高めに設定することができる。図１３Ｃ（Ａ）のような水平方向移動は検知しながら（Ｂ）のような垂直方向移動はブザーを鳴らさないようにするために、まず識別パラメータ設定ツールで垂直方向に移動した場合のデータを取得し、それが静止と識別できるよう［移動識別パラメータ］欄の調整を行う。
多くの場合フィルタを厳しめに設定することで水平方向は移動と識別しながら垂直方向を静止と識別することができる。しかし、特にアンテナ近くの裏面から表面への移動時、垂直方向も識別されてしまう場合は、［識別オプション］の「垂直方向移動フィルタ」を利用すること（チェックをＯＮ）により解決することが可能である。

さらに、倉庫内盗難品検知への応用例について説明する。図３０で示すようにＲＦＩＤアンテナ２０を出入口付近の天井等に設置することによって、倉庫内のタグつき物品の移動が発生した場合に検知を行うことができる。

図３０は倉庫４階で、“緑花瓶”の盗難が発生したことを地下の中央監視室で検知した例である。カメラや赤外線センサで監視する方法と違い、暗闇の中でも検知できたり、鞄の中などに隠された場合でもどの物品が盗まれたのかわかったりするのが大きな特徴である。位相は人の移動だけでも振れることがあり静止タグを移動と誤検知してしまう可能性があるが、通常人の出入りが発生することのない夜間等の監視時には位相振れ少フィルタのみを使用することによりＲＳＳＩ非連続上昇フィルタ等使用時よりも広範囲の進入検知が可能となる。進入検知後にＲＳＳＩ非連続上昇フィルタ等で物品の特定を行うことができる。

次に、倉庫内棚卸への応用例について説明する。図３１で示すように、倉庫２００は既存のものであり、物品（商品）の高さ等に応じた棚が設けられている。この棚にはＲＦＩＤタグからなるロケーションタグ２０１が付されており、ロケーションタグ２０１が示す位置情報が得られることになる。一方、棚に格納される物品（商品）にも、その性質や数量等を含めて商品情報が記憶されたＲＦＩＤタグからなる商品タグ２０２が付されている。アンテナ２０とＲＦＩＤリーダ１０が格納されている制御装置（コンピュータ装置Ｃ）を備えた自動読取装置をフォークリフト２１０で移動させることにより、静止したアンテナ前をタグが移動したのと同じ状態となる。

フォークリフト移動中にアンテナ前にあるタグのＲＳＳＩのみが大きな山型になることを利用しＲＳＳＩ非連続上昇フィルタのパラメータを適切に設定することによって丁度アンテナ前に位置するロケーションタグ及び商品タグのみを検知することができる。これにより、どの位置にどの商品があるという棚卸を安全・高速に行うことができる。無線回線を通して、読取った商品情報、位置情報、紐付けた商品情報（有無含む）と位置情報をリアルタイム転送する。長距離読取に優れたＵＨＦ帯ＲＦＩＤを使用することによって同種の装置でバーコード技術を使う場合に比べタグの貼付位置を自由に配置できる。また、箱の内側の非可視位置に貼付することも可能となる。

１０ ＲＦＩＤリーダ
２０ アンテナ
３０ ＲＦＩＤタグ
Ｆ 静止ＲＦＩＤタグフィルタ
Ｃ コンピュータ装置

Claims (6)

1. 複数のＲＦＩＤタグと、該複数のＲＦＩＤタグからの電波を受信するＲＦＩＤリーダを用いて、複数のＲＦＩＤタグのうち移動ＲＦＩＤタグを識別するＲＦＩＤタグ移動識別方法において、
   前記ＲＦＩＤリーダによって移動ＲＦＩＤタグ及び静止ＲＦＩＤタグを含む複数のＲＦＩＤタグの時系列読取データを取得し、一定のパラメータに基づく静止ＲＦＩＤタグフィルタによって、前記読取った複数のＲＦＩＤタグの情報から前記静止ＲＦＩＤタグを特定し、前記複数のＲＦＩＤタグから前記移動ＲＦＩＤタグを識別可能とするもので、
   前記静止ＲＦＩＤタグフィルタは、以下の（１）乃至（６）の工程の何れかからなる個別静止ＲＦＩＤタグフィルタを備え、
   前記静止ＲＦＩＤタグフィルタは、複数の前記個別静止ＲＦＩＤタグフィルタのうち少なくとも一つ以上の個別静止ＲＦＩＤタグフィルタを用いて静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別し、移動ＲＦＩＤタグを識別することを特徴とするＲＦＩＤタグ移動識別方法；
   （１）．読取中のタグ枚数を変数とした関数を用いて、秒間読取可能回数を計算し、前記ＲＦＩＤリーダによって読み取った秒間読取回数との比較に基づいて識別する工程；
   （２）．電波受信強度の時系列分析によって、静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する工程；
   （３）．前記ＲＦＩＤリーダで読みとられたＲＦＩＤタグの電波受信強度時系列データが、谷型となる場合、すなわち電波受信強度の高い点発生時に所定の時系列データまで遡り、別の電波受信強度が高い点があるかを判定し、高い点があった場合に静止ＲＦＩＤタグと判別する工程；
   （４）．前記ＲＦＩＤタグの電波受信強度時系列データが示す山型全体、山型の登り部分を線形・非線形曲線として適合する度合い、山型の各強度最大値と最小値の差、山型の時系列で前の点から上昇した点の割合のうち少なくとも一つを用いて静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する工程；
   （５）．前記ＲＦＩＤリーダの送信搬送波とＲＦＩＤタグからの反射で返却された搬送波の位相の時系列分析によって、静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する工程；
   （６）．ＲＦＩＤタグ応答の位相変曲点またはドップラー周波数０時点と電波受信強度山頂発生時点を比較することによって移動方向を識別する工程。
2. 前記（１）における読取中タグ枚数を変数とした関数は、読取中タグ枚数を独立変数とした累乗曲線関数または多項式関数であり、この関数を用い計算した計算値に、所定の割合を乗算した閾値を判断基準として、前記判断基準より秒間読取回数が少ない場合に静止タグと識別することを特徴とする請求項１記載のＲＦＩＤタグ移動識別方法。
3. 前記各種ＲＦＩＤタグフィルタを使用しアンテナ近くの水平方向移動のみを検知し遠くの水平方向移動や垂直方向移動は検知しないことを特徴とする請求項１記載のＲＦＩＤタグ移動識別方法。
4. 前記（５）における時系列分析において、位相の標準偏差の値に閾値を設けること、二項検定、カイ二乗検定、Ｇ検定の少なくとも一つを用いて静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する工程を備えたことを特徴とする請求項１記載のＲＦＩＤタグ移動識別方法。
5. 前記一定のパラメータは、予め試験的ＲＦＩＤタグを移動した後に、移動時の開始点及び終了点を選択することでＲＦＩＤリーダやアンテナ、タグ、移動速度といった環境の違いを吸収した前記パラメータを自動で計算することを特徴とする請求項１記載のＲＦＩＤタグ移動識別方法。
6. コンピュータ装置に、移動ＲＦＩＤタグ及び静止ＲＦＩＤタグを含む複数のＲＦＩＤタグをＲＦＩＤリーダで読取らせた時系列読取データを、一定のパラメータに基づく静止ＲＦＩＤタグフィルタよって、複数の前記ＲＦＩＤタグの情報から前記静止ＲＦＩＤタグを特定させ、前記複数のＲＦＩＤタグから前記移動ＲＦＩＤタグを識別させるものであって、静止ＲＦＩＤタグフィルタとして、
   （１）．読取中のタグ枚数を変数とした関数を用いて、秒間読取可能回数を計算し、前記ＲＦＩＤリーダによって読み取った秒間読取回数との比較に基づいて識別する手順、
   （２）．電波受信強度の時系列分析によって、静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する手順、
   （３）．前記ＲＦＩＤリーダで読みとられたＲＦＩＤタグの電波受信強度時系列データが谷型となる場合、すなわち電波受信強度の高い点発生時に所定の時系列データまで遡り、別の電波受信強度が高い点があるかを判定し、高い点があった場合に静止ＲＦＩＤタグと判別する手順、
   （４）．前記ＲＦＩＤタグの電波受信強度時系列データが示す山型全体、山型の登り部分を線形・非線形曲線として適合する度合い、山型の各強度最大値と最小値の差、山型の時系列で前の点から上昇した点の割合のうち少なくとも一つを用いて静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する手順、
   （５）．前記ＲＦＩＤリーダの送信搬送波とＲＦＩＤタグからの反射で返却された搬送波の位相の時系列分析によって、静止ＲＦＩＤタグと移動ＲＦＩＤタグを識別する手順、
   （６）．ＲＦＩＤタグ応答の位相変曲点またはドップラー周波数０時点と電波受信強度山頂発生時点を比較することによって移動方向を識別する手順、
   のうち少なくとも一つ以上の手順からなる個別静止ＲＦＩＤタグフィルタを任意の順でコンピュータ装置に実行させるためのＲＦＩＤタグ移動識別プログラム。