FI80527B - Foerfarande foer optisk faststaellning av en ytas hastighet. - Google Patents
Foerfarande foer optisk faststaellning av en ytas hastighet. Download PDFInfo
- Publication number
- FI80527B FI80527B FI885758A FI885758A FI80527B FI 80527 B FI80527 B FI 80527B FI 885758 A FI885758 A FI 885758A FI 885758 A FI885758 A FI 885758A FI 80527 B FI80527 B FI 80527B
- Authority
- FI
- Finland
- Prior art keywords
- image
- determining
- value
- calculated
- pixel
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
- G01P3/64—Devices characterised by the determination of the time taken to traverse a fixed distance
- G01P3/68—Devices characterised by the determination of the time taken to traverse a fixed distance using optical means, i.e. using infrared, visible, or ultraviolet light
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Image Analysis (AREA)
Description
1 80527
Menetelmä pinnan nopeuden optiseksi määrittämiseksi Tämän keksinnön kohteena on menetelmä pinnan nopeuden optiseksi määrittämiseksi liikkeen suuntaiselta kuva-5 detektoririviltä saatujen, ajallisesti peräkkäisten kuva-rivien avulla.
Aikaisemmin tunnetussa menetelmässä on määritetty kahden ajallisesti peräkkäisen kuvarivin välinen korrelaatio, josta voidaan päätellä kuinka pitkän matkaa pinta on 10 siirtynyt näiden kuvarivien kuvaushetkien välillä. Näin on pystytty määrittämään pinnan nopeus. Tällainen korre-lointimenetelmä edellyttää suhteellisen selvien muusta ympäristöstään erottuvien pisteiden olemassaoloa kuvattavalla pinnalla, jotta nopeuden määritystä voitaisiin pitää 15 luotettavana. Sellaisten pintojen, joista saatavien kuva-rivien kuvapisteiden intensiteetit vaihtelevat näennäisen satunnaisesti, nopeutta ei ole juurikaan voitu tällaisella korrelointimenetelmällä määrittää.
Esillä olevan keksinnön tavoitteena onkin tuoda 20 esiin menetelmä lähes minkä tahansa ainakin jossain määrin kuviollisen pinnan nopeuden määrittämiseksi. Tällaisena kuviollisena pintana voidaan mainita esimerkiksi puupinta, metallipinta, tienpinta, kankaanpinta tai vastaavat pinnat, joista saatavien kuvarivien kuvapisteiden intensitee-25 tit vaihtelevat ainakin riittävää suurennusta käytettäessä. Keksinnön mukaisesti pinnan nopeuden määritys käsittää kuvadetektoririviltä saatujen analogiasignaalien digitoinnin; digitoitujen ajallisesti peräkkäisten kuvarivien tallentamisen, jolloin rinnakkain sovitettujen kuvarivien 30 voidaan ajatella muodostavan paikka-aikatasossa kaksiulotteisen kuvan; ja kunkin kuvarivin pisteen mainittuun paikka-aikatasoon piirtämän viivan kulmakertoimen tai vastaavan suureen arvon määrittämisen pinnan nopeuteen verrannollisen arvon saamiseksi. Edullisesti viivan kulmakertoi-35 men tai vastaavan suureen arvon määrittäminen käsittää 2 80527 sekä paikka- että aika-akselin suuntaisten gradienttien määrittämisen kullekin kuvarivin pisteelle ja kunkin kuva-pisteen paikka- ja aika-akselin suuntaisten gradienttien suhteen määrittämisen. Mittaustarkkuus voidaan moninker-5 taistaa määrittämällä kunkin kuvarivin pisteille laskettujen arvojen keskiarvo. Tuloksen luotettavuuden parantamiseksi voidaan haluttaessa ennen keskiarvon määritystä suorittaa keskiarvoistettavien kuvapisteiden valinta esimerkiksi vertaamalla kullekin kuvapisteelle laskettua, sen 10 gradienteista riippuvaa arvoa ennalta määrättyyn asetusar-voon. Tällä tavoin sellaiset kuvapisteet, jotka eivät erotu juurikaan ympäristöstään saadaan poistettua laskentatuloksesta.
Keksinnön mukaisen menetelmän perusideana on täten, 15 että keräämällä muistiin detektoririviltä tulevia digitoituja kuvarivejä, saadaan muodostettua kuva, jossa muusta ympäristöstä intensiteetiltään poikkeavat kuvapisteet muodostavat viivoja, joiden kulmakerroin on verrannollinen pinnan nopeuteen. Keksinnön mukaisesti viivojen kulmaker-20 roin mitataan digitaalisesti esimerkiksi konvoloimalla kuva paikka-akselin ja vastaavasti aika-akselin suuntaisten gradienttien saamiseksi, jolloin kulmakerroin saadaan näiden gradienttien suhteena. Tällä tavoin saadaan nopeus tai siihen verrannollinen arvo kullekin kuvan pisteelle. 25 Kuten edellä mainittiin mittaustarkkuutta voidaan tällöin moninkertaistaa keskiarvoistamalla kunkin kuvarivin no-peusarvot. Keskiarvoon mukaan otettavien arvojen valintaa voidaan ohjata ottamalla summattaviin arvoihin mukaan ainoastaan ne, joiden neliösumma ylittää ennalta määrätyn 30 arvon. Tällä tavoin keskiarvoistettaviin arvoihin otetaan mukaan vain ne kuvapisteet, jotka ovat mukana jossakin kuvan viivassa.
Esillä olevan keksinnön mukainen menetelmä on käyttökelpoinen suurella nopeusalueella ja sen mittausresoluu-35 tio on riittävä menetelmän soveltamiseen myös kappaleen 3 80527 pituuden mittaukseen. Menetelmä ottaa huomioon myös mitattavan pinnan suunnan tai sen liikkeen mahdollisen pysähtymisen. Menetelmä soveltuu toteutettavaksi rinnakkaisesta kuvan käsittelyyn kehitettyjen laitteistoratkaisujen avul-5 la suurella prosessointinopeudella.
Seuraavassa keksinnön mukaista menetelmää kuvataan yksityiskohtaisemmin viitaten oheiseen piirustukseen, jossa kuvio 1 esittää kaaviollisesti kuvattavaa pintaa 10 ja sitä kuvaavaa fotodetektoririviä, kuvio 2 esittää keksinnön mukaisesti muodostetun kaksiulotteisen kuvan ja kuvio 3 esittää keksinnön mukaisen menetelmän vuo-kaavion.
15 Kuviossa 1 on esitetty kaaviollisesti liikkuva pin ta 1, jota kuvataan kuvadetektoririvillä 2. Tämä kuvade-tektoririvi voi olla käytännössä muodostettu esimerkiksi videokameran tai rivikameran avulla tai miksei myös erillisten kuvadetektorien avulla. Tämä kuvadetektoririvi 2 20 kuvaa liikkuvalta pinnalta 1 kuvaelementtirivin 3. Pinnan 1 liikkeen suuntaa on kuviossa 1 kuvattu edestakaisella nuolella V.
Kuviossa 2 on esitetty kuvion 1 mukaisen järjestelyn avulla muodostettavissa oleva kaksiulotteinen paikka-25 aikatasossa oleva kuva. Tätä kuvaa muodostettaessa kuvion 1 pinta 1 on liikkunut edestakaisin kuvion vaakatasossa. Kuvio 2 on muodostettu sovittamalla detektoririviltä 2 saatavat, ajallisesti peräkkäiset kuvarivit päällekkäin siten, että yläreunassa on hetkellä t saatu kuvarivi, sen 30 alapuolella on hivenen myöhemmin hetkellä t+ t saatu kuva-rivi jne. Kuten kuviosta 2 voidaan havaita muodostavat tällöin muusta ympäristöstään intensiteetiltään eroavat kuvapisteet kuvaan viivoja, kun ne etenevät aikaa myöten liikkeen suuntaan (kuvion 2 yläosassa vasemmalle). Keksin-35 nön mukaisesti pinnan 1 nopeus määritetään määrittämällä 4 80527 näiden viivojen kulmakerroin tai vastaava suure, josta voidaan helposti laskea myös kuvapisteen nopeus.
Kuviossa 3 on esitetty eräs esimerkinomainen vuo-kaavio keksinnön mukaisen menetelmän toteuttamiseksi. Ku-5 vion 3 ensimmäisessä vaiheessa 4 kuvataan kuvadetektoriri-villä 2 yksi kuvarivi, joka seuraavassa vaiheessa 5 digitoidaan, jolloin kuvarivin kullekin kuvaelementille annetaan intensiteettiä vastaava digitaalinen arvo. Seuraavis-sa vaiheissa 6 ja 7 nämä digitoidut ajallisesti peräkkäi-10 set kuvarivit tallennetaan siirtorekisterien avulla siten, että kulloinkin on käytettävissä kaksi halutun lukumäärän kuvarivejä käsittävää kokonaisuutta, joiden voidaan ajatella kummankin muodostavan saman osan kuviossa 2 esitetystä kuvasta. Tästä kuvasta sitten vaiheessa 6 määrite-15 tään aika-akselin, lyhyesti y-akselin suuntainen gradient-ti kullekin kuvapisteelle. Vastaavasti vaiheessa 7 muodostetaan paikka-akselin, lyhyesti x-akselin, suuntainen gra-dientti kullekin kuvapisteelle. Nämä gradientit voidaan määrittää esimerkiksi niin sanottujen konvoluutiomaskien 20 avulla. Tällaista kuvapisteiden gradienttien määrittämistä konvoluutiomaskien avulla on kuvattu esimerkiksi teoksessa "Handbook of Pattern Recognition and Image Processing", Academic Press, Inc., 1986, toimittajina Tzay Y. Young ja King-Sun Fu sivuilla 192-194, 216 ja 217.
25 Vaiheilta 6 ja 7 saatujen kunkin kuvapisteen aika- gradientin dy ja paikkagradientin dx suhteena dy/dx määritetään nyt vaiheessa 8 kuvapisteen kulmakerroin. Täten on onnistuttu määrittämään jo joihinkin sovellutuksiin riittävä tieto pinnan nopeudesta. Suurempaan tarkkuuteen kui-30 tenkin päästään, jos lasketaan saman kuvarivin kuvapisteiden kulmakertoimien keskiarvo. Nämä toimenpiteet on toteutettu summanlaskennan osalta vaiheelta 9 ja keskiarvon osalta vaiheessa 10. Jotta kuitenkin summattavien kulmakertoimien arvoista saataisiin poistettua sellaiset, jotka 35 eivät osu minkään kuviossa 2 esiintyvän viivan kohdalle 5 80527 ja joiden sisältämä informaatio on täten lähinnä satunnaista, on vaiheessa 11 määritetty gradienttien dy ja dx neliösumman neliöjuuri kullekin kuvapisteelle. Vaiheessa 12 tätä neliösumman neliöjuurta on sitten verrattu ennalta 5 määrättyyn asetusarvoon. Tämän vertailun tuloksen perusteella ohjataan vaiheessa 9 summattavien kulmakertoimien laskentaa. Tätä on kuviossa 3 kuvattu vaiheiden 12 ja 9 välisellä nuolella. Tällä tavoin vaiheessa 9 summataan vain niin sanottujen viivojen kohdalle osuvien kuvapistei-10 den kulmakertoimet. Vaiheelta 12 saadaan tieto hyväksyttyjen kuvapisteiden lukumäärästä n vaiheelle 10, joka samoin sai vaiheelta 9 tiedon hyväksyttyjen kulmakertoimien summasta. Tällä tavoin vaiheessa 10 voidaan jakaa valittujen kulmakertoimien summa Σ dy/dx tähän summaan sisälty-15 vien kulmakertoimien lukumäärällä n, jolloin saadaan rivin kulmakertoimien keskiarvo, joka erittäin luotettavalla tavalla määrittää nopeuden kyseisen kuvarivin kuvaushetkellä.
Keksinnön mukaisen menetelmän käytännön toteutuk-20 seen esimerkiksi kuvion 3 mukaisten menetelmävaiheiden avulla voidaan käyttää VLSI-tason digitaalisia piirejä. Täten esimerkiksi analogia/digitaalimuunnosvaiheessa 5 voidaan toteuttaa esimerkiksi piirillä TDC1048 valmistaja TRW Inc.. Vaiheiden 6 ja 7 siirtorekisteritoiminnot ja 25 gradienttien laskentatoiminnot voidaan toteuttaa esimerkiksi piirin PDSP16488 avulla, jonka valmistaja on Plessey Semiconductors. Vaiheiden 8 kulmakertoimen laskennan ja gradienttien 4 summan neliöjuuren laskenta voidaan puolestaan toteuttaa pythagorasprosessorilla tyypiltään 30 PDSP16330 valmistajana Plessey Semiconductors, jolloin toisena ulostulona kulmakertoimen dy/dx sijasta onkin vastaava kulma eli atan (dy/dx). Vaiheen 9 mukainen kulma-arvojen summanlaskenta voidaan toteuttaa piirillä PDSP1601 valmistajana Plessey Semiconductors. Vaiheen 12 vertailu 35 ja hyväksytyksi tulleiden arvojen lukumäärän laskenta voi- 6 80527 daan toteuttaa varsin tavanomaisten MSI-tason komponenteista muodostettujen komparaattorin ja laskurin avulla. Vaiheen 10 lopullisen keskiarvon laskenta on yksinkertaisimmin toteutettavissa tähän tehtävään ohjelmoidun mikro-5 tietokoneen avulla, joka tällöin siis lukee kuvarivin kulmakertoimien summan sekä summaan laskettujen kulma-arvojen lukumäärän ja laskee näiden osamäärän eli rivin kulmakertoimien keskiarvon. Tämä keskiarvohan puolestaan on suoraan verrannollinen tarkasteltavan pinnan nopeuteen.
10 Yllä keksinnön mukaista menetelmää on kuvattu vain yhden esimerkinomaisen suoritusmuodon avulla eli kunkin kuvarivin pisteen paikka-aikatason piirtämän viivan kulmakertoimen tai vastaavan suureen arvon määrittämiseksi on annettu vain yksi esimerkinomainen tapa. Tämä esimerk-15 kina annettu tapa on yksinkertaisuutensa johdosta edullinen, mutta on ymmärrettävää, että myös muiden matemaattisten ja kuvan käsittelyteknisten menetelmien avulla voidaan määrittää muodostetusta kaksiulotteisesta kuvasta siinä esiintyvien viivojen kulmakertoimet tai niihin verrannol-20 liset suureet.
Claims (5)
1. Menetelmä pinnan (1) nopeuden optiseksi määrittämiseksi liikkeen suuntaiselta kuvadetektoririviltä (2) 5 saatujen, ajallisesti peräkkäisten kuvarivien (3) avulla, tunnettu siitä, että se käsittää - kuvadetektoririviltä (2) saatujen analogiasignaa-lien digitoinnin (5), - digitoitujen ajallisesti peräkkäisten kuvarivien 10 tallentamisen, jolloin rinnakkain sovitettujen kuvarivien voidaan ajatella muodostavan paikka-aikatasossa kaksiulotteisen kuvan, ja - kunkin kuvarivin pisteen mainittuun paikka-aikatasoon piirtämän viivan kulmakertoimen tai vastaavan suu- 15 reen arvon määrittämisen (8) pinnan (1) nopeuteen verrannollisen arvon saamiseksi.
2. Patenttivaatimuksen 1 mukainen menetelmä, tunnettu siitä, että viivan kulmakertoimen tai vastaavan suureen arvon määrittäminen käsittää 20. sekä paikka- että aika-akselin suuntaisten gra- dienttien määrittämisen (6, 7) kullekin kuvarivin pisteelle ja - kunkin kuvapisteen paikka- ja aika-akselin suuntaisten gradienttien suhteen määrittämisen (8).
3. Patenttivaatimuksen 1 tai 2 mukainen menetelmä, tunnettu siitä, että se lisäksi käsittää saman kuvarivin pisteille laskettujen arvojen keskiarvon määrittämisen (10).
4. Patenttivaatimuksen 3 mukainen menetelmä, 30 tunnettu siitä, että se ennen keskiarvon määritystä (10) käsittää keskiarvoistettavien kuvapisteiden valinnan (12).
5. Patenttivaatimuksen 4 mukainen menetelmä, tunnettu siitä, että keskiarvoistettavien kuvapis- 35 teiden valinta suoritetaan vertaamalla (12) kullekin kuva-pisteelle laskettua (12), sen gradienteista riippuvaa arvoa ennalta määrättyyn asetusarvoon. 8 80527
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FI885758A FI80527C (fi) | 1988-12-13 | 1988-12-13 | Foerfarande foer optisk faststaellning av en ytas hastighet. |
| DE90900875T DE68913993D1 (de) | 1988-12-13 | 1989-12-11 | Verfahren zur optischen bestimmung der geschwindigkeit einer fläche. |
| PCT/FI1989/000231 WO1990007123A1 (en) | 1988-12-13 | 1989-12-11 | A method of optically determining the speed of a surface |
| AT90900875T ATE103074T1 (de) | 1988-12-13 | 1990-07-06 | Verfahren zur optischen bestimmung der geschwindigkeit einer flaeche. |
| EP19900900875 EP0448615B1 (en) | 1988-12-13 | 1990-07-06 | A method of optically determining the speed of a surface |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FI885758 | 1988-12-13 | ||
| FI885758A FI80527C (fi) | 1988-12-13 | 1988-12-13 | Foerfarande foer optisk faststaellning av en ytas hastighet. |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| FI885758A0 FI885758A0 (fi) | 1988-12-13 |
| FI80527B true FI80527B (fi) | 1990-02-28 |
| FI80527C FI80527C (fi) | 1990-06-11 |
Family
ID=8527559
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FI885758A FI80527C (fi) | 1988-12-13 | 1988-12-13 | Foerfarande foer optisk faststaellning av en ytas hastighet. |
Country Status (4)
| Country | Link |
|---|---|
| EP (1) | EP0448615B1 (fi) |
| DE (1) | DE68913993D1 (fi) |
| FI (1) | FI80527C (fi) |
| WO (1) | WO1990007123A1 (fi) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007085704A1 (en) | 2006-01-30 | 2007-08-02 | Visicamet Oy | Method and measuring device for measuring translation of surface |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102006030130B3 (de) * | 2006-06-28 | 2007-09-27 | Scansonic Gmbh | Verfahren und Vorrichtung zum Bearbeiten eines Werkstücks mittels eines Energiestrahls, insbesondere Laserstrahls |
| DE102010022592B4 (de) * | 2010-05-31 | 2013-08-01 | Siemens Aktiengesellschaft | Verfahren und Einrichtung zum Ermitteln der Geschwindigkeit eines Fahrzeugs |
| ITVR20100143A1 (it) * | 2010-07-15 | 2012-01-16 | Microtec Srl | Metodo ed apparecchiatura per determinare lunghezza e/o velocità di pezzi di legname che avanzano longitudinalmente su un trasportatore |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3890462A (en) * | 1974-04-17 | 1975-06-17 | Bell Telephone Labor Inc | Speed and direction indicator for video systems |
| US4096525A (en) * | 1976-03-08 | 1978-06-20 | William James Lathan | Video scanning change discriminator |
| US4162509A (en) * | 1978-06-21 | 1979-07-24 | The United States Of America As Represented By The Secretary Of The Army | Non-contact velocimeter using arrays |
| DE3110828C2 (de) * | 1981-03-19 | 1984-05-30 | Endress U. Hauser Gmbh U. Co, 7867 Maulburg | Verfahren und Anordnung zur Erzeugung von miteinander zu korrelierenden Signalen, insbesondere zur berührungslosen Geschwindigkeitsmessung |
| SE438091B (sv) * | 1983-06-08 | 1985-04-01 | Gote Palsgard | Anordning for koordinatregistrering |
| DE3539793A1 (de) * | 1985-11-09 | 1987-05-14 | Bbc Brown Boveri & Cie | Verfahren zur beruehrungslosen geschwindigkeits- und laengenmessung |
-
1988
- 1988-12-13 FI FI885758A patent/FI80527C/fi not_active IP Right Cessation
-
1989
- 1989-12-11 DE DE90900875T patent/DE68913993D1/de not_active Expired - Lifetime
- 1989-12-11 WO PCT/FI1989/000231 patent/WO1990007123A1/en not_active Ceased
-
1990
- 1990-07-06 EP EP19900900875 patent/EP0448615B1/en not_active Expired - Lifetime
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007085704A1 (en) | 2006-01-30 | 2007-08-02 | Visicamet Oy | Method and measuring device for measuring translation of surface |
| EA011851B1 (ru) * | 2006-01-30 | 2009-06-30 | Визикамет Оу | Способ и устройство для измерения перемещения поверхности |
| US8107089B2 (en) | 2006-01-30 | 2012-01-31 | Visicamet Oy | Method and measuring device for measuring translation of surface |
| AU2007209273B2 (en) * | 2006-01-30 | 2012-03-22 | Visicamet Oy | Method and measuring device for measuring translation of surface |
Also Published As
| Publication number | Publication date |
|---|---|
| DE68913993D1 (de) | 1994-04-21 |
| FI80527C (fi) | 1990-06-11 |
| WO1990007123A1 (en) | 1990-06-28 |
| EP0448615B1 (en) | 1994-03-16 |
| EP0448615A1 (en) | 1991-10-02 |
| FI885758A0 (fi) | 1988-12-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7750970B2 (en) | Distance-measuring device installed in camera | |
| US4842411A (en) | Method of automatically measuring the shape of a continuous surface | |
| US7084905B1 (en) | Method and apparatus for obtaining high dynamic range images | |
| JP2696044B2 (ja) | 合焦検出方法、これを用いた非接触変位測定方法及び装置 | |
| Seitz | Optical superresolution using solid-state cameras and digita; signal processing | |
| CN113191954B (zh) | 一种基于双目相机的全景图像拼接方法 | |
| US6876458B2 (en) | Method and device for determining the absolute coordinates of an object | |
| CN109883391B (zh) | 基于微透镜阵列数字成像的单目测距方法 | |
| Spirig et al. | The multitap lock-in CCD with offset subtraction | |
| US7315383B1 (en) | Scanning 3D measurement technique using structured lighting and high-speed CMOS imager | |
| JP3342588B2 (ja) | 不鮮明ビデオラインの補正方法およびカメライメージャ | |
| FI80527B (fi) | Foerfarande foer optisk faststaellning av en ytas hastighet. | |
| US5432331A (en) | Method and apparatus for detecting focus of moving images with tilted plane detector and time delay means | |
| EP0493105A2 (en) | Data processing method and apparatus | |
| JP2002286433A (ja) | 連続移動物体のリアルタイム形状計測方法及びシステム | |
| Leinss et al. | Measuring glacier velocity by autofocusing temporally multilooked SAR time series | |
| Konig et al. | Image-based motion compensation for structured light scanning of dynamic surfaces | |
| JPH0618779A (ja) | 焦点自動検出装置 | |
| KR20150111944A (ko) | 전기 회로의 3차원 이미지를 결정하기 위한 시스템 | |
| Hutber | Automatic inspection of 3D objects using stereo | |
| JP2000055651A (ja) | 測距演算装置およびそれを用いた測距演算方法 | |
| Pedersini et al. | Subpixel edge localization with statistical error compensation | |
| McIvor | An alternative interpretation of structured light system data | |
| Zhang et al. | A Fast Center of Mass Estimation Algorithm for Coordinates of IR Markers | |
| Edelman | Image Understanding and Computational Vision |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM | Patent lapsed | ||
| MM | Patent lapsed |
Owner name: VALTION TEKNILLINEN TUTKIMUSKESKUS |