FI113564B - Kappaleiden paikannus - Google Patents

Kappaleiden paikannus Download PDF

Info

Publication number
FI113564B
FI113564B FI20010558A FI20010558A FI113564B FI 113564 B FI113564 B FI 113564B FI 20010558 A FI20010558 A FI 20010558A FI 20010558 A FI20010558 A FI 20010558A FI 113564 B FI113564 B FI 113564B
Authority
FI
Finland
Prior art keywords
signal
signals
sources
amplitudes
signal sources
Prior art date
Application number
FI20010558A
Other languages
English (en)
Swedish (sv)
Other versions
FI20010558A (fi
FI20010558A0 (fi
Inventor
Samu Taulu
Matti Kajola
Lauri Parkkonen
Original Assignee
Elekta Neuromag Oy
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 Elekta Neuromag Oy filed Critical Elekta Neuromag Oy
Priority to FI20010558A priority Critical patent/FI113564B/fi
Publication of FI20010558A0 publication Critical patent/FI20010558A0/fi
Priority to JP2002575573A priority patent/JP4203592B2/ja
Priority to EP02708381.5A priority patent/EP1370826B1/en
Priority to CA2440988A priority patent/CA2440988C/en
Priority to PCT/FI2002/000225 priority patent/WO2002077565A1/en
Publication of FI20010558A publication Critical patent/FI20010558A/fi
Priority to US10/665,543 priority patent/US6876196B1/en
Application granted granted Critical
Publication of FI113564B publication Critical patent/FI113564B/fi

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B7/00Measuring arrangements characterised by the use of electric or magnetic techniques
    • G01B7/004Measuring arrangements characterised by the use of electric or magnetic techniques for measuring coordinates of points
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0247Determining attitude
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1113Local tracking of patients, e.g. in a hospital or private home
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/242Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents
    • A61B5/245Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents specially adapted for magnetoencephalographic [MEG] signals
    • A61B5/246Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents specially adapted for magnetoencephalographic [MEG] signals using evoked responses

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Measuring Magnetic Variables (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Description

1 113564
KAPPALEIDEN PAIKANNUS KEKSINNÖN ALA
Esillä oleva keksintö liittyy kappaleiden paikannukseen. Erityisesti esillä olevan keksinnön 5 kohteena on uusi ja parannettu menetelmä kappaleiden paikan ja asennon määrittämiseksi toistensa suhteen sähkömagneettisten signaalien avulla.
KEKSINNÖN TAUSTA
10 Sähkömagneettisiin signaaleihin perustuvaa paikannusmenetelmää on kuvattu hyvin yleisellä, sovellutuskohteesta riippumattomalla tavalla mm. patenttijulkaisuissa US5747996, US4346384 ja DE3326476. Eräässä sovelluksessa laitteeseen kuuluu joukko signaali-15 lähteitä, joukko vastaanottimia ja yksi tai useampi signaaligeneraattori, jolla generoidaan joukko aikamuodoltaan tunnettuja lähetinsignaaleita lähetettäväksi signaalilähteillä. Lisäksi mainituissa patenttijulkaisuissa on kuvattu analyysimenetelmä vastaanottimien 20 ulostulosignaalien käsittelemiseksi ja niiden käyttämiseksi kappaleen paikan laskemisessa suhteessa toi-V.’ seen kappaleeseen. Yhteistä patenttijulkaisuissa kuva- jtj ; tuille laitteille on, että signaalilähettimet on kiin- • nitetty kappaleeseen geometrisesti melko tiukasti ra- 25 joitetulla tavalla.
Patenttijulkaisussa US5747996 ja US4346384 ,··_ vaaditaan lisäksi, että signaalilähteet ovat keskenään ortogonaaliset. Ortogonaalisuuden ansiosta signaali-lähteiden lähettämien signaalien välillä ei ole korre-•**: 30 laatiota eli signaalit eivät vaikuta toisiinsa paikan- '...· nusta häiritsevästi. Patenttijulkaisussa US5747996 J vaaditaan lisäksi, että vastaanottimet ovat samaan ta- .···. soon sijoitettuja keloja. Geometrisilla vaatimuksilla ·’ pyritään helpottamaan ja nopeuttamaan signaalianalyy- 35 siä ja poistamaan mahdollisia paikannustulokseen vai- kuttavia virhelähteitä.
2 113564
Kohdekappaleen koordinaatistossa tunnettujen signaalilähteiden käyttöön perustuvaa paikannusmenetelmää käytetään esimerkiksi magnetoenkefalografiässä (MEG), jossa mitataan ihmisen tai muun eliön hermotoi-5 minnoista peräisin olevia heikkoja, ajasta ja paikasta riippuvia magneettikenttiä. Mitattujen magneettikent-täarvojen perusteella pyritään paikantamaan lähdealueet, jotka synnyttivät havaitun kentän. Magnetoenkefa-lografiassa koehenkilön pää on mahdollisimman lähellä 10 erittäin herkistä suprajohtavista antureista koostuvaa anturi- eli vastaanotinjoukkoa, jonka geometria tunnetaan. Pään paikka mittalaitteen suhteen määritetään käyttäen tunnettuina signaalilähteinä pieniä pään pinnalle kiinnitettyjä keloja, joiden tuottamaa magneet-15 tikenttää voidaan approksimoida magneettisen dipolin kentällä.
Vastaanottimina käytetään mittalaitteen mittausantureita, joita käytetään myös varsinaisten mitattavien aivosignaalien vastaanottamiseen ja mittaa-20 miseen. Menetelmän perusperiaatteet on selostettu esimerkiksi julkaisuissa SQUID'85: Superconducting Quan tum Interference Devices and their Applications, 1985, ss. 93 9-944 sekä Proceedings of the 7th International :: Conference on Biomagnetism, 1989, ss. 693-696.
: 25 Varsinaiset MEG-mittaukset toteutetaan yleen- sä toistomittauksina, joissa esimerkiksi tiettyä ärsy-kettä seuraavaa aivojen tuottamaa vastetta mitataan useita kertoja peräkkäin, ja lasketaan ärsykkeen suhteen aikalukittujen mittaustulosten keskiarvo. Käytet-30 täessä mittaustulosten keskiarvoa, voidaan kohinan ·>·: vaikutus vaimentaa tekijällä, joka on kääntäen verran- nollinen toistojen lukumäärän neliöjuureen. Toistomit- : tausten eräs ongelma on niiden pitkä kesto, minkä • » .**·. vuoksi koehenkilön pää saattaa liikkua mittauksen ai- ’·’ 35 kana. Tästä automaattisesti seuraa, että aivojen tuot- tämän vasteen lähteen paikka muuttuu mittalaitteen 113564 3 suhteen kesken mittauksen ja aiheuttaa siten virhettä lopulliseen analyysiin.
Perinteisesti pää on paikannettu ainoastaan mittauksen alussa siten, että kukin päänpaikannuskela 5 on aktivoitu ja syntynyt magneettikenttä mitattu yksitellen, jolloin paikannusmenetelmä on ollut verrattain hidas. Paikannuksen jälkeen koehenkilöä on pyydetty pitämään pää mahdollisimman liikkumatta toistomittauksen loppuun asti.
10 Mittauksen aikaisesta pään liikkeestä aiheu tuvat virheet voidaan välttää jatkuvalla paikanmitta-uksella. Tällöin mittalaitetta on voitava käyttää samanaikaisesti myös muiden kuin paikannuksessa tuotettavien lähetinsignaalien mittaamiseen. Eräs tapa 15 eliminoida lähetinsignaalien vaikutus mitattavaan hyö-tysignaaliin eli aivojen tuottamaan vastesignaaliin on asettaa lähetinsignaalien taajuudet kauaksi tutkittavalta taajuuskaistalta ja suodattaa mittausdataa sopivasti taajuustasossa. Tällainen ratkaisu esitetään 20 julkaisussa Biomag2000, 12th International Conference on Biomagnetism, Book of Abstracts, s. 188, Peters, H. et ai. Toinen ratkaisu on lähetinsignaalien suodatta-, .· minen vastaanottimien ulostulosignaaleista vähentämällä | lä lähetinsignaaleja vastaavat osuudet mitatuista sig- • 25 naaleista, jolloin täytyy tuntea mitattavien lähetin- • t » » signaalien voimakkuudet ja aaltomuodot.
Pyrittäessä paikantamaan kappaletta jatkuvas-ti tai toistuvasti lyhyin väliajoin on signaali lähet- • · timet aktivoitava yhtä aikaa ja pystyttävä erottamaan 30 samanaikaiset, eri lähettimien synnyttämät komponentit ·>·: mittaussignaaleista. Menetelmän tulisi erottaa taa- juuskomponentit mahdollisimman tehokkaasti ja tarkasti käyttäen mahdollisimman lyhyttä tiedonkeruuaikaa. Yleisesti käytetyssä erotusmenetelmässä taajuudet ja •t 35 tiedonkeruuaika sovitetaan siten, että signaalikom- '·!·* ponentit ovat keskenään ortogonaaliset tarkasteltaval- la aikavälillä. Jos lähetinsignaalin vaihe tunnetaan, 113564 4 niin kunkin signaalikomponentin amplitudi saadaan suoraan laskemalla mittaustuloksista koostuvan signaali-vektorin projektio tutkittavaa signaalikomponenttia vastaavalle kantavektorille, joka koostuu taajuudel-5 taan tunnetun kantafunktion laskennallisista arvoista. Kantavektoreiden ortogonaalisuuteen perustuvia sovelluksia on kuvattu mm. julkaisussa "The use of an MEG device as a 3D digitizer and a motion correction system", de Munck et ai Proceedings of the 12th Interna-10 tional Conference on Biomagnetism, Helsinki, Finland. Tässä kuvauksessa epäortogonaalisuuden vaikutus on otettu huomioon periaatteellisella tasolla. Kuvatussa paikannusmenetelmässä lähetinsignaalien ortogonali-sointi kuitenkin vähentää oleellisesti paikannukseen 15 liittyvän laskennan määrää, joten käytännön toteutuksessa lähetinsignaalit on ortogonalisoitu.
Ortogonaalisuusvaatimus asettaa rajoituksia käytettäville taajuuksille sekä tiedonkeruuajalle ja lisäksi ortogonaalisuusoletus epäortogonaalisille sig-20 naaleille aiheuttaa suuria virheitä laskettuihin amp-litudikertoimiin ja täten myös paikannukseen. Edellä kuvatussa signaalianalyysissä pyritään käyttämään mahdollisimman lyhyeltä aikaväliltä kerättyjä signaaleja, jotta paikannus olisi mahdollisimman reaaliaikaista ja : 25 kappaleiden liike olisi mahdollisimman vähäistä pai- kannusmittauksen tiedonkeruun aikana. Luotettavia pai-kannusmittauksia on tehty magnetoenkef alograf iässä ... käyttäen ainoastaan 100 ms:n pituista tiedonkeruu- aikaa.
30 Tällaisellakin aikavälillä voivat kappaleet ·.: kuitenkin liikkua, mikä huonontaa paikannustulosta.
• < ·
Kappaleiden mahdollisesti suuren liikkeen vuoksi vas- : taanottimien mittaamien signaalien voimakkuudet saat- ,···, tavat vaihdella nopeasti havaittavissa olevan signaa- • « 35 Iin alarajalta aina vastaanottimien dynaamisen alueen ylärajoille asti. Vaihtelu voi olla merkittävää eri-tyisesti pienillä etäisyyksillä, sillä mitatun signaa- 113564 5
Iin voimakkuus on kääntäen verrannollinen kappaleiden etäisyyden kolmanteen potenssiin. Tämän lisäksi samoja lähettimiä saatetaan eri mittauksissa käyttää hyvin erikokoisilla ja vastaanottimiin nähden erilaisilla 5 etäisyyksillä sijaitsevilla kappaleilla. Erilaisissa tilanteissa tehtävien mittausten toistuva onnistuminen edellyttää, että vastaanottimien mittaamien lähetin-signaalien voimakkuudet pysyvät jatkuvasti tiettyjen rajojen sisällä. Ongelma on ratkaistu käyttämällä sää-10 töalgoritmia, joka ohjaa lähettimien tehoa siten, että kaikkien vastaanottimien mittaamien signaalien amplitudit pysyvät jatkuvasti jonkin tietyn alarajan yläpuolella ja jonkin tietyn ylärajan alapuolella. Lähe-tinsignaalien takaisinkytkentää on kuvattu esimerkiksi 15 patenttijulkaisussa US5747996.
KEKSINNÖN TARKOITUS
Keksinnön tarkoituksena on poistaa edellä mainitut epäkohdat tai ainakin merkittävästi lieventää 20 niitä. Erityisesti keksinnön tarkoituksena on tuoda esiin uudentyyppinen menetelmä paikannusmittauksen toteuttamiseksi mahdollisimman nopeasti ja tarkasti. Li-säksi keksinnön tarkoituksena on tuoda esiin mittaus- · menetelmä, joka on laskennallisesti yksinkertainen ja : 25 tehokas ja jossa signaalilähteiden epäortogonaalisuu- desta ei ole haittaa lopputuloksen tarkkuuteen. Edel-leen keksinnön tarkoituksena on tuoda esiin laskenta-*·», menetelmä, jolla voidaan lähes reaaliaikaisesti laskea kohdekappaleen paikka ja siten eliminoida kohdekappa-, 30 leen liikkeestä aiheutuneet virheet varsinaiseen mit- taukseen, esimerkiksi magnetoenkef alograf iässä. Edel-'···* leen keksinnön tarkoituksena on tuoda esiin uudenlai- nen laskennan korjausmenetelmä, jolla mitatut amplitu-dijakaumat voidaan korjata epäortogonaalisuudesta ja ·, 35 mahdollisista muista häiriöistä johtuvien virheiden t · · *·:·" poistamiseksi signaalille määritetystä amplitudista tai amplitudi jakaumasta riippuen vastaanottimien mää- 113564 6 rästä. Esillä olevalle keksinnölle on tunnusomaista se, mitä on esitetty patenttivaatimuksissa.
KEKSINNÖN KUVAUS
5 Keksinnön kohteena on menetelmä, jonka avulla kappaleen paikka ja asento voidaan määrittää toisen kappaleen suhteen sähkömagneettisten signaalien avulla. Keksinnön mukaisessa järjestelyssä on kaksi kappaletta, joista toiseen kappaleeseen on kiinnitetty sig-10 naalilähteitä eli lähettimiä, jotka tuottavat sähkömagneettisia signaaleja, ja toinen kappale sisältää yhden tai useampia vastaanottimia lähetinsignaalien mittaamista varten. Yleensä lähettimiä käsittävä kappale on se, jonka paikka tai asento kiinnostaa ja on 15 mittauksen kohteena. Esimerkiksi biomagneettisissa mittauksissa lähettimiin liittyvä kappale on ihmisen pää tai muu kehon rajattu osa, jonka pinnalle lähettimet sijoitetaan. Keksinnön mukaisella järjestelyllä voidaan pään paikka ja asento selvittää, jolloin aivo-20 jen tuottamien signaalien lähdealueet voidaan ilmaista pään koordinaatistossa. Samoilla vastaanottimina mitataan sekä aivojen että lähettimien tuottamia signaaleita.
• Lähetinkappaleen koordinaatistossa tunnetuis- I | j 25 sa paikoissa sijaitsevat signaalilähettimet voidaan ; aktivoida tuottamaan samanaikaisesti tai vuorotellen eri taajuisia signaaleja siten, että taajuudet ja aal-. ··. tomuodot ovat vapaasti valittavissa. Tämän ansiosta lähettimiin liittyvät järjestelyt, esimerkiksi geomet-, 30 rian ja käytettyjen signaalien suhteen, tulevat huo- ;j mattavasti aikaisempaa yksinkertaisemmiksi. Eri lähet timien tuottamien signaalien amplitudit mitataan toi-sen laitteen vastaanottimella, joiden keskinäinen geo-, \ metriä on joko ennalta tunnettu tai tulee määritettyä 35 paikannuksen aikana. Varsinaisten mittausten kannalta * ‘ riittää, että pystytään selvittämään lähettimien ja • » » vastaanottimien keskinäinen geometria, koska lähetti- 113564 7 mien sijoittelu kohdekappaleessa yleensä tunnetaan. Siten tunnettaessa vastaanottimien ja lähettimien geometria sekä lähettimistä lähetettyjen signaalien amplitudijakauma vastaanottimilla kohdekappaleesta mi-5 tatut varsinaiset kohdesignaalit ja niiden lähtöpiste voidaan paikantaa suhteessa lähettimiin ja siten myös suhteessa kohdekappaleseen.
Keksintö perustuu kappaleeseen kiinnitettyjen sähkömagneettisten tai akustisten lähettimien tuotta-10 mien ja toisen kappaleen vastaanottimien mittaamien signaalien käyttöön analyysissä, jonka tuloksena kappaleiden suhteellinen paikka tai asento tai molemmat voidaan laskea. Paikannus voi olla jatkuvaa käyttäen nopeaa amplitudien laskentamenetelmää ja peräkkäisissä 15 mittauksissa mahdollisesti osittain ajallisesti pääl lekkäin olevia mitattujen signaaliarvojen muodostamia vektoreita.
Keksinnön mukaisessa paikannusmenetelmässä sallitaan epäortogonaaliset lähetinsignaalien kanta-20 vektorit, jolloin lähetettävien signaalien taajuudet, aaltomuoto ja tiedonkeruuaika voidaan valita varsin vapaasti. Laskennassa käytettävä kanta- , vektorisisätuloilla painotettu epäortogonaalinen pro- : ; : jektiomenetelmä on laskennallisesti hyvin nopea ja • 25 tarkka operaatio verrattuna muihin käytettyihin sig- naalianalyysimenetelmiin. Toisin kuin edellä mainitus-sa de Munckin artikkelissa, tässä keksinnössä lähetin- i · ... signaalien epäortogonaalisuuden aiheuttaman lisälas- ’· kennan määrä on käytännössä merkityksetön. Menetelmäs- 30 sä voidaan käyttää lähetinsignaaleja, joiden vaiheet * t · ovat tuntemattomia. Tällöin vaiheet voidaan ratkaista * I ft *,.,: laskemalla amplitudit sopivasti valituille kantavekto- .·. : reille, joilla saa olla vaihe-ero varsinaisiin lähet- » »« timen kantavektoreihin nähden.
I I
'!* 35 Keksinnön mukaisella menetelmällä voidaan edelleen estimoida lähetinsignaalien lisäksi tunnettu-: : jen häiriölähteiden amplitudeja, jolloin niiden häi- 113564 8 ritsevä vaikutus voidaan poistaa. Tällaisia häiriöitä ovat erityisesti verkkotaajuuden perustaajuus ja sen harmoniset komponentit.
Jotta mitattujen ja laskettujen amplitudien 5 yhteensovitus olisi mahdollisimman tarkka, tarvitaan tietoa kanavilla esiintyvistä häiriöistä ja kohinasta. Koska nämä voivat vaihdella ajan funktiona, on edullista mitata kyseiset parametrit paikannusmittauksen yhteydessä. Tämä voidaan toteuttaa vähentämällä mita-10 tuista signaaleista estimoidut signaalit ja käyttämällä jäljelle jääneen signaalin tehoa jollakin taajuuskaistalla .
Keksinnön erään sovelluksen mukaisesti voidaan estimoida signaalikomponentteja, jotka poikkeavat 15 varsinaisista lähettimien tai tunnettujen häiriölähteiden synnyttämistä estimoitavista signaalimuodoista.
Vastaanottimia voidaan käyttää paikannusmittauksen aikana mittaamaan myös muita signaalilähteitä kuin lähetinsignaaleja. Tämä toteutetaan vähentämällä 20 laskettujen lähetinsignaaliamplitudien osuus kunkin vastaanottimen ulostulosignaalista kullakin ajanhet-kellä. Keksinnön ansiosta vähentäminen onnistuu aikai-:J sempaa paremmin, koska signaalien amplitudien ja vai- ; heen estimointi on tarkempaa kuin aiemmin tunnetuissa I 25 laitteissa. Suodatus mahdollistaa jatkuvan paikannuk- ,*··, sen, kun mitataan myös muita signaaleita kuin lähetin- signaaleita.
Vastaanottimissa havaittavien lähetinsignaa-· liamplitudien voimakkuuksia voidaan säätää takaisin- 30 kytkennällä, jossa otetaan huomioon koko vastaanotin- joukon mittaamat signaalit. Näin varmistetaan riittävä tll signaalikohinasuhde kaikissa mittaustilanteissa.
.*. : Keksinnön mukaisen paikannuksen tarkkuutta voidaan parantaa mittaamalla signaalin häiriötaso vä- i · 'I' 35 hentämällä mitatuista signaaleista tunnettujen kan- tasignaalimuotojen avulla määritetyt signaalit. Täl-: ( : löin jäljelle jäävä erosignaali kertoo siitä, kuinka 113564 9 luotettava kukin estimaatti on, ja tätä tietoa voidaan käyttää sovituksen tarkentamiseen ottamalla siinä ko-hinataso huomioon. Häiriötaso voidaan myös mitata estimoimalla signaali tai signaaleja, jotka eroavat lä-5 hetinsignaalien ja tuntemiemme, esimerkiksi verkkotaa-juushäiriön, häiriöiden muodoista. Estimoimalla tällaiselle signaalille amplitudin keksinnön mukaisella menetelmällä saadaan tietoa siitä, kuinka paljon oletettu signaaliavaruuden malli poikkeaa todellisuudes-10 ta.
Esillä olevan keksinnön etuna tunnettuun tekniikkaan verrattuna on, että keksinnön mukaisella järjestelyllä tietyn lähetinjoukon lähettämien ja tietyllä vastaanotinjoukolla, yhdellä tai useammalla vas-15 taanottimella vastaanotettujen signaalien amplitudija-kaumien selvittäminen tulee aikaisempaa tarkemmaksi ja tehokkaammaksi. Samaten keksinnön ansiosta lähetettävien signaalien valinta tulee aikaisempaa merkittävästi vapaammaksi erityisesti taajuuden ja vaiheen suh-20 teen. Signaalien ei keksinnön ansiosta tarvitse olla keskenään ortogonaalisia.
Edelleen keksinnön ansiosta voidaan ottaa huomioon esimerkiksi MEG-mittauksissa pään mahdolliset liikkeet varsinaisen hyötysignaalin mittauksen aikana.
j * · ,·, 25 Lisäksi keksintö mahdollistaa ulkoisten virhelähteiden i t eliminoinnin hyötysignaalin mittaustulosten tarkkuuden parantamiseksi.
I 1 *
Edelleen keksinnön mukainen menetelmä ja laite ovat helposti muunneltavia ja keksinnössä toteutet-30 tava korjauslaskenta voidaan toteuttaa kulloinkin par- ;·' haiten sopivassa mittaus- ja laskentavaiheessa.
» I I ·
/ , KEKSINNÖN YKSITYISKOHTAINEN SELOSTUS
% · ,**. Seuraavassa keksintöä selostetaan yksityis- 35 kohtaisen sovellusesimerkin avulla viittaamalla ohei- t i * seen piirustukseen, joka esittää kaaviomaisesti erästä » » * ·...* keksinnön mukaista mittauslaitejärjestelyä.
10 113564
Seuraavassa kuvataan yksityiskohtaisesti keksinnön eräs toteutustapa. Piirustuksessa esitetään periaatteellisella tasolla keksinnön mukainen mittausjärjestely, johon kuuluu lähetinosa glt ..., gn ja vas-5 taanotinrakenne, johon kuuluu joukko vastaanottimia 1, K. Piirustuksessa esitetään ainoastaan yksi vastaanotinantenni, mutta ammattimiehelle on selvää miten useamman antennin järjestely toteutetaan piirustuksen mukaisena. Vastaanottimella mitataan lähetinsignaalien 10 clk, c2k, ..., cNk ja mahdollisen häiriön i]k amplitudit. Signaalien merkinnöissä parametri k viittaa yhdestä lähetetystä signaalista saatavaan k mittaustulokseen johtuen siitä, että kukin signaali vastaanotetaan lähtökohtaisesti k.-11a vastaanottimella.
15 Lisäksi piirustuksessa esitetään lähettimiä ohjaava takaisinkytkentä AGC, jolle syötteenä on lopulliset mitatut signaaliamplitudit i^. Saatujen amplitudien perusteella takaisinkytkentä ohjaa lähettimien lähetystehoa niin, että signaali-kohinasuhde pysyy ko-20 ko ajan halutulla tasolla parhaan mahdollisen mittaustuloksen saavuttamiseksi. Edelleen piirustuksessa esitetään suodatin AF, joka on järjestetty vastaanottimi- * en 1, ..., k, ja mittauslaitteen lähdön väliin lähetti-millä lähetettyjen signaalien suodattamiseksi varsi-25 naisesta hyötysignaalista sk(t).
* t . ·«, Mittausjärjestelyssä on siis edellä esitetyn mukaisesti K vastaanotinta sisältävä mittalaite, jota ... eksitoidaan samanaikaisesti N:llä paikannettavaan kap- ·** paleeseen kiinnitetyllä signaalilähteellä g. Vastaan- 30 ottimella k kerätään signaalia aikavälin T yli, ja eri lähettimiin g liittyvien signaalikomponenttien e^t), » · · ..., eM(t) (M > N kun halutaan estimoida lähetinsignaa-.·. : leista muitakin kuin niiden peruskomponentteja 1, ..., ,·>, N) amplitudien mittaamista varten projisoidaan aluksi 35 eri ajanhetkiltä kerätyistä signaaliarvoista koostuva ·.:.* signaalivektori eri taajuuksia vastaaville
I · I
:it>: kantavektoreille sekä näiden kanssa samalla taajuudel- 113564 11 la oleville kantavektoreille, joilla on noin 90 asteen vaihe-ero edellisiin kantavektoreihin nähden.
Piirustuksessa projektioiden laskemista on kuvattu jatkuvien funktioiden tapaan mitatun signaalin 5 ja kantafunktioiden tulojen määrätyillä integraaleilla ajan T yli. Huomattakoon, että digitaalisessa toteutuksessa käytetään integroinnin tilalla numeerisena integrointina summausta, joka on oikeastaan jatkuva-muotoisen integroinnin estimaatti. Digitaalinen toteu-10 tus ei kuitenkaan mitenkään ole rajattu pois esillä olevan keksinnön sovellusmahdollisuuksista.
Haluttaessa integrandia voidaan painottaa jollakin ikkunafunktiolla w. Projisoinnin jälkeen saatu M-ulotteinen projektiovektori u korjataan matrii-15 silla F2, joka on (M x M) - ulotteinen kantavektoreista riippuva matriisi siten, että u = F2*u, jolloin u on lähetinsignaalien ja näiden kanssa 90 astetta eri vaiheessa olevien identtisten signaalien sekä muiden estimoitavien aaltomuotojen amplitudikertoimet sisältävä 20 vektori.
Korjauksen muodollinen matemaattinen perustelu on seuraavanlainen. Merkitään eri ajanhetkillä mi-tatuista signaaliarvoista koostuvaa signaalivektoria :.· : s:llä ja kantavektorit vaakavektoreina sisältävää mat- j 25 riisiä E:llä, jolloin mitattu signaali on muotoa s = ET*u, josta saadaan pseudoinverssiratkaisu u = inv (E*ET) *E*s = inv(E*ET)*u, missä T viittaa matriisin transpoosiin ja inv() matriisin käänteismatriisiin. Tällöin matriisi F2 = inv(E*ET) ja F2 on , 30 yksikkömatriisi. Epäortogonaalisuuden vaatima korjaus ;;; voidaan tehdä myös yksin kantavektoreille käyttämällä « t *·;·' matriiseja Fx = inv(E*ET)*E ja F2 = yksikkömatriisi tai yhdessä sekä kantavektoreille että projektiovektorille siten, että F1 ja F2 yhdessä muodostavat korjausoperaa-35 tion. Eri kanavilla mitattujen signaaliamplitudien pe- rusteella voidaan laskea kappaleiden suhteellinen paikka tai asento tai molemmat.
i » · 113564 12
Matriisit F1 ja F2 ovat samat kaikille K: lie vastaanottimelle, joten ainoastaan projektiovektori u täytyy laskea erikseen eri vastaanottimille, ja täten lähetinsignaalien epäortogonaalisuudesta johtuva lisä-5 laskennan määrä on erittäin pieni. Eri kanavilla mitattujen signaaliamplitudien perusteella voidaan laskea kappaleiden suhteellinen paikka tai asento tai molemmat .
Koska laitteella pystytään mittaamaan mieli-10 valtaisten signaalikomponenttien amplitudit, voidaan edellä mainitussa kuvauksessa käytetyt paikallaan pysyvän kappaleen synnyttämät signaalit korvata laajemmalla joukolla kantavektoreita, jotka pystyvät esittämään liikkuvan kappaleen synnyttämät aaltomuodot. Lii-15 kettä voidaan mallintaa esimerkiksi amplitudimoduloi- malla alkuperäistä signaalia vakiosta poikkeavilla aaltomuodoilla. Näin voidaan liikkuvakin kappale paikantaa tarkemmin ja siten ottaa huomioon mittauksen aikana kappaleen liikkeestä johtuva signaaliamplitudi-20 en muutos.
Suodatin AF vähentää kullakin ajanhetkellä mitatusta signaalista lähetinsignaalien cx, ..., cN ja ’.V näiden ajassa siirrettyjen vastineiden amplitudit, jo- 1,· f ten lähettimien signaaleihin perustuva paikannusmitta- : : : 25 us ei häiritse yhtäaikaista varsinaista mittausta.
Takaisinkytkentä AGC toteutetaan laskemalla kunkin mitatun lähetinsignaalin RMS-arvot aikaväliltä ···, T yli vastaanotinjoukon 1...K ja säätämällä tämän pe rusteella yksittäisten lähettimien lähetystehoa siten, 30 että estimoitujen signaalikomponenttien maksimi yli ;; vastaanotinjoukon pysyy suunnilleen vakiona.
'··' Keksintöä ei rajata pelkästään edellä esitet- ’·.· tyjä sovellusesimerkkejä koskevaksi, vaan monet muun- .·**. nokset ovat mahdollisia pysyttäessä patenttivaatimus- *, 35 ten määrittelemän keksinnöllisen ajatuksen puitteissa.

Claims (16)

1. Menetelmä kappaleen paikan ja/tai asennon määrittämiseen ennalta määrätyssä koordinaatistossa, jossa menetelmässä kappaleeseen on järjestetty joukko 5 signaalilähteitä kappaleen koordinaatiston suhteen tunnetusti ja jossa lähetetään signaalilähteistä ennalta määrättyä signaalia, vastaanotetaan vastaanottimella, johon kuuluu 10 ainakin yksi signaalivastaanotin, signaalilähteistä lähetetty signaali, ja lasketaan kappaleen paikka ja/tai asento vastaanotettujen signaalien amplitudien perusteella, tunnettu siitä, että 15 määritetään kustakin lähettimestä vastaan otettujen signaalien amplitudit käyttämällä hyväksi kunkin lähetteen tunnettua aaltomuotoa ja korjaamalla aaltomuotojen korrelaatiosta aiheutuva virhe ratkaisemalla amplitudit saatavasta lineaarisesta ongelmasta, 20 määritetään jokainen signaalilähde erillisenä riippumattomista amplitudeista, ja .. lasketaan kappaleen paikka tarkasteltavalla aikavälillä signaalilähteisiin liittyvien riippumatto-’** ' mien amplitudijakaumien perusteella sovittamalla sig- : 25 naalilähteiden laskennalliset amplitudit vastaanotti- :t(,! mella mitattuihin amplitudeihin. *:"i
2. Patenttivaatimuksen 1 mukainen menetelmä, ;‘tunnettu siitä, että • · · sovitetaan amplitudit siten, että asetetaan 30 signaalilähteiden ja/tai vastaanottimen geometriset * vapaat parametrit arvoihin, joilla laskettujen ja mi-T tattujen amplitudijakaumien välinen ero on pienimmil- * t lään, : : lasketaan signaalilähteiden paikka sekä kap- 35 paleen että mittalaitteen koordinaatistossa vapaille .‘li parametreille asetetuista arvoista, ja 113564 14 lasketaan kappaleen paikka ja/tai asento vastaanottimen suhteen käyttämällä signaalilähteiden tunnettuja paikkoja.
3. Patenttivaatimuksen 2 mukainen menetelmä, 5 tunnettu siitä, että yksittäisen signaalilähteen määrittämiseksi: muodostetaan signaalilähdekohtaisesti estimoitavan signaalin ja vastaanottimella vastaanotetun signaalin tulo, 10 integroidaan tulot ennalta määrätyn ajan T yli alustavan mittaustuloksen saamiseksi signaalilähteiden lähettämien signaalien mitatuille amplitudeille, ja muodostetaan alustavan mittaustuloksen ja 15 korjauskertoimen tulo, jossa korjauskerroin on eri signaalilähteistä lähetettyjen signaalien välistä korrelaatiota kuvaava suure, vastaanotetun signaalin amplitudin saamiseksi signaalilähdekohtaisesti.
4. Patenttivaatimuksen 2 mukainen menetelmä, 20 tunnettu siitä, että muodostetaan estimoitavan signaalin, korjauskertoimen ja vastaanotetun signaalin tulo, jossa korjauskerroin on eri signaalilähteistä lähetettyjen sig- ,* ; naalien välistä korrelaatiota kuvaava suure, ja t * · j 25 integroidaan tulot ennalta määrätyn ajan T yli mittaustuloksen saamiseksi signaalilähteiden lä- i * · hettämien signaalien mitatuille amplitudeille.
··*, 5. Patenttivaatimuksen 2 mukainen menetelmä, tunnettu siitä, että , 30 muodostetaan estimoitavan signaalin ja väli- • 4 I ;;; tun kertoimen signaalitulo, ···’ muodostetaan saadun signaalitulon ja vastaan- ;*·,· otetun signaalin tulo, integroidaan tulot ennalta määrätyn ajan T *, 35 yli alustavan mittaustuloksen saamiseksi signaaliläh- teiden lähettämien signaalien mitatuille amplitudeil-le, 113564 15 muodostetaan alustavan mittaustuloksen ja korjauskertoimen tulo, jossa korjauskerroin on eri signaalilähteistä lähetettyjen signaalien välistä korrelaatiota ja valitun kertoimen vaikutusta kuvaava 5 suure, vastaanotetun signaalin amplitudin saamiseksi signaalilähdekohtaisesti.
6. Jonkin patenttivaatimuksista 3-5 mukainen menetelmä, tunnettu siitä, että painotetaan tulot ikkunafunktiolla w.
7. Patenttivaatimuksen 1 mukainen menetelmä, tunnettu siitä, että lähetetään signaalilähteistä sinimuotoista signaalia; ja että käytetään laskennassa estimoituna signaalina 15 lähes samanmuotoista signaalia kuin lähetetty signaa li.
8. Patenttivaatimuksen 7 mukainen menetelmä, tunnettu siitä, että käytetään laskennassa toista lähetetyn signaalin kanssa samalla taajuudella olevaa 20 signaalia, jolla on vaihe-ero estimoituun signaaliin nähden.
9. Patenttivaatimuksen 1 mukainen menetelmä, : tunnettu siitä, että lähetetään signaalit saman- ; aikaisesti kultakin signaalilähteeltä. • 25
10. Patenttivaatimuksen 1 mukainen menetelmä, tunnettu siitä, että iit_; vastaanotetaan vastaanottimella hyötysignaa- ... lia, ja suodatetaan signaalilähteellä lähetetyt sig-30 naalit hyötysignaalista.
11. Patenttivaatimuksen 1 mukainen menetelmä, tunnettu siitä, että estimoidaan signaaleja, jot- : ka vastaavat ennalta määrätyllä tavalla liikkuvaan kappaleeseen kiinnitettyjen signaalilähteiden signaa-'!* 35 leita kappaleen liikkeen estimoimiseksi.
12. Jonkin patenttivaatimuksista 1-6 mukai- ^ · nen menetelmä, tunnettu siitä, että toistetaan 113564 16 kappaleen paikan ja/tai asennon määritys kappaleen suhteellisen paikan määrittämiseksi toistamalla ajallisesti päällekkäisiä mittausjaksoja.
13. Patenttivaatimuksen 6 mukainen menetelmä, 5 tunnettu siitä, että käytetään estimoituna signaalina tunnettujen häiriölähteiden signaalimuotoja.
13 113564
14. Patenttivaatimuksen 1 mukainen menetelmä, tunnettu siitä, että muodostetaan takaisinkytkentä saaduista amp-10 litudeista signaalilähteisiin, ja ohjataan signaalilähteiden lähetystehoa ta-kaisinkytkennällä.
15. Jonkin patenttivaatimuksista 1-14 mukainen menetelmä, tunnettu siitä, että 15 vähennetään mitatuista signaaleista lasketut signaalit, ja tarkennetaan mittaustulosta jäljelle jäävän signaalin avulla.
16. Jonkin patenttivaatimuksista 1-15 mu-20 kainen menetelmä, tunnettu siitä, että estimoidaan yksi tai useampi signaali, joka poikkeaa signaalilähteiden tai tunnettujen häiriöläh-teiden signaaleista, ja :,· · tarkennetaan paikannustulosta saadun mittaus- i 25 tuloksen avulla. » · • · » * # * · • * t ! * t * · * » » * * » * » i7 1 13564
FI20010558A 2001-03-19 2001-03-19 Kappaleiden paikannus FI113564B (fi)

Priority Applications (6)

Application Number Priority Date Filing Date Title
FI20010558A FI113564B (fi) 2001-03-19 2001-03-19 Kappaleiden paikannus
JP2002575573A JP4203592B2 (ja) 2001-03-19 2002-03-19 対象物の位置の決定
EP02708381.5A EP1370826B1 (en) 2001-03-19 2002-03-19 Determining a position of objects
CA2440988A CA2440988C (en) 2001-03-19 2002-03-19 Determining a position of objects
PCT/FI2002/000225 WO2002077565A1 (en) 2001-03-19 2002-03-19 Determining a position of objects
US10/665,543 US6876196B1 (en) 2001-03-19 2003-09-22 Determining a position of objects in a predetermined coordinate system

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FI20010558A FI113564B (fi) 2001-03-19 2001-03-19 Kappaleiden paikannus
FI20010558 2001-03-19

Publications (3)

Publication Number Publication Date
FI20010558A0 FI20010558A0 (fi) 2001-03-19
FI20010558A FI20010558A (fi) 2002-09-20
FI113564B true FI113564B (fi) 2004-05-14

Family

ID=8560785

Family Applications (1)

Application Number Title Priority Date Filing Date
FI20010558A FI113564B (fi) 2001-03-19 2001-03-19 Kappaleiden paikannus

Country Status (6)

Country Link
US (1) US6876196B1 (fi)
EP (1) EP1370826B1 (fi)
JP (1) JP4203592B2 (fi)
CA (1) CA2440988C (fi)
FI (1) FI113564B (fi)
WO (1) WO2002077565A1 (fi)

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020115925A1 (en) * 1996-06-25 2002-08-22 Avrin William F. Ferromagnetic foreign body detection utilizing eye movement
US6702789B1 (en) 1997-03-11 2004-03-09 Alcove Medical, Inc. Catheter having insertion control mechanism and anti-bunching mechanism
US7976518B2 (en) 2005-01-13 2011-07-12 Corpak Medsystems, Inc. Tubing assembly and signal generator placement control device and method for use with catheter guidance systems
US20110098559A1 (en) * 2005-01-14 2011-04-28 William John Besz Guiding insert assembly for a catheter used with a catheter position guidance system
JP4868382B2 (ja) * 2005-05-17 2012-02-01 公立大学法人広島市立大学 磁気刺激における刺激部位の特定あるいはターゲッティングを行うための装置
US8197494B2 (en) * 2006-09-08 2012-06-12 Corpak Medsystems, Inc. Medical device position guidance system with wireless connectivity between a noninvasive device and an invasive device
WO2008059636A1 (fr) * 2006-11-13 2008-05-22 Olympus Medical Systems Corp. Système d'analyse de forme lors de l'insertion d'un endoscope et système d'observation d'un corps vivant
ATE545039T1 (de) 2008-09-05 2012-02-15 Fiat Ricerche Bestimmung der relativen position von zwei relativ beweglichen elementen
US20110137102A1 (en) * 2009-06-04 2011-06-09 Mayo Foundation For Medical Education And Research Stereotactic intracranial target localization guidance systems and methods
DE102009049978A1 (de) * 2009-10-20 2011-05-05 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur relativen Lageerkennung von Funkortungsstationen
US9028441B2 (en) 2011-09-08 2015-05-12 Corpak Medsystems, Inc. Apparatus and method used with guidance system for feeding and suctioning
CN104169738B (zh) * 2012-04-26 2016-11-09 英特尔公司 确定相对定位信息
WO2019060298A1 (en) 2017-09-19 2019-03-28 Neuroenhancement Lab, LLC METHOD AND APPARATUS FOR NEURO-ACTIVATION
US11717686B2 (en) 2017-12-04 2023-08-08 Neuroenhancement Lab, LLC Method and apparatus for neuroenhancement to facilitate learning and performance
WO2019133997A1 (en) 2017-12-31 2019-07-04 Neuroenhancement Lab, LLC System and method for neuroenhancement to enhance emotional response
US11364361B2 (en) 2018-04-20 2022-06-21 Neuroenhancement Lab, LLC System and method for inducing sleep by transplanting mental states
EP3849410A4 (en) 2018-09-14 2022-11-02 Neuroenhancement Lab, LLC SLEEP ENHANCEMENT SYSTEM AND METHOD
US11786694B2 (en) 2019-05-24 2023-10-17 NeuroLight, Inc. Device, method, and app for facilitating sleep
CN113552536B (zh) * 2021-07-30 2022-08-09 重庆大学 一种含圆孔结构的声发射/微震事件定位方法、系统、终端及可读储存介质

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4407785A1 (de) 1994-03-09 1995-09-14 Philips Patentverwaltung Anordnung zur Bestimmung der räumlichen Position eines gegenüber einem Bezugselement verschiebbaren Abtastelements
US5729129A (en) 1995-06-07 1998-03-17 Biosense, Inc. Magnetic location system with feedback adjustment of magnetic field generator
IL122578A (en) 1997-12-12 2000-08-13 Super Dimension Ltd Wireless six-degree-of-freedom locator
WO2000017603A1 (en) * 1998-09-23 2000-03-30 The Johns Hopkins University Magnetic sensor system for fast-response, high resolution, high accuracy, three-dimensional position measurements
AU4579600A (en) * 1999-04-23 2000-11-21 Sentec Limited Catheter Device

Also Published As

Publication number Publication date
JP2004534209A (ja) 2004-11-11
CA2440988A1 (en) 2002-10-03
EP1370826A1 (en) 2003-12-17
US6876196B1 (en) 2005-04-05
CA2440988C (en) 2010-05-18
FI20010558A (fi) 2002-09-20
EP1370826B1 (en) 2013-04-24
FI20010558A0 (fi) 2001-03-19
WO2002077565A1 (en) 2002-10-03
JP4203592B2 (ja) 2009-01-07

Similar Documents

Publication Publication Date Title
FI113564B (fi) Kappaleiden paikannus
US7373271B1 (en) System and method for measuring position and orientation using distortion-compensated magnetic fields
US8121812B2 (en) AC magnetic tracking with phase disambiguation
Queiros et al. Cross-correlation and sine-fitting techniques for high-resolution ultrasonic ranging
US8611991B2 (en) System and method for conducting multiplexed electrical impedance tomography
US20180313908A1 (en) System for continuously calibrating a magnetic imaging array
AU776644B2 (en) Metal immune system for tracking spatial coordinates of an object in the presence of a perturbed energy field
WO2001033162A1 (en) Eddy current error-reduced ac magnetic position measurement system
JPH03180689A (ja) ボーリング機械の位置を特定する装置
CA2388328A1 (en) Method and apparatus for electromagnetic position and orientation tracking with distortion compensation
JP2008508070A (ja) 位相マッピングと、位相基準として用いる基準媒体が関係するmri温度測定
CN105548917B (zh) 一种非屏蔽环境下磁传感器阵列的校准方法
DE60130968D1 (de) Apparat und verfahren zur oberflächenkonturmessung
JPH10179551A (ja) 少なくとも1つの局所アンテナの位置を検出する方法
WO2018137977A1 (en) A method and apparatus for measuring a physiological characteristic of a subject
JPH04215084A (ja) 電界又は磁界の界パターンを測定する方法
US3249911A (en) Method and device for determining the position of sound sources in water
US5257625A (en) Method of noninvasive motion analysis by using forced closure of phase contrast MRI maps of velocity
KR100598146B1 (ko) 체임피던스 측정용 벨트 전극 장치
Chen et al. Multiple-target vital signs sensing using 77GHz FMCW radar
ATE494832T1 (de) Ekg-system und verfahren zur grossflächigen messung von ekg-signalen
SU1840066A1 (ru) Способ определения дистанции до источника излучения в режиме шумопеленгования
SU987671A1 (ru) Способ измерени интенсивности звука локального источника на дискретной частоте
WO2018138171A1 (en) A method and apparatus for measuring a physiological characteristic of a subject
JPH0947439A (ja) Mri装置による流速イメージング方法

Legal Events

Date Code Title Description
PC Transfer of assignment of patent

Owner name: ELEKTA AB (PUBL).

Free format text: ELEKTA AB (PUBL).

PC Transfer of assignment of patent

Owner name: MEGIN OY

MA Patent expired