JP4468447B2 - 自給式呼吸装置システム内で発生した音声信号を等化するための方法および装置 - Google Patents

自給式呼吸装置システム内で発生した音声信号を等化するための方法および装置 Download PDF

Info

Publication number
JP4468447B2
JP4468447B2 JP2007519235A JP2007519235A JP4468447B2 JP 4468447 B2 JP4468447 B2 JP 4468447B2 JP 2007519235 A JP2007519235 A JP 2007519235A JP 2007519235 A JP2007519235 A JP 2007519235A JP 4468447 B2 JP4468447 B2 JP 4468447B2
Authority
JP
Japan
Prior art keywords
noise
audio signal
mask
noise model
lpc
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.)
Expired - Fee Related
Application number
JP2007519235A
Other languages
English (en)
Other versions
JP2008505538A (ja
Inventor
エム. カシュナー、ウィリアム
エム. ハートン、サラ
エイ. ジャシューク、マーク
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Motorola Solutions Inc
Original Assignee
Motorola Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Motorola Inc filed Critical Motorola Inc
Publication of JP2008505538A publication Critical patent/JP2008505538A/ja
Application granted granted Critical
Publication of JP4468447B2 publication Critical patent/JP4468447B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63CLAUNCHING, HAULING-OUT, OR DRY-DOCKING OF VESSELS; LIFE-SAVING IN WATER; EQUIPMENT FOR DWELLING OR WORKING UNDER WATER; MEANS FOR SALVAGING OR SEARCHING FOR UNDERWATER OBJECTS
    • B63C11/00Equipment for dwelling or working underwater; Means for searching for underwater objects
    • B63C11/02Divers' equipment
    • B63C11/26Communication means, e.g. means for signalling the presence of divers
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L15/00Speech recognition
    • G10L15/20Speech recognition techniques specially adapted for robustness in adverse environments, e.g. in noise, of stress induced speech
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0316Speech enhancement, e.g. noise reduction or echo cancellation by changing the amplitude
    • G10L21/0364Speech enhancement, e.g. noise reduction or echo cancellation by changing the amplitude for improving intelligibility
    • G10L2021/03643Diver speech
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208Noise filtering
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/78Detection of presence or absence of voice signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Ocean & Marine Engineering (AREA)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)
  • Percussion Or Vibration Massage (AREA)
  • Transmitters (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Description

本発明は、概して、通信システムと結合している加圧空気供給システムに関する。
消防のような危険な環境活動に従事しているスタッフ間での明瞭で信頼性の高い通信は、自分自身の健康および安全を維持しながらその使命を達成するために不可欠なものである。作業状況により、例えば、自給式呼吸装置(SCBA)マスクおよび空気供給システム、自給式水中呼吸装置(SCUBA)マスクおよび空気供給システム、または航空機用酸素マスク・システムのような加圧空気供給システムを使用しなければならない場合が出てくる。
米国特許出願第10/882,452号 米国特許出願第10/882,715号
スタッフがこのような加圧空気供給システムを使用している場合でも、明瞭で信頼性の高い通信を維持し、スタッフの健康および安全を効果的に監視することが望ましい。
図1は、通信システム130と結合している加圧空気供給システム110を含む従来技術によるシステム100の簡単なブロック図である。加圧空気供給システムは、通常、SCBAマスクのような呼吸マスク112と、空気シリンダ(図示せず)と、レギュレータ118と、レギュレータ118を空気シリンダに接続している高圧ホース120とを含む。使用する空気供給システム110のタイプにより、システム110は、例えば、ユーザに呼吸のためのクリーンな空気を供給し、有害な毒素がユーザの肺に達するのを防止し、ユーザの肺が燃焼している構造の内部の過熱した空気により火傷することから保護し、ユーザの肺を水から保護し、さらにユーザを顔面および呼吸器の火傷から保護することにより、ユーザを保護することができる。さらに、通常、マスクは圧力を必要とする呼吸システムであるとみなされる。何故なら、通常、マスクの着用者が息を吸い込むときだけ空気が供給されるからである。
通信システム130は、通常、マスクの着用者の音声を記録するように設計されていて、マスクの内側および外側に装着、マスクに取り付け、またはマスク112上の音声送信ポート上で手の中に保持することができる従来のマイクロホン132を含む。通信システム130は、さらに、マスクの着用者が、例えば他の通信装置に自分の音声を送るために使用することができる双方向無線機のような通信装置134を含む。マスクのマイクロホン・デバイス132は、直接、または中間電子処理装置138を通して無線機134に接続することができる。この接続は、従来のワイヤ・ケーブル(例えば136)を通して、または従来のRF、赤外線または超音波短距離送信機/受信機システムにより無線で行うことができる。中間電子処理装置138は、例えば、デジタル信号プロセッサとして実施することができ、インタフェース電子回路、音声増幅器、およびデバイス用およびマスク・マイクロホン用のバッテリー電力を含むことができる。
システム100のようなシステムは、その使用に関連するいくつかの問題点がある。これらの問題点については、説明を分かりやすくするために、図1のシステム100のマスク−無線機音声パスを示す図2のブロック図を参照しながら説明する。唇からの音声入力210(例えば、S(f))は、音響共鳴およびヌルを特徴とする音響伝達関数220(例えば、MSK(f))を有する、マスク(例えば、SCBAマスク)に入る。これらの共鳴およびヌルは、マスクの空洞容積およびマスクの内面からの音の反射によるものである。伝達関数MSK(f)を特徴とするこれらの効果は、入力音声波形S(f)を歪ませ、そのスペクトル内容を変化させる。もう1つの音源は、同様にマスクに入り、MSK(f)の影響を受ける呼吸装置(例えば、レギュレータ吸気ノイズ)が発生するノイズ230である。もう1つの伝達関数240(例えば、NP(f))が、ノイズが、音声のノイズではなく、マスク内の若干異なる位置から発生するノイズによるものであることを示している。音声およびノイズS(f)は、自分自身の伝達関数250(例えば、MIC(f))を有するマイクロホンにより音響エネルギーから電子信号に変換される。次に、マイクロホン信号は、通常、音声増幅器、およびまた、伝達関数260(例えば、MAA(f))を有する他の回路を通る。次に、MAA(f)からの出力信号270(例えば、S(f))は、無線機内に入力され、さらに処理され送信される。
システム100のようなシステムの問題点に戻って説明すると、このような問題点の一例は、その動作の一部としての大きな音響ノイズのこれらシステムによる発生に関連する。より詳細に説明すると、これらのノイズは、特に無線機のような電子システムと一緒に使用した場合、通信の品質を有意に劣化させる恐れがある。SCBAシステムのような加圧空気供給システムが導入する顕著な音声アーティファクトであるこのような1つのノイズとしては、図2にボックス230として示すレギュレータ吸気ノイズがある。
レギュレータ吸気ノイズは、マスクの着用者が息を吸い込む度に発生する広帯域のノイズ・バーストとして発生する。マスク内の圧力が負になると、エア・レギュレータの弁が開き、高圧の空気がマスク内に入ることができるようになり、大きなヒス音を発生する。このノイズは、後続の音声と一緒にマスク通信システムのマイクロホンによりピックアップされ、音声とほぼ同じエネルギーを有する。吸気ノイズは、通常、音声を遮断しない。何故なら、吸気ノイズは、通常、息を吸い込む時だけ発生するからである。しかし、このことがいくつかの問題を引き起こす恐れがあるが、そのいくつかの例について以下に説明する。例えば、吸気ノイズはVOX(音声作動スイッチ)回路をトリガする恐れがあり、それにより無線チャネルが開きまた占拠され、同じ無線チャネル上の他のスピーカと干渉を起こさないとも限らない。さらに、デジタル無線機を使用する通信システムにおいては、吸気ノイズがVAD(音声活動検出装置)アルゴリズムをトリガして、ノイズ抑制アルゴリズムでノイズ推定を無線信号処理チェーンのさらに下まで混乱させる恐れがある。さらに、吸気ノイズは、通常、聞き手にとって不快なものである。
システム100のようなシステムの第2の問題点について以下に説明する。これらシステムは、通常、鼻および口、または顔全体を覆うマスクを使用する。空気システムのマスクは、マスクの容積および内部反射面の幾何学的形状の関数である特定の組の音響共鳴および反共鳴(ヌル)を有し、マスク内で発生する音声のスペクトル特性を変化させる一定の幾何学的形状の閉空気空洞を形成する。より詳細に説明すると、エアマスク音声パス(図2)を特徴付ける際に、システムの最も難しい部分は、話す人の唇からマスクのマイクロホンまでの音響伝達関数(220)である。これらのスペクトルの歪みは、取り付けられている音声通信システム、特に、劣化した音声を処理するために最適化されないパラメータ的デジタル・コーデックを使用するシステムの性能を有意に劣化する恐れがある。音響マスク歪みは、パラメータ的デジタル・コーデックが関連する場合に、通信システムの品質および明瞭度に影響することが分かっている。一般的に、吸気ノイズの他に、音声品質を最も劣化させる空気システムの影響は、マスクの音響特性がよくないためによるものと思われる。
図3は、マスク(320)内およびマスク・マイクロホン出力(310)のところで測定したスペクトル振幅応答、およびマスク、マイクロホンおよびマイクロホン増幅器に対する結合伝達関数(330)の一例を示す。これらの特定のデータは、頭部および胴体のシミュレータ上に装着したSCBAマスクにより入手したものである。音響励起は、人工口部シミュレータを駆動する3Hz〜10KHzの掃引正弦波からなる。図3に示すように、スペクトルは、ほとんどの場合、マイクロホン内のプリアンプ帯域フィルタにより、500Hz未満および4.0KHzを超える周波数のところで有意に減衰していて、50〜4.0KHz間の有意の音声通過帯域領域内に多数の大きなスペクトル・ピークおよびノッチを含む。これらのスペクトル・ピークおよびノッチは、通常、櫛フィルタリングを行うマスク内の反射、および空洞共鳴条件によるものである。有意のスペクトル・ピークおよびノッチは、通過帯域を前後に移動するにつれて音声ピッチ成分およびフォルマントを変調し、その結果、音声の品質が劣化し、音声が歪む。音声の歪みを低減する目的で等化システムを定義するために使用するこのような伝達関数を含むこのようなシステムを特徴付ける1つまたは複数の伝達関数を決定するのが望ましい場合がある。
システムの伝達関数を適合できるように決定し、送信チャネルを等化するための多数の実績のある技術が存在する。システム伝達関数を決定するための1つの効果的な方法は、システムを励起し、システム・パラメータを決定するために、広帯域基準信号を使用する方法である。多くの音声送信環境で伝達関数を推定する際の1つの問題は、適当な広帯域の励起信号が容易に入手できないことである。1つの普通のアプローチは、長期平均音声スペクトルを基準として使用するアプローチである。しかし、この基準を使用する適応時間は、特に音声入力がまばらである場合には、長い時間がかかる場合がある。さらに、長期音声スペクトルは、音声スペクトルをかなり変える恐れがある叫び声および感情的ストレスを多くの場合含んでいる公共サービス活動の場合には、個人によりまた個人間でかなり異なる恐れがある。
システム100のようなシステムに関連するもう1つの問題点は、例えばバイオメトリック・パラメータを含むマスク着用者のいくつかのパラメータを測定するためのもっと効率的な方法および装置がないことである。システム100のようなシステムを使用しているかもしれない、危険な環境内で作業している個人のこのようなパラメータの測定は、これらの個人の安全および性能を監視するために重要なものである。例えば、個人の呼吸数および空気消費量の測定は、その人の作業量、生理的適性、ストレス・レベル、貯蔵空気供給の消費(すなわち、使用できる作業時間)を特徴付ける重要なパラメータである。呼吸を測定する従来の方法は、胸部インピーダンス体積記録器の使用またはサーミスタ・センサを使用する空気の流れの温度の測定を含む。しかし、消防のような物理的に要求の厳しい環境で作業している個人から、これらの従来の方法により、信頼性の高い測定値を入手することは、通常測定を行うために使用する身体に装着しているセンサおよびアーティファクトのズレを起こす恐れがある激しい物理的な運動によりさらに難しくなる。
それ故、吸気ノイズを効果的に検出し、減衰し、音声を等化し(すなわち歪みの影響を除去し)、通信システムと結合している加圧空気供給システムを含むシステムで、ユーザに関連するパラメータを測定するための方法および装置が求められている。
以下に添付の図面を参照しながら本発明の好ましい実施形態について説明するが、これは単に例示としてのものに過ぎない。
本発明を、多くの異なる形の実施形態で表すことができるが、図面および本明細書で特定の実施形態について詳細に説明する。この開示は本発明の原理の一例であると見なすべきであり、本発明を図に示すおよび本明細書に記載する特定の実施形態に限定するものではないことを理解されたい。さらに、本明細書で使用する用語および語は、本発明を制限するものではなく、単に説明のためのものであると見なすべきである。また、図面を簡単にし、分かりやすくするために図の素子の縮尺は必ずしも正確なものでないことも理解されたい。例えば、素子のうちのあるものの寸法は他の素子に対して誇張されている。さらに、そうするのがよいと考えた場合には、対応する素子を示すために、いくつかの図面で参照番号を反復して使用している。
本発明の種々の態様を詳細に説明する前に、本発明の理解を助けるために、すでに簡単に説明したエア・レギュレータの吸気ノイズをより詳細に説明することにする。吸気ノイズは、誰かが息を吸い込んでレギュレータ弁が開いた場合に、高圧の空気がSCBAまたは他の加圧空気供給システム・マスクに入るために発生するノイズである。SCBAマスク内に直接結合している弁のところで空気の乱流が起こると、非常に大きな広帯域のヒス・ノイズが発生するが、このノイズは音声信号によるマイクロホンのところの振幅にほぼ等しい。図4および図5は、それぞれSCBAマスク内で記録した通常の吸気ノイズ400の一例とその広帯域の分光写真500を示す。
図5を見れば分かるように、ノイズ・スペクトルは、約500、1700、2700および6000Hzのところに高いスペクトル・ピークを有する広帯域のスペクトルである。ピークは、マスク内の共鳴および内部マスク反射による櫛フィルタリングによるものであり、マスク・モデル、サイズおよび構成が変化すると、周波数および振幅が変化する。ノイズ・スペクトルの色は、通常、特定のマスク/着用者の組合せに対して一定である。何故なら、マスクを顔に着用すると内部幾何学的形状全体がほぼ一定になるからである。図6はそのことを示す。この図においては、所与のSCBAマスクを着用している同じ話し手に対するSCBAマスク・マイクロホンからの3つの別々の吸気ノイズ(610、620および630)のスペクトルは重畳している。この一致は、何人かの話し手についても、いくつかのメーカーからのマスクについても観察された。さらに、同じマスクを着用している何人かの話し手からのエア・レギュレータ・ノイズの高いスペクトルの類似性も観察された。
最後に、図7は、SCBAシステムから記録した音声710の例を示す。図7に示すように、吸気ノイズ720の影響は、音声自身の上には認められない。何故なら、人々は息を吸っている場合、通常話そうとしないからである。しかし、ノイズのエネルギーおよびスペクトルは十分大きく、無線機内の音声検出装置およびノイズ抑制回路といくつかの問題を起こし、聞き手にとってうるさいものである。
第1の態様においては、本発明は、図1に示すシステム100のような通信システムと結合している加圧空気供給システム内の吸気ノイズを検出および除去するための方法および装置である。本明細書においては、本発明のこの実施形態によるこの方法をARINA(エア・レギュレータ吸気ノイズ減衰器)方法とも呼ぶ。エア・レギュレータの吸気ノイズを識別し、除去するためのARINA方法の基礎は、音声と比較した場合、また例えば種々の環境ノイズのような他のタイプのノイズと比較した場合、ノイズが比較的一定であることである。図8は、ARINA方法800のブロック図を示すが、このブロック図は4つの部分、すなわち、ノイズ・モデル・マッチング810、ノイズ検出830、ノイズ減衰850、およびノイズ・モデル更新870に分割することができる。
ARINA方法800の基本的方法は、下記のように要約することができる。方法800は、好ましくはデジタル・フィルタ(例えば、全極線形予測符号化(LPC)デジタル・フィルタ)を使用する吸気ノイズをモデル化する。次に、方法800は、ノイズ・モデル・フィルタの倒置により、音声入力信号(すなわち、マスク・マイクロホンがピックアップした音声およびノイズ)をフィルタリングし、逆ノイズ・モデル・フィルタの出力のエネルギーを、入力信号または他のエネルギー基準のエネルギーと比較する。入力信号とモデルとの間で、よく似たスペクトルの一致が見られる信号周期中に、入力信号を含むレギュレータ吸気ノイズを任意の所望のレベルに減衰することができる。
ここで、図8に示すARINA方法800の特異性について説明すると、処理の際の第1のステップは、許容することができる実施態様の複雑さに依存して、図9または図10により好ましい実施形態で実施することができる方法800のノイズ・モデル・マッチング・セクション810により、入力信号802を基準ノイズ・モデルと連続的に比較することにより吸気ノイズの発生を検出することである。しかし、当業者であれば、他のスペクトル・マッチング方法も使用することができることを理解することができるだろう。本明細書においては、図9および図10に示す上記2つの好ましいマッチング方法は、正規化モデル誤差(またはNME)方法、および板倉−斉藤(またはI−S)歪み方法と呼ばれる。両方の方法の場合、基準ノイズ・モデルは、吸気ノイズのスペクトル特性を近似するデジタル・フィルタ(912、1012)により表される。好ましい実施形態の場合には、このモデルは、一組のLPC係数により指定される全極(自己回帰)フィルタとして表される。しかし、当業者であれば、全極モデルの代わりに、例えば、周知のARMA(自己回帰移動平均)モデルのような他のフィルタ・モデルも使用することができることを理解することができるだろう。
基準ノイズ・モデル・フィルタ係数は、吸気ノイズの少なくとも1つのデジタル化したサンプルから得られる、一組の自己相関係数から入手することができる。最初のノイズ・サンプルおよび対応する最初の自己相関係数(872)は、任意の数のノイズ予備記録からオフラインで入手することができ、本発明の実施態様にとって重要なものではない。さらに、実験により、例えば、1つのSCBAマスクからの最初のノイズ・サンプルも、同じ設計の他のマスクおよび場合によっては異なる設計のマスクに対してうまく動作することが分かっている。自己相関係数は、生のサンプリングしたノイズ・データから直接計算することができるし、または当業者にとって周知の普通の方法を使用して、LPCまたは反射係数のような他の通常使用するスペクトル・パラメータ表示から入手することができる。
好ましい実施形態の場合には、ノイズ・モデル自己相関係数は、下記の標準的な式により計算することができる。
Figure 0004468447
ここで、Rはp個の自己相関係数の最大のi番目の係数であり、xは、N個のサンプルの最大を含む通常の吸気ノイズ信号サンプル・セグメントのn番目のサンプルであり、Rは全セグメントのエネルギーを表す。自己相関関数pの次数は、通常、14である好ましい実施形態に対する値を含む10〜20である。さらに、理想的には、N個の信号サンプルは、スペクトル推定を平滑にするために自己相関を実行する前に、ハミング窓により窓処理することが好ましい。ハミング窓は下式により表される。
Figure 0004468447
通常の当業者であれば、他の窓処理方法も使用することができることを理解することができるだろう。
下式のz領域表示伝達関数を含む全極線形予測モデル・フィルタを表す一組の10番目の次数のノイズ・モデルLPC係数a,a,...,aを決定するために、ノイズ・モデル自己相関係数が次に使用される。
Figure 0004468447
ここで、z=e−jnωTはz変換変数である。この例の場合、10番目の次数のLPC係数が決定された。しかし、異なる次数のLPC係数も特定の実施態様に基づいて選択することができる。自己相関からLCPパラメータへの変換(ステップ912、1012)は、当業者にとって周知の任意の数のパラメータ変換技術により行うことができる。好ましい実施形態の場合には、LPCパラメータは、当業者であれば周知のダービン方法により自己相関パラメータから入手することができる。
ここで、図9のNMEスペクトル・マッチング方法の特異性について説明すると、得られた全極LPCノイズ・フィルタ・モデルは、逆LPCフィルタを形成するために倒置される(ステップ914)。
Figure 0004468447
理想的には、マスク・マイクロホンから入手し、音声および吸気ノイズS(z)を含む、低域フィルタリングし、サンプリングした音声入力信号802は、出力信号を入手するために逆フィルタ
Figure 0004468447
(ステップ914)を通過することが好ましい。
Figure 0004468447
次に、逆フィルタ入力および出力信号のエネルギーEin、Eoutが計算され(それぞれステップ918および916で)、歪み測定値Dが、ステップ920で計算され、ノイズ・モデルと入力信号との間の類似性の測定値として使用される。Dの理論的なもっと低い境界は無限の次数に対してゼロであるが、実際には、もっと低い境界および有限の次数のノイズ・モデルにどれくらいよく一致するかは入力信号により決まる。この実施態様の場合には、歪み測定値は、下式によりステップ920で計算される正規化したモデル誤差(NME)と呼ばれるEinに対するEoutの比率により定義される。
Figure 0004468447
次に、入力信号のエネルギーを、ノイズ・モデルにどのくらいよく似ているかにより除去することができる。好ましい実施形態の場合には、上記信号フィルタリングは、時間領域内のコンボルーションにより行われるが、この信号フィルタリングは上式に示すように周波数領域内でも行うことができる。
ARINA方法800に対する信号処理は、通常、セグメント分割したフレームをベースとして行われる。好ましい実施形態の場合には、入力信号802は、8.0KHzで低域フィルタリングされ、サンプリングされ、80のサンプル(10ミリ秒)のブロック内にバッファされ、逆ノイズ・モデル・フィルタを通る(式5)。それ故、理想的には、すべてのフィルタリングは、連続的に入力信号802の80サンプル・セグメントで行うことが望ましい。次に、逆ノイズ・モデル・フィルタの正規化したモデル誤差(NME)が、フィルタ出力フレーム・エネルギーを入力信号フレーム・エネルギーで割ることにより計算される(式6)。しかし、理想的には、この計算は時間の解像度をよりよくするためにサブフレームに基づいて行うことが望ましい。それ故、各80点フレームは、サブフレーム、例えば、4、20点サブフレームに分割されるが、必要な精度の程度により、他のサブフレーム分割も使用することができる。次に、全正規化モデル誤差信号(NME)を、最後の16のサブフレームの出力フィルタ・エネルギーEoutを平均し、その値を対応する時間整合した16のサブフレーム入力フィルタ・エネルギーEinの平均で割ることにより平滑にすることができる。このようにしても、分析は遅れないが、遷移ドロップアウト、およびレギュレータ・ノイズ・スペクトルを変える恐れがある他の大きな背景ノイズの影響を容易に除去することができる。それにより、本発明のこの実施態様において、平均NME値が、入力信号スペクトルの類似性へのノイズ・モデルの測定値として使用される。
好ましい実施形態の場合には、図10の第2のもっと複雑な、しかしもっと正確なノイズ・モデル・マッチング方法810は、板倉−斉藤歪み方法を修正したものである。2つの信号間のスペクトルの類似性を決定するためのI−S方法は、当業者にとって周知のものである。この方法の場合、残留ノイズ・モデル逆フィルタ・エネルギーが、上記NME方法の場合のように、入力信号エネルギーとではなく、「最適」信号フィルタの残存エネルギーと比較される。フィルタは、現在の信号のセグメントのスペクトルと最もよく一致するという意味で「最適」なものである。
最適にフィルタリングした信号に対応する残存エネルギーはステップ1018〜1024で計算される。I−S方法の場合には、ステップ1018において、理想的には、入力信号802の2つの連続している80のサンプル・バッファを1つの160のサンプル・セグメントに結合することが望ましい。好ましくは、下式による160点ハミング窓により、160のサンプル・セグメントに窓処理する。
Figure 0004468447
次に、窓処理した信号データは、式1のところで説明した方法により自己相関がつけられる。ステップ1018で生成したこれらの自己相関係数は
Figure 0004468447
で表される。対応する一組のLPC係数は、ステップ1012において基準ノイズ・モデル・パラメータを生成するために使用したのと同じ方法で、ステップ1020において、好ましくはダービン・アルゴリズムにより、自己相関係数から入手する。ステップ1020において生成した信号モデルLPC係数は
Figure 0004468447
で表される。ステップ1022において、これらのLPC係数(ステップ1020)は下記の式9により自己相関がつけられ、
Figure 0004468447
となる。これらのパラメータを使用して、このフィルタを通過する信号Eの残存エネルギーが下式によりステップ1024において計算される。
Figure 0004468447
ノイズ・モデルを通過する入力信号のエネルギーは、ステップ1012〜1016において計算される。ステップ1012において、ノイズ・モデルLPC係数が、すでに説明したように、ノイズ・モデル自己相関係数(874)から計算される。ステップ1012において生成したこれらLPC係数は、a(i=1,2,...,p)で表される。ステップ1014において、(ステップ1012からの)LPC係数が、下式11により自己相関付けられ、bとなる。ステップ1018において計算したこれらのパラメータおよび自己相関シーケンスを使用して、基準ノイズ・モデルを通過する信号のエネルギー
Figure 0004468447
が、式10によりステップ1016において計算される。
Figure 0004468447
基準ノイズ・モデルに対する「最適」信号モデルのスペクトル歪みDの測定値が、下式に定義するように、ステップ1028において計算される。
Figure 0004468447
信号モデルが基準ノイズ・モデルに似れば似るほど、歪み測定値は下限である1.0に近づく。この歪み測定値は、吸気ノイズの有無を判定するために、ARINA方法800のノイズ検出セクション830により使用される。好ましい実施形態の場合には、I−S歪み測定値は、160のサンプルにより計算される。I−S歪み測定値により決定した吸気ノイズ分類は、160のサンプル・セグメントの各80のサンプル・フレームと関連する。さらに、最初のノイズ・モデル(例えば、最初の自己相関係数872に基づいて)を生成したり、すでに説明し、以下に詳細に説明するノイズ・モデル更新セクション870によりノイズ・モデルを更新するときだけステップ1012および1014を実行するだけでよい。
次に、ARINA方法800のノイズ検出部830において、スペクトル・マッチング810から入手した値(すなわち、入力信号とノイズ・モデル間の類似性の測定値を表すNMEまたはI−S歪み測定値)が、経験的に入手した閾値(例えば、Dmin1)と比較される(ステップ832)。この検出閾値は、音声または他のタイプのノイズを吸気ノイズと間違って分類しないようにしながら、吸気ノイズの有無を検出するために選択される。
さらに、ノイズ・フィルタ・モデルの特異性により、吸気ノイズのスペクトルの変動、および例えば、間違った検出のようなノイズ・モデルに対するいくつかの音声の音の類似性が発生する場合がある。それ故、本物のエア・レギュレータの吸気ノイズの持続時間が、かなり長い間音声アーティファクトと比較されるので、理想的には、ノイズ持続時間閾値試験も適用することが好ましい(ステップ834)。それ故、検出が確認される前に、検出閾値は、所定の数の連続しているフレーム「K」(例えば、4フレーム)に適合しなければならない。相対的信号エネルギー、波形0交差、および他の特徴パラメータ情報は、音声/吸気ノイズ識別を改善するために検出スキームに内蔵させることができる。それ故、両方の閾値基準が満たされた場合(ステップ832および834から)、スペクトルの一致は許容できるほどに近いと見なされ、吸気ノイズが現在存在するものと見なされる。
ARINA方法800のノイズ減衰部850においては、ノイズ検出部830の出力が、入力信号802が通過する出力信号乗算器(852)をゲート制御するために使用される。吸気ノイズを検出した場合には、ステップ854において、乗算器の利得Gがある所望の減衰値「Gmin」に設定される。この減衰の利得値は、ノイズを完全に除去するためには0.0にすることができ、吸気ノイズを完全に除去しないでこのノイズを抑制するだけでよいならもっと高い値に設定することができる。聞き手がエア・レギュレータが機能していることを確実にすることができるように、このノイズを完全に抑制しないことが望ましい。好ましい実施形態の場合には、Gminの値は0.05である。そうでない場合で、吸気ノイズが検出されない場合には、理想的には、利得Gは、音声信号を減衰しないように1.0に設定することが望ましい。このゲート制御/乗算スキームの変動を使用することもできる。例えば、ゲート制御の攻撃および劣化をもっと緩やかなものにすることができる変動を使用することができ、吸気ノイズの前後に直接発生する場合がある音声の減衰の可能性を低減し、それにより耳で聞いた音声の品質を改善することができる。さらに、方法800から容易に理解することができるように、本発明の重要な利点は、従来の連続ノイズ・フィルタリング方法とは異なり、レギュレータ・ノイズを検出した場合を除いてもとの信号が変化しないことである。
ARINA方法800の重要な構成要素は、検出のためにノイズ・モデルを周期的に更新することができることである。例えば、時間の経過中に顔の上でのエアマスクが移動すると、音響伝達関数に対するその影響が変化する。また、異なる人々が着用しているエアマスク、または異なるマスクを使用するということは、最初の基準ノイズ・モデルのスペクトルが実際の吸気ノイズのスペクトルからズレることを意味する。もとの基準ノイズ・モデルを周期的に更新することにより、正確な現在の基準ノイズ・モデルを維持することができる。それ故、ARINA方法800のノイズ・モデル更新セクション870は、ノイズ・モデルを更新するために使用される。
ノイズ・モデル更新セクション870は、レギュレータ吸気ノイズの基準LPCフィルタ・モデルをいつ更新すべきかを決定するためにノイズ検出セクション830の出力を使用する。例えば、ノイズ・モデルを更新するのかしないのかを決定するために、ステップ876において、ノイズ検出セクション830からの出力を、第2の経験的に決定した閾値(例えば、Dmin2)と比較することができる。閾値が一致した場合には、吸気ノイズとして検出した多数の連続しているサブフレームをカウントすることができ(ステップ878)、バッファ内に格納している各サブフレーム内の信号サンプルをカウントすることができる。連続しているノイズのサブフレームの数が、ある閾値「K」(例えば、好ましい実施形態の場合の8サブフレーム、160サンプル)を超えると、ステップ880においてノイズ・モデルを更新する決定が行われる。ノイズでないサブフレームが検出された場合には(例えば、ステップ832、834および876のいずれかのところで)、ノイズ・フレーム・カウントは、ステップ884においてゼロにリセットされ、ノイズ・フレーム・カウントがステップ878において更新される。次に、現在検出した吸気ノイズを表す「K」連続している信号サブフレームに対する自己相関係数を、上記式1および式2によりステップ882のところで計算することができる。
これらの新しい自己相関係数は、ステップ874において、ノイズ・モデル、自己相関係数を更新するために使用される。理想的には、ステップ882において計算した自己相関係数を、例えば、下式のような簡単な加重式により、ステップ874において前のノイズ・モデル自己相関係数で平均することが望ましい。
Figure 0004468447
ここで、
Figure 0004468447
は、現在の基準ノイズ・モデルの自己相関係数であり、
Figure 0004468447
は、現在検出している吸気ノイズ・サンプルの自己相関係数であり、αは、最初の基準モデルの更新速度を決定する1.0〜0.0の加重係数である。この加重係数は、すでに説明した通常は遅い吸気ノイズのスペクトル特性の変化速度により調整することができる。次に、ノイズ・モデル逆フィルタに対する新しい一組のLPC係数が、ステップ912および1012において更新したモデルの自己相関から再計算される。間違った検出によりノイズ・モデルから大きくズレることがないように、ノイズ・モデルに対する調整を制限することができる。さらに、必要な場合には、システムを最初のモデル状態にリセットすることができるように、最初の基準ノイズ・モデル係数(872)が格納される。ノイズ・モデル更新セクション870のところですでに説明した方法800の適応機能により、システムは、特定のマスクおよびレギュレータの特性に適合することができ、最適の検出性能にすることができる。
ARINA方法800の利点としては、従来の連続的フィルタリングを使用するアルゴリズムの場合のように、処理アルゴリズムにより音声信号自身が非可逆的な影響を受けないこと等がある。もう1つの利点は、ここで使用するLPCモデル化を簡単にリアルタイムで容易に適合させることができることであり、簡単で効率的に計算することができることである。当業者であれば、上記利点が本発明のARINA実施形態に関連するすべての利点を含んでいることを意味するものではなく、代表的な利点を表すためのものであることを理解することができるだろう。
第2の態様においては、本発明は、図1のシステム100のような通信システムと結合している加圧空気供給システムで音声信号を等化するための方法および装置である。本明細書においては、本発明の実施形態によるこの方法をAMSE(エアマスク音声イコライザ)方法とも呼ぶ。等化のためのAMSE方法の基礎は、音声と比較した場合、また例えば、種々の環境ノイズのような他のタイプのノイズと比較した場合、比較的一定していることである。同じマスクの共鳴条件は、レギュレータ・ノイズおよび音声信号の両方に影響を与えるので、ノイズを等化すると、音声信号の等化に適しているイコライザができる。しかし、音の反射によるピークおよびヌルは、音声とノイズとの間のもとの位置の違いによりノイズと音声との間で若干異なる。
AMSE方法は、マスク空洞および構造が発生する音響共鳴スペクトルのピークおよびヌル(すなわち、スペクトル振幅音響伝達関数)を推定するために、すべてのマスク・タイプの加圧空気呼吸システム(例えば、SCBA)が含む広帯域エア・レギュレータ吸気ノイズを使用する。次に、このスペクトル情報は、リアルタイムで補償デジタル逆フィルタを形成するために使用される。このフィルタは、スペクトルが歪んでいる音声信号を等化し、マスクを使用しない場合に発生する歪みのない音声を近似する出力信号を生成するために適用される。この適用により、マスク・マイクロホンからの音声の品質が改善し、通信の明瞭度が改善する。
ここで、AMSE方法の特異性について説明する。図11はこの方法1100のブロック図であり、このブロック図は、4つのセクション、すなわち、ノイズ・モデル・マッチング・セクション1110、ノイズ検出セクション1130、マスク音声等化セクション1150、およびノイズ・モデル更新セクション1170に分割することができる。理想的には、AMSE方法のノイズ・モデル・マッチング・セクション、ノイズ検出セクション、およびノイズ・モデル更新セクションは、すでに詳細に説明したARINA方法の対応するセクションと同じものであることが望ましい。それ故、説明を簡単にするために、これら3つのセクションの詳細な説明はここでは省略する。しかし、AMSE方法1100のマスク音声等化セクション1150(鎖線エリア内の)については以下に詳細に説明する。
AMSE方法1100の音声等化セクション1150により、吸気ノイズ基準自己相関係数が、上式3によりステップ1152において、ノイズのn番目の次数のLPCモデルを生成するために使用される。ステップ1152において生成されたLPCモデルは、例えば図2のMSK(f)のようなマスクの伝達関数を特徴とし、また、吸気ノイズに対してノイズ経路伝達関数NP(f)を含む。好ましくは、14番目の次数のモデルが適しているが、他の任意の次数のモデルも使用することができる。当業者であれば、例えば、周知のARMA(自己回帰移動平均)モデルのような全極モデルの代わりに、他のフィルタ・モデルも使用することができることを理解することができるだろう。さらに、本発明の好ましい実施形態のところですでに説明した時間領域フィルタリング動作の代わりに、周波数領域内でフィルタリング動作を行うこともできる。
次に、好ましくは、ステップ1156において音声信号が通過する(式4による)逆フィルタ内でLPCモデル係数を使用する。逆フィルタを音声信号が通過すると、入力信号が効果的に等化され、それにより図2のマスク伝達関数MSK(f)によるスペクトル歪み(ピークおよびノッチ)が除去される。適切な固定の後フィルタを使用するステップ1158におけるポストフィルタリングは、理想的には、吸気ノイズの任意の白でない部分を修正し、または後置の特定のコーデックまたは無線機の要件を最適に満足させる目的で音声信号を特定のトーン品質にするために、等化した信号上で行うことが望ましい。また、ポストフィルタリングは、図2のノイズ経路伝達関数NP(f)を補償するためにも使用することができる。
図12は、2つの異なる次数の等化フィルタの場合の、エア・レギュレータ・ノイズに対するAMSE方法800のイコライザの影響を示す。より詳細に説明すると、図12は、等化前の吸気ノイズ・バーストのスペクトル表示1210である。さらに、この図は、14番目の次数の等化フィルタ(1220)および20番目の次数の等化フィルタ(1230)による等化後の吸気ノイズのスペクトルも示す。図を見れば分かるように、スペクトルのピークは、20番目の次数の等化フィルタにより非常にうまく平らになっていて、14番目の次数の等化フィルタによりかなりうまく平らになっている。さらに、これらのフィルタにより等化したマスク音声の聴取試験により、等化しなかった音声と比較した場合、音声の品質が等化フィルタの使用により有意に改善されたことが分かった。さらに、2つのフィルタの次数の間に聴取した音声の品質の僅かな違いが認められた。
AMSEアルゴリズム・アプローチの利点としては、(1)マスク音響共鳴特性を判定するために、励起源としてエアマスク・システムに固有の一定のスペクトル的に安定な広帯域のレギュレータ・ノイズを使用していること、(2)簡単で、しっかり確立された効率的な技術により、リアルタイムでシステム伝達関数がモデル化されていること、(3)同じ効率的な技術で等化がリアルタイムで行われていること、(4)変化する状態にシステムの伝達関数モデルがリアルタイムで連続的に適合することができること等がある。通常の当業者であれば、上記利点は、本発明のAMSE実施形態に関連するすべての利点を含んでいるわけではなく、代表的なものを示すものであることを理解することができるだろう。
第3の態様においては、本発明は、吸気ノイズの持続時間および周波数を決定し、図1のシステム100のような通信システムと結合している加圧空気供給システムでの呼吸数および空気使用量を決定するための方法および装置である。本明細書においては、本発明のこの実施形態による方法はINRRA(吸気ノイズ呼吸器数分析装置)方法とも呼ばれる。このINRRA方法は、本質的には、人の呼吸音を測定する代わりに、エア・レギュレータが発生する音を監視することにより呼吸を測定する間接的な方法である。INRRA方法の基本は、SCBAのような加圧空気呼吸システムの空気は1つの方向にしか流れないことである。空気は、空気源およびレギュレータからだけシステム内に入り、排気弁を通してだけ外部にでることができる。吸気弁および排気弁を同時に開くことはできない。それ故、レギュレータの吸気弁の動作は、ユーザの呼吸サイクルに直接関連している。
レギュレータの吸気弁が開いたことの1つのインジケータは、レギュレータの吸気ノイズである。吸気ノイズは、SCBAまたは他の加圧空気供給システムのマスクにより高圧の空気が入った結果である。マスクは気密に維持されているため、人が空気を吸い込むと、マスク内に僅かな負の圧力が発生し、そのためレギュレータの弁が開き、加圧タンクの空気が入る。SCBAマスク内と直接結合している弁の前後の空気の乱流は、大きな広帯域のヒス・ノイズを発生し、マイクロホンでピックアップすることができ、また空気を吸い込む毎に発生する。すでに説明したように、ノイズは急激なもので、空気を吸い込んでいる間、非常に安定している振幅を有し、開始時間および終了時間が非常にはっきり分かる。所与のマスク・タイプおよび着用者の場合、開いた口の大きさ、発声管の状態、肺の空気の流れのようないくつかの要因によりかなり変化する人間の直接の呼吸音とは対照的に、吸気ノイズのスペクトル特性は非常に安定している。INRRAは、呼吸数の測定基準として、エア・レギュレータの吸気ノイズの安定性を利用している。
INRRAは、その全スペクトル特性により吸気ノイズの有無を識別するためにマッチド・フィルタリング・スキームを使用する。さらに、INRRAは、万一発生した場合、ノイズのスペクトル特性の変化に適合することができるので、吸気ノイズと他の音とを最もよく区別する。各吸気のスタートを計算することにより、発生した吸気ノイズから瞬間的な呼吸数およびその時間平均を容易に計算することができる。さらに、終わりを測定し、各吸気ノイズの持続時間を計算することにより、また予測できるマスク・レギュレータの流速についてのいくつかの情報を供給することにより、システムは空気の流れの容積を推定することができる。この推定は、吸気ノイズを記録しているマイクロホンからの信号だけにより行うことができる。
図13は、INRRA方法1300のブロック図である。このブロック図は、5つのセクション、すなわち、ノイズ・モデル・マッチング・セクション1310、ノイズ検出セクション1330、吸気定義装置セクション1350、パラメータ推定装置セクション1370、およびノイズ・モデル更新セクション1390に分割することができる。理想的には、INRRA方法のノイズ・モデル・マッチング・セクション、ノイズ検出セクションおよびノイズ・モデル更新セクションは、すでに詳細に説明したARINA方法の対応するセクションと同じものであることが望ましい。それ故、説明を簡単にするために、これら3つのセクションの詳細な説明はここでは省略する。しかし、INRRA方法1300の吸気定義装置セクション1350およびパラメータ推定装置セクション1370については以下に詳細に説明する。
最初に、吸気定義装置セクション1350について説明する。INRRA方法1300のセクション1350の目的は、少なくとも1つの要因に基づいて、例えば、この場合には、吸気に対応する1つまたは複数の完全な吸気ノイズ・バーストの一組の終点および持続時間に基づいて吸気ノイズを特徴付けることである。吸気ノイズ検出セクション1330の決定は、ステップ1352において、1および0の値を使用する時間指数mの関数としての吸気ノイズの有無を表す、好ましくは2進信号、INM(m=0,1,2,...,M−1)を発生するために使用される。この2進信号は、長さMのサンプルの回転バッファに格納される。Mは、少なくとも2つの吸気ノイズ・バースト、または予想される最も遅い呼吸数での息の時間周期を収容するために、2進信号の十分なサンプルを格納するのに十分な大きさを有する。好ましい実施形態の場合には、Mは約15秒である。この2進信号およびMの値の時間分解能は、吸気ノイズ・モデル・マッチング・セクションに依存するすでに説明した吸気ノイズ検出セクション1330で使用する最も小さいサブフレーム時間により決まり、20サンプル(2.5ミリ秒)または80サンプル(10ミリ秒)であり、これによりステップ1310で使用するスペクトル・マッチング方法が決まる。
1330からの吸気ノイズ検出装置の出力は、いつでも完全なものであるというわけではないので、吸気ノイズの検出中に検出ミスが起こる場合があり、本当のスタート、およびノイズの持続時間がある程度曖昧になる場合がある。それ故、ステップ1352において発生した2進吸気ノイズ信号は、ステップ1354において、周知の移動平均タイプまたは他の適当なフィルタにより統合される。このフィルタは、任意の短い持続時間の検出誤差を平滑化し、呼吸に対応する完全な吸気ノイズ・バーストを定義するもっと正確な信号を発生する。ステップ1354で発生したこの信号から、各ノイズ・バーストに対する正確なスタート時間S、終了時間E、および呼吸持続時間Dを含む少なくとも1つの要因を、ステップ1356において、処理フレーム持続時間精度内で決定することができる。2進信号INMで表す吸気ノイズ・バーストのスタート時間および終了時間は、信号バッファ内のその相対的指数を記録することにより入手される。持続時間Dは1つの吸気ノイズ・バーストに対して下式により定義される。
Figure 0004468447
ここで、iは、長さMおよび時間周期T秒の2進信号バッファに含まれているi番目のI吸気ノイズ・バーストを示す。理想的には、これらの吸気ノイズ・バースト要因の値は、ノイズ・バースト/呼吸毎に一組のパラメータで、回転有限長さバッファに格納することが望ましい。図14〜図17は、INRRAアルゴリズムのセクション1310、1330、1352および1354により処理したSCBAマスク・マイクロホン音声の結果を示す。この結果はSCBAを着用している男の話し手からの音声に基づくもので、静かな室内で記録したものである。図14は、ノイズ・バースト1410を含む入力音声1420を示す。図15は、吸気ノイズ・モデル・マッチング・セクション1310のスペクトル歪み測定値D出力の時間−振幅図1500である。図16は、吸気ノイズ検出装置1330の2進出力の時間−振幅図1600である。図17は、生の検出装置出力を統合し、各吸気の持続時間を正確に定義する呼吸定義装置アルゴリズム1350の移動平均フィルタ構成要素1354の出力の時間−振幅図1700である。
パラメータ推定装置1370セクションは、吸気定義装置セクション1350による吸気ノイズの特徴付け要因に基づいて推定することができるパラメータのいくつかの例を示す。決定することができるパラメータの2つのこのような例は、ユーザの呼吸数であり、大体の吸気の流れの量である。呼吸数は、吸気定義装置セクション内で判定することができる連続している吸気ノイズ・バーストのシーケンシャルなスタート時間情報Sにより容易に決定することができる。例えば、1分間の「瞬間的」呼吸数は、下式により計算することができる。
Figure 0004468447
ここで、Sは秒単位による2つの連続しているノイズ・バースト(吸気)である。それ故、平均呼吸数は、下式により計算することができる。
Figure 0004468447
ここで、Iは、指定の時間T中に検出した連続している呼吸(吸気ノイズ・バースト)の数である。
吸気中の大体の空気の流れの量は、例えば、吸気定義装置セクションにより決定することができる呼吸の持続時間から、初期エア・タンク充填圧力に関するいくつかの追加情報から、およびオフラインで決定することができるレギュレータの平均流量から推定することができる。吸気弁が開くと、空気供給タンク圧力が、エア・レギュレータの最低入力圧力レベル以上である場合には、顔マスクにほぼ一定の圧力で(周囲の空気/水圧の関数)ある量の空気がエア・レギュレータ内に流入することができる。さらに、マスク・レギュレータの吸気弁が開いている間は、マスク内への空気の流入速度はほぼ一定である。それ故、タンク供給源から出て呼吸装置に供給される空気の量は、吸気弁が開いている時間に比例する。弁が開いている時間は、各吸気ノイズの持続時間により測定することができる。
充填した場合、供給タンク内の最初の空気の量は、タンク容量V、充填圧力P、ガス温度T、および汎用ガス定数R、モル内のガスの質量Nの関数であり、周知の理想ガス式、PV=NRTにより計算することができる。初期充填圧力およびタンク・シリンダ容積は分かっているので、タンク・ガスの温度およびマスク・ガスの温度が同じであると仮定した場合、マスク圧力で呼吸するために使用することができる空気の容量は、下式で表すことができる。
Figure 0004468447
それ故、吸気イベントi中にユーザに供給される空気の大体の容量は、下式で表される。
Figure 0004468447
ここで、IVは空気の容積、Dは吸気ノイズから判定した吸気イベントの持続時間、Kは特定のエア・レギュレータに対する空気の流速に関連する校正係数である。Kは、個々のシステムに対して経験的に入手することができ、おそらくメーカーのデータから決定することができる。個々の吸気容積IVから、時間Tまでに使用される空気の大体の全量Vは下式のように定義することができる。
Figure 0004468447
ここで、Iは時間Tまでの吸気の全数である。それ故、残りのタンクの供給空気は下式のようになる。
Figure 0004468447
INRRA方法の利点としては、最小の音声帯域幅上で呼吸ノイズをピックアップするために任意のマイクロホン信号を使用することができること、および特殊なセンサを必要としないこと等がある。もう1つの利点は、呼吸検出装置が、安定しているスペクトル特性を有し、特性が変化する人間の呼吸ノイズを含んでいないエア・レギュレータが発生するノイズの検出をベースとしていることである。さらにもう1つの利点は、呼吸検出装置が、他のタイプの音響呼吸分析装置のように、特定のチェック周波数にロックされないことである。さらに、システムは、環境の変化、および異なるユーザ、および加圧空気呼吸装置マスク・システムに自動的に適合する。それ故、INRRA方法は、連続的、瞬間的、または平均の呼吸数、および例えば、モニタへの無線データ・チャネルを介してシステム100の外部に自動的に送信することができる貴重な情報である大体の空気使用量データを供給することができる。通常の当業者であれば、上記利点は、本発明のINRRA実施形態に関連するすべての利点を含んでいるのではなくて、利点の代表的なものに過ぎないことを理解することができるだろう。
本発明による3つのすべての方法(ARINA、AMSEおよびINRRA)は、好適には、(上記システム100によるシステムが内蔵している)メモリ・デバイス上に格納しているソフトウェア・アルゴリズムとして実施し、そのステップは、例えば、システム100のDSP138のような適切な処理装置で実施することが好ましい。本発明の自己相関およびLPCフィルタリング方法に対応するアルゴリズムは、処理時間のほとんどを使用する。しかし、別の方法としては、これらのアルゴリズム、またはARINA、AMSEおよびINRRA方法に対応するアルゴリズムのエンティティは、小さいハードウェア領域内で効率的に実施することができる。さらに、AMSE方法は、ARINA方法のように多くの方法を使用しているので、本発明の他の実施形態の場合には、これらの方法を効率的に結合することができる。
特定の実施形態により本発明を説明してきたが、当業者であれば他の利点および修正を容易に思い付くだろう。それ故、広義の意味で、本発明は、特定の詳細、代表的な装置、および図に示し、本明細書に記載する例示としての例に限定されない。上記説明を読めば、当業者であれば種々の変更、修正および変更を思い付くだろう。例えば、吸気ノイズを識別し、減衰するための方法についてはすでに説明したが、本発明に関連する方法は、呼気ノイズのような他のタイプのノイズ、または上記方法を使用する効率的な検出に役に立つ疑似定常スペクトル特性を有する他のタイプのノイズに適用することができる。それ故、本発明が上記説明により制限されないで、添付の特許請求の範囲の精神および範囲によるこのようなすべての変更、修正および変形を含むことを理解されたい。
通信システムと結合している呼吸のための加圧空気供給システムを含む従来技術のシステムの簡単なブロック図。 図1のシステムのマスク−無線機音声パス。 マスク内部およびマスク・マイクロホン出力のところで測定したスペクトル振幅応答の一例、およびマスク、マイクロホンおよびマイクロホン増幅器に対する計算した結合伝達関数。 SCBAエア・レギュレータが発生する吸気ノイズの一例。 図4の吸気ノイズの長期振幅スペクトル。 所与のSCBAマスクを着用している一人の話し手が発生する吸気ノイズの4つの重なっているスペクトル。 音声と混合している吸気ノイズ・バーストを示すSCBAマイクロホンからの音声出力。 本発明の一実施形態による吸気ノイズを検出し、除去するための方法の簡単なブロック図。 図8の方法で使用するスペクトル・マッチャの一実施形態の簡単なブロック図。 図8の方法で使用するスペクトル・マッチャの他の実施形態の簡単なブロック図。 本発明の他の実施形態による音声信号を等化するための方法の簡単なブロック図。 本発明による14番目の次数のおよび20番目の次数のLPC逆フィルタ等化後のスペクトルと比較した等化前の吸気ノイズ・スペクトル。 バイオメトリック・パラメータを測定する際に使用するための本発明の他の実施形態による、吸気ノイズの周波数の持続時間を決定し、呼吸数および空気使用量を決定するための方法の簡単なブロック図。 音声および空気規制吸気ノイズを含むマイクロホン入力からの信号を示す。 図13の方法により測定した図14の信号の平均正規化モデル誤差。 図13の方法により発生した吸気ノイズ検出装置の出力信号。 図13の方法により発生した統合吸気検出装置の出力。

Claims (9)

  1. 音声信号が入るマスクを含み、音声に影響を与える伝達関数を有する加圧空気供給システム内で発生する前記音声信号を等化するための方法であって、
    吸気ノイズに基づいて吸気ノイズ・モデルを生成するステップと、
    前記音声信号を含む入力信号を受信するステップと、
    前記ノイズ・モデルに基づいて前記音声信号を等化するステップであって、前記ノイズ・モデルが、前記音声信号の等化により、前記伝達関数の影響の少なくとも一部が除去されるように、前記マスクの伝達関数を特徴付けるステップと、を含む方法。
  2. 前記吸気ノイズ・モデルを生成する前記ステップが、
    前記吸気ノイズの少なくとも1つのデジタル化したサンプルを生成するために前記吸気ノイズをサンプリングするステップと、
    前記少なくとも1つのデジタル化したサンプルに窓処理するステップと、
    前記少なくとも1つの窓処理したデジタル化したサンプルから一組の自己相関係数を決定するステップと、
    前記一組の自己相関係数に基づいて一組の線形予測符号化(LPC)係数を生成するステップと、
    前記一組のLPC係数からLPCフィルタを生成するステップと、を含む請求項1に記載の方法。
  3. 前記ノイズ・モデルがデジタル信号として表示され、
    デジタル・フィルタが、一組の自己相関係数から生成する一組のLPC係数に基づく線形予測符号化(LPC)フィルタであり、
    前記音声信号が前記LPCフィルタの倒置により等化される請求項1に記載の方法。
  4. 前記等化した音声信号をポストフィルタ処理するステップをさらに含む請求項1に記載の方法。
  5. 前記入力信号が吸気ノイズをさらに含み、前記方法が前記ノイズ・モデルを更新するステップをさらに含む請求項1に記載の方法。
  6. 前記ノイズ・モデルを更新する前記ステップが、
    類似性の測定値を入手するために前記入力信号を前記ノイズ・モデルと比較するステップと、
    前記入力信号の前記吸気ノイズを検出するために前記類似性の測定値を少なくとも1つの閾値と比較するステップであって、前記ノイズ・モデルが、前記検出した吸気ノイズに基づいて更新されるステップとを含む請求項5に記載の方法。
  7. 前記ノイズ・モデルが、第1の組の自己相関係数に基づいて生成する一組のLPC係数に基づく線形予測符号化(LPC)フィルタであり、前記ノイズ・モデルを更新する前記ステップが、
    前記検出した吸気ノイズの少なくとも1つのデジタル化したサンプルを生成するために前記検出した吸気ノイズをサンプリングするステップと、
    前記少なくとも1つのデジタル化したサンプルから第2の組の自己相関係数を決定するステップと、
    前記第1および第2の組の自己相関係数の関数として前記第1の組の自己相関係数を更新するステップと、
    前記更新した一組の自己相関係数に基づいて前記一組のLPC係数を更新するステップと、
    前記更新した一組のLPC係数に基づいて前記LPCフィルタを更新するステップと、をさらに含む請求項6に記載の方法。
  8. 音声信号が入るマスクを含み、音声に影響を与える伝達関数を有する加圧空気供給システムで前記音声信号を等化するためのデバイスであって、
    処理素子と、
    前記処理素子と結合していて、
    吸気ノイズに基づいて吸気ノイズ・モデルを生成するステップと、
    前記音声信号を含む入力信号を受信するステップと、
    前記ノイズ・モデルに基づいて前記音声信号を等化するステップであって、前記ノイズ・モデルが、前記音声信号の等化により前記伝達関数の影響の少なくとも一部が除去されるように、前記マスクの伝達関数を特徴付けるステップと
    を実行するように、処理装置に命令するためのコンピュータ・プログラムを格納するためのメモリ素子と
    を備えるデバイス。
  9. 音声信号を等化するためのシステムであって、
    前記音声信号が入るマスクを含み、音声に影響を与える伝達関数を有する加圧空気供給システムと、
    前記加圧空気供給システムと結合している通信システムとを備え、該システムが、
    処理素子と、
    前記処理素子と結合していて、
    吸気ノイズに基づいて吸気ノイズ・モデルを生成するステップと、
    前記音声信号を含む入力信号を受信するステップと、
    前記ノイズ・モデルに基づいて前記音声信号を等化するステップであって、前記ノイズ・モデルが、前記音声信号の等化により、前記伝達関数の影響の少なくとも一部が除去されるように、前記マスクの伝達関数を特徴付けるステップと、
    を実行するように、処理装置に命令するためのコンピュータ・プログラムを格納するためのメモリ素子と、
    を備えるシステム。
JP2007519235A 2004-06-30 2005-06-06 自給式呼吸装置システム内で発生した音声信号を等化するための方法および装置 Expired - Fee Related JP4468447B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US10/882,450 US7254535B2 (en) 2004-06-30 2004-06-30 Method and apparatus for equalizing a speech signal generated within a pressurized air delivery system
PCT/US2005/019827 WO2006007290A2 (en) 2004-06-30 2005-06-06 Method and apparatus for equalizing a speech signal generated within a self-contained breathing apparatus system

Publications (2)

Publication Number Publication Date
JP2008505538A JP2008505538A (ja) 2008-02-21
JP4468447B2 true JP4468447B2 (ja) 2010-05-26

Family

ID=35658374

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007519235A Expired - Fee Related JP4468447B2 (ja) 2004-06-30 2005-06-06 自給式呼吸装置システム内で発生した音声信号を等化するための方法および装置

Country Status (9)

Country Link
US (1) US7254535B2 (ja)
EP (1) EP1769493B1 (ja)
JP (1) JP4468447B2 (ja)
CN (1) CN101010731A (ja)
AU (1) AU2005262623B2 (ja)
CA (1) CA2572715C (ja)
DE (1) DE602005008974D1 (ja)
MX (1) MXPA06015235A (ja)
WO (1) WO2006007290A2 (ja)

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7748380B1 (en) 2004-04-06 2010-07-06 Sti Licensing Corporation Combined air-supplying/air-purifying system
US7658190B1 (en) 2004-04-06 2010-02-09 Sti Licensing Corp. Portable air-purifying system utilizing enclosed filters
US7505902B2 (en) * 2004-07-28 2009-03-17 University Of Maryland Discrimination of components of audio signals based on multiscale spectro-temporal modulations
US9185487B2 (en) 2006-01-30 2015-11-10 Audience, Inc. System and method for providing noise suppression utilizing null processing noise subtraction
JP4866958B2 (ja) * 2006-05-04 2012-02-01 株式会社ソニー・コンピュータエンタテインメント コンソール上にファーフィールドマイクロフォンを有する電子装置におけるノイズ除去
WO2008067834A1 (en) * 2006-12-07 2008-06-12 Akg Acoustics Gmbh Dropout concealment for a multi-channel arrangement
US8265937B2 (en) * 2008-01-29 2012-09-11 Digital Voice Systems, Inc. Breathing apparatus speech enhancement using reference sensor
US8245706B2 (en) * 2009-06-22 2012-08-21 Motorola Solutions, Inc. Pressure activated remote microphone
US9192773B2 (en) 2009-07-17 2015-11-24 Peter Forsell System for voice control of a medical implant
JP5294085B2 (ja) * 2009-11-06 2013-09-18 日本電気株式会社 情報処理装置、その付属装置、情報処理システム、その制御方法並びに制御プログラム
US8473287B2 (en) 2010-04-19 2013-06-25 Audience, Inc. Method for jointly optimizing noise reduction and voice quality in a mono or multi-microphone system
US8538035B2 (en) 2010-04-29 2013-09-17 Audience, Inc. Multi-microphone robust noise suppression
US8781137B1 (en) 2010-04-27 2014-07-15 Audience, Inc. Wind noise detection and suppression
US9558755B1 (en) 2010-05-20 2017-01-31 Knowles Electronics, Llc Noise suppression assisted automatic speech recognition
US8447596B2 (en) * 2010-07-12 2013-05-21 Audience, Inc. Monaural noise suppression based on computational auditory scene analysis
US8965756B2 (en) 2011-03-14 2015-02-24 Adobe Systems Incorporated Automatic equalization of coloration in speech recordings
US9183844B2 (en) * 2012-05-22 2015-11-10 Harris Corporation Near-field noise cancellation
CN104937580B (zh) * 2012-09-24 2018-03-27 约翰·汉密尔顿 通信和语音增强系统
US9943712B2 (en) 2012-09-24 2018-04-17 Dolores Speech Products Llc Communication and speech enhancement system
US9640194B1 (en) 2012-10-04 2017-05-02 Knowles Electronics, Llc Noise suppression for speech processing based on machine-learning mask estimation
US9498658B2 (en) * 2013-02-01 2016-11-22 3M Innovative Properties Company Respirator mask speech enhancement apparatus and method
US9517366B2 (en) 2013-02-01 2016-12-13 3M Innovative Properties Company Respirator mask speech enhancement apparatus and method
FR3005823B1 (fr) * 2013-05-14 2016-10-14 Elno Microphone comportant un interrupteur de sourdine, et masque de respiration comprenant un tel microphone
GB2519117A (en) * 2013-10-10 2015-04-15 Nokia Corp Speech processing
DE112015003945T5 (de) 2014-08-28 2017-05-11 Knowles Electronics, Llc Mehrquellen-Rauschunterdrückung
US10163453B2 (en) * 2014-10-24 2018-12-25 Staton Techiya, Llc Robust voice activity detector system for use with an earphone
US9843859B2 (en) 2015-05-28 2017-12-12 Motorola Solutions, Inc. Method for preprocessing speech for digital audio quality improvement
US11631421B2 (en) * 2015-10-18 2023-04-18 Solos Technology Limited Apparatuses and methods for enhanced speech recognition in variable environments
CN105405447B (zh) * 2015-10-27 2019-05-24 航宇救生装备有限公司 一种送话呼吸噪声屏蔽方法
CN106935247A (zh) * 2017-03-08 2017-07-07 珠海中安科技有限公司 一种用于正压式空气呼吸器和狭小密闭空间的语音识别控制装置及方法
IT201700090078A1 (it) * 2017-08-03 2019-02-03 Mestel Safety S R L Maschera per uso subacqueo, in particolare di tipo granfacciale dotata di dispositivo di comunicazione.
US11404061B1 (en) * 2021-01-11 2022-08-02 Ford Global Technologies, Llc Speech filtering for masks

Family Cites Families (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3347229A (en) * 1965-02-23 1967-10-17 Sierra Eng Co Latch operated microphone switch for breathing mask
US3415245A (en) * 1965-03-08 1968-12-10 Y2 Associates Noise-suppression diving apparatus
US3850168A (en) * 1971-09-21 1974-11-26 Puritan Bennett Corp Oxygen mask apparatus
US4154981A (en) * 1977-12-16 1979-05-15 The United States Of America As Represented By The Secretary Of The Navy Telephone system for diver communication
US4376916A (en) * 1980-05-29 1983-03-15 Cbs Inc. Signal compression and expansion system
US4720802A (en) * 1983-07-26 1988-01-19 Lear Siegler Noise compensation arrangement
DE3609097A1 (de) * 1986-03-19 1987-09-24 Draegerwerk Ag Sprech- und hoereinrichtung fuer atemgeraete
US5143078A (en) * 1987-08-04 1992-09-01 Colin Electronics Co., Ltd. Respiration rate monitor
US4862503A (en) * 1988-01-19 1989-08-29 Syracuse University Voice parameter extractor using oral airflow
US4958638A (en) * 1988-06-30 1990-09-25 Georgia Tech Research Corporation Non-contact vital signs monitor
US5222190A (en) * 1991-06-11 1993-06-22 Texas Instruments Incorporated Apparatus and method for identifying a speech pattern
US5307793A (en) * 1992-06-29 1994-05-03 Puritan-Bennett Corporation Microphone signal attenuating apparatus for oxygen masks
US5444786A (en) * 1993-02-09 1995-08-22 Snap Laboratories L.L.C. Snoring suppression system
KR0121126B1 (ko) * 1994-06-09 1997-12-04 구자홍 전자악기의 자동반주시 코드변화처리방법
JPH10501951A (ja) * 1995-04-03 1998-02-17 フィリップス エレクトロニクス ネムローゼ フェンノートシャップ 自動イコライザを有する信号増幅システム
FI951918A (fi) * 1995-04-21 1996-10-22 Nokia Mobile Phones Ltd Lähetin-vastaanotinlaite ja menetelmä kompleksien I/Q-signaalin synnyttämiseksi ja käsittelemiseksi
US5730140A (en) * 1995-04-28 1998-03-24 Fitch; William Tecumseh S. Sonification system using synthesized realistic body sounds modified by other medically-important variables for physiological monitoring
US5605145A (en) 1995-07-18 1997-02-25 Puritan-Bennett Corporation Microphone attenuation device for use in oxygen breathing masks
US5579284A (en) * 1995-07-21 1996-11-26 May; David F. Scuba diving voice and communication system using bone conducted sound
US5822360A (en) * 1995-09-06 1998-10-13 Solana Technology Development Corporation Method and apparatus for transporting auxiliary data in audio signals
EP0776114A3 (fr) * 1995-11-22 1997-06-04 Laboratoires D'electronique Philips S.A.S. Appareil téléphonique avec volume réglable en fonction de bruit ambiant
US5727074A (en) * 1996-03-25 1998-03-10 Harold A. Hildebrand Method and apparatus for digital filtering of audio signals
US5734090A (en) * 1996-03-27 1998-03-31 Alcohol Sensors International, Ltd. Method and apparatus for sonic breath determination
GB9609321D0 (en) * 1996-05-03 1996-07-10 British Telecomm Automatic speech recognition
US6470315B1 (en) * 1996-09-11 2002-10-22 Texas Instruments Incorporated Enrollment and modeling method and apparatus for robust speaker dependent speech models
AUPO247496A0 (en) * 1996-09-23 1996-10-17 Resmed Limited Assisted ventilation to match patient respiratory need
US5890111A (en) * 1996-12-24 1999-03-30 Technology Research Association Of Medical Welfare Apparatus Enhancement of esophageal speech by injection noise rejection
US6295364B1 (en) * 1998-03-30 2001-09-25 Digisonix, Llc Simplified communication system
US8023665B2 (en) * 1998-05-04 2011-09-20 Schwartz Stephen R Microphone-tailored equalizing system
US6199550B1 (en) * 1998-08-14 2001-03-13 Bioasyst, L.L.C. Integrated physiologic sensor system
FR2786308B1 (fr) * 1998-11-20 2001-02-09 Sextant Avionique Procede de reconnaissance vocale dans un signal acoustique bruite et systeme mettant en oeuvre ce procede
US6324499B1 (en) * 1999-03-08 2001-11-27 International Business Machines Corp. Noise recognizer for speech recognition systems
US6304844B1 (en) * 2000-03-30 2001-10-16 Verbaltek, Inc. Spelling speech recognition apparatus and method for communications
US20020116186A1 (en) * 2000-09-09 2002-08-22 Adam Strauss Voice activity detector for integrated telecommunications processing
US6892175B1 (en) * 2000-11-02 2005-05-10 International Business Machines Corporation Spread spectrum signaling for speech watermarking
KR200226168Y1 (ko) * 2000-12-28 2001-06-01 엘지전자주식회사 이퀄라이저 기능을 구비한 휴대 통신 장치
JP3609748B2 (ja) * 2001-05-14 2005-01-12 株式会社バーテックススタンダード マイクロホンの特性調整装置
GB2380644A (en) * 2001-06-07 2003-04-09 Canon Kk Speech detection
US6702752B2 (en) * 2002-02-22 2004-03-09 Datex-Ohmeda, Inc. Monitoring respiration based on plethysmographic heart rate signal

Also Published As

Publication number Publication date
MXPA06015235A (es) 2007-03-26
EP1769493A4 (en) 2007-08-22
CA2572715A1 (en) 2006-01-19
WO2006007290A2 (en) 2006-01-19
AU2005262623A1 (en) 2006-01-19
AU2005262623B2 (en) 2008-07-03
EP1769493B1 (en) 2008-08-13
US7254535B2 (en) 2007-08-07
WO2006007290B1 (en) 2006-08-24
DE602005008974D1 (de) 2008-09-25
CN101010731A (zh) 2007-08-01
EP1769493A2 (en) 2007-04-04
US20060020451A1 (en) 2006-01-26
JP2008505538A (ja) 2008-02-21
WO2006007290A3 (en) 2006-06-01
CA2572715C (en) 2012-04-03

Similar Documents

Publication Publication Date Title
JP4468447B2 (ja) 自給式呼吸装置システム内で発生した音声信号を等化するための方法および装置
AU2005262624B2 (en) Method and apparatus for detecting and attenuating inhalation noise in a communication system
AU2005262586B2 (en) Method and apparatus for characterizing inhalation noise and calculating parameters based on the characterization
US8606572B2 (en) Noise cancellation device for communications in high noise environments
JP6220989B2 (ja) 吸気流を検出するための2つのマイクを備える吸入器
US9418675B2 (en) Wearable communication system with noise cancellation
CN110868921A (zh) 用于确定心律和/或呼吸速率的系统和方法
WO2012117393A1 (en) A survival necklace
US11295759B1 (en) Method and apparatus for measuring distortion and muffling of speech by a face mask
WO2016154139A1 (en) Sound-based spirometric devices, systems, and methods using audio data transmitted over a voice communication channel
RU2625274C1 (ru) Способ контроля физиологических параметров дыхательной системы водолазов
Radha et al. A Study on Alternative Speech Sensor
JP5136983B2 (ja) 気道系検査装置
Kushner et al. The acoustic properties of SCBA equipment and its effects on speech communication
Gerritsen et al. A source impedance compensated artificial voice for mask speech intelligibility testing: Thevenin model

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: 20100126

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100224

R150 Certificate of patent or registration of utility model

Ref document number: 4468447

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: 20130305

Year of fee payment: 3

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130305

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130305

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140305

Year of fee payment: 4

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

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees