JP2005124198A - デジタル信号中のノイズ誤りを訂正する方法 - Google Patents
デジタル信号中のノイズ誤りを訂正する方法 Download PDFInfo
- Publication number
- JP2005124198A JP2005124198A JP2004298490A JP2004298490A JP2005124198A JP 2005124198 A JP2005124198 A JP 2005124198A JP 2004298490 A JP2004298490 A JP 2004298490A JP 2004298490 A JP2004298490 A JP 2004298490A JP 2005124198 A JP2005124198 A JP 2005124198A
- Authority
- JP
- Japan
- Prior art keywords
- symbol
- signal
- value
- symbols
- sequence
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
Abstract
【解決手段】チャネル[20]で破損した受信信号[22]を処理して、受信信号[22]よりも信号破損の少ない処理済み信号[23]を生成する装置および方法が開示される。本装置[10]は、受信信号[22]、値Iのシンボルが値Jのシンボルにチャネル[20]で変換される確率を指定する情報、および、値Iのシンボルが値Jのシンボルで置換される場合に起こる信号劣化を測る信号劣化関数を指定する情報を記憶する。コントローラ[11]は、受信信号[22]中のシンボルのコンテクスト[24]において値Iの各シンボルを値Jのシンボルで置き換える。この値Jは、前記コンテクスト[24]で観察された統計を介して推定される、基礎となるノイズのない信号[21]に関連する処理済み信号[23]の全信号劣化を最小にする。
【選択図】図3
Description
電話での会話を考える。
この会話では、人はマイクロフォンに対して話しかけ、この出力は、デジタル化された後、さまざまな電話回線と切替システムとを介して、もうひとりの人に対して送られる。
もうひとりの人の場所で話す人は、一連のデジタル値を受信し、その後、このデジタル値は、もうひとりの人に再生される。
一般的に、この受信されたシーケンスは、伝送システム、デジタル−アナログ変換器、およびアナログ−デジタル変換器により導入された誤りのために、送信されたシーケンスとは異なる。
例えば、伝送システムのノイズにより、送信されたシーケンス中のいくつかのデジタル値は、変換されている。
ノイズ除去システムの1つの目的は、このようなノイズ誤りを、可能な限り多く除去することである。
一般的に、一連の「シンボル」を含む入力デジタル信号は、「伝達リンク」を通って伝送され、伝達リンクの出力で、出力デジタル信号として受信される。
この出力デジタル信号もまた、一連の「シンボル」を含む。
これらのシンボルそれぞれは、アルファベットとして触れられる所定のシンボルの集合から選択される。
出力信号は、出力信号と同じアルファベットで書かれていると仮定される。
この場合、入力信号および出力信号は、一連の0と1とを含む。
しかしながら、他のアルファベットが、通常使われる。
例えば、シンボルそれぞれが0〜M−1の整数で表されたデジタル信号は、デジタル加入回線(DSL)を介して、ユーザをインターネットにつなぐために、ブロードバンドデータ伝達システムで通常使われる。
例えば、磁気ディスクドライブにデータファイルを記憶することは、通信リンク、すなわちディスクドライブを通じてのデジタル信号の伝送とみなすことができる。
入力信号は、所定のアルファベットから選択されるシンボル(例えばバイトデータ)のシーケンスである。
バイトデータの場合、各シンボルは集合[0,1,...,255]から選択される整数値を取る。
ディスクドライブから取り出されるファイルもまた、この集合から選択されるシンボルのシーケンスからなる。
入力信号シンボルは、ディスクドライブの電子機器によって処理され、局所化された磁界の形態で記憶され、この磁界が読み出されて出力信号が生成される。
シンボルと磁界の間の変換を行うデジタル−アナログ回路中のノイズが、出力信号に誤りを導入する。
また、記憶されている間に、ランダム事象によって、磁界が変化することがあり、これがさらなる誤りを導入する。
この場合、信号は、光検出器中のノイズで破損を受ける画像である。
本装置は、第1、第2、および第3メモリならびにコントローラを含む。
第1メモリは、受信デジタル信号を記憶する。
受信デジタル信号は、チャネルで破損したシンボルのシーケンスを含む。
第2メモリは、値Iのシンボルが値Jのシンボルにチャネルで変換される確率を指定する情報を記憶する。
第3メモリは、値Iのシンボルが値Jのシンボルで置換される場合に起こる信号劣化を測る信号劣化関数を指定する情報を記憶する。
コントローラは、受信デジタル信号中のシンボルを置換することによって、受信デジタル信号から処理済みデジタル信号を生成する。
コントローラは、受信デジタル信号中のシンボルのコンテクストにおいて値Iの各シンボルを値Jのシンボルで置き換えれば、信号劣化関数で測った全信号劣化が受信デジタル信号に比べて処理済みデジタル信号において減少すると推定される場合に、この置換を行う。
コンテクストは、一定値のシンボルのシーケンスであって、置換が考慮されているシンボルに関して所定の位置にあるシンボルのシーケンスを含む。
コントローラは、第2メモリ内の情報を利用して、受信信号中で置換が考慮されているシンボルと同じ値およびコンテクストを有するシンボルにチャネルで変換された、ノイズのないもとの信号中の各シンボル値の出現回数を推定する。
以下の説明において、ノイズを導入する処理システムを「チャネル」と呼ぶ。
というのは、このようなシステムは、信号が通る伝送チャネルに類似しているからである。
次に図1を参照する。
図1は、本発明の一実施形態によるノイズ除去システム10のブロック図である。
ノイズ除去システム10は、チャネル20で破損した出力信号に作用する。
チャネルは、既知のアルファベットからのシンボルのシーケンスを含む入力信号21に作用して、同じくそのアルファベットからのシンボルのシーケンスを含む出力信号22を生成する。
また、本発明は、チャネルがメモリを有さない場合、すなわち、シンボルAがシンボルBに変換される確率が、入力信号中でAに先行または後続するシンボルとは独立であり、かつ先行または後続するシンボルがどのように変換されるかとも独立である場合に、最適化される。
一般的に、本発明の基礎となるアルゴリズムは、上記の仮定が成り立つことを仮定したデータに作用するが、アルゴリズムの作用は、仮定の妥当性とは無関係に明確に定義される。
具体的には、チャネル(変換)確率情報は、アルゴリズムのパラメータ付きクラスをもたらす調整可能なパラメータの集合とみなされることができる。
それらのアルゴリズムでは、データおよびノイズの実際の特性に応じて、一部のアルゴリズムが、他のアルゴリズムよりも良好に働く。
上記の仮定が成り立ち、かつ仮定された変換確率の結果として実際のノイズが生じる場合、アルゴリズムは、明確に定義された数学的意味での定理として表現されうるいくつかの興味深い性質を有する。
この状況は、望ましいが、本発明の実際の適用に必須ではない。
劣化は、AおよびBの値が異なれば異なってもよい。
以下の説明では、この劣化尺度は、単に「劣化」と呼ばれ、D(A,B)で表される。
例えば、8ビットアナログ−デジタル変換器を利用して、アナログ時変信号をデジタル値のシーケンスに変換することによって生成されるデジタル信号を考える。
結果として得られるデジタル信号は、デジタル値0〜255に対応する256個のシンボルを有するアルファベットから選択されたシンボルのシーケンスである。
出力信号がアナログ信号に逆変換されて、人間の観測者に対して再生されると仮定する。
1だけ変化したシンボルから生じる出力信号中の誤りは、2だけ変化したシンボル等から生じる誤りよりも、通常ははるかに小さい。
したがって、この場合、劣化関数は、シンボルが変化される量によって決まる。
そして、これらのシンボルを変更することにより、チャネルで導入された全信号劣化(劣化関数で測る)を減少させたノイズ除去出力信号23を生成する。
出力信号中のシンボルを一度に1つずつ考慮する場合、シンボルがほぼ等しい頻度で入力信号に現れるならば、特定のシンボルが誤っている確率に関して得ることができる情報は非常に少ない。
種々のシンボル値が入力信号中で等しい頻度を有すると仮定すると、シンボルが出力信号に現れる頻度もまた、対象となるほとんどの場合にほぼ等しくなる。
1を0に変換する確率が0を1に変換する確率と等しいチャネルを通る2値信号を考える。
入力信号中にほぼ等しい個数の1および0がある場合、出力信号もまたほぼ等しい個数の1および0を有する。
というのは、1が0に変換されるごとに、0が1に変換されることもあるからである。
したがって、他のシンボルに比べて、ノイズの結果である可能性が高いような特定のシンボルを識別することはできない。
例えば、入力信号中でのシンボルAの出現がシンボルBに比べて非常にまれであり、チャネル誤りによりAがBに変換されBがAに変換されると仮定する。
また、AがBに変換される確率は、BがAに変換される確率と等しいと仮定し、この確率をPで表す。
Pは既知であると仮定する。
入力信号中の値AおよびBのシンボルの個数を、それぞれNAおよびNBで表す。
出力信号中の値AおよびBのシンボルの個数を、それぞれN'AおよびN'Bで表す。
ノイズ除去器にとって利用可能であるのはN'AおよびN'Bだけであることに留意すべきである。
出力信号中の値Aのシンボルの個数は、入力信号中の値Aのシンボルの個数から、別のシンボルに変換された値Aのシンボルの個数を差し引き、Aに変換された他の値のシンボルの個数を加えたものに等しい。
この説明を簡単にするため、アルファベット内にあるのがAおよびBの2個のシンボルのみである場合を考える。
その場合、
N'A≒NA−P*NA+P*NB
および
N'B≒NB−P*NB+P*NA
である。
ここで、総カウント数で除算/正規化した後の近似における誤差は、高い確率で、特にシンボル数が大きくなるにつれて小さくなる。
本明細書全体を通じて、正規化されていないシンボル数の近似値または推定値における誤差が「小さい」という場合、または正規化されていないシンボル数に関わる量が互いに「近い」という場合には、関連する誤差は、必ずしも絶対的尺度で小さいのではなく、シンボルの総数に比べて小さい、すなわち、信号中のシンボルの総数で除算/正規化した後で小さいと理解するものとする。
したがって、与えられたP、N'A、およびN'Bに対して、上記の方程式系(≒を=で置き換えたもの)をNAおよびNBについて解くことができ、それによって、入力シーケンス中にAおよびBが現れる回数の推定値が、それぞれ得られる。
これらのAおよびBのカウントの推定値を、それぞれnAおよびnBで表すことにする。
なお、それらに対応するNAおよびNBは、入力シーケンス中の実際のカウントを表しているが、想定されるアプリケーションでは、ノイズ除去器にとって利用可能でないことに留意されたい。
この例では、NA≪NBおよびP≪1と仮定する。
したがって、N'AはNA+P*NBにほぼ等しく、N'BはNBにほぼ等しい。
この場合、チャネル誤りが、まれなシンボルの頻度を増大させることが分かる。
出力信号中のN'A個のシンボルは、シンボルの2つのクラスからなるとみなされることができる。
第1のクラスは、チャネルによって正しく伝送されたシンボルからなり、第2のクラスは、値Bのシンボルを値Aのシンボルに変換するという形で誤って伝送されたシンボルからなる。
残念ながら、この判定を行う方法はない。
本発明は、他の2つの可能なアクション、すなわち、何もしないか、すべてのAをBに逆変換するかについて調べる。
本発明は、D(A,B)で測った推定信号劣化が最も小さくなるアクションを実行する。
出力信号は、チャネルによるBからAへの変換の結果として、ほぼPNB≒N'A−NA個の誤りを有する。
これらの誤りに関して何もしない場合、(N'A−NA)Dという全信号劣化が、出力信号に存在する。
推定される信号劣化はPnBDとなり、これもまた、この例ではN'A−nAでよく近似される。
これは、BがAに変わったことから生じたN'A−NA個の誤りを訂正する。
この変更はまた、出力信号中におよそ(1−P)NA≒NA個の新たな誤りも導入する。
というのは、出力信号中でAとして正しく表現されているNA個のシンボルがBに変換されるからである。
したがって、出力信号中ですべてのAをBに変更することから生じる全信号劣化は、およそNA*Dであり、本発明は、これを(1−P)nAD(この例ではnADでよく近似される)であると推定する。
なお、(N'A−NA)D>NA*Dである場合、出力信号中のすべてのAをBとおくと、全信号劣化が減少することに留意されたい。
したがって、この例では、N'A>2NAの場合、すべてのAをBに変換すべきである。
実際にはNAおよびNBは未知であるので、本発明のアルゴリズムは、(N'A−NA)およびNAをそれぞれの推定値PnB(≒N'A−nA)および(1−P)nA(≒nA)で置き換えることによって得られるこの決定ルールの近似を適用する。
本発明のアルゴリズムは、シンボルがほぼ等しい頻度を有する場合にも適用可能であるが、その場合にアルゴリズムを適用すると、出力信号のシンボルは、ほとんど常に不変なままとなる。
入力信号および出力信号のシンボルがM個の値を有するアルファベットから選択されるシステムを考える。
図2を参照する。
図2は、出力信号の一部を示しており、シンボルa1,a2,...,ak(aiは一定)のシーケンスが、出力信号中のいくつかの位置24に現れるように示されている。
シーケンスa1,a2,...,ak,x(ここでxはシーケンスに続く次のシンボルである)は、xの値によって決定されるM個の値を取り得る「スーパーシンボル」とみなされることができる。
xおよびx'の或る値に対して、a1,a2,...,ak,xがシーケンスa1,a2,...,ak,x'よりもはるかに小さい頻度で現れるようなシーケンスa1,a2,...,akが見つけられれば、上記の説明で用いたアルゴリズムの一般化が適用され、全信号劣化が改善されることができる。
もう一度述べると、チャネルへの入力信号をシーケンス{yi}で表し、チャネルからの出力信号をシーケンス{zi}で表す。
ここで、i=1〜Nである。
本発明の目的は、{zi}中のシンボルを変化させて新たな出力信号{z'i}を得る際に、{yi}と{z'i}との差から生じる全劣化が{yi}と{zi}との差から生じる全劣化よりも小さくなるようにすることである。
一定値を取り、当該シンボルに関して所定の位置にあるシンボルのシーケンスを、当該シンボルの「コンテクスト」と呼ぶことにする。
例えば、シーケンスa1,a2,...,ak,xが出力信号中に現れる場合、上記のシーケンスa1,a2,...,akは、シンボルxのコンテクストである。
シンボルxは、多数の異なるコンテクスト中に現れる場合があることに留意すべきである。
また、コンテクストは、当該シンボルのちょうど左側にあるシンボル以外のシンボルに対して、他の何らかの所定の関係にあるシーケンスであってもよい。
このような場合については以下でさらに詳細に説明するが、以下の説明を簡単にするため、当面、出力信号中のシンボルzjのコンテクストは、出力シンボルのシーケンス中で、当該シンボルの左側にあるk個のシンボルであると仮定する。
また、kは一定であることも仮定する。
kを選択する方法についても、以下でさらに詳細に説明する。
そのシンボルのコンテクストは、シーケンス[zj−k,...,zj−2,zj−1]である。
すなわち、a1=zj−k、a2=zj−k+1等である。
以下の説明を簡単にするため、値a1,a2,...,akのシーケンスは、ベクトルaで表される。
本発明は、0からM−1までのxのそれぞれの値について、出力信号中でコンテクストaの後に値xが続くようなコンテクストaが見出される回数をカウントする。
観測されたカウントから、出力信号中のコンテクストの後にもとの信号{yi}中のそれぞれの可能なxの値が続くコンテクストが現れる回数が、そのチャネルに対する既知の誤り率を用いて、推定されることができる。
これらの推定値が求められた後、上記の方法に類似した方法で、既知の劣化関数に基づいて、a,xの形のすべてのシーケンスをa,x'に置き直す決定がされうる。
チャネルの出力が、2つの状態を取るシンボルのシーケンスであるような伝送、すなわち2値符号化伝送を考える。
チャネルが或るビットを0から1に変える確率は、そのビットを1から0に変える確率に等しいと仮定する。
また、0から1への切替から生じる信号劣化は、1が0に切り替わることから生じる信号劣化に等しいこと、すなわちD(0,1)=D(1,0)=Dであることも仮定する。
2よりも多くのシンボル値を取るアルファベットを用いて、より一般的な劣化関数を用いるという、より一般的な場合の説明は、後で行われる。
この例の目的上、ビットa1,a2,...,akの或る選択に対して、シーケンスa1,a2,...,ak,0は、シーケンスa1,a2,...,ak,1に比べて、伝送中にまれにしか出現しないと期待されることが既知であると仮定する。
上記の例で言えば、a,0は「スーパーシンボル」値Aに対応し、a,1はスーパーシンボル値Bに対応する。
本発明によるノイズ除去システムは、すべてのシーケンスa1,a2,...,ak,0の最終ビットを0から1に(またはこの逆に)置き直した結果として得られる信号について、劣化関数で測った全劣化が、チャネルから受信されるもとの信号よりも小さいと推定される場合に、この置き直しを行うことによって作用する。
zjに先行するkビットをベクトルa=[a1,a2,...,ak]で表す。
本発明の方法は、まず、出力信号中にシーケンスa,0およびa,1が現れる回数をカウントする。
これらの2つのカウントを、それぞれN'0およびN'1で表す。
この例の目的上、N'0≪N'1と仮定する。
すなわち、出力信号中にシーケンスa,0は、シーケンスa,1よりも、はるかに小さい頻度で出現する。
シーケンスa,0が出力信号中に生じる可能性がある理由として、入力信号がそのシーケンスを含んでいた場合がある。
また、このシーケンスは、チャネルがシーケンスa,1の一部のインスタンスをa,0に変換したために現れる場合もある。
同様に、シーケンスa,1が出力信号中に生じ得る理由として、入力信号がこのシーケンスを含んでおり、チャネルがそのシーケンスを正しく伝えた場合がある。
あるいは、a,1は、チャネルが誤ってシーケンスa,0をa,1に変換したために生じることもある。
出力シーケンス中のaと、対応する入力シーケンス位置における0および1とにより、シーケンスa,0およびa,1が現れる回数を、それぞれN0およびN1で表す。
具体的には、N0は、[zj−k,zj−k+1,...,zj−1]=aおよびyj=0となるようなインデックスjの数をカウントし、N1は、[zj−k,zj−k+1,...,zj−1]=aおよびyj=1となるようなインデックスjの数をカウントする。
なお、N0およびN1は、破損していない入力信号に関する観測されない情報に関わるので、ノイズ除去アルゴリズムにとって利用可能でないことに留意されたい。
或るビットが0から1へ、または1から0へ反転する確率をPで表す。
すると、N0およびN1は、n0およびn1についての次の方程式系を解くことによって推定されることができる。
N'0=n0−n0*P+n1*P (a)
N'1=n1−n1*P+n0*P (b)
ここで、N'0およびN'1はそれぞれ、出力シーケンス中にa,0およびa,1が現れる回数である。
本発明により上記の方程式を解くことに基づいて得られるN0およびN1の推定値をそれぞれ、n0およびn1で表すことにする。
チャネルにより導入される誤りの結果である出力信号中のシーケンスa,0のインスタンスの数は、高い確率でPN1に近づく。
本発明は、これをPn1により推定する。
したがって、何もしない場合に受ける信号劣化は、Pn1Dと推定される。
a,0のすべてのインスタンスがa,1に変換される場合、およそ(1−P)N0個の誤りが生じる。
というのは、これらのインスタンスのうちの(1−P)N0個は正しいと期待されるからである。
本発明は、この量を(1−P)n0により推定する。
したがって、a,0のすべてのインスタンスをa,1とおくことから生じる信号劣化は(1−P)n0Dと推定される。
よって、本発明は、(1−P)n0<Pn1の場合、a,0のすべてのインスタンスをa,1に置き直す。
P≪1およびN0≪N1の場合、高い確率で(1−P)n0はN0に近づき、Pn1はN'0−N0に近づく。
すなわち、この範囲のP、N0およびN1に対し、高い確率で、本発明によるノイズ除去は、コンテクストaとともに出現する0のすべてのインスタンスを、そのような0の出現の半数より多くに誤りがある場合に1とおくという合理的ルールに対応する。
チャネル誤り率Pが小さいと仮定すると、すべてのa,1シーケンスをa,0に置き直すコストは、ほとんど常に、何もしないよりもはるかに大きくなることを示すことができる。
というのは、これらの値のほとんどが正しいからである。
したがって、a,0をa,1に置き直す場合のみを詳細に考察した。
しかし、上記のアルゴリズムは、この場合にも適用可能である。
その結果は単に、シンボルを不変のままにすることとなる。
状況によっては、これが当てはまらない場合がある。
D(0,1)がD(1,0)とは異なる場合、ビットを置き直す判断は、D(0,1)とD(1,0)との相対値によって決まる。
その場合、上記の例において、何もしないことに伴う劣化は、Pn1D(1,0)と推定される。
a,0のすべてのインスタンスをa,1に変更することに伴う劣化は、同様にして(1−P)n0D(0,1)と推定される。
本例では、(1−P)n0D(0,1)<Pn1D(1,0)の場合に、a,0のすべてのインスタンスがa,1に変更される。
本発明は、入力信号が、統計的にかなりの回数繰り返されるコンテクストを有する場合、すなわち、入力信号中のxの前に出力シーケンス中のaが現れるようなxの出現回数を統計的に高い精度で推定することができるほど繰り返してシーケンスaが見出される場合に、最大の利益を提供する。
上記の分析は、チャネル誤りから生じる観測シーケンスの一部を推定することによって決まる。
この推定は、チャネルを出るシーケンス中で、各シーケンスの観測される出現頻度に基づく。
このデータの精度は、出力シーケンス中で観測されるシーケンスの頻度の統計精度によって制限される。
出力シーケンス中で観測されるシーケンスの出現回数が小さい場合、上記の式(a)および(b)の精度が低いため、これらのさらに精度の低い方程式を解くことから導かれる推定値n0およびn1の精度も同様に低くなる。
これらのカウントの精度が十分に低い場合、出力信号中のシンボルを置き直すことに関して、誤った判断がなされるであろう。
さらに、Q≫1として、pの任意の値に対して、シーケンスaはQ回現れるが、それよりも大きいコンテクストp,aは頻繁には現れないと仮定する。
すると、kよりも大きいコンテクストの出現回数ははるかに小さくなり、小さいほうのコンテクストが用いられた場合に比べて、カウントの統計精度が劣化することになる。
したがって、kに対して大きすぎる値を選択すると、決定誤りを引き起こす可能性がある。
kが小さすぎる場合、入力信号中のa,xに対する推定値の統計精度は高くなるが、xのそれぞれの値に対するa,xの相対出現回数はより同程度になる。
上記で留意したように、本ノイズ除去方法は、たとえノイズの結果であるxの値が相当数存在しても、xの値を変更しない。
というのは、種々のシーケンス値の相対頻度は、互いに類似するからである。
この場合、或るシンボルをすべて別のシンボルに置き直しても、信号の全劣化の改善は得られない。
この最適値は、経験的に決定されることができる。
入力信号中で相関しているシーケンスの長さが長時間にわたって著しく変化しない場合、kの最適値は、例示的な入力信号を利用し、kの種々の値に対するノイズ除去の結果を比較することによって、実験的に求められることができる。
このようなシステムでは、kの値は、満足な統計精度を与える値が見出される値まで、或る上限から減少させられることができる。
また、ノイズ除去アルゴリズムは、置き直すべき他のシンボルを発見するために、異なるkの値を用いて繰り返されることができる。
このような推定について、以下でさらに詳細に説明する。
しかし、他の形式のコンテクストも利用可能である。
例えば、シンボルxの両側にある一定シーケンスからなるコンテクスト、すなわちa,x,bが利用可能である。
同様に、シンボルの右側にある一定シーケンスからなるコンテクスト、すなわちx,bが利用可能である。
また、ベクトルaおよび/またはbが「ワイルドカード」を有するようなシーケンスも利用可能である。
すなわち、aは、a1,a2,...,aw,...,akの形式で書かれることができる。
ここで、awは、ストリング中のシンボルが任意の値を取り得るようなシンボルのストリングであることが可能である。
同様に、ベクトルは、処理されているシンボルに対して所定の位置にある限り、当該シンボルに隣接している必要はない。
なお、aとbとは異なる長さであってもよく、aまたはbのいずれかが欠けていても(すなわち長さ0であっても)よいことに留意すべきである。
しかし、本発明は、シンボルが2よりも多くの値を取るシーケンスにも適用可能である。
一般的な場合、各シンボルは、M個のシンボルを有するアルファベットから選択される。
すなわち、yj、zj、z'j、aj、およびbjは、値0,1,2,...,M−1を取り得る。
この場合、チャネルは、値Iのシンボルが値Jのシンボルにチャネルにより変換される確率を与える行列P(I,J)によって、特徴づけられる。
ここで、IおよびJはそれぞれ、0からM−1までになる。
同様に、シンボルIがシンボルJで置換される場合に信号に加わる劣化を測る所定の行列D(I,J)が存在すると仮定される。
処理されるべきシンボルzjが与えられた場合、当該シンボルと所定の関係にあるシンボルのシーケンスを調べることによって、当該シンボルのコンテクストを求める。
このコンテクストのシーケンスをaおよびbで表す。
0からM−1までのxのそれぞれの値について、出力信号中にa,x,bが現れる回数N'(x)をカウントする。
そして、確率行列Pを用いて、上記と類似の方法で、0からM−1までのxについて、入力信号中にa,x,bが現れる回数N(x)が推定される。
上記と同様に、これらの推定値で、次の方程式系(行列記法で表す)を満たすものをn(x)で表す。
N'=nP (c)
ここで、N'およびnは、M次元行ベクトルであり、その成分はそれぞれx=0〜M−1に対するカウントN'(x)およびn(x)からなる。
Pは、上記で定義したチャネル確率行列である。
したがって、n(x)は次のように表される。
n=N'P−1 (d)
ここで、P−1はPの逆行列である。
q=zjとして、xのそれぞれの可能な値について、a,q,bのすべてのインスタンスをa,x,bで置き換えた場合に生じる推定全信号劣化を、D(I,J)を用いて計算する。
x=qの場合は、シーケンスが不変のままである場合に対応する。
推定全信号劣化が最小となるxの値をpとして、zjをpで置き換える。
a,x,bからの寄与は、誤りの個数と、それぞれの誤りから生じる劣化との積である。
a,x,bから生じる誤ったa,q,bの個数は、高い確率で、N(x)P(x,q)に近く、これは、n(x)P(x,q)であると推定される。
それぞれの誤りから生じる劣化は、D(x,q)である。
したがって、何もしない場合の信号の推定劣化は、qと異なるすべてのxの値にわたるn(x)P(x,q)D(x,q)の総和である。
a,q,bのすべてのインスタンスがa,p,bに変換される場合に生じる劣化を考える。
a,q,bのインスタンスの一部は、出力信号コンテクストがa,bであるpのインスタンスを、チャネルがqに変換した結果である。
したがって、a,q,bのすべてのインスタンスをa,p,bに変換すれば、これらは補正される。
a,p,bに変換されたa,q,bの残りのインスタンスは、誤りとなる。
したがって、a,q,bをa,p,bに変換したことから生じる全劣化は、pとは異なるすべてのxの値にわたるn(x)P(x,q)D(x,p)の総和として推定される。
ここでも、n(x)P(x,q)は、N(x)P(x,q)の推定値であり、これは、入力信号中ではもともとxであった、コンテクストa,bを有するqの出現回数に、高い確率で近いことが期待される。
推定される劣化を最小にするpの値が、qを置き換えるために選択される。
最小にする値がqであることも十分あり得る。
その場合、a,q,bは、結局不変のままとなる。
本発明の好ましい一実施形態では、2パスシステムが用いられて、信号{zi}をノイズ除去する。
この例の目的上、各シンボルのコンテクストは、そのシンボルの左側のk個のシンボルと、そのシンボルの右側のk'個のシンボルであると仮定する。
第1パスで、コントローラ11は、シンボルが受信されるごとに、受信された{zi}を、メモリ13に記憶する。
また、コントローラ11は、{zi}中の長さK=k+k'+1のすべてのシーケンスのリストを作成する。
これらのシーケンスおよび各シーケンスに遭遇した回数が、コンテクストメモリ14に記憶される。
なお、受信信号は、必ずしも高速なメモリに記憶される必要はないことに留意すべきである。
任意の所与の時刻に、受信信号からK個のシンボルが、第1パス中に、コントローラ11によって必要とされる。
したがって、受信信号は、ディスクドライブに記憶されてもよい。
コンテクストメモリ14のみが、高速なメモリであればよい。
また、コントローラ11は、Kより小さい長さの各シーケンスについても、その識別およびカウントを記憶してもよい。
この情報は、好ましくは、行列Pとして、および行列Pの逆行列としての両方で記憶される。
というのは、チャネルによって、破損したシーケンス中のコンテクストa,bを有するシンボルに変換された、もとのシーケンス中のxの出現回数を推定するために、逆行列が必要となるからである。
本発明の好ましい実施形態では、この情報は、シンボルIがシンボルJで置換される場合に生じる劣化をI,J成分とする行列の形式で記憶される。
図3は、第2パスで用いられる信号処理アルゴリズムのフローチャートである。
受信信号がメモリ13に記憶された後、コントローラ11は、受信シンボルを順次調べて、シンボルを別の値に置き直すべきかどうかを判定する。
51に示すように、コントローラ11は、zjにある時、zjの左側のk個のシンボルと、zjの右側のk'個のシンボルとを読み出し、zjの最大のコンテクスト(そのカウントはメモリ14に記憶されている)を求める。
コントローラ11は、52に示すように、このコンテクストに関連するカウントをメモリ14から取り出し、53に示すように、記憶されているカウントが、上記の推定手続きによって、処理を進めるのに十分な統計精度を有する種々の入力シンボルカウントの推定をもたらすかどうかを判定する。
カウントが十分な精度を有する場合、54に示すように、コントローラ11は、破損した信号中の現在のコンテクストを有するもとの信号中の各シンボル値の出現回数を推定する。
次に、55に示すように、zjがそれぞれの可能なシンボル値で置換される場合に生じる信号劣化が推定される。
推定される信号劣化が最小となるxの値をxminで表す。
x=zjである場合は、zjが不変な場合に対応する。
次に、56に示すように、j番目の出力シンボルz'jをxminとする。
次に、57に示すように、jの値がインクリメントされ、コントローラ11は、同様にして次のシンボルの処理に進む。
例えば、現在のコンテクストが、zjの右側にk個のシンボルを含み、左側にk'個のシンボルを含む場合、右側に(k−1)個のシンボルおよび左側にk'個のシンボルのあるコンテクストが、テストされることができる。
メモリ14が長さK−1のストリングに対するカウントを含む場合、必要なデータが存在する。
より小さいコンテクストが見出された場合、コントローラ11は、60に示すように、まだ調べていないこの長さのコンテクストを取り出し、ステップ52から上記のプロセスを繰り返す。
より小さいコンテクストのデータがない場合、59に示すように、z'jは、zjとおかれ、コントローラ11は、次のシンボルに進む。
行列Pは、アルゴリズムの調整可能なパラメータとみなされることができ、Pの良好な選択は、いくつかの方法で得られることができる。
1つの選択肢は、結果として得られるノイズ除去された信号が十分な品質であると判断されるまで、Pを調整することによって、対話的にPを導出することである。
別法として、Pの選択は、チャネルの事前知識に基づくことも可能である。
場合によっては、入力信号中の対応するシンボルの真の値を強く示唆する、ある特定のコンテクストに現れるシンボルの分布に基づいて、破損した信号自体から直接Pを推定することが可能なことがある。
例えば、入力データが2値画像からなり、チャネルが2値チャネルである場合、すべて0のコンテクスト中に生じる1の割合は、チャネル変換確率P(0,1)の良好な推定値となる可能性が高い。
というのは、ほとんどすべてのインスタンスにおいて、対応する入力シンボルは0となるからである。
より大きいアルファベットの信号にこの技法を一般化することは、変動が最小の出力信号の領域を探索してから、これらの領域内の入力信号がほとんど一定であると仮定することに帰する。
所与のこのような領域に対して、この一定入力値の推定値をIとする。
これは、さまざまな方法(メジアン、平均等)で導出されることができる。
すると、Iに対応する領域に制限された出力信号の経験分布は、J=0,...,M−1に対するP(I,J)(これはPの第I行である)の妥当な推定値となる。
最終的に、上記の技法の組合せ(例えば、推定+事前知識+対話的実験)が最良の結果をもたらすであろう。
しかし、通常は、ノイズ除去問題は、自然な、広く受け入れられている劣化の尺度を有し、それを適切なD行列に反映させることができる。
しかし、これらの項目の2つ以上を記憶するために単一のメモリを利用する実施形態もまた、本発明の教示から逸脱することなく構成されうる。
しかし、本発明は、適当なデータ処理プログラムを汎用データ処理システムにロードすることによって、チャネルからの受信信号のコピーが転送された当該データ処理システム上で実施されることも可能である。
よって、本発明は、添付の特許請求の範囲のみによって限定されるべきである。
12、13・・・メモリ
14・・・コンテクストメモリ
20・・・チャネル
Claims (17)
- ノイズを導入するチャネル[20]によりノイズのない信号[21]から生成されたシンボルのシーケンスを含む受信信号[22]と、
前記チャネル[20]により、値Iのシンボルが値Jのシンボルに変換される確率を指定する情報と、
値Iのシンボルが値Jのシンボルに置換される結果として起こる信号劣化を測る信号劣化関数を指定する情報と
をそれぞれ記憶する第1、第2、および第3のメモリ[12、13]と、
前記受信信号[22]のシンボルを置換することにより、前記受信信号[22]から処理済み信号[23]を生成するコントローラ[11]と
を含み、
前記コントローラ[11]は、前記受信信号[22]の中の前記シンボルのコンテクスト[24]において値Iのシンボルそれぞれを値Jのシンボルに置換すると、前記ノイズのない信号[21]に関連する前記処理済み信号[23]において、前記信号劣化関数で測定される全信号劣化が減少すると推定される場合に、この置換を行い、
前記コンテクスト[24]は、一定値のシンボルのシーケンスであって、置換される前記シンボルに関して所定の位置にあるシンボルのシーケンスを含む
装置。 - 前記受信信号[22]の中の第1の長さのシーケンスそれぞれの識別と出現数とを有する第1のリストを記憶する第4のメモリ[14]をさらに含む
請求項1に記載の装置。 - 前記コントローラ[11]は、前記第2のメモリ[12]の情報を利用して、前記チャネル[20]により、前記受信信号[22]において、前記値Iを有し、前記コンテクスト[24]を有するシンボルに変換された前記ノイズのない信号[21]のシンボル値それぞれの出現数を推定する
請求項1に記載の装置。 - 前記コントローラ[11]は、
前記コンテクスト[24]のシンボルの値すべてをJに設定した場合の前記ノイズのない信号[21]に関する劣化を推定し、
この推定された劣化を最小化するJを選択する
請求項3に記載の装置。 - 前記第4のメモリ[14]は、
前記受信信号[22]の中の第2の長さのシーケンスそれぞれの識別と出現数とを有する第2のリストをさらに記憶し、
前記第2の長さは、前記第1の長さより短い
請求項2に記載の装置。 - 前記コンテクスト[24]は、前記第1のリストおよび前記第2のリストの1つから選ばれ、
この選択は、前記第1のリストおよび第2のリストの中のシーケンスの前記出現数に依存する
請求項5に記載の装置。 - チャネル[20]によりノイズのない信号[21]から生成されたシンボルのシーケンスを含む受信信号[22]を処理して、処理済み信号[23]を生成する方法であって、
前記受信信号[22]を記憶すること、
前記チャネル[20]により、値Iのシンボルが値Jのシンボルに変換される確率を指定するチャネル[20]破損情報を記憶すること、
値Iのシンボルが値Jのシンボルに変換される場合に起こる信号劣化を測る信号劣化関数を指定する情報を記憶すること、および、
前記受信信号の中の前記シンボルのコンテクスト[24]において値Iのシンボルそれぞれを値Jのシンボルに置換すると、前記ノイズのない信号[21]に関連する前記処理済み信号[23]において、前記信号劣化関数で測定される全信号劣化が減少すると推定される場合に、この置換を行うことによって、前記処理済み信号[23]を生成すること
を含み、
前記コンテクスト[24]は、一定値のシンボルのシーケンスであって、置換される前記シンボルに関して所定の位置にあるシンボルのシーケンスを含む
方法。 - 前記受信信号[22]の中の第1の長さのシーケンスそれぞれの識別と出現数とを有する第1のリストを記憶することをさらに含む
請求項7に記載の方法。 - 前記チャネル[20]破損情報は、チャネル[20]により、前記受信信号[22]において、前記値Iを有し、前記コンテクスト[24]を有するシンボルに変換された前記ノイズのない信号[21]のシンボル値それぞれの出現数を推定するために用いられる
請求項7に記載の方法。 - 前記劣化は、前記コンテクスト[24]のシンボルの値すべてをJに設定することによって、場合の前記ノイズのない信号[21]に関して推定され、
この推定された劣化を最小化するJを選択する
請求項9に記載の方法。 - 前記受信信号[22]の中の第2の長さのシーケンスそれぞれの識別と出現数とを有する第2のリストを記憶することをさらに含み、
前記第2の長さは、前記第1の長さより短い
請求項8に記載の方法。 - 前記コンテクスト[24]は、前記第1のリストおよび前記第2のリストの1つから選ばれ、
この選択は、前記第1のリストおよび第2のリストの中のシーケンスの前記出現数に依存する
請求項11に記載の方法。 - チャネル[20]により破損されたシンボルのシーケンスを含む受信信号[22]を処理して、処理済み信号[23]を生成するデータ処理プログラムで符号化されたコンピュータ可読媒体であって、
前記データ処理プログラムは、
前記受信信号を記憶すること、
前記チャネル[20]により、値Iのシンボルが値Jのシンボルに変換される確率を指定するチャネル[20]破損情報を記憶すること、
値Iのシンボルが値Jのシンボルに変換される場合に起こる信号劣化を測る信号劣化関数を指定する情報を記憶すること、および、
前記受信信号の中の前記シンボルのコンテクスト[24]において値Iのシンボルそれぞれを値Jのシンボルに置換すると、前記受信信号[21]に関連する前記処理済み信号[23]において、前記信号劣化関数で測定される全信号劣化が減少すると推定される場合に、この置換を行うことによって、前記処理済み信号[23]を生成すること
をデータ処理システムに実行させ、
前記コンテクスト[24]は、一定値のシンボルのシーケンスであって、置換される前記シンボルに関して所定の位置にあるシンボルのシーケンスを含む
コンピュータ可読媒体。 - 前記データ処理プログラムは、前記受信信号[22]の中の第1の長さのシーケンスそれぞれの識別と出現数とを有する第1のリストを記憶することを前記データ処理システムにさらに実行させる
請求項13に記載のコンピュータ可読媒体。 - 前記チャネル[20]破損情報は、チャネル[20]により、前記受信信号[22]において、前記値Iを有し、前記コンテクスト[24]を有するシンボルに変換された前記ノイズのない信号[21]のシンボル値それぞれの出現数を推定するために用いられる
請求項13に記載のコンピュータ可読媒体。 - 前記データ処理プログラムは、
前記受信信号[22]の中の第2の長さのシーケンスそれぞれの識別と出現数とを有する第2のリストを記憶することを前記データ処理システムにさらに実行させ、
前記第2の長さは、前記第1の長さより短い
請求項14に記載のコンピュータ可読媒体。 - 前記コンテクスト[24]は、前記第1のリストおよび前記第2のリストの1つから選ばれ、
この選択は、前記第1のリストおよび第2のリストの中のシーケンスの前記出現数に依存する
請求項16に記載のコンピュータ可読媒体。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/688,520 US7047472B2 (en) | 2003-10-17 | 2003-10-17 | Method for correcting noise errors in a digital signal |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2005124198A true JP2005124198A (ja) | 2005-05-12 |
JP3903052B2 JP3903052B2 (ja) | 2007-04-11 |
Family
ID=34465596
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004298490A Expired - Fee Related JP3903052B2 (ja) | 2003-10-17 | 2004-10-13 | デジタル信号中のノイズ誤りを訂正する方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US7047472B2 (ja) |
JP (1) | JP3903052B2 (ja) |
DE (1) | DE102004037857A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008504749A (ja) * | 2004-06-25 | 2008-02-14 | ヒューレット−パッカード デベロップメント カンパニー エル.ピー. | 信頼性情報による離散普遍雑音除去(discreteuniversaldenoising) |
US7912159B2 (en) | 2004-01-26 | 2011-03-22 | Hewlett-Packard Development Company, L.P. | Enhanced denoising system |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050162336A1 (en) * | 2004-01-27 | 2005-07-28 | Dell Products L.P. | System and apparatus for providing dual independent displays |
US7271749B2 (en) * | 2004-09-14 | 2007-09-18 | Hewlett-Packard Development Company, L.P. | Context-based denoiser that simultaneously updates probabilities for multiple contexts |
US7433427B2 (en) * | 2004-11-29 | 2008-10-07 | Hewlett-Packard Development Company, L.P. | Enhanced denoising system utilizing incremental parsing |
US7971131B1 (en) | 2005-05-06 | 2011-06-28 | Hewlett-Packard Development Company, L.P. | System and method for iterative denoising and error correction decoding |
US7434146B1 (en) * | 2005-05-06 | 2008-10-07 | Helwett-Packard Development Company, L.P. | Denoising and error correction for finite input, general output channel |
US7592936B2 (en) * | 2005-07-12 | 2009-09-22 | Hewlett-Packard Development Company, L.P. | Input distribution determination for denoising |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE69535027T2 (de) * | 1994-03-19 | 2007-01-11 | Sony Corp. | Optische Platte, Verfahren und Gerät zur Aufzeichnung und Wiedergabe von Informationen |
US5533033A (en) * | 1994-12-19 | 1996-07-02 | The United States Of America As Represented By The Director, National Security Agency | Device for and method of correcting errors in formatted modem transmissions |
US6161209A (en) * | 1997-03-28 | 2000-12-12 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of Industry Through The Communications Research Centre | Joint detector for multiple coded digital signals |
US6105158A (en) * | 1998-04-03 | 2000-08-15 | Lucent Technologies, Inc. | Screening for undetected errors in data transmission systems |
-
2003
- 2003-10-17 US US10/688,520 patent/US7047472B2/en not_active Expired - Fee Related
-
2004
- 2004-08-04 DE DE102004037857A patent/DE102004037857A1/de not_active Withdrawn
- 2004-10-13 JP JP2004298490A patent/JP3903052B2/ja not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7912159B2 (en) | 2004-01-26 | 2011-03-22 | Hewlett-Packard Development Company, L.P. | Enhanced denoising system |
JP2008504749A (ja) * | 2004-06-25 | 2008-02-14 | ヒューレット−パッカード デベロップメント カンパニー エル.ピー. | 信頼性情報による離散普遍雑音除去(discreteuniversaldenoising) |
Also Published As
Publication number | Publication date |
---|---|
JP3903052B2 (ja) | 2007-04-11 |
DE102004037857A1 (de) | 2005-05-19 |
US7047472B2 (en) | 2006-05-16 |
US20050097421A1 (en) | 2005-05-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AlSabhany et al. | Digital audio steganography: Systematic review, classification, and analysis of the current state of the art | |
US7271749B2 (en) | Context-based denoiser that simultaneously updates probabilities for multiple contexts | |
Wang et al. | An informed watermarking scheme using hidden Markov model in the wavelet domain | |
US7434146B1 (en) | Denoising and error correction for finite input, general output channel | |
JP2008504749A (ja) | 信頼性情報による離散普遍雑音除去(discreteuniversaldenoising) | |
US7498961B2 (en) | Context identification using a denoised signal | |
JP3903052B2 (ja) | デジタル信号中のノイズ誤りを訂正する方法 | |
JP4104552B2 (ja) | ウォーターマークデータの埋込み及び抽出 | |
US10666490B1 (en) | Method and apparatus for estimation of a bit rate of a digital signal | |
KR100737070B1 (ko) | 신호 디코딩 방법 및 디코딩 시스템과, 머신 액세스 가능한 기록 매체 | |
JP2003046395A (ja) | 積符号の復号方法および積符号の復号装置 | |
Han et al. | Erasure recovery matrices for encoder protection | |
US20090285410A1 (en) | Restoring audio signals | |
JP2008504748A (ja) | 誤り訂正符号化による離散普遍雑音除去(discreteuniversaldenoising) | |
Obidin et al. | Signal denoising with the use of the wavelet transform and the Kalman filter | |
JP3933668B2 (ja) | 信号処理装置および信号処理方法およびコンピュータ読取可能記憶媒体 | |
Hussain et al. | A novel wavelet thresholding method for adaptive image denoising | |
Sasaki et al. | Multiple matrix rank minimization approach to audio declipping | |
US7433427B2 (en) | Enhanced denoising system utilizing incremental parsing | |
Ulacha et al. | Highly effective predictor blending method for lossless image coding | |
CN111462765B (zh) | 一种基于一维卷积核的自适应音频复杂度表征方法 | |
Rashidha et al. | An adaptive-size median filter for impulse noise removal using neural network-based detector | |
CN117014034A (zh) | 脉冲干扰抑制方法及相关设备 | |
Abid et al. | Consistent reconstruction of the input of an oversampled filter bank from noisy subbands | |
Levy et al. | An image watermarking scheme based on information theoretic principles |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20061207 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070105 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 3903052 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110112 Year of fee payment: 4 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110112 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120112 Year of fee payment: 5 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130112 Year of fee payment: 6 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R360 | Written notification for declining of transfer of rights |
Free format text: JAPANESE INTERMEDIATE CODE: R360 |
|
R360 | Written notification for declining of transfer of rights |
Free format text: JAPANESE INTERMEDIATE CODE: R360 |
|
R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |