腹部大動脈瘤検出装置
【課題】腹部大動脈瘤が生じているか否かを容易に判定可能にする検出装置を提供する。
【解決手段】人体の臍周辺における正中を挟んだ少なくとも左右2箇所のそれぞれにエアクッション10を当接し、各エアクッション10を当接した部位における脈波によって生じる各空気圧変動を検出するセンサ14からの出力信号の信号波形間のずれ度合いを求め、各信号波形間のずれ度合いから腹部大動脈瘤の有無を判別する構成である。腹部大動脈瘤が生じている場合、その部分の血管はうねっており、正常状態と比較して左右いずれかに突出している。このため、少なくとも左右2箇所のエアクッション10の空気圧変動の信号波形は実質的に同じとはならず、一方の振幅が大きくなったり、乱れが生じたりする。この信号波形間のずれ度合いから腹部大動脈瘤の有無を簡易に判別することができる。
【解決手段】人体の臍周辺における正中を挟んだ少なくとも左右2箇所のそれぞれにエアクッション10を当接し、各エアクッション10を当接した部位における脈波によって生じる各空気圧変動を検出するセンサ14からの出力信号の信号波形間のずれ度合いを求め、各信号波形間のずれ度合いから腹部大動脈瘤の有無を判別する構成である。腹部大動脈瘤が生じている場合、その部分の血管はうねっており、正常状態と比較して左右いずれかに突出している。このため、少なくとも左右2箇所のエアクッション10の空気圧変動の信号波形は実質的に同じとはならず、一方の振幅が大きくなったり、乱れが生じたりする。この信号波形間のずれ度合いから腹部大動脈瘤の有無を簡易に判別することができる。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、脈波によって生じるエアクッションの空気圧変動の信号波形から腹部大動脈瘤の有無を簡易的に検出する腹部大動脈瘤検出装置に関する。
【背景技術】
【0002】
腹部大動脈瘤は、超音波診断装置やCTを用いれば容易に発見することができる。しかし、超音波診断装置やCT自体、相応の導入コストがかかるため、医療機関における保有台数にも限りがある。従って、自覚症状を訴えた患者、医師が触診により疑わしいと判定した患者に対して超音波診断装置やCTによる検査が行われていることが一般的である。しかし、より安価で簡易な装置が開発できれば、例えば、健康診断などにおいて、健診対象者の全てに対して腹部大動脈瘤に関する検査を行って、自覚症状の有無に拘わらず腹部大動脈瘤の早期発見を行うことが可能である。このため、超音波診断装置やCTなどよりも診断精度が劣るとしても、安価に入手でき、かつ簡易に操作可能な装置であって、腹部大動脈瘤の保有者をある程度の精度でスクリーニングできる装置の開発が望まれていた。
【0003】
一方、本出願人は、例えば、特許文献1として、内部に三次元立体編物を挿入したエアクッションを備え、このエアクッションを例えば人の腰部に対応する部位に配置し、エアクッションの空気圧変動を測定し、得られた空気圧変動の時系列データから人の生体信号を検出し、人の生体の状態を分析するシステムを開示している。また、非特許文献1及び2においても、腰腸肋筋に沿うようにエアクッションを配置して人の生体信号を検出する試みを報告している。
【先行技術文献】
【特許文献】
【0004】
【特許文献1】特開2007−90032号公報
【非特許文献】
【0005】
【非特許文献1】「非侵襲型センサによって測定された生体ゆらぎ信号の疲労と入眠予知への応用」、落合直輝(外6名)、第39回日本人間工学会 中国・四国支部大会 講演論文集、平成18年11月25日発行、発行所:日本人間工学会 中国・四国支部事務局
【非特許文献2】「非侵襲生体信号センシング機能を有する車両用シートの試作」、前田慎一郎(外4名)、第39回日本人間工学会 中国・四国支部大会 講演論文集、平成18年11月25日発行、発行所:日本人間工学会 中国・四国支部事務局
【発明の概要】
【発明が解決しようとする課題】
【0006】
特許文献1及び非特許文献1、2によれば、腰部付近の動脈(大動脈)の脈波を検知し、得られた脈波の時系列信号データを用い、例えば、本出願人が特開2004−344612において提案した手法により入眠予兆を判定することができる。そして、内部に三次元立体編物を挿入したエアクッションを備えた生体信号測定手段は、低コストで製作できるという利点がある。
【0007】
そこで、本発明は、上記エアクッションを備えた生体信号測定手段を用いて、腹部大動脈瘤の簡易的な検出を行う技術を提供することを課題とする。
【課題を解決するための手段】
【0008】
上記課題を解決するため、本発明の腹部大動脈瘤検出装置は、少なくとも、人体の臍周辺における正中を挟んだ左右2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかに当接可能な2つのエアクッションと、前記各エアクッションを当接した部位における脈波によって生じる各空気圧変動を検出する少なくとも2つのセンサとを備えた生体信号測定手段と、前記生体信号測定手段の各センサの出力信号を受信する受信部と、各出力信号の信号波形間のずれ度合いを求める演算部と、前記演算部により求められた信号波形間のずれ度合いから腹部大動脈瘤の有無を判別する判別部とが設定された脈波検出手段とを備えてなることを特徴とする。
【0009】
前記生体信号測定手段は、人体の臍周辺における正中を挟んだ左右2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかに当接可能な前記2つのエアクッションに加えて、さらに、左右2箇所の腸骨動脈に対応する部位及び上腹部正中のうちの1箇所以上にも当接可能な1以上のエアクッションを備えると共に、これらの各エアクッションのそれぞれにも前記センサが設けられていることが好ましい。また、前記生体信号測定手段は、前記各エアクッションに積層されるビーズ発泡樹脂弾性部材を備える構成とすることができる。前記ビーズ発泡樹脂弾性部材は、ビーズ発泡体と該ビーズ発泡体の外面を被覆する被覆材とからなることが好ましい。
【0010】
前記生体信号測定手段は、前記各エアクッションが三次元立体編物の一面に保持されており、該三次元立体編物を腹部に巻き付けて脈波を測定する構成とすることができる。また、前記生体信号測定手段は、椅子のシートバック部に張設された三次元立体編物の一面に前記各エアクッションが保持されており、該椅子に着席して脈波を測定する構成とすることができる。
【発明の効果】
【0011】
本発明によれば、人体の臍周辺における正中を挟んだ左右2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかに前記2つのエアクッションを当接し、各エアクッションを当接した部位における脈波によって生じる各空気圧変動を検出するセンサからの出力信号の信号波形間のずれ度合いを求め、各信号波形間のずれ度合いから腹部大動脈瘤の有無を時系列解析や周波数解析などによって判別する構成である。腹部大動脈瘤が生じている場合、その部分の血管には部分的な変形が生じており、正常状態と比較して上下、左右、前後のいずれかの方向に膨出している。このため、腹部大動脈瘤が生じている付近における少なくとも左右2箇所、上下2箇所又は前後(腹側と背側)2箇所のエアクッションの空気圧変動の信号波形間に変位が生じ、一方の振幅が大きくなったり、波形に極端な乱れが生じたりする。従って、信号波形間のずれ度合いから腹部大動脈瘤の有無を簡易に判別することができる。生体信号測定手段はエアクッションから構成されるため、低コストで製造可能であることから、装置全体の導入コストも低く抑えられる。このため、健康診断等において、超音波診断装置やCTでの精密検査を要するか否かを判定する簡易的なスクリーニング手段として有益である。
【図面の簡単な説明】
【0012】
【図1】図1は、本発明の一の実施形態に係る腹部大動脈瘤検出装置の概略構成を示した図である。
【図2】図2(a)は、エアクッションにビーズ発泡樹脂弾性部材を積層した状態の一例を示し、図2(b)は、ビーズ発泡樹脂弾性部材の構成を説明するための図である。
【図3】図3は、脈波検出手段の構成を説明するための図である。
【図4】図4は、エアクッションのセンサから検出した信号波形の一例を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形である。
【図5】図5は、指尖容積脈波計から検出した信号波形の一例を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形である。
【図6】図6は、エアクッションのセンサから検出した被験者Aの信号波形の周波数解析結果を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図7】図7は、被験者Aの上記信号波形のカオス解析により得られるアトラクタを示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図8】図8(a)は、被験者Aの上記信号波形のカオス解析により得られる原波形、1階微分波形、2階微分波形のリアプノフ指数の時系列波形であり、(b)は原波形、1階微分波形、2階微分波形のリアプノフ指数の周波数スペクトルを示した図である。
【図9】図9は、エアクッションのセンサから検出した被験者Bの信号波形の周波数解析結果を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図10】図10は、被験者Bの上記信号波形のカオス解析により得られるアトラクタを示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図11】図11(a)は、被験者Bの上記信号波形のカオス解析により得られる原波形、1階微分波形、2階微分波形のリアプノフ指数の時系列波形であり、(b)は原波形、1階微分波形、2階微分波形のリアプノフ指数の周波数スペクトルを示した図である。
【図12】図12は、エアクッションのセンサから検出した被験者Cの信号波形の周波数解析結果を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図13】図13は、被験者Cの上記信号波形のカオス解析により得られるアトラクタを示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図14】図14(a)は、被験者Cの上記信号波形のカオス解析により得られる原波形、1階微分波形、2階微分波形のリアプノフ指数の時系列波形であり、(b)は原波形、1階微分波形、2階微分波形のリアプノフ指数の周波数スペクトルを示した図である。
【図15】図15は、エアクッションのセンサから検出した被験者Dの信号波形の周波数解析結果を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図16】図16は、被験者Dの上記信号波形のカオス解析により得られるアトラクタを示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図17】図17(a)は、被験者Dの上記信号波形のカオス解析により得られる原波形、1階微分波形、2階微分波形のリアプノフ指数の時系列波形であり、(b)は原波形、1階微分波形、2階微分波形のリアプノフ指数の周波数スペクトルを示した図である。
【図18】図18は、エアクッションのセンサから検出した被験者Eの信号波形の周波数解析結果を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図19】図19は、被験者Eの上記信号波形のカオス解析により得られるアトラクタを示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図20】図20(a)は、被験者Eの上記信号波形のカオス解析により得られる原波形、1階微分波形、2階微分波形のリアプノフ指数の時系列波形であり、(b)は原波形、1階微分波形、2階微分波形のリアプノフ指数の周波数スペクトルを示した図である。
【図21】図21は、被験者F,Gの臍を挟んだ上下2カ所にエアクッションを当てて胸部脈波と腹部脈波を検出する試験の試験方法を説明するための図である。
【図22】図22(a)は、エアクッションのセンサから検出した被験者Fの胸部脈波の波形(胸部脈波形)を示した図であり、図22(b)は、エアクッションのセンサから検出した被験者Fの腹部脈波の波形(腹部脈波形)を示した図であり、図22(c)は、エアクッションのセンサから検出した被験者Gの胸部脈波の波形(胸部脈波形)を示した図であり、図22(d)は、エアクッションのセンサから検出した被験者Gの腹部脈波の波形(腹部脈波形)を示した図である。
【図23】図23(a)は、図22(a),(b)に示した被験者Fの胸部脈波形と腹部脈波形の一部拡大図であり、図23(b)は、図22(c),(d)に示した被験者Gの胸部脈波形と腹部脈波形の一部拡大図である。
【図24】図24は、図22(a),(b)に示した被験者Fの胸部脈波形と腹部脈波形の位相差と、図22(c),(d)に示した被験者Gの胸部脈波形と腹部脈波形の位相差を示した図である。
【図25】図25(a)は、図22(a),(b)に示した被験者Fの胸部脈波形と腹部脈波形の周波数解析を示した図であり、図22(c),(d)に示した被験者Gの胸部脈波形と腹部脈波形の周波数解析結果を示した図である。
【図26】図26は、エアクッションを保持部材に装填した状態の一例を示した図である。
【発明を実施するための形態】
【0013】
以下、図面に示した本発明の実施形態に基づき、本発明をさらに詳細に説明する。図1は、本実施形態に係る腹部大動脈瘤検出装置1を示した構成図であり、生体信号測定手段2と脈波検出手段3とを備えて構成される。
【0014】
生体信号測定手段2は、2つのエアクッション10を備えて構成されている。各エアクッション10は、各エアクッション10は、図1及び図2に示したように、表側エアクッション11と裏側エアクッション12とが積層されて構成される。表側エアクッション11及び裏側エアクッション12とは、それぞれ空気袋11a,12aとその内部に配置される復元力付与部材としての三次元立体編物11b,12bとを備えている。表側エアクッション11と裏側エアクッション12とは、長手方向に沿った一方の側縁同士が接合され、接合された側縁を中心にして2つ折りにされて、相互に重ね合わせられて用いられる。
【0015】
表側エアクッション11にはセンサ取付チューブ13が設けられ、その内側に空気圧変動を測定するセンサ14が固定されている。なお、センサ取付チューブ13は密閉されている。また、空気袋11a,12aの素材は限定されるものではないが、例えば、ポリウレタンエラストマー(例えば、シーダム株式会社製、品番「DUS605−CDR」)からなるシートを用いて形成することができる。センサ14としては、空気袋11a,12a内の空気圧を測定できるものであればよく、例えば、コンデンサ型マイクロフォンセンサを用いることができる。また、空気袋11a,12a内に装填される三次元立体編物11b,12bは、例えば、特開2002−331603号公報に開示されているように、互いに離間して配置された一対のグランド編地と、該一対のグランド編地間を往復して両者を結合する多数の連結糸とを有する立体的な三次元構造となった編地である。
【0016】
一方のグランド編地は、例えば、単繊維を撚った糸から、ウェール方向及びコース方向のいずれの方向にも連続したフラットな編地組織(細目)によって形成され、他方のグランド編地は、例えば、短繊維を撚った糸から、ハニカム状(六角形)のメッシュを有する編み目構造に形成されている。もちろん、この編地組織は任意であり、細目組織やハニカム状以外の編地組織を採用することもできるし、両者とも細目組織を採用するなど、その組み合わせも任意である。連結糸は、一方のグランド編地と他方のグランド編地とが所定の間隔を保持するように、2つのグランド編地間に編み込んだものである。このような三次元立体編物としては、例えば、以下のようなものを用いることができる。なお、各三次元立体編物は、必要に応じて複数枚積層して用いることもできる。
【0017】
(1)製品番号:49076D(住江織物(株)製)
材質:
表側のグランド編地・・・300デシテックス/288fのポリエチレンテレフタレート繊維仮撚加工糸と700デシテックス/192fのポリエチレンテレフタレート繊維仮撚加工糸との撚り糸
裏側のグランド編地・・・450デシテックス/108fのポリエチレンテレフタレート繊維仮撚加工糸と350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメントとの組み合わせ
連結糸・・・・・・・・・350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメント
【0018】
(2)製品番号:49013D(住江織物(株)製)
材質:
表側のグランド編地・・・450デシテックス/108fのポリエチレンテレフタレート繊維仮撚加工糸の2本の撚り糸
裏側のグランド編地・・・450デシテックス/108fのポリエチレンテレフタレート繊維仮撚加工糸の2本の撚り糸
連結糸・・・・・・・・・350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメント
【0019】
(3)製品番号:69030D(住江織物(株)製)
材質:
表側のグランド編地・・・450デシテックス/144fのポリエチレンテレフタレート繊維仮撚加工糸の2本の撚り糸
裏側のグランド編地・・・450デシテックス/144fのポリエチレンテレフタレート繊維仮撚加工糸と350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメントとの組み合わせ
連結糸・・・・・・・・・350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメント
【0020】
(4)旭化成せんい(株)製の製品番号:T24053AY5−1S
【0021】
生体信号測定手段2は、図2(a)に示したように、上記したエアクッション10に積層されるビーズ発泡樹脂弾性部材15,16を備えていることが好ましい。ビーズ発泡樹脂弾性部材15,16を人体とエアクッション10との間に設けることにより、測定時において、人がエアクッション10の異物感を感じにくくなる。その一方、ビーズ発泡樹脂弾性部材15,16として、好ましくは、直径30mmの加圧板によってたわみ量1mmまで加圧した際の荷重−たわみ特性から得られるバネ定数が、エアクッション10よりも高いものを用いる。これにより、生体信号を増幅してエアクッション10に伝達できる。
【0022】
ビーズ発泡樹脂弾性部材15,16は、図2(b)に示したように、略平板状に形成されたビーズ発泡体15a,16aと、その外面に貼着される被覆材15b,16bとから構成されている。ビーズ発泡体15a,16aとしては、ポリスチレン、ポリプロピレン及びポリエチレンのいずれか少なくとも一つを含む樹脂のビーズ法による発泡成形体が用いられる。なお、発泡倍率は任意であり限定されるものではない。被覆材15b,16bは、ビーズ発泡体15a,16aの外面に接着により貼着され、高い伸度と回復率を有する素材とすることが好ましい。例えば、伸度200%以上、100%伸長時の回復率が80%以上である弾性繊維不織布を用いることができる。このような弾性繊維不織布としては、例えば、特開2007−92217号公報に開示された熱可塑性エラストマー弾性繊維が相互に溶融接着された不織布を用いることができる。具体的には、KBセーレン(株)製、商品名「エスパンシオーネ」を用いることができる。
【0023】
また、ビーズ発泡樹脂弾性部材15,16は、図2(a)に示したように、2枚を重ねて用いるだけでなく、エアクッション10を挟むようにして用いることもできる。この場合、2枚のビーズ発泡樹脂弾性部材は、バネ定数の異なるものを用いることが好ましい。例えば、上記の商品名「エスパンシオーネ」で被覆したビーズ発泡樹脂弾性部材を用いた場合、これと共に使用する他のビーズ発泡樹脂弾性部材は、被覆材として伸縮性の小さい素材、例えば、熱可塑性ポリエステルからなる不織布を用いた構成とすることができる。このような被覆材としては、例えば、帝人(株)製のポリエチレンナフタレート(PEN)繊維(1100dtex)から形成した2軸織物(縦:20本/inch、横:20本/inch)を用いることができる。
【0024】
次に、脈波検出手段3の構成について図3に基づき説明する。脈波検出手段3は、受信部31、演算部32、判別部33とが設定されたコンピュータから構成される。
【0025】
受信部31は、生体信号測定手段2を構成する2つのエアクッション10の空気圧変動に応じた出力信号を上記した2つのセンサ14からそれぞれ受信する機能を備えている。
【0026】
演算部32は、受信部31により受信した各出力信号の信号波形を演算し、それらを重ね合わせて2つの信号波形間のずれ度合いを求める。信号波形は、出力信号の原波形、原波形を微分した波形(1階微分波形、2階微分波形等)のいずれかを採用することができる。
【0027】
ここで、図4(a)は、一つのエアクッション10が臍部に当たるように被験者がうつぶせになって得られたセンサの出力信号の原波形であり、図4(b)は、図4(a)の1階微分波形、図4(c)は、図4(a)の2階微分波形である。また、図5(a)は、図4と同時に測定した被験者の指尖容積脈波(指尖容積脈波計((株)アムコ製、フィンガークリッププローブ SR−5C)により測定)の原波形であり、図5(b)は、図5(a)の1階微分波形、図5(c)は、図5(a)の2階微分波形である。これらの図から明らかなように、原波形、1階微分波形、2階微分波形のいずれも、エアクッションのセンサの出力信号は、指尖容積脈波の出力信号とよく一致しており、エアクッションによって腹部大動脈の脈波が捉えられていることがわかる。
【0028】
腹部大動脈溜が生じていると、上記のように血管がうねるなどの部分的な変形が生じていたり、あるいは、動脈硬化を併発していることが多いため血管雑音が伴いやすい。従って、2つのセンサから得られる信号波形は、一方の振幅が他方よりも大きく、きれいな波形にならなかったり、一方の信号波形が他方の信号波形では見られない複数の凹凸部を有していたり、あるいは所定の位相差があったりするなど、信号波形間にずれがある。そこで、演算部32はこのずれ度合いを求める。ずれ度合いの求め方としては、例えば、一方の振幅を基準として、他方の振幅の比率を求めることで算出することができる。
【0029】
判別部33は、演算部32により求められた信号波形間のずれ度合いに基づき、腹部大動脈瘤の有無を判別する。例えば、上記の振幅の違いからずれ度合いを求める場合には、一方の振幅を基準とした他方の振幅の比率が所定の閾値以上の場合に、「腹部大動脈瘤:有」と判定する。もちろん、本発明は、腹部大動脈瘤を生じている可能性のある人をスクリーニングすることが目的であり、ここでいう「腹部大動脈瘤:有」とは、腹部大動脈瘤が生じていることが信号波形の比較から疑われる人の意味であり、「腹部大動脈瘤:有」と判定された場合でも、実際には腹部大動脈瘤を有していない人が含まれる場合もある。閾値を低めに設定すれば、腹部大動脈瘤を生じていない者が含まれる比率が高くなるが、腹部大動脈瘤を有する人を大雑把に絞り込むという目的に照らすと、閾値はそれほど高く設定しない方が好ましい。
【0030】
なお、原波形を微分処理すると、原波形が強調されるが、1階微分と2階微分とでは強調の程度が異なるため、原波形の形やそれに含有される高周波成分のレベルにより、いずれの結果を採用すべきかは一概には言えず個人差がある。従って、原波形、1階微分波形、2階微分波形のいずれを採用するかは、得られた各信号波形から、比較しやすいものを操作者が選択するように設定することができる。また、原波形、1階微分波形、2階微分波形の全ての信号波形をそれぞれ比較して、そのうちの1つだけでも、ずれ度合いが閾値を超えている場合に腹部大動脈瘤「有」と判定するようにすれば、見落とし率が小さくなる。
【0031】
本実施形態に係る腹部大動脈瘤検出装置1を用いて腹部大動脈瘤の検査を行う場合には、検査対象者を寝具にうつぶせにして腹側における臍周辺の正中を挟んだ2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかにエアクッション10を当接するだけでよい。これにより、エアクッション10のセンサの出力信号が、脈波検出手段3の受信部31により受信され、演算部32、判別部33によって処理されて、2つの信号波形が比較される。その結果、信号波形間のずれ度合いが所定の閾値以上か否かにより、腹部大動脈溜の疑いがあるかどうかが上記のように出力される。
【0032】
なお、信号波形間のずれ度合いは、振幅のほか、周波数、位相、低周波及び高周波成分の割合、切痕の位置などから判断することもできる。また、振幅のずれ度合いは、時系列の信号波形を比較して求めることができ、低周波及び高周波成分の割合のずれ度合いは、信号波形を周波数解析することによって求めることができる。
【0033】
図6〜図20は、被験者A〜Eが、一つのエアクッション10を臍部に当ててうつぶせになってセンサから得られた脈波信号のうち、ノイズ及び体動のない200秒間を抽出して分析した結果である。図6、図9、図12、図15、図18は各被験者の脈波信号の周波数解析結果であり、図7、図10、図13、図16、図19は各被験者の脈波信号をカオス解析して描いたアトラクタを示した図であり、これらの図において、(a)は脈波信号の原波形の解析結果を、(b)は1階微分波形の解析結果を、(c)は2階微分波形の解析結果をそれぞれ示す。また、図8、図11、図14、図17、図20の(a)はカオス解析により得られるリアプノフ指数の時系列波形であり、(b)はリアプノフ指数の周波数スペクトルを示した図である。この周波数解析結果、アトラクタ図形、リアプノフ指数の波形は、出力結果に乱れがなく安定しているほど、被験者が健康であることを示し、出力結果の乱れが大きいほど、体調不良であることを示している。
【0034】
なお、各被験者の性別、年齢、体格等は次の通りであった。
・被験者A:男性、33歳、身長168cm、体重65.0kg、体脂肪率13.1%
・被験者B:男性、31歳、身長171cm、体重65.2kg、体脂肪率16.4%
・被験者C:女性、25歳、身長159cm、体重44.8kg、体脂肪率23.4%
・被験者D:男性、36歳、身長167cm、体重70.9kg、体脂肪率24.2%
・被験者E:女性、31歳、身長159cm、体重55.0kg、体脂肪率25.5%
【0035】
被験者Aは、図6から、原波形、1階微分波形、2階微分波形のいずれの周波数解析も、200秒間安定した周波数になっており、健康であることがわかる。図7のアトラクタを考察した場合も、全体的には安定しているが、2階微分波形のアトラクタにおいては、「00:02:13〜00:02:43」の時間帯において乱れがあり、基本的には健康であるが、若干の不調を読み取ることができる。図8のリアプノフ指数の解析においても図7と同様の傾向を読み取ることができる。被験者Aのデータの場合、この若干の不調を示すデータは、2階微分波形において読み取れることから、被験者Aの場合には、特に、2階微分波形による解析が適しているということもわかる。
【0036】
被験者Bは、図9の周波解析がいずれも安定していることから健康であることがわかる。図10及び図11においても、アトラクタ及び波形に乱れがほとんどなく、健康であると読み取れるが、「00:02:47〜00:03:17」の時間帯において若干の乱れがある。これは、試験後半において、疲労が生じたことによるものと思われる。
【0037】
被験者Cは、図12の周波数解析が安定していると共に、図13及び図14のアトラクタ及び波形が共に乱れがなく、健康であるとこがわかる。
【0038】
被験者Dは、図15の周波数解析において、周波数が安定していない部分があり、図16のアトラクタ図形においても乱れの生じた図形が多い。また、図17のリアプノフ指数の波形では、特に原波形の乱れが大きいことが読み取れる。これらのことから、被験者Dは体力があまり高くなく、しかも体調が万全でないと判断できる。
【0039】
被験者Eは、図18の周波数解析において、周波数が安定していない部分があり、図19及び図20のアトラクタ図形及び波形に乱れが多いことが読み取れる。従って、被験者Eも体調が万全でないと判断できる。実際、被験者Eの場合、本試験後に循環器系疾患が発見されており、本試験による解析が体調の健康度合いを示す指標として有用であることがわかった。
【0040】
図6〜図20の各解析手法は、検出された信号波形から、被験者の健康状態を判定できることを示したものであるが、この解析手法を複数のエアクッション10から得られる信号波形の解析に応用した場合には、例えば、腹部大動脈瘤の生じていない範囲に当接されたエアクッション10の検出信号の解析結果と、腹部大動脈瘤が生じている範囲に当接されたエアクッション10の検出信号の解析結果とは明確に異なることになると予測できる。従って、それらの解析結果のずれ度合い(周波数解析の差、アトラクタ図形の形の差、リアプノフ指数波形の差等)を求め、所定の閾値を設定することにより、腹部大動脈瘤の有無を判定できる。
【0041】
このように、本実施形態によれば、極めて簡易な構成のエアクッション10を用いて測定することができるため、導入コストが低く抑えられると共に、検査方法も簡易であることから、健康診断などにおいて、腹部大動脈瘤の有無(疑いがあるか否か)を容易に判定することができる。
【0042】
図21〜図25は、腹部大動脈瘤を患ってはいないが動脈疾患を指摘されている50歳代の肥満体型の被験者F(性別:男性、身長:174cm、体重:94kg)と、健康な20歳代の標準体型の被験者G(性別:女性、身長:167cm、体重:57kg)についての生体信号の解析結果を示した図である。試験は、被験者F,Gを寝具(マットレス)にうつぶせに寝かせ、図21に示したように、臍を挟んだ上下2カ所に相当する、臍よりやや下側の腹部とその位置から約30cm上側に離間した胸部に、一方のエアクッション10と他方のエアクッション10を当てて行った。
【0043】
図22(a)〜(d)は各エアクッション10,10から得られた被験者F,Gの検出信号の原波形であるが、その一部を拡大した波形が図23(a),(b)である。臍を挟んだ下側に位置する一方のエアクッション10から得られた検出信号の波形を「腹部脈波形」として、臍を挟んだ上側に位置する他方のエアクッション10から得られた検出信号の波形を「胸部脈波形」として示している。図23の(a)と(b)を比較すると、腹部脈波形と胸部脈波形の位相差は、被験者Gのデータ(図23(b))に対し、被験者Fのデータ(図23(a))の方が大きいことがわかる。すなわち、2つのエアクッション10,10の離間距離が30cm程度に過ぎないため、動脈疾患がなければ、被験者Gのように2つの脈波形間に位相差はほとんど生じないはずである。しかし、動脈疾患を有する場合には、被験者Fのように2つの脈波形間に明らかな位相差が生じることがわかる。図22の900秒間の原波形について腹部脈波形と胸部脈波形の位相差を計算により求めたところ、図24に示したように、被験者Fのデータの位相差が約60度であったのに対し、被験者Gのデータの位相差が約5度以下であり、動脈疾患を有する被験者Fの位相差が顕著に大きかった。
【0044】
そこで、脈波検出手段3を構成するコンピュータの記憶部に、腹部大動脈瘤を有する場合の位相差の大きさを閾値として予め設定しておく。閾値は、腹部大動脈瘤を有する場合の被験者の位相差をできるたけ多く測定して定めておくことが好ましい。そして、実際の判定に当たっては、演算部32において、2つの脈波形間の位相差を求めた後、判別部33において、求めた位相差を記憶部に設定されている上記閾値と比較して行う。その結果、予め記憶されている所定の閾値以上であったならば、腹部大動脈瘤を有すると判定する。これにより腹部大動脈瘤の有無を容易に判定できる。なお、閾値を調整すれば、腹部大動脈瘤に至る前の動脈疾患を判定することができる。
【0045】
また、図23(a)では、492〜495秒の間のように腹部脈波形の振幅が胸部脈波形の振幅よりも2倍前後大きいといったように、2つの脈波形間の振幅差が大きい。これに対し、図23(b)では、2つの脈波間の振幅差は図23(a)よりも小さい。そこで、この振幅差について所定の閾値を設定することにより、腹部大動脈瘤の有無を含む動脈疾患の判定に利用できる。
【0046】
また、被験者Gの場合には、胸部脈波形において切痕が捉えられている。切痕は、拍出期の終わりに大動脈が急に閉鎖することに対応して現れる信号であり、健康な人であればこれが捉えられる。これに対し、被験者Fの場合には、胸部脈波形及び腹部脈波形のいずれも切痕が明確には捉えられていない。従って、この切痕の有無も、腹部大動脈瘤を含む動脈疾患の有無の判定に利用できる。
【0047】
一方、図25の周波数解析結果を見ると、(b)の被験者Gのデータは、0.4Hz以上の周波数帯域の全てにおいて、胸部脈波よりも腹部脈波のパワースペクトルが大きくなっている。これに対し、(a)の被験者Fのデータは、0.4Hz以上の周波数帯域の全てにおいて、腹部脈波のパワースペクトルの方が大きくなっているというわけではなく、2.5〜3.0Hz付近では逆転し、胸部脈波のパワースペクトルの方が大きくなっている。このことからも、被験者Fの場合には、胸部脈波形と腹部脈波形との間で所定のずれが生じていることがわかる。従って、この周波数解析を利用して、2つの脈波形間のパワースペクトルの比較により、例えば、このようなパワースペクトルの逆転現象が生じている場合には、腹部大動脈瘤を含む動脈疾患有りと判定するように設定することもできる。
【0048】
なお、上記した説明では、エアクッション10とセンサとの組み合わせを2つ用いているだけであるが、3つ以上であってもよい。この場合、上記の2つのエアクッションを当接する位置に加えて、左右の腸骨動脈に対応する部位及び上腹部正中に対応する部位にもエアクッションを当接できるようにし、腹部大動脈瘤が形成される位置の周辺に合計で4〜5箇所のエアクッションを当接できる構成であることが好ましい。このようにすると腹部大動脈瘤の判定がより正確になる。
【0049】
また、エアクッション10は、それぞれ個別のものにするのではなく、図26に示したように、各エアクッション10を所定の位置で保持可能な保持部材40に装填した構成とすることもできる。このようにすれば、各エアクッション10を、臍周辺の左右2箇所、左右の腸骨動脈対応箇所等に容易にセッティングすることができる。保持部材40としては三次元立体編物を用いることが好ましい。三次元立体編物は、人の筋肉の荷重特性に近似した荷重特性、すなわち、直径30mmの小さな圧縮板で加圧して得られる荷重−たわみ特性と、それよりも大きな直径98mmの圧縮板で加圧した場合の荷重−たわみ特性とが異なるという特性を有している。従って、三次元立体編物を人体に押し付けられても、腹部大動脈瘤を強く圧迫せずに測定することができる。
【0050】
また、上記した説明では、エアクッション10をいずれも腹側に当接しているが、背側に当接するようにしてもよい。この場合、図26に示したようにマットレス型として構成することもできるが、椅子のシートバック部に三次元立体編物を張設し、その表面又は裏面にエアクッションを保持した構成とすることもできる。さらに、エアクッションを三次元立体編物の一面に保持したベルト状あるいは腹巻き状のものを製作し、この三次元立体編物を腹部に巻き付けるようにしても、容易に脈波を測定できる。いずれも三次元立体編物が所定の張力を備えているため、脈波の伝達性に優れていると共に、脈波によって変動するエアクッションの確実な支えとして機能する。
【符号の説明】
【0051】
1 腹部大動脈瘤検出装置
2 生体信号測定手段
3 脈波検出手段
10 エアクッション
11 表側エアクッション
11a 空気袋
11b 三次元立体編物
12 裏側エアクッション
12a 空気袋
12b 三次元立体編物
15,16 ビーズ発泡樹脂弾性部材
31 受信部
32 演算部
33 判別部
40 保持部材
【技術分野】
【0001】
本発明は、脈波によって生じるエアクッションの空気圧変動の信号波形から腹部大動脈瘤の有無を簡易的に検出する腹部大動脈瘤検出装置に関する。
【背景技術】
【0002】
腹部大動脈瘤は、超音波診断装置やCTを用いれば容易に発見することができる。しかし、超音波診断装置やCT自体、相応の導入コストがかかるため、医療機関における保有台数にも限りがある。従って、自覚症状を訴えた患者、医師が触診により疑わしいと判定した患者に対して超音波診断装置やCTによる検査が行われていることが一般的である。しかし、より安価で簡易な装置が開発できれば、例えば、健康診断などにおいて、健診対象者の全てに対して腹部大動脈瘤に関する検査を行って、自覚症状の有無に拘わらず腹部大動脈瘤の早期発見を行うことが可能である。このため、超音波診断装置やCTなどよりも診断精度が劣るとしても、安価に入手でき、かつ簡易に操作可能な装置であって、腹部大動脈瘤の保有者をある程度の精度でスクリーニングできる装置の開発が望まれていた。
【0003】
一方、本出願人は、例えば、特許文献1として、内部に三次元立体編物を挿入したエアクッションを備え、このエアクッションを例えば人の腰部に対応する部位に配置し、エアクッションの空気圧変動を測定し、得られた空気圧変動の時系列データから人の生体信号を検出し、人の生体の状態を分析するシステムを開示している。また、非特許文献1及び2においても、腰腸肋筋に沿うようにエアクッションを配置して人の生体信号を検出する試みを報告している。
【先行技術文献】
【特許文献】
【0004】
【特許文献1】特開2007−90032号公報
【非特許文献】
【0005】
【非特許文献1】「非侵襲型センサによって測定された生体ゆらぎ信号の疲労と入眠予知への応用」、落合直輝(外6名)、第39回日本人間工学会 中国・四国支部大会 講演論文集、平成18年11月25日発行、発行所:日本人間工学会 中国・四国支部事務局
【非特許文献2】「非侵襲生体信号センシング機能を有する車両用シートの試作」、前田慎一郎(外4名)、第39回日本人間工学会 中国・四国支部大会 講演論文集、平成18年11月25日発行、発行所:日本人間工学会 中国・四国支部事務局
【発明の概要】
【発明が解決しようとする課題】
【0006】
特許文献1及び非特許文献1、2によれば、腰部付近の動脈(大動脈)の脈波を検知し、得られた脈波の時系列信号データを用い、例えば、本出願人が特開2004−344612において提案した手法により入眠予兆を判定することができる。そして、内部に三次元立体編物を挿入したエアクッションを備えた生体信号測定手段は、低コストで製作できるという利点がある。
【0007】
そこで、本発明は、上記エアクッションを備えた生体信号測定手段を用いて、腹部大動脈瘤の簡易的な検出を行う技術を提供することを課題とする。
【課題を解決するための手段】
【0008】
上記課題を解決するため、本発明の腹部大動脈瘤検出装置は、少なくとも、人体の臍周辺における正中を挟んだ左右2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかに当接可能な2つのエアクッションと、前記各エアクッションを当接した部位における脈波によって生じる各空気圧変動を検出する少なくとも2つのセンサとを備えた生体信号測定手段と、前記生体信号測定手段の各センサの出力信号を受信する受信部と、各出力信号の信号波形間のずれ度合いを求める演算部と、前記演算部により求められた信号波形間のずれ度合いから腹部大動脈瘤の有無を判別する判別部とが設定された脈波検出手段とを備えてなることを特徴とする。
【0009】
前記生体信号測定手段は、人体の臍周辺における正中を挟んだ左右2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかに当接可能な前記2つのエアクッションに加えて、さらに、左右2箇所の腸骨動脈に対応する部位及び上腹部正中のうちの1箇所以上にも当接可能な1以上のエアクッションを備えると共に、これらの各エアクッションのそれぞれにも前記センサが設けられていることが好ましい。また、前記生体信号測定手段は、前記各エアクッションに積層されるビーズ発泡樹脂弾性部材を備える構成とすることができる。前記ビーズ発泡樹脂弾性部材は、ビーズ発泡体と該ビーズ発泡体の外面を被覆する被覆材とからなることが好ましい。
【0010】
前記生体信号測定手段は、前記各エアクッションが三次元立体編物の一面に保持されており、該三次元立体編物を腹部に巻き付けて脈波を測定する構成とすることができる。また、前記生体信号測定手段は、椅子のシートバック部に張設された三次元立体編物の一面に前記各エアクッションが保持されており、該椅子に着席して脈波を測定する構成とすることができる。
【発明の効果】
【0011】
本発明によれば、人体の臍周辺における正中を挟んだ左右2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかに前記2つのエアクッションを当接し、各エアクッションを当接した部位における脈波によって生じる各空気圧変動を検出するセンサからの出力信号の信号波形間のずれ度合いを求め、各信号波形間のずれ度合いから腹部大動脈瘤の有無を時系列解析や周波数解析などによって判別する構成である。腹部大動脈瘤が生じている場合、その部分の血管には部分的な変形が生じており、正常状態と比較して上下、左右、前後のいずれかの方向に膨出している。このため、腹部大動脈瘤が生じている付近における少なくとも左右2箇所、上下2箇所又は前後(腹側と背側)2箇所のエアクッションの空気圧変動の信号波形間に変位が生じ、一方の振幅が大きくなったり、波形に極端な乱れが生じたりする。従って、信号波形間のずれ度合いから腹部大動脈瘤の有無を簡易に判別することができる。生体信号測定手段はエアクッションから構成されるため、低コストで製造可能であることから、装置全体の導入コストも低く抑えられる。このため、健康診断等において、超音波診断装置やCTでの精密検査を要するか否かを判定する簡易的なスクリーニング手段として有益である。
【図面の簡単な説明】
【0012】
【図1】図1は、本発明の一の実施形態に係る腹部大動脈瘤検出装置の概略構成を示した図である。
【図2】図2(a)は、エアクッションにビーズ発泡樹脂弾性部材を積層した状態の一例を示し、図2(b)は、ビーズ発泡樹脂弾性部材の構成を説明するための図である。
【図3】図3は、脈波検出手段の構成を説明するための図である。
【図4】図4は、エアクッションのセンサから検出した信号波形の一例を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形である。
【図5】図5は、指尖容積脈波計から検出した信号波形の一例を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形である。
【図6】図6は、エアクッションのセンサから検出した被験者Aの信号波形の周波数解析結果を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図7】図7は、被験者Aの上記信号波形のカオス解析により得られるアトラクタを示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図8】図8(a)は、被験者Aの上記信号波形のカオス解析により得られる原波形、1階微分波形、2階微分波形のリアプノフ指数の時系列波形であり、(b)は原波形、1階微分波形、2階微分波形のリアプノフ指数の周波数スペクトルを示した図である。
【図9】図9は、エアクッションのセンサから検出した被験者Bの信号波形の周波数解析結果を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図10】図10は、被験者Bの上記信号波形のカオス解析により得られるアトラクタを示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図11】図11(a)は、被験者Bの上記信号波形のカオス解析により得られる原波形、1階微分波形、2階微分波形のリアプノフ指数の時系列波形であり、(b)は原波形、1階微分波形、2階微分波形のリアプノフ指数の周波数スペクトルを示した図である。
【図12】図12は、エアクッションのセンサから検出した被験者Cの信号波形の周波数解析結果を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図13】図13は、被験者Cの上記信号波形のカオス解析により得られるアトラクタを示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図14】図14(a)は、被験者Cの上記信号波形のカオス解析により得られる原波形、1階微分波形、2階微分波形のリアプノフ指数の時系列波形であり、(b)は原波形、1階微分波形、2階微分波形のリアプノフ指数の周波数スペクトルを示した図である。
【図15】図15は、エアクッションのセンサから検出した被験者Dの信号波形の周波数解析結果を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図16】図16は、被験者Dの上記信号波形のカオス解析により得られるアトラクタを示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図17】図17(a)は、被験者Dの上記信号波形のカオス解析により得られる原波形、1階微分波形、2階微分波形のリアプノフ指数の時系列波形であり、(b)は原波形、1階微分波形、2階微分波形のリアプノフ指数の周波数スペクトルを示した図である。
【図18】図18は、エアクッションのセンサから検出した被験者Eの信号波形の周波数解析結果を示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図19】図19は、被験者Eの上記信号波形のカオス解析により得られるアトラクタを示した図であり、(a)は原波形、(b)は1階微分波形、(c)は2階微分波形に関するものである。
【図20】図20(a)は、被験者Eの上記信号波形のカオス解析により得られる原波形、1階微分波形、2階微分波形のリアプノフ指数の時系列波形であり、(b)は原波形、1階微分波形、2階微分波形のリアプノフ指数の周波数スペクトルを示した図である。
【図21】図21は、被験者F,Gの臍を挟んだ上下2カ所にエアクッションを当てて胸部脈波と腹部脈波を検出する試験の試験方法を説明するための図である。
【図22】図22(a)は、エアクッションのセンサから検出した被験者Fの胸部脈波の波形(胸部脈波形)を示した図であり、図22(b)は、エアクッションのセンサから検出した被験者Fの腹部脈波の波形(腹部脈波形)を示した図であり、図22(c)は、エアクッションのセンサから検出した被験者Gの胸部脈波の波形(胸部脈波形)を示した図であり、図22(d)は、エアクッションのセンサから検出した被験者Gの腹部脈波の波形(腹部脈波形)を示した図である。
【図23】図23(a)は、図22(a),(b)に示した被験者Fの胸部脈波形と腹部脈波形の一部拡大図であり、図23(b)は、図22(c),(d)に示した被験者Gの胸部脈波形と腹部脈波形の一部拡大図である。
【図24】図24は、図22(a),(b)に示した被験者Fの胸部脈波形と腹部脈波形の位相差と、図22(c),(d)に示した被験者Gの胸部脈波形と腹部脈波形の位相差を示した図である。
【図25】図25(a)は、図22(a),(b)に示した被験者Fの胸部脈波形と腹部脈波形の周波数解析を示した図であり、図22(c),(d)に示した被験者Gの胸部脈波形と腹部脈波形の周波数解析結果を示した図である。
【図26】図26は、エアクッションを保持部材に装填した状態の一例を示した図である。
【発明を実施するための形態】
【0013】
以下、図面に示した本発明の実施形態に基づき、本発明をさらに詳細に説明する。図1は、本実施形態に係る腹部大動脈瘤検出装置1を示した構成図であり、生体信号測定手段2と脈波検出手段3とを備えて構成される。
【0014】
生体信号測定手段2は、2つのエアクッション10を備えて構成されている。各エアクッション10は、各エアクッション10は、図1及び図2に示したように、表側エアクッション11と裏側エアクッション12とが積層されて構成される。表側エアクッション11及び裏側エアクッション12とは、それぞれ空気袋11a,12aとその内部に配置される復元力付与部材としての三次元立体編物11b,12bとを備えている。表側エアクッション11と裏側エアクッション12とは、長手方向に沿った一方の側縁同士が接合され、接合された側縁を中心にして2つ折りにされて、相互に重ね合わせられて用いられる。
【0015】
表側エアクッション11にはセンサ取付チューブ13が設けられ、その内側に空気圧変動を測定するセンサ14が固定されている。なお、センサ取付チューブ13は密閉されている。また、空気袋11a,12aの素材は限定されるものではないが、例えば、ポリウレタンエラストマー(例えば、シーダム株式会社製、品番「DUS605−CDR」)からなるシートを用いて形成することができる。センサ14としては、空気袋11a,12a内の空気圧を測定できるものであればよく、例えば、コンデンサ型マイクロフォンセンサを用いることができる。また、空気袋11a,12a内に装填される三次元立体編物11b,12bは、例えば、特開2002−331603号公報に開示されているように、互いに離間して配置された一対のグランド編地と、該一対のグランド編地間を往復して両者を結合する多数の連結糸とを有する立体的な三次元構造となった編地である。
【0016】
一方のグランド編地は、例えば、単繊維を撚った糸から、ウェール方向及びコース方向のいずれの方向にも連続したフラットな編地組織(細目)によって形成され、他方のグランド編地は、例えば、短繊維を撚った糸から、ハニカム状(六角形)のメッシュを有する編み目構造に形成されている。もちろん、この編地組織は任意であり、細目組織やハニカム状以外の編地組織を採用することもできるし、両者とも細目組織を採用するなど、その組み合わせも任意である。連結糸は、一方のグランド編地と他方のグランド編地とが所定の間隔を保持するように、2つのグランド編地間に編み込んだものである。このような三次元立体編物としては、例えば、以下のようなものを用いることができる。なお、各三次元立体編物は、必要に応じて複数枚積層して用いることもできる。
【0017】
(1)製品番号:49076D(住江織物(株)製)
材質:
表側のグランド編地・・・300デシテックス/288fのポリエチレンテレフタレート繊維仮撚加工糸と700デシテックス/192fのポリエチレンテレフタレート繊維仮撚加工糸との撚り糸
裏側のグランド編地・・・450デシテックス/108fのポリエチレンテレフタレート繊維仮撚加工糸と350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメントとの組み合わせ
連結糸・・・・・・・・・350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメント
【0018】
(2)製品番号:49013D(住江織物(株)製)
材質:
表側のグランド編地・・・450デシテックス/108fのポリエチレンテレフタレート繊維仮撚加工糸の2本の撚り糸
裏側のグランド編地・・・450デシテックス/108fのポリエチレンテレフタレート繊維仮撚加工糸の2本の撚り糸
連結糸・・・・・・・・・350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメント
【0019】
(3)製品番号:69030D(住江織物(株)製)
材質:
表側のグランド編地・・・450デシテックス/144fのポリエチレンテレフタレート繊維仮撚加工糸の2本の撚り糸
裏側のグランド編地・・・450デシテックス/144fのポリエチレンテレフタレート繊維仮撚加工糸と350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメントとの組み合わせ
連結糸・・・・・・・・・350デシテックス/1fのポリトリメチレンテレフタレートモノフィラメント
【0020】
(4)旭化成せんい(株)製の製品番号:T24053AY5−1S
【0021】
生体信号測定手段2は、図2(a)に示したように、上記したエアクッション10に積層されるビーズ発泡樹脂弾性部材15,16を備えていることが好ましい。ビーズ発泡樹脂弾性部材15,16を人体とエアクッション10との間に設けることにより、測定時において、人がエアクッション10の異物感を感じにくくなる。その一方、ビーズ発泡樹脂弾性部材15,16として、好ましくは、直径30mmの加圧板によってたわみ量1mmまで加圧した際の荷重−たわみ特性から得られるバネ定数が、エアクッション10よりも高いものを用いる。これにより、生体信号を増幅してエアクッション10に伝達できる。
【0022】
ビーズ発泡樹脂弾性部材15,16は、図2(b)に示したように、略平板状に形成されたビーズ発泡体15a,16aと、その外面に貼着される被覆材15b,16bとから構成されている。ビーズ発泡体15a,16aとしては、ポリスチレン、ポリプロピレン及びポリエチレンのいずれか少なくとも一つを含む樹脂のビーズ法による発泡成形体が用いられる。なお、発泡倍率は任意であり限定されるものではない。被覆材15b,16bは、ビーズ発泡体15a,16aの外面に接着により貼着され、高い伸度と回復率を有する素材とすることが好ましい。例えば、伸度200%以上、100%伸長時の回復率が80%以上である弾性繊維不織布を用いることができる。このような弾性繊維不織布としては、例えば、特開2007−92217号公報に開示された熱可塑性エラストマー弾性繊維が相互に溶融接着された不織布を用いることができる。具体的には、KBセーレン(株)製、商品名「エスパンシオーネ」を用いることができる。
【0023】
また、ビーズ発泡樹脂弾性部材15,16は、図2(a)に示したように、2枚を重ねて用いるだけでなく、エアクッション10を挟むようにして用いることもできる。この場合、2枚のビーズ発泡樹脂弾性部材は、バネ定数の異なるものを用いることが好ましい。例えば、上記の商品名「エスパンシオーネ」で被覆したビーズ発泡樹脂弾性部材を用いた場合、これと共に使用する他のビーズ発泡樹脂弾性部材は、被覆材として伸縮性の小さい素材、例えば、熱可塑性ポリエステルからなる不織布を用いた構成とすることができる。このような被覆材としては、例えば、帝人(株)製のポリエチレンナフタレート(PEN)繊維(1100dtex)から形成した2軸織物(縦:20本/inch、横:20本/inch)を用いることができる。
【0024】
次に、脈波検出手段3の構成について図3に基づき説明する。脈波検出手段3は、受信部31、演算部32、判別部33とが設定されたコンピュータから構成される。
【0025】
受信部31は、生体信号測定手段2を構成する2つのエアクッション10の空気圧変動に応じた出力信号を上記した2つのセンサ14からそれぞれ受信する機能を備えている。
【0026】
演算部32は、受信部31により受信した各出力信号の信号波形を演算し、それらを重ね合わせて2つの信号波形間のずれ度合いを求める。信号波形は、出力信号の原波形、原波形を微分した波形(1階微分波形、2階微分波形等)のいずれかを採用することができる。
【0027】
ここで、図4(a)は、一つのエアクッション10が臍部に当たるように被験者がうつぶせになって得られたセンサの出力信号の原波形であり、図4(b)は、図4(a)の1階微分波形、図4(c)は、図4(a)の2階微分波形である。また、図5(a)は、図4と同時に測定した被験者の指尖容積脈波(指尖容積脈波計((株)アムコ製、フィンガークリッププローブ SR−5C)により測定)の原波形であり、図5(b)は、図5(a)の1階微分波形、図5(c)は、図5(a)の2階微分波形である。これらの図から明らかなように、原波形、1階微分波形、2階微分波形のいずれも、エアクッションのセンサの出力信号は、指尖容積脈波の出力信号とよく一致しており、エアクッションによって腹部大動脈の脈波が捉えられていることがわかる。
【0028】
腹部大動脈溜が生じていると、上記のように血管がうねるなどの部分的な変形が生じていたり、あるいは、動脈硬化を併発していることが多いため血管雑音が伴いやすい。従って、2つのセンサから得られる信号波形は、一方の振幅が他方よりも大きく、きれいな波形にならなかったり、一方の信号波形が他方の信号波形では見られない複数の凹凸部を有していたり、あるいは所定の位相差があったりするなど、信号波形間にずれがある。そこで、演算部32はこのずれ度合いを求める。ずれ度合いの求め方としては、例えば、一方の振幅を基準として、他方の振幅の比率を求めることで算出することができる。
【0029】
判別部33は、演算部32により求められた信号波形間のずれ度合いに基づき、腹部大動脈瘤の有無を判別する。例えば、上記の振幅の違いからずれ度合いを求める場合には、一方の振幅を基準とした他方の振幅の比率が所定の閾値以上の場合に、「腹部大動脈瘤:有」と判定する。もちろん、本発明は、腹部大動脈瘤を生じている可能性のある人をスクリーニングすることが目的であり、ここでいう「腹部大動脈瘤:有」とは、腹部大動脈瘤が生じていることが信号波形の比較から疑われる人の意味であり、「腹部大動脈瘤:有」と判定された場合でも、実際には腹部大動脈瘤を有していない人が含まれる場合もある。閾値を低めに設定すれば、腹部大動脈瘤を生じていない者が含まれる比率が高くなるが、腹部大動脈瘤を有する人を大雑把に絞り込むという目的に照らすと、閾値はそれほど高く設定しない方が好ましい。
【0030】
なお、原波形を微分処理すると、原波形が強調されるが、1階微分と2階微分とでは強調の程度が異なるため、原波形の形やそれに含有される高周波成分のレベルにより、いずれの結果を採用すべきかは一概には言えず個人差がある。従って、原波形、1階微分波形、2階微分波形のいずれを採用するかは、得られた各信号波形から、比較しやすいものを操作者が選択するように設定することができる。また、原波形、1階微分波形、2階微分波形の全ての信号波形をそれぞれ比較して、そのうちの1つだけでも、ずれ度合いが閾値を超えている場合に腹部大動脈瘤「有」と判定するようにすれば、見落とし率が小さくなる。
【0031】
本実施形態に係る腹部大動脈瘤検出装置1を用いて腹部大動脈瘤の検査を行う場合には、検査対象者を寝具にうつぶせにして腹側における臍周辺の正中を挟んだ2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかにエアクッション10を当接するだけでよい。これにより、エアクッション10のセンサの出力信号が、脈波検出手段3の受信部31により受信され、演算部32、判別部33によって処理されて、2つの信号波形が比較される。その結果、信号波形間のずれ度合いが所定の閾値以上か否かにより、腹部大動脈溜の疑いがあるかどうかが上記のように出力される。
【0032】
なお、信号波形間のずれ度合いは、振幅のほか、周波数、位相、低周波及び高周波成分の割合、切痕の位置などから判断することもできる。また、振幅のずれ度合いは、時系列の信号波形を比較して求めることができ、低周波及び高周波成分の割合のずれ度合いは、信号波形を周波数解析することによって求めることができる。
【0033】
図6〜図20は、被験者A〜Eが、一つのエアクッション10を臍部に当ててうつぶせになってセンサから得られた脈波信号のうち、ノイズ及び体動のない200秒間を抽出して分析した結果である。図6、図9、図12、図15、図18は各被験者の脈波信号の周波数解析結果であり、図7、図10、図13、図16、図19は各被験者の脈波信号をカオス解析して描いたアトラクタを示した図であり、これらの図において、(a)は脈波信号の原波形の解析結果を、(b)は1階微分波形の解析結果を、(c)は2階微分波形の解析結果をそれぞれ示す。また、図8、図11、図14、図17、図20の(a)はカオス解析により得られるリアプノフ指数の時系列波形であり、(b)はリアプノフ指数の周波数スペクトルを示した図である。この周波数解析結果、アトラクタ図形、リアプノフ指数の波形は、出力結果に乱れがなく安定しているほど、被験者が健康であることを示し、出力結果の乱れが大きいほど、体調不良であることを示している。
【0034】
なお、各被験者の性別、年齢、体格等は次の通りであった。
・被験者A:男性、33歳、身長168cm、体重65.0kg、体脂肪率13.1%
・被験者B:男性、31歳、身長171cm、体重65.2kg、体脂肪率16.4%
・被験者C:女性、25歳、身長159cm、体重44.8kg、体脂肪率23.4%
・被験者D:男性、36歳、身長167cm、体重70.9kg、体脂肪率24.2%
・被験者E:女性、31歳、身長159cm、体重55.0kg、体脂肪率25.5%
【0035】
被験者Aは、図6から、原波形、1階微分波形、2階微分波形のいずれの周波数解析も、200秒間安定した周波数になっており、健康であることがわかる。図7のアトラクタを考察した場合も、全体的には安定しているが、2階微分波形のアトラクタにおいては、「00:02:13〜00:02:43」の時間帯において乱れがあり、基本的には健康であるが、若干の不調を読み取ることができる。図8のリアプノフ指数の解析においても図7と同様の傾向を読み取ることができる。被験者Aのデータの場合、この若干の不調を示すデータは、2階微分波形において読み取れることから、被験者Aの場合には、特に、2階微分波形による解析が適しているということもわかる。
【0036】
被験者Bは、図9の周波解析がいずれも安定していることから健康であることがわかる。図10及び図11においても、アトラクタ及び波形に乱れがほとんどなく、健康であると読み取れるが、「00:02:47〜00:03:17」の時間帯において若干の乱れがある。これは、試験後半において、疲労が生じたことによるものと思われる。
【0037】
被験者Cは、図12の周波数解析が安定していると共に、図13及び図14のアトラクタ及び波形が共に乱れがなく、健康であるとこがわかる。
【0038】
被験者Dは、図15の周波数解析において、周波数が安定していない部分があり、図16のアトラクタ図形においても乱れの生じた図形が多い。また、図17のリアプノフ指数の波形では、特に原波形の乱れが大きいことが読み取れる。これらのことから、被験者Dは体力があまり高くなく、しかも体調が万全でないと判断できる。
【0039】
被験者Eは、図18の周波数解析において、周波数が安定していない部分があり、図19及び図20のアトラクタ図形及び波形に乱れが多いことが読み取れる。従って、被験者Eも体調が万全でないと判断できる。実際、被験者Eの場合、本試験後に循環器系疾患が発見されており、本試験による解析が体調の健康度合いを示す指標として有用であることがわかった。
【0040】
図6〜図20の各解析手法は、検出された信号波形から、被験者の健康状態を判定できることを示したものであるが、この解析手法を複数のエアクッション10から得られる信号波形の解析に応用した場合には、例えば、腹部大動脈瘤の生じていない範囲に当接されたエアクッション10の検出信号の解析結果と、腹部大動脈瘤が生じている範囲に当接されたエアクッション10の検出信号の解析結果とは明確に異なることになると予測できる。従って、それらの解析結果のずれ度合い(周波数解析の差、アトラクタ図形の形の差、リアプノフ指数波形の差等)を求め、所定の閾値を設定することにより、腹部大動脈瘤の有無を判定できる。
【0041】
このように、本実施形態によれば、極めて簡易な構成のエアクッション10を用いて測定することができるため、導入コストが低く抑えられると共に、検査方法も簡易であることから、健康診断などにおいて、腹部大動脈瘤の有無(疑いがあるか否か)を容易に判定することができる。
【0042】
図21〜図25は、腹部大動脈瘤を患ってはいないが動脈疾患を指摘されている50歳代の肥満体型の被験者F(性別:男性、身長:174cm、体重:94kg)と、健康な20歳代の標準体型の被験者G(性別:女性、身長:167cm、体重:57kg)についての生体信号の解析結果を示した図である。試験は、被験者F,Gを寝具(マットレス)にうつぶせに寝かせ、図21に示したように、臍を挟んだ上下2カ所に相当する、臍よりやや下側の腹部とその位置から約30cm上側に離間した胸部に、一方のエアクッション10と他方のエアクッション10を当てて行った。
【0043】
図22(a)〜(d)は各エアクッション10,10から得られた被験者F,Gの検出信号の原波形であるが、その一部を拡大した波形が図23(a),(b)である。臍を挟んだ下側に位置する一方のエアクッション10から得られた検出信号の波形を「腹部脈波形」として、臍を挟んだ上側に位置する他方のエアクッション10から得られた検出信号の波形を「胸部脈波形」として示している。図23の(a)と(b)を比較すると、腹部脈波形と胸部脈波形の位相差は、被験者Gのデータ(図23(b))に対し、被験者Fのデータ(図23(a))の方が大きいことがわかる。すなわち、2つのエアクッション10,10の離間距離が30cm程度に過ぎないため、動脈疾患がなければ、被験者Gのように2つの脈波形間に位相差はほとんど生じないはずである。しかし、動脈疾患を有する場合には、被験者Fのように2つの脈波形間に明らかな位相差が生じることがわかる。図22の900秒間の原波形について腹部脈波形と胸部脈波形の位相差を計算により求めたところ、図24に示したように、被験者Fのデータの位相差が約60度であったのに対し、被験者Gのデータの位相差が約5度以下であり、動脈疾患を有する被験者Fの位相差が顕著に大きかった。
【0044】
そこで、脈波検出手段3を構成するコンピュータの記憶部に、腹部大動脈瘤を有する場合の位相差の大きさを閾値として予め設定しておく。閾値は、腹部大動脈瘤を有する場合の被験者の位相差をできるたけ多く測定して定めておくことが好ましい。そして、実際の判定に当たっては、演算部32において、2つの脈波形間の位相差を求めた後、判別部33において、求めた位相差を記憶部に設定されている上記閾値と比較して行う。その結果、予め記憶されている所定の閾値以上であったならば、腹部大動脈瘤を有すると判定する。これにより腹部大動脈瘤の有無を容易に判定できる。なお、閾値を調整すれば、腹部大動脈瘤に至る前の動脈疾患を判定することができる。
【0045】
また、図23(a)では、492〜495秒の間のように腹部脈波形の振幅が胸部脈波形の振幅よりも2倍前後大きいといったように、2つの脈波形間の振幅差が大きい。これに対し、図23(b)では、2つの脈波間の振幅差は図23(a)よりも小さい。そこで、この振幅差について所定の閾値を設定することにより、腹部大動脈瘤の有無を含む動脈疾患の判定に利用できる。
【0046】
また、被験者Gの場合には、胸部脈波形において切痕が捉えられている。切痕は、拍出期の終わりに大動脈が急に閉鎖することに対応して現れる信号であり、健康な人であればこれが捉えられる。これに対し、被験者Fの場合には、胸部脈波形及び腹部脈波形のいずれも切痕が明確には捉えられていない。従って、この切痕の有無も、腹部大動脈瘤を含む動脈疾患の有無の判定に利用できる。
【0047】
一方、図25の周波数解析結果を見ると、(b)の被験者Gのデータは、0.4Hz以上の周波数帯域の全てにおいて、胸部脈波よりも腹部脈波のパワースペクトルが大きくなっている。これに対し、(a)の被験者Fのデータは、0.4Hz以上の周波数帯域の全てにおいて、腹部脈波のパワースペクトルの方が大きくなっているというわけではなく、2.5〜3.0Hz付近では逆転し、胸部脈波のパワースペクトルの方が大きくなっている。このことからも、被験者Fの場合には、胸部脈波形と腹部脈波形との間で所定のずれが生じていることがわかる。従って、この周波数解析を利用して、2つの脈波形間のパワースペクトルの比較により、例えば、このようなパワースペクトルの逆転現象が生じている場合には、腹部大動脈瘤を含む動脈疾患有りと判定するように設定することもできる。
【0048】
なお、上記した説明では、エアクッション10とセンサとの組み合わせを2つ用いているだけであるが、3つ以上であってもよい。この場合、上記の2つのエアクッションを当接する位置に加えて、左右の腸骨動脈に対応する部位及び上腹部正中に対応する部位にもエアクッションを当接できるようにし、腹部大動脈瘤が形成される位置の周辺に合計で4〜5箇所のエアクッションを当接できる構成であることが好ましい。このようにすると腹部大動脈瘤の判定がより正確になる。
【0049】
また、エアクッション10は、それぞれ個別のものにするのではなく、図26に示したように、各エアクッション10を所定の位置で保持可能な保持部材40に装填した構成とすることもできる。このようにすれば、各エアクッション10を、臍周辺の左右2箇所、左右の腸骨動脈対応箇所等に容易にセッティングすることができる。保持部材40としては三次元立体編物を用いることが好ましい。三次元立体編物は、人の筋肉の荷重特性に近似した荷重特性、すなわち、直径30mmの小さな圧縮板で加圧して得られる荷重−たわみ特性と、それよりも大きな直径98mmの圧縮板で加圧した場合の荷重−たわみ特性とが異なるという特性を有している。従って、三次元立体編物を人体に押し付けられても、腹部大動脈瘤を強く圧迫せずに測定することができる。
【0050】
また、上記した説明では、エアクッション10をいずれも腹側に当接しているが、背側に当接するようにしてもよい。この場合、図26に示したようにマットレス型として構成することもできるが、椅子のシートバック部に三次元立体編物を張設し、その表面又は裏面にエアクッションを保持した構成とすることもできる。さらに、エアクッションを三次元立体編物の一面に保持したベルト状あるいは腹巻き状のものを製作し、この三次元立体編物を腹部に巻き付けるようにしても、容易に脈波を測定できる。いずれも三次元立体編物が所定の張力を備えているため、脈波の伝達性に優れていると共に、脈波によって変動するエアクッションの確実な支えとして機能する。
【符号の説明】
【0051】
1 腹部大動脈瘤検出装置
2 生体信号測定手段
3 脈波検出手段
10 エアクッション
11 表側エアクッション
11a 空気袋
11b 三次元立体編物
12 裏側エアクッション
12a 空気袋
12b 三次元立体編物
15,16 ビーズ発泡樹脂弾性部材
31 受信部
32 演算部
33 判別部
40 保持部材
【特許請求の範囲】
【請求項1】
少なくとも、人体の臍周辺における正中を挟んだ左右2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかに当接可能な2つのエアクッションと、前記各エアクッションを当接した部位における脈波によって生じる各空気圧変動を検出する少なくとも2つのセンサとを備えた生体信号測定手段と、
前記生体信号測定手段の各センサの出力信号を受信する受信部と、各出力信号の信号波形間のずれ度合いを求める演算部と、前記演算部により求められた信号波形間のずれ度合いから腹部大動脈瘤の有無を判別する判別部とが設定された脈波検出手段と
を備えてなることを特徴とする腹部大動脈瘤検出装置。
【請求項2】
前記生体信号測定手段は、人体の臍周辺における正中を挟んだ左右2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかに当接可能な前記2つのエアクッションに加えて、さらに、左右2箇所の腸骨動脈に対応する部位及び上腹部正中のうちの1箇所以上にも当接可能な1以上のエアクッションを備えると共に、これらの各エアクッションのそれぞれにも前記センサが設けられていることを特徴とする請求項1記載の腹部大動脈瘤検出装置。
【請求項3】
前記生体信号測定手段は、前記各エアクッションに積層されるビーズ発泡樹脂弾性部材を備えることを特徴とする請求項1又は2記載の腹部大動脈瘤検出装置。
【請求項4】
前記ビーズ発泡樹脂弾性部材は、ビーズ発泡体と該ビーズ発泡体の外面を被覆する被覆材とからなることを特徴とする請求項3記載の腹部大動脈瘤検出装置。
【請求項5】
前記生体信号測定手段は、前記各エアクッションが三次元立体編物の一面に保持されており、該三次元立体編物を腹部に巻き付けて脈波を測定する構成であることを特徴とする請求項1〜4のいずれか1に記載の腹部大動脈瘤検出装置。
【請求項6】
前記生体信号測定手段は、椅子のシートバック部に張設された三次元立体編物の一面に前記各エアクッションが保持されており、該椅子に着席して脈波を測定する構成であることを特徴とする請求項1〜4のいずれか1に記載の腹部大動脈瘤検出装置。
【請求項1】
少なくとも、人体の臍周辺における正中を挟んだ左右2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかに当接可能な2つのエアクッションと、前記各エアクッションを当接した部位における脈波によって生じる各空気圧変動を検出する少なくとも2つのセンサとを備えた生体信号測定手段と、
前記生体信号測定手段の各センサの出力信号を受信する受信部と、各出力信号の信号波形間のずれ度合いを求める演算部と、前記演算部により求められた信号波形間のずれ度合いから腹部大動脈瘤の有無を判別する判別部とが設定された脈波検出手段と
を備えてなることを特徴とする腹部大動脈瘤検出装置。
【請求項2】
前記生体信号測定手段は、人体の臍周辺における正中を挟んだ左右2箇所、臍を挟んだ上下2箇所、又は、臍周辺における腹側と背側の2箇所のいずれかに当接可能な前記2つのエアクッションに加えて、さらに、左右2箇所の腸骨動脈に対応する部位及び上腹部正中のうちの1箇所以上にも当接可能な1以上のエアクッションを備えると共に、これらの各エアクッションのそれぞれにも前記センサが設けられていることを特徴とする請求項1記載の腹部大動脈瘤検出装置。
【請求項3】
前記生体信号測定手段は、前記各エアクッションに積層されるビーズ発泡樹脂弾性部材を備えることを特徴とする請求項1又は2記載の腹部大動脈瘤検出装置。
【請求項4】
前記ビーズ発泡樹脂弾性部材は、ビーズ発泡体と該ビーズ発泡体の外面を被覆する被覆材とからなることを特徴とする請求項3記載の腹部大動脈瘤検出装置。
【請求項5】
前記生体信号測定手段は、前記各エアクッションが三次元立体編物の一面に保持されており、該三次元立体編物を腹部に巻き付けて脈波を測定する構成であることを特徴とする請求項1〜4のいずれか1に記載の腹部大動脈瘤検出装置。
【請求項6】
前記生体信号測定手段は、椅子のシートバック部に張設された三次元立体編物の一面に前記各エアクッションが保持されており、該椅子に着席して脈波を測定する構成であることを特徴とする請求項1〜4のいずれか1に記載の腹部大動脈瘤検出装置。
【図1】
【図2】
【図3】
【図4】
【図5】
【図7】
【図8】
【図10】
【図11】
【図13】
【図14】
【図16】
【図17】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図6】
【図9】
【図12】
【図15】
【図18】
【図2】
【図3】
【図4】
【図5】
【図7】
【図8】
【図10】
【図11】
【図13】
【図14】
【図16】
【図17】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図6】
【図9】
【図12】
【図15】
【図18】
【公開番号】特開2010−284517(P2010−284517A)
【公開日】平成22年12月24日(2010.12.24)
【国際特許分類】
【出願番号】特願2010−112372(P2010−112372)
【出願日】平成22年5月14日(2010.5.14)
【出願人】(594176202)株式会社デルタツーリング (111)
【出願人】(504136568)国立大学法人広島大学 (924)
【Fターム(参考)】
【公開日】平成22年12月24日(2010.12.24)
【国際特許分類】
【出願日】平成22年5月14日(2010.5.14)
【出願人】(594176202)株式会社デルタツーリング (111)
【出願人】(504136568)国立大学法人広島大学 (924)
【Fターム(参考)】
[ Back to top ]