KR100570202B1 - 도전율 및 전류밀도를 영상화하는 시스템 및 방법 - Google Patents
도전율 및 전류밀도를 영상화하는 시스템 및 방법 Download PDFInfo
- Publication number
- KR100570202B1 KR100570202B1 KR1020040001993A KR20040001993A KR100570202B1 KR 100570202 B1 KR100570202 B1 KR 100570202B1 KR 1020040001993 A KR1020040001993 A KR 1020040001993A KR 20040001993 A KR20040001993 A KR 20040001993A KR 100570202 B1 KR100570202 B1 KR 100570202B1
- Authority
- KR
- South Korea
- Prior art keywords
- conductivity
- magnetic flux
- flux density
- measurement object
- measured
- Prior art date
Links
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
- E01D—CONSTRUCTION OF BRIDGES, ELEVATED ROADWAYS OR VIADUCTS; ASSEMBLY OF BRIDGES
- E01D19/00—Structural or constructional details of bridges
-
- E—FIXED CONSTRUCTIONS
- E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
- E01D—CONSTRUCTION OF BRIDGES, ELEVATED ROADWAYS OR VIADUCTS; ASSEMBLY OF BRIDGES
- E01D19/00—Structural or constructional details of bridges
- E01D19/12—Grating or flooring for bridges; Fastening railway sleepers or tracks to bridges
- E01D19/125—Grating or flooring for bridges
Landscapes
- Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
본 발명은 측정대상의 도전율 및 전류밀도를 보다 정확하게 영상화할 수 있는 도전율 및 전류밀도를 영상화하는 시스템 및 방법에 관한 것으로, 측정대상 내에 서로 다른 방향의 전류들을 주입하는 전류주입수단들과, 상기 전류들에 의해 상기 측정대상 내에서 발생되는 자속밀도의 한 방향 성분을 측정하는 MRI 장치와, 상기 측정된 자속밀도의 한 방향 성분을 이용하여 상기 측정대상의 내부 전압을 계산하며, 상기 자속밀도 성분을 한번 미분하는 연산과정을 통해 상기 측정대상 내부의 도전율을 계산하는 연산부와, 상기 계산된 도전율을 이용하여 상기 측정대상의 내부를 영상화하는 디스플레이 수단을 포함한다.
도전율, 전류밀도, MRI, MREIT
Description
도 1은 본 발명에 따른 도전율 및 전류밀도 영상화 시스템을 나타낸 도면
도 2는 도 1의 전류주입수단의 일 실시예를 나타낸 도면
도 3a 및 3b는 도 1의 전류주입수단의 다른 실시예를 나타낸 도면
도 4는 측정대상 내부의 전류 흐름을 나타낸 도면
도면의 주요부분에 대한 부호의 설명
100 : 전류주입수단 200 : MRI 장치
300 : 연산부 400 : 디스플레이 수단
500 : 전류원 21,31 : 전극
23,33 : 도선 25,35 : 절연용기
37 : 전해물질
본 발명은 인체나 물체 등과 같은 측정대상 내의 도전율(또는 저항률)과 전류밀도를 영상화하는 방법 및 그를 위한 시스템에 관한 것이다.
일반적으로 인체 또는 물체의 내부 구조를 영상화하는데는 X-선, MRI, 초음파 등이 사용되고 있다. 그러나 이들 방법들은 인체 또는 물체의 전기적 특성을 영상화하지는 못한다.
이러한 문제를 해결하기 위한 노력으로, 1989년에는 MRI 기술을 활용하여 측정대상 내부의 전류밀도를 영상화하는 기술이 캐나다 토론토대학의 연구팀에 의해 최초로 제안되었고, 그 이후에도 관련된 연구가 활발히 진행되어 왔다. 이것이 CDI(Current Density Imaging) 기술이다.
그러나 이 CDI 기술은 3가지 성분(Bx,By,Bz)을 갖는 자속밀도 B를 얻기 위해서는 측정대상을 MRI 장비 내에서 회전시켜야 한다는 단점을 가진다. 이는 MRI 장비 내에 측정대상이 위치하고 있을 때, MRI 장비의 특성상 주자장(main magnetic field)의 방향과 일치되는 방향인 z-방향 성분의 자속밀도 즉 Bz 만이 측정 가능하기 때문이다.
즉 MRI 장비가 한 번에 측정 가능한 성분은 BZ 성분 하나이기 때문에, 기존의 CDI 기술은 세 방향성분의 자속밀도 벡터 값을 모두 획득하기 위하여 MRI 장비 내에서 측정대상(인체 또는 물체)을 회전시켜야 한다는 심각한 문제점을 가지고 있다.
또한 인체 또는 물체의 전기적 특성을 영상화하기 위한 기존 기술로는, 1970년대 말부터 활발히 연구되어온 전기임피던스단층촬영(EIT : Electrical Impedance Tomography) 기술이 있다. 이 EIT 기술은 측정대상의 전기적 특성을 나타내는 저항률 (또는 도전율) 영상을 제공한다.
그 EIT 기술은 주로 인체를 측정대상으로 하여, 그 인체의 표면에 여러 개의 전극을 부착하여 인체 내부의 저항률을 영상화한다. 즉 인체의 표면에 부착된 여러 전극들을 통해 전류를 인가한 후에 다시 표면에 부착된 전극을 통해 전압을 측정하여 인체 내부의 저항률을 영상화하는 기술이다.
이렇게 인체 내부를 저항률에 따라 영상화할 수 있는 것은 혈액, 뼈, 허파, 심장, 뼈 등의 인체 내부의 생체조직들이 서로 다른 전기적 특성을 가지기 때문이다.
그러나 EIT 기술이 가지고 있는 근본적인 결점 때문에, EIT 기술에 의해 복원된 영상은 부정확하고 저해상도만을 지원하였다. 즉 EIT 기술에 의해 측정된 전류-전압 데이터는 인체 내부의 저항률 변화에 지극히 둔감하다는 근본적인 결점을 갖는다. 따라서 EIT 기술은 임상적인 이용이 어려운 실정이다.
EIT 기술에서 기존의 CDI 기술을 적용한다 하더라도, 측정대상(인체 또는 물체)을 회전시켜야 하기 때문에 그 기술적 문제점을 해결할 수 없다.
본 발명의 목적은 상기한 기존의 기술들의 단점을 보완하기 위해 안출한 것으로써, 측정대상을 회전시키기 않고도 간단한 연산과정을 통해 측정대상의 전기적 특성들(도전율 및 전류밀도)을 정확한 고해상도로 영상화할 수 있는 도전율 및 전류밀도의 영상화 방법 및 그를 위한 시스템을 제공하는 데 있다.
상기와 같은 목적을 달성하기 위한 본 발명의 도전율 및 전류밀도의 영상화 시스템은, 측정대상 내에 서로 다른 방향의 전류들을 주입하는 전류주입수단들과, 상기 전류들에 의해 상기 측정대상 내에서 발생되는 자속밀도의 한 방향 성분을 측정하는 MRI 장치와, 상기 측정된 자속밀도의 한 방향 성분을 이용하여 상기 측정대상의 내부 전압을 계산하며, 상기 자속밀도 성분을 한번 미분하는 연산과정을 통해 상기 측정대상 내부의 도전율을 계산하는 연산부와, 상기 계산된 도전율을 이용하여 상기 측정대상의 내부를 영상화하는 디스플레이 수단을 포함한다.
본 발명의 도전율 및 전류밀도 영상화 방법은, 측정대상 표면에 부착된 전류주입수단들을 통해 상기 측정대상 내부에 서로 다른 방향의 전류들을 공급하는 단계와; 상기 전류들에 의해 상기 측정대상 내에서 발생되는 자속밀도의 한 방향 성분을 측정하는 단계와; 상기 측정된 자속밀도의 한 방향 성분을 한번 미분하는 연산과정을 통해 상기 측정대상 내부의 도전율을 계산하는 단계와; 상기 계산된 도전율을 이용하여 상기 측정대상 내부를 영상화하는 단계를 포함한다.
상기 도전율을 계산하는 단계는, 임의의 도전율에 상응하는 상기 측정대상의 내부전압을 계산하는 단계와; 상기 자속밀도의 성분과 상기 계산된 내부전압을 이용하여 도전율을 계산하는 단계를 포함한다.
상기 도전율을 계산하는 단계에서, 아래의 수학식들 중 하나에 상기 자속밀 도 성분 B
j
z과 상기 내부 전압 V
j 을 대입하여 상기 측정대상 내부의 도전율을 계산한다.
이하, 첨부된 도면을 참조하여 본 발명의 바람직한 일실시예를 상세히 설명하면 다음과 같다.
도 1은 본 발명에 따른 영상화 시스템의 구성을 나타낸 도면이다. 도 1에 도시한 바와 같이, 본 발명에 따른 영상화 시스템은, 인체나 물체 등과 같은 측정대상(S)의 내부에 여러 방향(j=1,2,...N)의 전류들을 차례로 주입하는 다수의 전류주입수단(100)들과, 상기 측정대상(S) 내부에 주입된 전류들에 따른 상기 측정대상(S) 내의 자속밀도들을 측정하는 MRI 장치(200)와, 상기 주입된 전류들에 따른 상기 측정대상(S) 표면의 전압과 상기 측정된 자속밀도를 이용하여 상기 측정 대상(S) 내부의 도전율 분포(conductivity distribution)와 전류밀도 분포(current density distribution)를 계산하는 연산부(300)와, 상기 측정대상(S)의 내부를 도전율(또는 저항률) 영상 및 전류밀도 영상으로 나타내는 디스플레이 수단(400)과, 상기 전류주입수단(100)에 전류를 공급하는 전류원(500)을 포함한다.
상기 전류주입수단(100)은, 구리와 같은 비자성의 금속이나 도전율이 매우 큰 재료로 이루어진 전극(21)과, 상기 전극(21)에 전류를 공급하는 도선(wire)(23)를 포함한다. 상기 측정대상(S)의 표면에 상기 전극(21)을 직접적으로 부착시키는 것이 가능하지만, 상기 측정대상(S)의 표면(boundary)에서 자속밀도의 왜곡을 방지하기 위해 상기 전극(21)을 상기 측정대상(S)의 표면으로부터 이격(separate)시키는 것이 바람직하다. 따라서 본 발명의 전류주입수단(100)의 일실시예는, 도 2에 도시한 바와 같이 상기 전극(21)을 상기 측정대상(S)으로부터 이격시키기 위해 상기 전극(21)과 부착되는 절연용기(25)를 더 포함한다. 상기 절연용기(25)의 한쪽 면은 상기 전극(21)과 부착되고, 다른 면은 개방된 상태로 상기 측정대상(S)과 접촉한다. 그리고 상기 절연용기(25)의 내부에는 전해질 젤이나 전해질 용액을 흡수한 스폰지가 삽입되어져 있다. 전해질 젤이나 전해질 용액을 흡수한 스폰지는 전극(21)으로부터 공급되는 전류가 절연용기(25) 내에 균일하게 흐르게 한다. 이처럼 절연용기(25)를 이용하여 상기 전극(21)을 측정대상(S)의 표면으로부터 일정 간격만큼 이격시킴과 동시에 전해질 젤이나 전해질 용액을 흡수한 스폰지를 이용하여 상기 측정대상(S)에 전류를 균일하게 공급함으로써, 전극(21)과 접촉하는 측정대상(S)의 표면 근처에서 MRI 영상의 왜곡을 줄일 수 있다.
본 발명의 전류주입수단(100)의 다른 실시예는, 도 3a 및 3b에 도시한 바와 같이 전극(31)과, 상기 전극(31)에 연결된 도선(33)와, 상기 전극(31)을 부착한 반구형(半球形)의 절연용기(35)를 포함한다. 상기 절연용기(35)의 한쪽 면은 상기 전극(31)과 부착되고, 다른 면은 개방된 상태로 상기 측정대상(S)과 접촉한다. 그리고 상기 절연용기(35)의 내부에는 전해질 젤이나 전해질 용액을 흡수한 스폰지(37)가 삽입되어져 있다.
상기 전류주입수단(100)들은 측정대상(S)의 둘레에 부착되고, 한번에 두 개의 전류주입수단(100)에만 전류가 흐르게 된다. 따라서 두 개의 전류주입수단(100)에 의해 상기 측정대상(S)내에 전류가 공급된다. 전류주입수단(100)의 개수를 E라고 가정했을 때, 상기 측정대상(S) 내에 전류를 주입하기 위한 전류주입수단(100) 쌍에 대한 경우의 수 N은 E(E-1)/2 개이다. 이는 상기 측정대상(S)내에서 흐르는 전류 경로의 개수와 같다. 상기 전류주입수단(100)의 쌍 사이의 주입 전류를 I j (j=1,2,...N)라고 한다. 상기 전류 I j는 상기 측정대상(S) 내에서 전류 밀도 를 발생시킨다.
상기 MRI 장치(200)는 상기 전류 I
j에 따른 상기 측정대상(S) 내의 자속 밀도 B
j를 측정한다. 특히, 상기 자속 밀도 B
j 중 상기 MRI 장치(200)의 주자장(main magnetic field)과 평행한 z방향에 해당하는 자속 밀도의 성분 B
j
z를 측정한다. 즉, 상기 MRI 장치(200)는 상기 측정대상(S) 내에 주입된 각 전류(I
1,I
2,
I
3,....I
N)에 의해 발생하는 자속 밀도 를 측정한다.
상기 연산부(300)는 전류주입수단(100) 쌍을 차례로 선택하며 상기 선택된 전류주입수단(100) 쌍으로 하여금 상기 측정대상(S) 내에 전류를 주입하도록 상기 전류주입수단(100)들을 제어한다. 또한 상기 연산부(300)는 상기 MRI 장치(200)에 의해 측정된 자속 밀도 B
j
z를 설정된 알고리즘에 적용하여 상기 측정대상(S) 내부의 도전율 과 전류밀도 J
j를 계산한다. 여기서, 상기 도전율 은 저항률 의 역수이므로 1/로 표현될 수 있다. 그리고, 상기 디스플레이 수단(400)은 상기 도전율 및 전류밀도 J
j에 따라 상기 측정대상(S)의 내부를 각각 영상화한다.
본 발명에 따른 영상화 방법과 그에 이용되는 알고리즘을 구체적으로 설명하면 다음과 같다.
먼저, 도 1에 도시한 바와 같이, 측정대상(S)의 둘레 즉, 신체의 둘레에 다수의 전류주입수단(100)이 부착된다. 도전율 를 구하기 위해서는 측정대상(S) 내부에 차례로 주입되는 전류들은 2가지 이상의 방향으로 주입되어야 하기 때문에 상기 전류주입수단(100)은 적어도 3개 이상 필요하며 적절한 간격으로 배치되어야 한다.
상기 측정대상(S)의 표면에 상기 전류주입수단(100)이 부착된 상태로 상기 측정대상(S)을 MRI 장치(200) 내에 위치시킨 후에, 상기 전류주입수단(100) 쌍을 하나씩 차례로 선택하고, 차례로 선택되는 상기 전류주입수단(100) 쌍을 통해 상기 측정대상(S) 내부에 전류 I
j (j=1,2,...N)를 차례로 공급한다. 이때, MRI 장치(200)는 상기 측정대상(S) 내에 주입된 각 전류들에 의해 각각 발생하는 상기 측정대상(S) 내의 자속 밀도와 상기 측정대상(S)의 표면 전압 을 측정한다. 상기 MRI 장치(200)는 z 방향에 해당하는 자속 밀도의 성분만을 측정한다. 여기서, 상기 표면 전압 은 상기 선택된 전류주입수단(100) 쌍을 통해 상기 측정대상(S)에 전류가 공급될 때, 선택되지 않은 전류주입수단(100)들과 측정대상(S) 사이의 경계면(boundary)에서의 전압이다. 즉, 상기 표면 전압 은 상기 선택되지 않은 전류주입수단(100)들에 의해 측정된다.
이어, 상기 연산부(300)는, 상기 전류 에 의해 발생하는 상기 측정대상(S)의 표면 전압 과, 상기 MRI 장치(200)에 의해 측정된 자속 밀도 를 설정된 알고리즘에 적용하여 상기 측정대상(S) 내부의 도전율와 전류밀도를 계산한다.
g j : 측정대상 표면에서 주입전류에 따른 전류밀도 성분
도전율 의 계산에 자속밀도 와 전류밀도 사이의 막스웰 관계식 을 이용한다. 상기 관계식의 양변에 벡터함수 를 내적한 후 양변을 측정영역 S에서 삼중적분하고, 오옴의 법칙 을 이용하면 다음의 수학식 1을 얻게 된다. 여기서, 는 자유공간에서의 투자율(magnetic permeability)이다.
여기에서 이다. 측정대상을 회전시키지 않고 측정 가능한 자속밀도 성분은 만이므로 수학식 1에서 와 를 제거하는 것이 필요하다. 이를 위해 서 다음의 수학식 2와 같은 함수의 집단을 설정한다.
상기 MRI 장치(200)에 의해 측정된 자속밀도 Bjz와 상기 측정대상(S) 내부의 전압 Vj (j=1,2,...N)를 상기 수학식 3에 대입하면 상기 도전율 를 구할 수 있다. 그러나 상기 측정대상(S) 내부의 전압 Vj 를 알지 못하면 상기 측정대상(S) 내의 실제 도전율 를 구할 수 없다. 상기 측정대상(S) 내부의 전압 Vj 는 측정이 불가능하므로 이를 구하기 위해 수학식 4을 이용한 반복알고리즘을 사용한다.
여기에서 은 새로 계산하는 도전율이고 는 이전단계의 도전율 을 이용하여 계산한 측정대상(S) 내부의 전압이다. 이때, 는 수학식 5를 이용하여 FEM(finite element method)와 같은 방법으로 계산할 수 있다.
여기서, 상기 전류밀도 성분 g
j 은 선택된 전류주입수단(100) 쌍과 접촉하지 않는 측정대상(S)의 표면에서 제로(0)값을 갖는다.
상기의 수학식들(수학식4, 수학식5)을 이용하여 상기 연산부(300)가 도전율 및 전압 를 계산하는 과정은 다음과 같다. 임의로 초기 도전율 (m=0)을 설정하고, 수학식 5에 상기 초기 도전율 과 상기 측정대상(S) 표면에서의 전류밀도 성분 g
j 을 대입하여 상기 측정대상(S) 내부 및 표면에서의 전압 과 을 계산한다. 여기서, 상기 초기 도전율 은 매질이 균질하다는 가정 하에 동일 값으로 설정되거나 또는 인체의 조직별로 알려져 있는 도전율의 값으로 설정 될 수 있다. 그리고 상기 전압들은 임의로 설정된 초기 도전율 을 이용하여 계산된 전압들이다.
상기 계산된 도전율 은 이전 단계의 도전율 에 비하여 실제의 도전율 값에 더 가까워진다. 왜냐하면, 상기 계산된 도전율 은 상기 측정대상(S) 내부의 실제 자속밀도 B
j
z가 고려된 값이기 때문이다.
본 발명에서는 한층 더 실제 도전율에 가까운 도전율을 얻기 위해 반복(iteration) 형태의 알고리즘을 이용한다. 즉, 이전 단계의 도전율 과 상기 개선된 도전율 의 차이를 적분값 로 계산하고, 상기 적분 값이 설정된 값 보다 작은지를 판단한다. 이때 상기 적분 값이 설정된 값 보다 크다면 작아질 때까지 m 값을 증가시켜 상술한 바와 같이 상기 수학식 4과 수학식 5를 반복 계산한다. 만약, 상기 이전 단계의 도전율 의 성분과 상기 계산된 도전율 의 성분간의 차가 설정된 값 보다 작으면 상기 도전율 의 성분들을 실제 도전율로 여기고, 상기 도전율 을 이용하여 상기 측정대상(S) 내부를 영상화한다.
그러나 상기 계산된 도전율 이 실제 자속밀도 B
j
z를 고려해서 얻어진 것이라 해도 실제 도전율과 다를 수 있다. 즉, 상기 계산된 도전율 은 상기 자속밀도 B
j
z에 따라 다른 값을 갖지만 실제 도전율에 비해 상수 배만큼 크거나 작을 수 있다. 실제 도전율에 대한 상기 계산된 도전율 의 상수 배 차이를 제거하기 위해 본 발명에서는 상기 측정된 표면 전압 와 상기 계산된 표면 전압 을 이용하여 상기 도전율 를 스케일(scale)하는 과정을 추가할 수 있다. 다시 말해서, 상기 측정된 표면 전압 와 상기 계산된 표면 전압 의 비율에 따라 상기 도전율 을 상수 배로 곱하거나 나눈다. 이로써 보다 더 정확한 도전율 값을 얻을 수 있다.
이후, 상기 도전율 을 상기 수학식 5에 대입하여 그에 대응하는 전압, 즉 상기 측정대상(S) 내부의 전압 V
j 을 계산하고, 상기 도전율 과 상기 전압 V
j 을 이용하여 전류밀도 J
j를 계산한다. 그리고, 상기 전류밀도
J
j를 이용하여 상기 측정대상(S) 내부를 영상화한다.
이상에서는 수학식 1-5를 이용하는 방법을 기술하였다. 이제 또 다른 수학식들을 이용하여 상기 자속밀도 B
j
z를 한번만 미분하는 방법을 기술한다.
수학식 7을 수학식 6에 적용하면 다음 수학식 8을 얻게 된다.
상기 수학식 8에서 기호 를 사용한 이유는 상기 두 함수 값이 와 측정데이터 값 에 의해 결정되기 때문이다. 상기 수학식 8과 오옴의 법칙 를 이용하면 다음 수학식 9를 얻을 수 있다.
상기 MRI 장치(200)에 의해 측정된 자속밀도 와 상기 측정대상(S) 내부의 전압 V
j (j=1,2,...N)를 상기 수학식 9에 대입하면 상기 도전율 를 구할 수 있다. 상기 수학식 9를 이용하면 상기 자속밀도 와 상기 내부 전압 V
j 를 한번만 미분하여 상기 도전율 를 구할 수 있다. 그러나 상기 측정대상(S) 내부의 전압 V
j 를 알지 못하면 상기 측정대상(S) 내의 실제 도전율을 구할 수 없다. 상기 측정대상(S) 내부의 전압 V
j 는 측정이 불가능하므로 이를 구하기 위해 수학식 10을 이용한 반복알고리즘을 사용한다.
여기에서 은 새로 계산하는 도전율이고 는 이전단계의 도전율 을 이용하여 계산한 측정대상(S) 내부의 전압이다. 이때, 는 수학식 11을 이용하여 FEM(finite element method)와 같은 방법으로 계산할 수 있다.
여기서, 상기 전류밀도 성분 g
j 은 선택된 전류주입수단(100) 쌍과 접촉하지 않는 측정대상(S)의 표면에서 제로(0)값을 갖는다.
상기의 수학식들(수학식10, 수학식11)을 이용하여 상기 연산부(300)가 도전율 및 전압 를 계산하는 과정은 다음과 같다. 임의로 초기 도전율 (m=0)을 설정하고, 수학식 11에 상기 초기 도전율 과 상기 측정대상(S) 표면에서의 전류밀도 성분 g
j 을 대입하여 상기 측정대상(S) 내부 및 표면에서의 전압
과 을 계산한다. 여기서, 상기 초기 도전율 은 매질이 균질하다는 가정 하에 동일 값으로 설정되거나 또는 인체의 조직별로 알려져 있는 도전율의 값으로 설정될 수 있다. 그리고 상기 전압들은 임의로 설정된 초기 도전율 이용하여 계산된 전압들이다.
상기 계산된 도전율 은 이전 단계의 도전율 에 비하여 실제의 도전율 값에 더 가까워진다. 왜냐하면, 상기 계산된 도전율 은 상기 측정대상(S) 내부의 실제 자속밀도 가 고려된 값이기 때문이다.
본 발명에서는 한층 더 실제 도전율에 가까운 도전율을 얻기 위해 반복(iteration) 형태의 알고리즘을 이용한다. 즉, 이전 단계의 도전율 과 상기 개선된 도전율 의 차이를 적분값 로 계산하고, 상기 적분 값이 설정된 값 보다 작은지를 판단한다. 이때 상기 적분 값이 설정된 값 보다 크다면 작아질 때까지 m 값을 증가시켜 상술한 바와 같이 상기 수학식 10과 수학식 11을 반복 계산한다. 만약, 상기 이전 단계의 도전율 의 성분과 상기 계산된 도전율 의 성분간의 차가 설정된 값 보다 작으면 상기 도전율 의 성분들을 실제 도전율로 여기고, 상기 도전율 을 이용하여 상기 측정대상(S) 내부를 영상화한다.
그러나 상기 계산된 도전율 이 실제 자속밀도 를 고려해서 얻어진 것이라 해도 실제 도전율과 다를 수 있다. 즉, 상기 계산된 도전율 은 상기 자속밀도 에 따라 다른 값을 갖지만 실제 도전율에 비해 상수 배만큼 크거나 작을 수 있다. 실제 도전율에 대한 상기 계산된 도전율 의 상수 배 차이를 제거하기 위해 본 발명에서는 상기 측정된 표면 전압 와 상기 계산된 표면 전압 을 이용하여 상기 도전율 를 스케일(scale)하는 과정을 추가할 수 있다. 다시 말해서, 상기 측정된 표면 전압 와 상기 계산된 표면 전압 의 비율에 따라 상기 도전율 을 상수 배로 곱하거나 나눈다. 이로써 보다 더 정확한 도전율 값을 얻을 수 있다.
상술한 바와 같이, 본 발명은 인체 또는 물체 등과 같은 측정대상을 MRI 장 치 내에서 회전시키지 않고도 자속밀도의 한 방향 성분만을 이용하여 측정대상 내부의 도전율 영상 및 전류밀도 영상을 얻을 수 있다. 또한, 상기 자속밀도의 한 방향 성분을 한번만 미분하여 도전율을 계산하므로 간단한 연산과정을 통해 도전율을 얻을 수 있다.
또한, 전류주입수단을 기존의 MRI 장치와 함께 사용함으로써, 기존의 MRI 장치만으로는 얻을 수 없었던 측정대상 내부의 도전율 및 전류밀도를 얻을 수 있으며, 이를 이용하여 고화질의 영상을 재현할 수 있다. 게다가 인체 내의 각 지점에서 도전율의 방향성분을 알 수 있기 때문에 근육과 같은 인체조직에서의 도전율을 보다 정확하게 측정할 수 있다.
또한, 본 발명의 도전율 영상 및 전류밀도 영상을 기존의 의료장비에 적용시킬 수 있다. 따라서 생리작용에 따라 생체조직의 저항률이 변하는 장기들을 대상으로 한 기능영상(functional imaging)을 제공할 수 있으며 이로 인해 새로운 의료 검진 기술을 확보할 수 있다. 예를 들어, 호흡에 따른 폐의 용적 측정(폐에 공기가 흡입될 때와 배출될 때의 폐 내부의 저항률이 달라짐을 이용한 측정)이나, 심장의 심박출량 측정, 위장기능의 측정, 방광 잔료량 측정, 골절 치유상태의 측정, 또는 뇌기능 측정 등에 필요한 영상을 제공할 수 있다.
그리고, 인체조직의 전기적인 특성을 정량적으로 알 수 있으므로 각종 전자기 신경/근육 자극기(electric/magnetic stimulator), 세동 제거기(defibrillator), 심박 조율기(cardiac pacemaker) 등 전자기 에너지를 이용하는 진단기기나 치료기기의 최적화가 가능해진다.
이상에서 설명한 내용을 통해 당업자라면 본 발명의 기술 사상을 이탈하지 아니하는 범위에서 다양한 변경 및 수정이 가능함을 알 수 있을 것이다. 따라서 본 발명의 기술적 범위는 실시예에 기재된 내용으로 한정되는 것이 아니라 특허 청구의 범위에 의해 정해져야 한다.
Claims (7)
- 삭제
- 측정대상 내에 서로 다른 방향의 전류들을 주입하는 전류주입수단들과,상기 전류들에 의해 상기 측정대상 내에서 발생되는 자속밀도의 한 방향 성분 Bjz을 측정하는 MRI 장치와,상기 측정된 자속밀도의 한 방향 성분 Bjz을 이용하여 상기 측정대상의 내부 전압 Vj 을 계산하고, 아래의 연산식상기 계산된 도전율을 이용하여 상기 측정대상의 내부를 영상화하는 디스플레이 수단을 포함하는 도전율 및 전류밀도 영상화 시스템.
- 측정대상 표면에 부착된 전류주입수단들을 통해 상기 측정대상 내부에 서로 다른 방향의 전류들을 공급하는 단계와;상기 전류들에 의해 상기 측정대상 내에서 발생되는 자속밀도의 한 방향 성분 Bjz을 측정하는 단계와;상기 측정된 자속밀도의 한 방향 성분 Bjz을 이용하여 상기 측정대상의 내부 전압 Vj 을 계산하는 단계와;상기 계산된 도전율을 이용하여 상기 측정대상 내부를 영상화하는 단계를 포함하는 도전율 및 전류밀도 영상화 방법.
- 제 4 항에 있어서,상기 전류주입수단을 통해 상기 측정대상의 표면전압을 측정하는 단계와;임의의 도전율에 상응하는 상기 측정대상의 표면전압을 계산하는 단계와;상기 측정된 표면전압 및 상기 계산된 표면전압의 비율을 근거로 하여 상기 계산된 도전율을 보정하는 단계를 더 포함하는 것을 특징으로 하는 도전율 및 전류밀도 영상화 방법.
- 삭제
Priority Applications (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP08006022A EP1952764B1 (en) | 2003-11-25 | 2004-11-23 | System and method for visualizing conductivity and current density distribution in object |
DE602004027096T DE602004027096D1 (de) | 2003-11-25 | 2004-11-23 | System und Verfahren zum Anzeigen der Leitfähigkeit und Stromdichteverteilung eines Objekts |
EP04800117A EP1689295B1 (en) | 2003-11-25 | 2004-11-23 | System and method for visualizing conductivity and current density distribution in object |
US10/595,980 US7514921B2 (en) | 2003-11-25 | 2004-11-23 | System and method for visualizing conductive and current density distribution in object |
DE602004021625T DE602004021625D1 (de) | 2003-11-25 | 2004-11-23 | System und verfahren zur visualisierung der leitfähigkeit und stromdichteverteilung in einem objekt |
PCT/KR2004/003026 WO2005051163A2 (en) | 2003-11-25 | 2004-11-23 | System and method for visualizing conductivity and current density distribution in object |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR20030041569 | 2003-06-25 | ||
KR1020030041569 | 2003-06-25 |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20050005740A KR20050005740A (ko) | 2005-01-14 |
KR100570202B1 true KR100570202B1 (ko) | 2006-04-12 |
Family
ID=37219910
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020040001993A KR100570202B1 (ko) | 2003-06-25 | 2004-01-12 | 도전율 및 전류밀도를 영상화하는 시스템 및 방법 |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR100570202B1 (ko) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7204908B2 (ja) * | 2019-05-31 | 2023-01-16 | 旭化成株式会社 | 計測装置、計測方法、およびプログラム |
-
2004
- 2004-01-12 KR KR1020040001993A patent/KR100570202B1/ko not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
KR20050005740A (ko) | 2005-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1581110B1 (en) | System and method for three-dimensional visualization of conductivity and current density distribution in electrically conducting object | |
Davalos et al. | A feasibility study for electrical impedance tomography as a means to monitor tissue electroporation for molecular medicine | |
US7514921B2 (en) | System and method for visualizing conductive and current density distribution in object | |
US6397095B1 (en) | Magnetic resonance—electrical impedance tomography | |
US20020106681A1 (en) | High definition electrical impedance tomography methods for the detection and diagnosis of early stages of breast cancer | |
US7603158B2 (en) | Current density impedance imaging (CDII) | |
Kranjc et al. | Magnetic resonance electrical impedance tomography for monitoring electric field distribution during tissue electroporation | |
Seo et al. | Magnetic resonance electrical impedance tomography (MREIT): conductivity and current density imaging | |
Joy | MR current density and conductivity imaging: the state of the Aart | |
Li et al. | A novel design of implantable esophageal stent to reduce the MRI RF-induced heating | |
Wenger et al. | Water-content electrical property tomography (wEPT) for mapping brain tissue conductivity in the 200–1000 kHz range: Results of an animal study | |
KR100570202B1 (ko) | 도전율 및 전류밀도를 영상화하는 시스템 및 방법 | |
KR101787902B1 (ko) | 전기 임피던스 단층 촬영 영상 생성을 위한 전도도 산출 방법 및 장치 | |
KR100575210B1 (ko) | 이방성 도전율 및 전류밀도를 영상화하는 시스템 및 방법 | |
Nakane et al. | Forward electrocardiogram modeling by small dipoles based on whole-body electric field analysis | |
Kuratko et al. | Forward model of rat electroencephalogram: comparative study of numerical simulations with measurements on rat head phantoms | |
Gibson | Electrical impedance tomography of human brain function | |
KR100642837B1 (ko) | 도전율 및 전류밀도의 3차원 영상화 방법 및 그를 위한시스템 | |
Šemrov et al. | Numerical modeling for in vivo electroporation | |
KR100660557B1 (ko) | 자속밀도 처리 장치 및 방법 | |
KR100400978B1 (ko) | 인체내의 임피던스를 영상화하는 방법 및 장치 | |
Oh et al. | A single current density component imaging by MRCDI without subject rotations | |
KR100746274B1 (ko) | 측정대상 내부를 영상화하는 장치 및 방법 | |
Woo et al. | Conductivity images of biological tissue phantoms using a 3.0 Tesla MREIT system | |
Jog | Visualizing Transcranial Direct Current Stimulation (tDCS) in vivo using Magnetic Resonance Imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A201 | Request for examination | ||
E902 | Notification of reason for refusal | ||
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant | ||
FPAY | Annual fee payment |
Payment date: 20130404 Year of fee payment: 8 |
|
FPAY | Annual fee payment |
Payment date: 20140324 Year of fee payment: 9 |
|
FPAY | Annual fee payment |
Payment date: 20160317 Year of fee payment: 11 |
|
FPAY | Annual fee payment |
Payment date: 20170328 Year of fee payment: 12 |
|
FPAY | Annual fee payment |
Payment date: 20180404 Year of fee payment: 13 |
|
LAPS | Lapse due to unpaid annual fee |