KR101287426B1 - 극저자기장 핵자기 공명 물체 식별 방법 및 극저자기장 핵자기 공명 물체 식별 장치 - Google Patents

극저자기장 핵자기 공명 물체 식별 방법 및 극저자기장 핵자기 공명 물체 식별 장치 Download PDF

Info

Publication number
KR101287426B1
KR101287426B1 KR1020110109971A KR20110109971A KR101287426B1 KR 101287426 B1 KR101287426 B1 KR 101287426B1 KR 1020110109971 A KR1020110109971 A KR 1020110109971A KR 20110109971 A KR20110109971 A KR 20110109971A KR 101287426 B1 KR101287426 B1 KR 101287426B1
Authority
KR
South Korea
Prior art keywords
magnetic field
signal
parameter vector
spin
measuring
Prior art date
Application number
KR1020110109971A
Other languages
English (en)
Other versions
KR20130045637A (ko
Inventor
김기웅
강찬석
이성주
이용호
유권규
Original Assignee
한국표준과학연구원
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 한국표준과학연구원 filed Critical 한국표준과학연구원
Priority to KR1020110109971A priority Critical patent/KR101287426B1/ko
Priority to PCT/KR2012/008426 priority patent/WO2013062260A1/ko
Priority to CN201280052725.4A priority patent/CN103890573B/zh
Publication of KR20130045637A publication Critical patent/KR20130045637A/ko
Application granted granted Critical
Publication of KR101287426B1 publication Critical patent/KR101287426B1/ko
Priority to US14/262,251 priority patent/US9759793B2/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/448Relaxometry, i.e. quantification of relaxation times or spin density
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/445MR involving a non-standard magnetic field B0, e.g. of low magnitude as in the earth's magnetic field or in nanoTesla spectroscopy, comprising a polarizing magnetic field for pre-polarisation, B0 with a temporal variation of its magnitude or direction such as field cycling of B0 or rotation of the direction of B0, or spatially inhomogeneous B0 like in fringe-field MR or in stray-field imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/323Detection of MR without the use of RF or microwaves, e.g. force-detected MR, thermally detected MR, MR detection via electrical conductivity, optically detected MR
    • G01R33/326Detection of MR without the use of RF or microwaves, e.g. force-detected MR, thermally detected MR, MR detection via electrical conductivity, optically detected MR involving a SQUID

Abstract

본 발며은 극저자기장 핵자기 공명 물체 식별 방법 및 극저자기장 핵자기 공명 물체 식별 장치를 제공한다. 이 방법은 측정 대상에 사전 자화 자기장의 세기 또는 사전 자화 자기장의 인가 시간을 변경하면서 스핀-격자 완화 시간(T1)을 측정하는 단계, 및 스핀-격자 완화 시간을 이용하여 상기 측정 대상를 분류하는 단계를 포함한다.

Description

극저자기장 핵자기 공명 물체 식별 방법 및 극저자기장 핵자기 공명 물체 식별 장치{Target Discrimination Method by using Ultra-Low Field NMR and Ultra-Low Field NMR Target Discrimination Apparatus}
본 발명은 극저자기장 물체 식별 방법에 관한 것으로, 더 구체적으로 스핀-격자 완화시간이 포화되지 않는 낮은 사전 자화 자기장 또는 검출 자기장을 변경하면서 측정한 스핀-격자 완화 시간을 이용하여 측정 대상를 구분하는 극저자기장 물체 식별 방법에 관한 것이다.
NMR은 모든 물질을 구성하는 원자핵에서 자기장을 인가할 때 원자핵의 자기 스핀이 상기 자기장에 공명하여 세차하는 현상이다. 기존의 NMR의 주자석은 0.1 테슬라 내지 수 테슬라의 자기장을 공간적으로 균일하게 형성해야 한다. 따라서, 상기 자기장을 형성하는 초전도 주자석은 부피가 크고, 비싸다.
저 자기장/ 극저 자기장 NMR는 기존의 NMR의 주자석에서 생성되는 자기장을 사전자화 자기장과 검출 자기장의 두 부분으로 분리한다. 이에 따라, 상기 검출 자기장의 크기는 수십 마이크로 테슬라 내지 수 마이크로 테슬라 정도로 낮출 수 있다.
저 자기장/ 극저 자기장 NMR은 상기 사전 자화 자기장을 형성하는 사전 자화 코일과 상기 검출 자기장을 형성하는 검출 코일을 포함한다. 상기 사전 자화 코일은 사전 자화 자기장을 생성하여 측정하고자 하는 대상 물체를 사전 자화시킨다. 이어서, 상기 사전자화 자기장을 제거하고, 상기 검출 자기장을 상기 대상 물체에 인가하여, 상기 대상 물체의 자성이 이완되는 시간 동안 나타나는 핵자기 공명 신호를 검출기로 측정한다. 이에 따라, 상기 사전 자화 코일은 상기 사전 자화 자기장의 균일도가 떨어지더라도 강한 자기장만 만들면 된다. 또한, 상기 검출 코일은 규일하지만 약한 검출 자기장만 만들면 된다. 따라서 극저 자기장/ 저 자기장 NMR은 그 구조가 간단하고 가격도 저렴할 수 있다. 상기 검출 자기장의 크기가 감소함에 따라 자기장의 크기에 비례하는 라모어(Larmor) 공명 주파수에 해당하는 검출 신호의 주파수는 수 kHz 내지 수백 Hz 정도로 감소한다.
이에 따라, 종래의 고 자기장 NMR에서 볼 수 없었던 현상의 측정이 가능하다. 저 자기장/ 극저 자기장 NMR은 금속에 의한 왜곡이 현저히 적다. 따라서, 금속 캔 내부의 물질 또는 포장된 물질은 비파괴적으로 측정될 수 있다.
본 발명의 해결하고자 하는 일 기술적 과제는 극저자기장 NMR을 이용한 물체 식별 방법을 제공하는 것이다.
본 발명의 해결하고자 하는 일 기술적 과제는 극저자기장 NMR을 이용한 물체 식별 장치를 제공하는 것이다.
본 발명의 일 실시예에 따른 극저자기장 핵자기 공명 물체 식별 방법은 측정 대상에 복수의 사전 자화 자기장의 세기에서 각각의 스핀-격자 완화 시간(T1)을 측정하는 단계; 및 상기 스핀-격자 완화 시간을 이용하여 상기 측정 대상을 분류하는 단계를 포함한다.
본 발명의 일 실시예에 있어서, 각각의 스핀-격자 완화 시간(T1)들을 측정하는 단계는 상기 측정 대상의 개별적 완화시간 특성에 적합한 펄스 시퀀스를 설정하는 단계; 상기 사전 자화 자기장 및 검출 자기장을 인가하는 단계; 자유 유도 감쇄(FID) 신호를 측정하는 단계; 상기 FID 신호를 푸리어 변환하여 소정의 원자핵의 공명 주파수에서의 신호 크기 또는 선폭을 측정하는 단계; 상기 사전 자화 자기장의 인가 시간을 변경하는 단계; 상기 사전 자화 자기장의 인가 시간에 따른 신호 크기 변화를 피팅하여 T1을 구하는 단계; 및 변경된 복수의 상기 사전 자화 자기장 크기에 대해, 위의 단계들을 반복하여 각각의 T1을 구하여 제1 파라미터 벡터를 구성하는 단계를 포함한다. 상기 제1 파라미터 벡터를 이용하여 상기 측정 대상의 품종이 확인된다.
본 발명의 일 실시예에 있어서, 그레디언트 에코 신호를 이용하여 상기 측정 대상의 스핀-스핀 완화 시간(T2)을 측정하는 단계; 및 상기 스핀-스핀 완화 시간을 이용하여 상기 측정 대상을 분류하는 단계를 더 포함할 수 있다.
본 발명의 일 실시예에 있어서, 상기 측정 대상의 스핀-스핀 완화 시간(T2)을 측정하는 단계는 상기 측정 대상의 개별적 완화시간 특성에 적합한 펄스 시퀀스를 설정하는 단계; 상기 사전 자화 자기장 및 검출 자기장을 인가하는 단계; 양의 경사 자기장 및 음의 경사 자기장을 연속적으로 인가하는 단계; 그레디언트 에코 신호를 측정하는 단계; 상기 그레디언트 에코 신호를 푸리어 변환하여 소정의 원자핵의 공명 주파수에서의 신호 선폭을 측정하는 단계; 및 소정의 원자핵의 공명 주파수에 따른 선폭 또는 선폭의 역수를 제2 파라미터 벡터로 구성하는 단계를 포함한다. 상기 제2 파라미터 벡터를 이용하여 품종이 하여 품종이 확인될 수 있다.
본 발명의 일 실시예에 따른 극저자기장 핵자기 공명 물체 식별 방법은 측정 대상에 검출 자기장의 세기를 변경하면서 자유 유도 감쇄 신호를 측정하고 푸리어 변환하여 소저의 원자핵의 공명 주파수에서의 신호 선폭을 측정하는 단계; 및 상기 원자핵의 공명 주파수에서의 신호 선폭을 이용하여 상기 측정 대상을 분류하는 단계를 포함한다.
본 발명의 일 실시예에 있어서, 측정 대상에 검출 자기장의 세기를 변경하면서 자유 유도 감쇄 신호를 측정하고 푸리어 변환하여 원자핵의 공명 주파수의 선폭을 측정하는 단계는 상기 측정 대상의 개별적 완화시간 특성에 적합한 펄스 시퀀스를 설정하는 단계; 사전 자화 자기장을 인가하는 단계; 상기 검출 자기장을 인가하는 단계; 자유 유도 감쇄(FID) 신호를 측정하는 단계; 상기 FID 신호를 푸리어 변환하여 소정의 원자핵의 공명 주파수에서의 신호 선폭을 측정하는 단계; 상기 검출 자기장의 세기를 변경하는 단계; 및 위의 단계의 반복적인 측정으로 상기 검출 자기장의 세기에 따른 소정의 원자핵의 공명 주파수에서의 신호 선폭으로 제3 파라미터 벡터를 구성하는 단계를 포함할 수 있다. 상기 제3 파라미터 벡터를 이용하여 상기 측정 대상의 품종이 확인될 수 있다.
본 발명의 일 실시예에 따른 극저자기장 핵자기 공명 물체 식별 방법은 복수의 사전 자화 자기장 크기에 대해 각각의 스핀-격자 완화 시간(T1)을 구하여 제1 파라미터 벡터를 구성하는 단계; 그레디언트 에코 신호를 이용하여 소정의 원자핵의 공명 주파수에 따른 선폭 혹은 선폭의 역수(T2)을 제2 파라미터 벡터로 구성하는 단계; 검출 자기장의 세기에 따른 소정의 원자핵의 공명 주파수에서의 신호 선폭의 값들로 제3 파라미터 벡터를 구성하는 단계; 상기 제1 내지 제3 파라미터 벡터들의 성분들을 조합하여 새로운 파라미터 벡터를 구성하는 단계; 및 상기 새로운 파라미터 벡터를 이용하여 측정 대상의 품종을 확인하는 단계를 포함한다.
본 발명의 일 실시예에 따른 극저자기장 핵자기 공명 물체 식별 방법은 복수의 사전 자화 자기장 크기에 대해 각각의 스핀-격자 완화 시간(T1)을 구하여 제1 파라미터 벡터를 구성하는 단계; 그레디언트 에코 신호를 이용하여 소정의 원자핵의 공명 주파수에 따른 선폭 혹은 선폭의 역수(T2)을 제2 파라미터 벡터로 구성하는 단계; 검출 자기장의 세기에 따른 소정의 원자핵의 공명 주파수에서의 신호 선폭의 값들로 제3 파라미터 벡터를 구성하는 단계; 상기 제1 내지 제3 파라미터 벡터들의 성분들을 조합하여 새로운 파라미터 벡터를 구성하는 단계; 구별하려는 여러 가지 물질의 반복 측정으로 각각의 물질에 대한 새로운 파라미터 벡터들을 수집하는 단계; 수집된 새로운 파라미터 벡터들에 서포트벡터머신(support vector machine;SVM) 또는 주요소분석법(Principal component analysis;PCA)을 적용하여 구분에 효과적인 주요 벡터 성분을 찾아내고, 새로운 파라미터 벡터를 주요 벡터 성분으로 사영시킴으로써 파라미터 벡터의 차원을 줄이는 단계; 및 차원 감소된 파라미터 벡터 성분을 이용하여 고속으로 물질의 품종을 확인하는 단계를 포함한다.
본 발명의 일 실시예에 따른 극저자기장 핵자기 공명 물체 식별 장치는 저온 액체 냉매를 담고 있는 듀아; 측정 대상을 사전 자화시키는 사전자화 코일; 상기 사전자화 코일에 전류를 단속적으로 인가하여 사전자화 자기장(Bp)을 형성하는 사전자화 코일 구동부; 상기 사전자화 코일에 의하여 상기 사전자화 자기장(Bp)이 인가되는 측정 대상에 의한 핵자기 공명 신호를 측정하는 센서부; 측정 대상에 검출 자기장을 인가하는 검출 코일; 및 상기 검출 코일에 전류를 인가하여 상기 측정 대상에 검출 자기장(readout magnetic field)을 인가하는 검출 자기장 전원부를 포함한다. 상기 센서부는 상기 측정 대상에 복수의 사전 자화 자기장의 세기에서 각각의 스핀-격자 완화 시간(T1)을 측정하고, 상기 스핀-격자 완화 시간을 이용하여 상기 측정 대상을 분류한다.
본 발명의 일 실시예에 있어서, 상기 측정 대상에 경사 자기장을 인가하는 경사 자기장 코일; 상기 경사 자기장 코일에 전류를 공급하는 경사 자기장 구동부; 상기 경사 자기장 구동부에 전원을 제공하는 경사 자기장 전원부; 및 상기 사전 자화 코일 구동부, 상기 검출 자기장 전원부, 및 검출 자기장 전원부에 펄스 스퀀스를 제공하는 펄스 시퀀스 발생부를 더 포함할 수 있다.
본 발명의 일 실시예에 따른 물체 식별 방법은 종래의 DNA 검사를 하지 않고 비파괴적으로 농축산물의 품종을 알아낼 수 있다.
도 1은 본 발명의 일 실시예에 따른 품종 판별을 위한 저자기장 핵자기 공명 장치를 설명하는 도면이다.
도 2는 도 1의 극저자기장 핵자기 공명 장치를 이용한 물체 식별 방법을 설명하는 흐름도이다.
도 3a 및 도 3b는 본 발명의 일 실시예에 따른 펄스 시퀀스를 설명하는 도면이다.
도 4는 도 3a의 펄스 시퀀스에 의하여 획득된 FID 신호의 푸리어 변환 결과이다.
도 5는 사전 자화 자기장의 세기 또는 인가 시간에 따른 라모어 공명 주파수에서의 진폭을 나타내는 도면이다.
도 6은 본 발명의 다른 실시예에 따른 펄스 시퀀스를 설명하는 도면이다.
도 7은 본 발명의 일 실시예에 따른 Bp에 따른 T1의 의존성을 설명하는 도면이다.
도 8은 본 발명의 또 다른 실시예에 따른 T2를 측정하는 펄스 시퀀스를 나타내는 도면이다.
도 9는 도 8의 T2를 측정하기 위한 펄스 시퀀스를 설명하는 도면이다.
도 10은 본 발명의 일 실시예에 따른 물체 식별 방법을 설명하는 도면이다.
도 11은 도 10의 FID 신호를 푸리어 변환하여 얻은 신호의 선폭을 측정하기 위한 펄스 시퀀스를 설명하는 도면이다.
도 12는 도 10의 물체 식별 방법을 이용하여 획득한 가상 데이터이다.
종래의 고 자기장 NMR/MRI에서, MR Signal은 수소 원자핵에서 나온 신호이며, 몸속의 물 분자의 H 원자핵 농도에 따라 영상이 나타난다. 따라서, 자기 공명 영상은 공간적인 수소 원자핵 농도를 나타낸다. 고 자기장 NMR/MRI에서, T1은 스핀-격자 완화 시간(Spin-lattice relaxation time; T1)이고, T2는 스핀-스핀 완화 시간(spin-spin relaxation time; T2)이다. 생체 조직(tissue)은 고정된 주 자기장하에서 서로 다른 T1 및/또는 T2를 가질 수 있다. MRI(Magnetic Resonance Imagnig)는 공간적인 T1 또는 T2의 공간적 분포이다. 그러나, 고 자기장 NMR/MRI는 주자기장의 세기를 변화시키기 어려워 농축산물의 품종을 알아내기 어렵다.
극저 자기장 NMR에서 스핀-격자 완화 시간(Spin-lattice relaxation time; T1)은 사전 자화 자기장(Bp) 및/또는 검출 자기장(Bm)의 함수일 수 있다. 극저 자기장 NMR는 사전 자화 자기장(Bp) 및/또는 검출 자기장(Bm)을 용이하게 변경할 수 있다. 따라서, 사전 자화 자기장(Bp)에 따라 측정된 T1, 검출 자기장(Bm)에 따라 측정된 FID 신호의 선폭, 또는 그레디언트 에코 신호를 처리하여 얻은 T2는 측정 대상의 분별에 사용될 수 있다.
예를 들어, 한 품종의 특정한 생체 조직은 다른 품종의 특정한 생체 조직과 서로 다른 T1 의존성을 가질 수 있다. 측정된 T1 의존성은 육안으로 식별 불가능한 측정 대상의 품종을 구분하기 위하여 사용될 수 있다. 예를 들어, 측정 대상은 과일, 채소, 주류, 어류, 및 육류일 수 있다.
한국산 소(한우)와 미국산 소는 서로 다른 품종으로, 현재 DNA 검사에 의해서만 판별될 수 있다. 그러나, DNA 검사는 많은 시간 및 측정 대상의 파괴를 요구한다. 따라서, 축산물의 품종을 구별하는 새로운 방법이 요구된다. 따라서, 극저 자기장 NMR은 축산물의 품종을 구별에 사용될 수 있다.
본 발명의 일 실시예에 따르면, 측정 대상의 품종을 판별하기 위하여 서로 다른 조건에서의 T1이 측정될 수 있다.
이하, 첨부한 도면들을 참조하여 본 발명의 바람직한 실시예들을 상세히 설명하기로 한다. 그러나, 본 발명은 여기서 설명되어지는 실시예들에 한정되지 않고 다른 형태로 구체화될 수도 있다. 오히려, 여기서 소개되는 실시예는 개시된 내용이 철저하고 완전해질 수 있도록 그리고 당업자에게 본 발명의 사상이 충분히 전달될 수 있도록 하기 위해 제공되어지는 것이다. 도면들에 있어서, 구성요소는 명확성을 기하기 위하여 과장되어진 것이다. 명세서 전체에 걸쳐서 동일한 참조번호로 표시된 부분들은 동일한 구성요소들을 나타낸다.
도 1은 본 발명의 일 실시예에 따른 품종 판별을 위한 극저자기장 핵자기 공명 장치를 설명하는 도면이다.
도 1을 참조하면, 상기 극저자기장 핵자기 공명 물체 식별 장치(100)는 저온 액체 냉매를 담고 있는 듀아(112), 측정 대상을 사전 자화시키는 사전자화 코일(122), 상기 사전자화 코일(122)에 전류를 단속적으로 인가하여 사전자화 자기장(Bp)을 형성하는 사전자화 코일 구동부(124), 상기 사전자화 코일(122)에 의하여 상기 사전자화 자기장(Bp)이 인가되는 상기 측정 대상(130)에 의한 핵자기 공명 신호를 측정하는 센서부(160), 상기 측정 대상(130)에 검출 자기장(Bm)을 인가하는 검출 코일(152), 및 상기 검출 코일(152)에 전류를 인가하여 상기 측정 대상(130)에 상기 검출 자기장(readout magnetic field, Bm)을 인가하는 검출 자기장 전원부(154)를 포함한다.
상기 센서부(160)는 상기 측정 대상(130)에 사전 자화 자기장(Bp)의 세기 또는 사전 자화 자기장(Bp)의 인가 시간을 변경하면서 획득된 핵자기 공명 신호를 처리하여 추출된 스핀-격자 완화 시간(T1)을 이용하여 상기 측정 대상의 품종을 식별한다.
상기 센서부(160)는 상기 측정 대상에 복수의 사전 자화 자기장의 세기에서 각각의 스핀-격자 완화 시간(T1)을 측정하고, 상기 스핀-격자 완화 시간을 이용하여 상기 측정 대상을 분류한다.
경사 자기장 코일(142)은 상기 측정 대상(130)에 경사 자기장(BG)을 인가할 수 있다. 경사 자기장 구동부(144)는 상기 경사 자기장 코일(142)에 전류를 공급하여 상기 측정 대상(130)에 상기 경사 자기장(BG)을 인가한다. 경사 자기장 전원부는 상기 경사 자기장 구동부에 전원을 제공할 수 있다. 상기 경사 자기장 전원부(146)는 펄스 시퀸스 발생부(172)로부터 펄스 시퀀스를 제공받아 상기 경사 자기장 구동부(144)에 전력을 공급한다.
상기 펄스 시퀀스 발생부(172)는 제어부(165)의 제어 신호를 입력받아 상기 사전 자화 코일 구동부(124), 상기 검출 자기장 전원부(154), 및 경사 자기장 전원부(146)에 펄스 스퀀스를 제공한다.
상기 센서부(160)는 자속을 감지 및/또는 감쇠/증폭하는 자속변환부(161), 상기 자속변환부(161)의 출력 신호를 입력받아 자기장을 검출해 전압신호로 변환하는 SQUID(163), 상기 SQUID(163)의 전압 신호를 선형화하여 검출된 자기장에 비례하는 전압 신호로 제공하는 FLL부(168), 상기 FLL부(168)의 상기 선형화된 전압 신호를 처리하여 노이즈를 제거하고 증폭하는 센서 신호 처리부(166), 및 상기 FLL부(168)에 제어 신호를 제공하는 센서 제어부(167)를 포함할 수 있다.
상기 SQUID(163)는 초전도체만이 지닌 조셉슨 효과와 자속의 양자화 현상을 결합하여 외부 자속의 변화를 전압으로 변환하는 일종의 트랜스듀서(transducer)이다. 상기 SQUID(163)는 한 개의 초전도 루프에 한 개 또는 두 개의 조셉슨 접합을 끼워 결합한 자기센서이다. RF SQUID는 초전도루프에 한 개의 조셉슨 접합을 삽입하여 결합할 수 있다. DC SQUID는 초전도루프에 두 개의 조셉슨 접합을 삽입하여 결합할 수 있다. 상기 RF SQUID는 RF 주파수 대역의 교류 전압을 출력하고 그 주파수가 가해진 자속에 따라 바뀌는 방식으로 작동하며, 상기 DC SQUID는 가해진 자속에 대한 함수로 직류 전압을 발생하는 방식으로 작동하며 이 함수는 자속의 양자 값인 Φ0(= 2.07 x 10^(-15) Wb)를 주기로 하여 진동하는 형태로 주어진다. 상기 자속/전압 변환 함수의 구체적인 형태는 상기 DC SQUID의 구체적인 구조에 따라 결정될 수 있다.
상기 자속변환부(161)는 자속을 감지하여 초전도 전류로 변환하는 픽업(pick-up) 코일 및/또는 이를 상기 SQUID(163)에 다시 자속 형태로 증폭 또는 감쇠시켜 전달해 주는 입력(input)코일을 포함할 수 있다. 상기 자속변환부(161)는 초전도체로 구성될 수 있다. 상기 픽업 코일은 많은 자속을 감지하기 위하여 넓은 면적을 지닐 수 있다. 상기 입력 코일은 SQUID에 집속시키기 위하여 SQUID와 비슷한 면적을 지니고 있고 그 증폭율 또는 감쇠율을 변화시키기 위하여 여러 번 감길 수 있다. 상기 자속변환부(161)는 픽업 코일이 하나의 루프로 구성된 마그네토미터(Magnetometer) 또는 픽업 코일이 서로 반대 방향으로 감긴 루프 쌍이 하나 또는 그 이상의 수 로 구성된 그레디오미터 (Gradiometer)를 포함할 수 있다.
상기 SQUID(163)는 도선을 통하여 상기 FLL부(168)와 연결될 수 있다. 상기SQUID(163)를 사전자호 자기장(Bp)와 같이 매우 큰 자기장하에서 안정적으로 동작시키기 위해서 보호가 필요하다. 따라서 극저자장-MRI 시스템에서는 SQUID의 보호를 위해 초전도 차폐(164)를 이용한다. 그러나 SQUID 센서 전체를 초전도 차폐할 경우 SQUID가 자기장 감지 센서로서의 기능을 할 수 없다. 따라서 초전도체를 이용한 차폐를 할 때는 SQUID 부분과 자속 변환부의 입력 코일 부분만을 초전도 차폐하고, 검출코일은 초전도 차폐 외부에 놓인다. 이 경우, SQUID 자체는 초전도 차폐(164)를 함으로써 강한 자기장으로부터 보호되나, 검출 코일로부터 유도된 전류가 SQUID로 인가되는 것은 막을 수가 없다. 그러므로 극저자장-NMR 시스템에서는 이러한 검출코일에서 유도된 과전류가 SQUID로 인가되는 것을 막기 위한 전류 제한부(162)가 배치된다.
FLL(flux locked loop)부(168)는 상기 SQUID(163)의 출력신호를 받아들이는 입력단, 적분기, 피드백(feedback) 방식의 선형화 회로, 및 피드백 코일 등을 포함할 수 있다. 상기 FLL부(168)는 자속의 변화량을 자속 양자값 Φ0 보다 훨씬 넓은 범위를 가지는 전압 신호로 변환하여 출력할 수 있다.
검출 자기장 생성부(155)는 공간적으로 균일하고 낮은 검출 자기장(Bm)을 생성할 수 있다. 상기 검출 자기장 생성부(155)는 검출 자기장 전원부(154), 검출 자기장 게이트부(156), 및 검출 코일(152)을 포함할 수 있다. 상기 검출 자기장 게이트부(156)는 단속적으로 상기 검출 자기장(Bm)을 생성하도록 상기 검출 코일(152)에 인가되는 전류를 조절할 수 있다.
경사 자기장 생성부(145)는 y 성분의 경사 자기장(dBy/dy)을 생성하도록 경사 자기장 전원부(146), 경사 자기장 구동부(144), 및 경사 자기장 코일(142)을 포함할 수 있다.
펄스 시퀀스 발생부(172)는 펄스 시퀀스를 생성하여 T1, T2, FID 신호의 선폭을 얻을 수 있도록 사전 자화 코일 구동부(124), 검출 자기장 전원부(154), 및 상기 경사 자기장 전원부(146)에 제공할 수 있다.
제어부(165)는 상기 센서 신호처리부(166)의 신호를 처리하고 상기 펄스 시퀀스 발생부(172), 및 센서 센서부(167)를 제어할 수 있다.
한편, 사전 자화 자기장(Bp), 검출 자기장(Bm), 및 경사 자기장(BG)을 켜고 끄는 스위치로 광학 에스에스알(Optical Solid State Relay; SSR)이 사용될 수 있다. 상기 SSR이 꺼져있는 동안에는 사전 자화 코일(122), 검출 자기장 코일(152), 및 경사 자기장 코일(142)은 전류원과 완전히 단락될 수 있다. 상기 SSR을 구동하기 위한 TTL 신호는 광통신을 통해 인가될 수 있다. 이에 따라, 상기 센서부(160)에 악영향을 줄 수 있는 모든 전기적인 연결은 제거될 수 있다.
상기 사전 자화 코일(122)은 측정 대상(130)을 사전 자화시키기 위한 사전 자화 자기장(Bp)을 생성할 수 있다. 상기 사전 자화 코일(Bp)은 외부 지름이 30 mm 이고 길이가 60 mm 인 보어에 지름 1 mm 인 구리선이 총 240 회 감겨있는 솔레노이드 형태일 수 있다. 상기 사전 자화 코일(122)의 총 저항은 0.75 오옴이고, 인덕턴스는 1 kHz에서 0.79 mH이다.
상기 센서부(160)에 미치는 잡음의 영향을 최소화하기 위하여 전류원으로 12 V Pb 배터리가 사용될 수 있다. 분해능이 0.01 mT 인 가우스 미터로 측정한 사전 자화 코일 내부 중심에서의 자기장은 약 50 mT일 수 있다.
검출 코일(152)은 측정 대상에 검출 자기장(Bm)을 인가할 수 있다. 상기 검출 코일은 한 변의 길이가 1340 mm인 사각형 헬름홀츠 코일(Helmholtz coil)일 수 있다. 상기 검출 코일의 중심에서의 자기장은 플럭스 케이트를 이용하여 측정하였을 때 약 2.82 Guass이다. 이 크기의 자기장에 해당하는 수소원자핵의 라모어(Lamour) 주파수는 약 120 Hz 이다.
상기 검출 자기장의 세기에 따라 수소 원자핵, 탄소 원자핵, 또는 질소 원자핵의 공명(Lamour) 주파수는 변경될 수 있다. 따라서, 소정의 검출 자기장에서 수소 원자핵, 탄소 원자핵, 또는 질소 원자핵에 대한 라모어(Lamour) 주파수는 고정될 수 있다. 자유 유도 감쇄 신호(Free induction decay; FID Singnal)를 푸리어 변환한 경우, 소정의 라모어 주파수에서 진폭은 사전 자화 자기장의 세기 또는 고정된 사전 자화 자기장의 세기에서 인가 시간에 따라 변한다. 따라서, T1이 획득될 수 있다.
경사 자기장 코일(142)은 경사 자기장을 인가할 수 있다. 경사 자기장 코일(142)은 마주보는 두 코일에 극성이 반대인 자장이 인가되는 맥스웰 페어(Maxwell pair) 형태일 수 있다. 상기 경사 자기장 코일(142)의 한 변의 길이는 1420 mm 이며, 중심에서의 자기장 기울기는 플럭스 게이트를 이용하여 측정한 결과 약 31.3 pT/cm 이다. 상기 자기장 기울기(dBy/dy)는 주파수 성분으로 표시하면 약 1.33 Hz/cm 에 해당한다.
도 2는 도 1의 극저 자기장 핵자기 공명 장치를 이용한 물체 식별 방법을 설명하는 흐름도이다.
도 3a 및 도 3b는 본 발명의 일 실시예에 따른 펄스 시퀀스를 설명하는 도면이다.
도 2, 도 3a 및 도 3b를 참조하면, 극저 자기장 핵자기 공명 물체 식별 방법은 측정 대상에 복수의 사전 자화 자기장의 세기에서 각각의 스핀-격자 완화 시간(T1)을 측정하는 단계(S111), 및 상기 스핀-격자 완화 시간을 이용하여 상기 측정 대상을 분류하는 단계(S144)를 포함한다.
측정 대상은 극저자기장 핵자기 공명 장치에 설치된다. 상기 극저자기장 핵자기 공명 장치의 펄스 시퀀스는 결정된다. T1을 측정하기 위한 펄스 시퀀스는 2 가지 방식을 가질 수 있다.
하나는 고정된 사전 자화 자기장의 세기에서 사전 자화 자기장의 인가 시간을 변경하면서 자유 유도 감쇄(FID) 신호를 측정하는 것이다. 이에 따라, 사전 자화 자기장의 인가 시간에 따라 측정 대상의 자화 정도가 결정된다. 측정된 FID 신호는 푸리어 변환(Fourier Transformation)될 수 있다(S116). 이어서, 주파수에 따른 푸리어 변환에 의하여 진폭이 구해진다. 수소 원자핵의 공명 주파수에서, 진폭이 구해진다(S118).
수소 원자핵의 공명 주파수에서 진폭은 사전 자화 자기장의 인가 시간에 따라 측정될 수 있다(S124). 상기 인가 시간에 따른 수소 원자핵의 공명 주파수에서 진폭은 상기 인가 시간이 증가함에 따라 포화된다. 상기 진폭은 인가 시간에 따라 지수 함수 의존성을 가지고, T1은 지수 함수로 피팅되어 구해질 수 있다(S142).
각각의 스핀-격자 완화 시간(T1)을 측정하는 단계(S111)는 상기 측정 대상의 개별적 완화시간 특성에 적합한 펄스 시퀀스를 설정하는 단계(S110), 상기 사전 자화 자기장 및 검출 자기장을 인가하는 단계(S112), 자유 유도 감쇄(FID) 신호를 측정하는 단계(S114), 상기 FID 신호를 푸리어 변환하여 소정의 원자핵의 공명 주파수에서의 신호 크기 또는 선폭을 측정하는 단계(S118), 상기 사전 자화 자기장의 인가 시간을 변경하는 단계(S124), 상기 사전 자화 자기장의 인가 시간에 따른 신호 크기 변화를 피팅하여 T1을 구하는 단계(142), 및 변경된 복수의 상기 사전 자화 자기장 크기에 대해, 위의 단계들을 반복하여 각각의 T1을 구하여 제1 파라미터 벡터를 구성하는 단계(S134)를 포함한다. 상기 제1 파라미터 벡터를 이용하여 상기 측정 대상의 품종이 확인된다.
도 3a를 참조하면, T1 측정을 위한 펄스 시퀸스는 사전 자화 자기장(Bp1) 및 검출 자기장(Bm)을 순차적으로 또는 동시에 인가할 수 있다. 상기 검출 자기장(Bm)은 상기 사전 자화 자기장(Bp1)이 꺼지 전에 인가되는 것이 바람직하다. 상기 사전 자화 자기장의 세기(Bp1)는 일정하다. 상기 사전 자화 자기장의 인가 시간(t_P1)은 펄스 시퀸스에 따라 변경된다. 상기 사전 자화 자기장(Bp1)이 꺼진 후 소정의 FLL 구동 시간(t_FLL)이 지난 시간부터 FID 신호를 측정 시간(t_measure) 동안 측정할 수 있다.
도 3b를 참조하면, 사전 자화 자기장의 인가 시간(t_P2)은 도 3a의 사전 자화 자기장의 인가 시간(t_P1)과 다르다(S124). 따라서, 측정 대상의 자화 정도에 따라, FID 신호의 초기 진폭은 변한다.
도 4는 도 3a의 펄스 시퀀스에 의하여 획득된 FID 신호의 푸리어 변환 결과이다.
도 4를 참조하면, 사전 자화 자기장의 인가 시간에 따른 수소 원자핵의 공명 주파수에서의 진폭을 나타내는 도면이다.
도 4를 참조하면, 측정 대상의 1H 양성자의 자유유도감쇠(Free Induction Decay : FID) 신호의 FFT(Fast Fourier Transforamtion)를 보여준다. T2 *는 약 2.5 초 이상으로 나타났고, 신호 대 잡음비는 약 14.5 로 나타났다. FFT는 제로 필링(zero filling) 후에 수행되었다. 약 120 Hz 근방에서 주파수 피크가 나타난다. 이 피크의 위치는 인가한 2.82 uT (micro Tesla)크기의 검출 자기장에 해당하는 물의 공명 주파수인 120 Hz와 잘 일치한다. NMR 시스템의 공간 분해능을 의미하는 스펙트랄 레졸류션(spectral resolution)의 반치폭 (Full-Width-Half-Maximum : FWHM)는 약 0.5 Hz이다.
도 5는 사전 자화 자기장의 세기 또는 인가 시간에 따른 라모어 공명 주파수에서의 진폭을 나타내는 도면이다.
도 5를 참조하면, 수소 원자핵의 공명 주파수에서 진폭은 사전 자화 자기장(Bp1)의 인가 시간(t_P) 또는 세기에 따라 변한다. 사전 자화 자기장의 세기 또는 인가 시간이 증가함에 따K라, 측정 대상의 자화는 포화된다. 따라서, 측정 대상의 자화는 상기 진폭에 의존한다. 따라서, 상기 진폭은 인가 시간에 따라 지수 함수 의존성을 가지고, T1은 지수 함수로 피팅되어 구해질 수 있다(S142).
도 6은 본 발명의 다른 실시예에 따른 펄스 시퀀스를 설명하는 도면이다.
도 2 및 도 6을 참조하면, T1은 고정된 사전 자화 자기장(Bp)의 인가시간에서 사전 자화 자기장(Bp)의 세기를 변경하면서 FID 신호를 측정하여 구해질 수 있다(S134). 이에 따라, 사전 자화 자기장의 세기에 따라 측정 대상의 자화 정도가 결정된다. 측정된 FID 신호는 푸리어 변환(Fourier Transformation)될 수 있다(S116). 이어서, 주파수에 따른 푸리어 변환된 진폭이 구해진다(S118).
수소 원자핵의 공명 주파수에서 진폭은 사전 자화 자기장의 세기에 따라 측정될 수 있다. 사전 자화 자기장의 세기에 따른 수소 원자핵의 공명 주파수에서 진폭은 상기 세기가 증가함에 따라 포화된다. 상기 진폭은 세기에 따라 지수 함수의존성을 가지고, T1은 지수 함수로 피팅되어 구해질 수 있다(S142).
사전 자화 자기장의 인가 시간(t_P)은 일정하고, 사전 자화 자기장(Bp)의 세기가 펄스 시퀀스에 따라 증가 또는 감소한다. 사전 자화 자기장(Bp)이 꺼진 후 소정의 FLL 구동 시간(t_FLL)이 지난 시간부터 FID 신호를 측정 시간(t_measure) 동안 측정할 수 있다. 따라서, 고정된 사전 자화 자기장의 인가 시간에서 사전 자화 자기장의 세기에 따라, 수소핵의 공명 주파수에서의 진폭이 추출될 수 있다. 상기 진폭은 사전 자화 자기장의 세기에 따라 지수 함수 의존성을 가지고, 지수 함수 의존성에 의하여 T1이 추출될 수 있다.
수소 이외에 탄소, 질소 등의 다른 원자핵에 대하여, T1 및/또는 T2는 측정 대상의 품종에 따라 분류될 수 있다.
도 7은 본 발명의 일 실시예에 따른 사전 자화 자기장(Bp)에 따른 T1의 의존성을 설명하는 도면이다.
도 7을 참조하면, 수소 원자핵에 대하여 A 품종과 B 품종은 사전 자화 자기장(Bp)에 따라 서로 다른 T1 의존성을 가질 수 있다. 또한, 탄소 원자핵에 대하여 A 품종과 B 품종은 서로 다른 경향성의 T1 의존성을 가질 수 있다. 또한, 질소 원자핵에 대하여 A 품종과 B 품종은 서로 다른 경향성의 T1 의존성을 가질 수 있다. 따라서, 사전자화 자기장(Bp)에 따른 T1 의존성을 가지고 품종이 구별될 수 있다.
도 8은 본 발명의 다른 실시예에 따른 T2를 측정하는 펄스 시퀀스를 나타내는 도면이다.
도 9는 도 8의 T2를 측정하기 위한 펄스 시퀀스를 설명하는 도면이다.
도 8 및 도 9를 참조하면, 사전 자화 자기장(BP)이 인가되고, 상기 사전 자호 자기장이 꺼지지 전에 검출 자기장(Bm)의 인가 시점에 동시에 양의 경사 자기자계(B_G+)이 측정 대상에 인가된다. 또한, 양의 경사 자기장이 꺼지는 시간과 동시에 음의 경사 자기장(B_G-)가 측정 대상에 인가된다. 따라서, 측정 대상에서 탈위상화(dephasing)가 진행되던 스핀들이 극성이 반대인 음의 경사 자기장(B_G-)에 의해 재초점화된다. 재초점화는 그레디언트 에코(gradient echo) 신호를 발생시킬 수 있다. 양의 경사 자기장(B_G+)의 인가 시간은 t_G+이고, 음의 경사 자기장(B_G-)의 인가 시간은 t_G-이다.
사전 자화 자기장(Bp)이 꺼지기 이전에 검출 자기장(Bm)과 양의 경사자장(B_G+)이 켜지도록 펄스 시퀸스가 설계되었다. 양의 경사 자기장(B_G+)과 사전 자화 자기장(Bp)이 동시에 껴지는 시간(t_x)는 30 ms 이다. 양의 경사자장(B_G+)이 꺼짐과 동시에 음의 경사자장(B_G-)이 켜지도록 펄스열이 설계되었다. 음의 경사 자기장(B_G-)이 켜진 후 최대 전압에 도달하는 데에는 35 ms 가량의 시간 지연이 있었다. 음의 경사 자기장(B_G-)이 켜진 후에 FLL이 측정을 시작하는 시간(t_FLL)은 8 ms 이다.
그레디언트 에코(gradient echo) 신호는 푸리어 변환될 수 있다. 상기 수소 핵의 공명 주파수에서의 피크의 선폭의 역수는 T2를 나타낸다.
극저자기장 핵자기 공명 물체 식별 방법은 그레디언트 에코 신호를 이용하여 상기 측정 대상의 스핀-스핀 완화 시간(T2)을 측정하는 단계(S211), 및 상기 스핀-스핀 완화 시간을 이용하여 상기 측정 대상을 분류하는 단계(244)를 포함할 수 있다.
상기 측정 대상의 스핀-스핀 완화 시간(T2)을 측정하는 단계(S211)는 상기 측정 대상의 개별적 완화시간 특성에 적합한 펄스 시퀀스를 설정하는 단계(S210), 상기 사전 자화 자기장 및 검출 자기장을 인가하는 단계(S212), 양의 경사 자기장 및 음의 경사 자기장을 연속적으로 인가하는 단계(S213), 그레디언트 에코 신호를 측정하는 단계(S214), 상기 그레디언트 에코 신호를 푸리어 변환하여(S216) 소정의 원자핵의 공명 주파수에서의 신호 선폭을 측정하는 단계(S218), 및 소정의 원자핵의 공명 주파수에 따른 선폭 또는 선폭의 역수를 제2 파라미터 벡터로 구성하는 단계를 포함할 수 있다(S242). 상기 제2 파라미터 벡터를 이용하여 품종이 하여 품종이 확인될 수 있다(S244).
도 10은 본 발명의 또 다른 실시예에 따른 물체 식별 방법을 설명하는 도면이다.
도 11은 도 10의 FID 신호를 푸리어 변환하여 얻은 신호의 선폭을 측정하기 위한 펄스 시퀀스를 설명하는 도면이다.
도 12는 도 10의 물체 식별 방법을 이용하여 획득한 가상 데이터이다.
도 10 내지 도 12를 참조하면, 극저자기장 핵자기 공명 물체 식별 방법은 측정 대상에 검출 자기장의 세기를 변경하면서 자유 유도 감쇄 신호를 측정하고 푸리어 변환하여 소저의 원자핵의 공명 주파수에서의 신호 선폭을 측정하는 단계(S311); 및 상기 원자핵의 공명 주파수에서의 신호 선폭을 이용하여 상기 측정 대상을 분류하는 단계(S344)를 포함한다.
측정 대상에 검출 자기장의 세기를 변경하면서 자유 유도 감쇄 신호를 측정하고 푸리어 변환하여 원자핵의 공명 주파수의 선폭을 측정하는 단계(S311)는 상기 측정 대상의 개별적 완화시간 특성에 적합한 펄스 시퀀스를 설정하는 단계(S310), 사전 자화 자기장을 인가하는 단계(S312), 상기 검출 자기장을 인가하는 단계(S313), 자유 유도 감쇄(FID) 신호를 측정하는 단계(S314), 상기 FID 신호를 푸리어 변환하여 소정의 원자핵의 공명 주파수에서의 신호 선폭을 측정하는 단계(S316,S318), 상기 검출 자기장의 세기를 변경하는 단계(S334), 및 위의 단계의 반복적인 측정으로 상기 검출 자기장의 세기에 따른 소정의 원자핵의 공명 주파수에서의 신호 선폭으로 제3 파라미터 벡터를 구성하는 단계(S342)를 포함한다. 상기 제3 파라미터 벡터를 이용하여 상기 측정 대상의 품종이 확인된다(S344).
수소핵의 공명 주파수의 선폭을 측정하기 위한 펄스 시퀸스는 사전 자화 자기장(Bp1) 및 검출 자기장(Bm)을 순차적으로 또는 동시에 인가할 수 있다. 상기 검출 자기장(Bm)은 상기 사전 자화 자기장(Bp1)이 꺼지 전에 인가되는 것이 바람직하다. 상기 사전 자화 자기장의 세기(Bp)는 일정하다. 상기 사전 자화 자기장의 인가 시간(t_P)은 일정하다. 상기 사전 자화 자기장(Bp1)이 꺼진 후 소정의 FLL 구동 시간(t_FLL)이 지난 시간부터 FID 신호를 측정 시간(t_measure) 동안 측정할 수 있다. 검출 자기장(Bm)의 세기는 펄스 시퀀스에 따라 증가 또는 감소하도록 변경될 수 있다. 이에 따라, 검출 자기장(Bm)의 세기에 따른 소정의 원자핵의 공명 주파수에 따른 선폭은 품종 구분에 사용될 수 있다.
도 13은 본 발명의 또 다른 실시예에 따른 극저자기장 핵자기 공명 물체 식별 방법을 설명하는 흐름도이다.
도 13을 참조하면, 극저자기장 핵자기 공명 물체 식별 방법은 복수의 사전 자화 자기장 크기에 대해 각각의 스핀-격자 완화 시간(T1)을 구하여 제1 파라미터 벡터를 구성하는 단계(S410), 그레디언트 에코 신호를 이용하여 소정의 원자핵의 공명 주파수에 따른 선폭 혹은 선폭의 역수(T2)을 제2 파라미터 벡터로 구성하는 단계(S420), 검출 자기장의 세기에 따른 소정의 원자핵의 공명 주파수에서의 신호 선폭의 값들로 제3 파라미터 벡터를 구성하는 단계(S430), 상기 제1 내지 제3 파라미터 벡터들의 성분들을 조합하여 새로운 파라미터 벡터를 구성하는 단계(S440), 및 상기 새로운 파라미터 벡터를 이용하여 측정 대상의 품종을 확인하는 단계(S470)를 포함한다.
본 발명의 변형된 실시예에 따른 극저자기장 핵자기 공명 물체 식별 방법은 복수의 사전 자화 자기장 크기에 대해 각각의 스핀-격자 완화 시간(T1)을 구하여 제1 파라미터 벡터를 구성하는 단계(S410), 그레디언트 에코 신호를 이용하여 소정의 원자핵의 공명 주파수에 따른 선폭 혹은 선폭의 역수(T2)을 제2 파라미터 벡터로 구성하는 단계(S420), 검출 자기장의 세기에 따른 소정의 원자핵의 공명 주파수에서의 신호 선폭의 값들로 제3 파라미터 벡터를 구성하는 단계(S430), 상기 제1 내지 제3 파라미터 벡터들의 성분들을 조합하여 새로운 파라미터 벡터를 구성하는 단계(S440), 및 구별하려는 여러 가지 물질의 반복 측정으로 각각의 물질에 대한 새로운 파라미터 벡터들을 수집하는 단계(S450), 및 수집된 새로운 파라미터 벡터들에 서포트벡터머신(support vector machine;SVM) 또는 주요소분석법(Principal component analysis;PCA)을 적용하여 구분에 효과적인 주요 벡터 성분을 찾아내고, 새로운 파라미터 벡터를 주요 벡터 성분으로 사영시킴으로써 파라미터 벡터의 차원을 줄이는 단계(S460), 및 차원 감소된 파라미터 벡터 성분을 이용하여 고속으로 물질의 품종을 확인하는 단계(S470)를 포함한다.
복수의 사전 자화 자기장 세기에 의해 구해진 복수의 T1들, 복수의 검출 자기장 세기에 의해 구해진 복수의 신호 선폭(혹은 그의 역수)들, 그레디언트 에코에 의해 구해진 T2를 각각의 차원 값으로 조합하여 파라미터 벡터(g)를 구성할 수 있다. 예를 들어, 파라미터 벡터(g)는 다음과 같이 주어질 수 있다. g=(T1@10mT, T1@20mT, T1@30mT, ... ,T2)).
구분하려는 물질군에 대해 반복적으로 측정하여 구성된 파라미터 벡터는 다차원 파라미터 벡터 공간상에 도시(plot)될 수 있다. 이에 따라, 물질군에 따른 파라미터 벡터 공간상에서의 군집이 관찰될 수 있다. 파라미터 벡터 공간상에서의 물질군 구분은 패턴인식연구에서 잘 알려진 서포트 벡터 머신(support vector machine:SVM)이나 인공지능망(artificial neural network; ANN)등의 방법을 통해 수학적으로 구분 기준(criteria)을 사용하여 수행될 수 있다.
파라미터의 수가 많아지면 측정 및 파라미터 공간에서의 구분에 시간이 걸릴 수도 있다. 이를 위해 SVM이나 주요소분석(Principal component analysis)를 통해서 파라미터의 성분 중 중요한 차원만을 남기고 다른 차원을 그 중요차원에 사영시킴으로써, 파라미터 공간의 차원을 축소할 수 있다.
축소된 파라미터 벡터 구성요소를 이용하여 보다 고속으로 물질을 구분하는 것이 가능하다.
이상에서는 본 발명을 특정의 바람직한 실시예에 대하여 도시하고 설명하였으나, 본 발명은 이러한 실시예에 한정되지 않으며, 당해 발명이 속하는 기술분야에서 통상의 지식을 가진 자가 특허청구범위에서 청구하는 본 발명의 기술적 사상을 벗어나지 않는 범위 내에서 실시할 수 있는 다양한 형태의 실시예들을 모두 포함한다.
100: 극저자기장 핵자기 공명 물체 식별 장치
112: 듀어
122: 사전자화 코일
124: 사전자화 코일 구동부
130: 측정 대상
160: 센서부
152: 검출 코일
154: 검출 자기장 전원부

Claims (10)

  1. 측정 대상에 사전 자화 자기장의 세기를 변경하면서 복수의 스핀-격자 완화 시간(T1)을 측정하는 단계; 및
    상기 스핀-격자 완화 시간을 이용하여 상기 측정 대상을 분류하는 단계를 포함하고,
    복수의 상기 스핀-격자 완화 시간은 제1 파라미터 벡터를 구성하고, 상기 제1 파라미터 벡터에 의하여 상기 측정 대상이 분류되는 것을 특징으로 하는 극저자기장 핵자기 공명 물체 식별 방법.
  2. 제1 항에 있어서,
    복수의 스핀-격자 완화 시간(T1)을 측정하는 단계는:
    상기 측정 대상의 개별적 완화시간 특성을 측정하기 위한 펄스 시퀀스를 설정하는 단계;
    상기 사전 자화 자기장 및 검출 자기장을 인가하는 단계;
    자유 유도 감쇄(FID) 신호를 측정하는 단계;
    상기 FID 신호를 푸리어 변환하여 소정의 원자핵의 공명 주파수에서의 신호 크기 또는 선폭을 측정하는 단계;
    상기 사전 자화 자기장의 인가 시간을 변경하는 단계;
    상기 사전 자화 자기장의 인가 시간에 따른 신호 크기 변화를 피팅하여 T1을 구하는 단계; 및
    변경된 복수의 상기 사전 자화 자기장 크기에 대해, 위의 단계들을 반복하여 각각의 T1을 구하여 상기 제1 파라미터 벡터를 구성하는 단계를 포함하고,
    상기 제1 파라미터 벡터를 이용하여 상기 측정 대상의 품종이 확인되는 것을 특징으로 하는 극저자기장 핵자기 공명 물체 식별 방법.
  3. 제 1항에 있어서,
    상기 스핀-격자 완화 시간(T1)을 측정한 후에 그레디언트 에코 신호를 이용하여 상기 측정 대상의 스핀-스핀 완화 시간(T2)을 측정하는 단계; 및
    상기 스핀-스핀 완화 시간을 이용하여 상기 측정 대상을 분류하는 단계를 더 포함하는 것을 특징으로 하는 극저자기장 핵자기 공명 물체 식별 방법.
  4. 제 3항에 있어서,
    상기 측정 대상의 스핀-스핀 완화 시간(T2)을 측정하는 단계는:
    상기 측정 대상의 개별적 완화시간 특성을 측정하기 위한 펄스 시퀀스를 설정하는 단계;
    상기 사전 자화 자기장 및 검출 자기장을 인가하는 단계;
    양의 경사 자기장 및 음의 경사 자기장을 연속적으로 인가하는 단계;
    그레디언트 에코 신호를 측정하는 단계;
    상기 그레디언트 에코 신호를 푸리어 변환하여 소정의 원자핵의 공명 주파수에서의 신호 선폭을 측정하는 단계; 및
    소정의 원자핵의 공명 주파수에 따른 선폭 또는 선폭의 역수를 제2 파라미터 벡터로 구성하는 단계를 포함하고,
    상기 제2 파라미터 벡터를 이용하여 품종이 확인되는 것을 특징으로 하는 극저자기장 핵자기 공명 물체 식별 방법.
  5. 측정 대상에 검출 자기장의 세기를 변경하면서 자유 유도 감쇄 신호를 측정하고 푸리어 변환하여 소정의 원자핵의 공명 주파수에서의 신호 선폭을 측정하는 단계; 및
    상기 원자핵의 공명 주파수에서의 신호 선폭을 이용하여 상기 측정 대상을 분류하는 단계를 포함하고,
    측정 대상에 검출 자기장의 세기를 변경하면서 자유 유도 감쇄 신호를 측정하고 푸리어 변환하여 원자핵의 공명 주파수의 선폭을 측정하는 단계는:
    상기 측정 대상의 개별적 완화시간 특성을 측정하기 위한 펄스 시퀀스를 설정하는 단계;
    사전 자화 자기장을 인가하는 단계;
    상기 검출 자기장을 인가하는 단계;
    자유 유도 감쇄(FID) 신호를 측정하는 단계;
    상기 FID 신호를 푸리어 변환하여 소정의 원자핵의 공명 주파수에서의 신호 선폭을 측정하는 단계;
    상기 검출 자기장의 세기를 변경하는 단계; 및
    위의 단계의 반복적인 측정으로 상기 검출 자기장의 세기에 따른 소정의 원자핵의 공명 주파수에서의 신호 선폭으로 제3 파라미터 벡터를 구성하는 단계를 포함하고,
    상기 제3 파라미터 벡터를 이용하여 상기 측정 대상의 품종이 확인되는 것을 특징으로 하는 극저자기장 핵자기 공명 물체 식별 방법.
  6. 삭제
  7. 복수의 사전 자화 자기장 크기에 대해 각각의 스핀-격자 완화 시간(T1)을 구하여 제1 파라미터 벡터를 구성하는 단계;
    그레디언트 에코 신호를 이용하여 소정의 원자핵의 공명 주파수에 따른 선폭 혹은 선폭의 역수(T2)을 제2 파라미터 벡터로 구성하는 단계;
    검출 자기장의 세기에 따른 소정의 원자핵의 공명 주파수에서의 신호 선폭의 값들로 제3 파라미터 벡터를 구성하는 단계;
    상기 제1 내지 제3 파라미터 벡터들의 성분들을 조합하여 새로운 파라미터 벡터를 구성하는 단계; 및
    상기 새로운 파라미터 벡터를 이용하여 측정 대상의 품종을 확인하는 단계를 포함하는 것을 특징으로 하는 극저자기장 핵자기 공명 물체 식별 방법.
  8. 복수의 사전 자화 자기장 크기에 대해 각각의 스핀-격자 완화 시간(T1)을 구하여 제1 파라미터 벡터를 구성하는 단계;
    그레디언트 에코 신호를 이용하여 소정의 원자핵의 공명 주파수에 따른 선폭 혹은 선폭의 역수(T2)을 제2 파라미터 벡터로 구성하는 단계;
    검출 자기장의 세기에 따른 소정의 원자핵의 공명 주파수에서의 신호 선폭의 값들로 제3 파라미터 벡터를 구성하는 단계;
    상기 제1 내지 제3 파라미터 벡터들의 성분들을 조합하여 새로운 파라미터 벡터를 구성하는 단계;
    구별하려는 여러 가지 물질의 반복 측정으로 각각의 물질에 대한 새로운 파라미터 벡터들을 수집하는 단계;
    수집된 새로운 파라미터 벡터들에 서포트벡터머신(support vector machine;SVM) 또는 주요소분석법(Principal component analysis;PCA)을 적용하여 구분에 효과적인 주요 벡터 성분을 찾아내고, 새로운 파라미터 벡터를 주요 벡터 성분으로 사영시킴으로써 파라미터 벡터의 차원을 줄이는 단계; 및
    차원 감소된 파라미터 벡터 성분을 이용하여 고속으로 물질의 품종을 확인하는 단계를 포함하는 것을 특징으로 하는 극저자기장 핵자기 공명 물체 식별 방법.
  9. 저온 액체 냉매를 담고 있는 듀아;
    측정 대상을 사전 자화시키는 사전자화 코일;
    상기 사전자화 코일에 전류를 단속적으로 인가하고 시간에 따라 그 세기가 변경되는 사전자화 자기장(Bp)을 형성하는 사전자화 코일 구동부;
    상기 사전자화 코일에 의하여 상기 사전자화 자기장(Bp)이 인가되는 측정 대상에 의한 핵자기 공명 신호를 측정하는 센서부;
    측정 대상에 검출 자기장을 인가하는 검출 코일; 및
    상기 검출 코일에 전류를 인가하여 상기 측정 대상에 검출 자기장(readout magnetic field)을 인가하는 검출 자기장 전원부를 포함하고,
    상기 센서부는 상기 측정 대상에 사전 자화 자기장의 세기가 변경됨에 따라 복수의 스핀-격자 완화 시간(T1)을 측정하고, 상기 스핀-격자 완화 시간을 이용하여 제1 파라미터 벡터를 구성하고, 상기 제1 파라미터 벡터에 의하여 상기 측정 대상을 분류하는 것을 특징으로 하는 극저자기장 핵자기 공명 물체 식별 장치.
  10. 제9 항에 있어서,
    상기 측정 대상에 경사 자기장을 인가하는 경사 자기장 코일;
    상기 경사 자기장 코일에 전류를 공급하는 경사 자기장 구동부;
    상기 경사 자기장 구동부에 전원을 제공하는 경사 자기장 전원부; 및
    상기 사전 자화 코일 구동부, 상기 검출 자기장 전원부, 및 검출 자기장 전원부에 펄스 스퀀스를 제공하는 펄스 시퀀스 발생부를 더 포함하는 것을 특징으로 하는 극저자기장 핵자기 공명 물체 식별 장치.
KR1020110109971A 2011-10-26 2011-10-26 극저자기장 핵자기 공명 물체 식별 방법 및 극저자기장 핵자기 공명 물체 식별 장치 KR101287426B1 (ko)

Priority Applications (4)

Application Number Priority Date Filing Date Title
KR1020110109971A KR101287426B1 (ko) 2011-10-26 2011-10-26 극저자기장 핵자기 공명 물체 식별 방법 및 극저자기장 핵자기 공명 물체 식별 장치
PCT/KR2012/008426 WO2013062260A1 (ko) 2011-10-26 2012-10-16 극저자기장 핵자기 공명 물체 식별 방법 및 극저자기장 핵자기 공명 물체 식별 장치
CN201280052725.4A CN103890573B (zh) 2011-10-26 2012-10-16 使用超低磁场核磁共振的物体识别方法和物体识别装置
US14/262,251 US9759793B2 (en) 2011-10-26 2014-04-25 Object discrimination method using ultra-low magnetic field nuclear magnetic resonance and an object discrimination apparatus of the same

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020110109971A KR101287426B1 (ko) 2011-10-26 2011-10-26 극저자기장 핵자기 공명 물체 식별 방법 및 극저자기장 핵자기 공명 물체 식별 장치

Publications (2)

Publication Number Publication Date
KR20130045637A KR20130045637A (ko) 2013-05-06
KR101287426B1 true KR101287426B1 (ko) 2013-07-18

Family

ID=48168043

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020110109971A KR101287426B1 (ko) 2011-10-26 2011-10-26 극저자기장 핵자기 공명 물체 식별 방법 및 극저자기장 핵자기 공명 물체 식별 장치

Country Status (4)

Country Link
US (1) US9759793B2 (ko)
KR (1) KR101287426B1 (ko)
CN (1) CN103890573B (ko)
WO (1) WO2013062260A1 (ko)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8593141B1 (en) 2009-11-24 2013-11-26 Hypres, Inc. Magnetic resonance system and method employing a digital squid
US8970217B1 (en) 2010-04-14 2015-03-03 Hypres, Inc. System and method for noise reduction in magnetic resonance imaging
US9864033B1 (en) * 2013-08-26 2018-01-09 University Of New Brunswick Free induction decay based magnetic resonance imaging methods
KR101520814B1 (ko) * 2013-11-01 2015-05-29 한국표준과학연구원 원편광 펄스를 사용하는 극저자기장 핵자기공명 장치 및 극저자기장 핵자기 공명 측정 방법
CN105137374B (zh) * 2014-06-03 2018-09-25 中国科学院上海微系统与信息技术研究所 一种超高分辨率的磁共振成像方法及装置
JP6674958B2 (ja) * 2014-11-11 2020-04-01 ハイパーファイン リサーチ,インコーポレイテッド 低磁場磁気共鳴のためのパルス・シーケンス
CN105806872A (zh) * 2016-04-29 2016-07-27 大连工业大学 利用低场核磁技术鉴别不同鲟鱼卵的方法
CN106018452A (zh) * 2016-04-29 2016-10-12 大连工业大学 一种基于核磁共振技术的花生品种无损检测方法
CN107102369B (zh) * 2017-06-05 2018-09-28 吉林大学 机载低温超导核磁共振浅层油渗漏探测装置及探测方法
CN107543838B (zh) * 2017-09-20 2019-07-26 北京市食品安全监控和风险评估中心(北京市食品检验所) 一种稀奶油中掺假植脂奶油的核磁共振检测方法
CN107703174B (zh) * 2017-10-30 2020-05-08 中国石油天然气股份有限公司 一种核磁共振流体分析仪及其制备方法
JP2021530709A (ja) * 2018-07-20 2021-11-11 エヌピーエル マネージメント リミテッド 材料応答を検出する方法およびシステム
TW202012951A (zh) 2018-07-31 2020-04-01 美商超精細研究股份有限公司 低場漫射加權成像
GB2576579A (en) * 2018-08-24 2020-02-26 Magnetic Shields Ltd Magnetic shield device
CN110006938B (zh) * 2019-04-12 2022-06-17 上海理工大学 一种基于svm的用于现场快速筛选掺兑橄榄油的方法
US10914802B2 (en) * 2019-06-24 2021-02-09 Synaptive Medical Inc. Magnetic resonance imaging systems and methods
US11510588B2 (en) 2019-11-27 2022-11-29 Hyperfine Operations, Inc. Techniques for noise suppression in an environment of a magnetic resonance imaging system
CN112345815B (zh) * 2020-11-06 2021-06-29 中国科学院精密测量科学与技术创新研究院 测量动力电池组充放电电流分布的阵列式磁场计及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07198636A (ja) * 1993-12-29 1995-08-01 Yamamoto Mfg Co Ltd 果実の品質評価装置
KR20110031723A (ko) * 2009-09-21 2011-03-29 한국표준과학연구원 자기장 생성 장치, 자기장 생성 방법, 저 자기장 핵자기 공명 장치, 및 저 자기장 자기 공명 영상 장치

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10142175A (ja) * 1996-11-15 1998-05-29 Jeol Ltd 核磁気共鳴測定方法
US6268727B1 (en) * 1997-06-24 2001-07-31 Southwest Research Institute Measurement of flow fractions flow velocities and flow rates of a multiphase fluid using ESR sensing
JP2002090321A (ja) * 2000-09-21 2002-03-27 Matsushita Refrig Co Ltd 冷凍食品の物理的損傷を定量化する方法
EA006832B1 (ru) * 2002-02-06 2006-04-28 Те Риджентс Оф Те Юниверсити Оф Калифорния Ядерный магнитный резонанс, измеряемый с помощью сверхпроводящего квантового интерференционного датчика, и формирование изображения с помощью магнитного резонанса при сверхслабых полях
JP2008229048A (ja) * 2007-03-20 2008-10-02 Hamano Life Science Research Foundation 皮膚評価方法
US8390286B2 (en) * 2007-05-18 2013-03-05 Los Alamos National Security, Llc Ultra-low field nuclear magnetic resonance and magnetic resonance imaging to discriminate and identify materials
US7688069B2 (en) * 2007-05-18 2010-03-30 Los Alamos National Security, Llc Ultra-low field nuclear magnetic resonance and magnetic resonance imaging to discriminate and identify materials
WO2011112743A1 (en) * 2010-03-09 2011-09-15 Los Alamos National Security, Llc Ultra-low field nuclear magnetic resonance method to discriminate and identify materials
CN101975788B (zh) * 2010-09-01 2012-02-15 苏州纽迈电子科技有限公司 低场核磁共振鉴别食用油品质的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07198636A (ja) * 1993-12-29 1995-08-01 Yamamoto Mfg Co Ltd 果実の品質評価装置
KR20110031723A (ko) * 2009-09-21 2011-03-29 한국표준과학연구원 자기장 생성 장치, 자기장 생성 방법, 저 자기장 핵자기 공명 장치, 및 저 자기장 자기 공명 영상 장치

Also Published As

Publication number Publication date
US9759793B2 (en) 2017-09-12
CN103890573A (zh) 2014-06-25
US20140232400A1 (en) 2014-08-21
KR20130045637A (ko) 2013-05-06
CN103890573B (zh) 2017-02-15
WO2013062260A1 (ko) 2013-05-02

Similar Documents

Publication Publication Date Title
KR101287426B1 (ko) 극저자기장 핵자기 공명 물체 식별 방법 및 극저자기장 핵자기 공명 물체 식별 장치
US8179135B2 (en) Low field electron paramagnetic resonance imaging with SQUID detection
US6885192B2 (en) SQUID detected NMR and MRI at ultralow fields
JP5552583B2 (ja) 極微弱磁場におけるsquid検出核磁気共鳴およびイメージング
Bielecki et al. Zero‐field NMR and NQR spectrometer
US20150369887A1 (en) Methods and apparatuses related to instrumentation for magnetic relaxometry measurements
Augustine et al. Low field magnetic resonance images of polarized noble gases obtained with a dc superconducting quantum interference device
US9632151B2 (en) Magnetic resonance measuring equipment
Pannetier-Lecoeur et al. Magnetoresistive-superconducting mixed sensors for biomagnetic applications
TWI429935B (zh) 低磁場核磁共振系統
JPH04138130A (ja) 磁気共鳴像作成装置
Ropp et al. Electropermanent magnets for variable-field NMR
KR20150036941A (ko) 하모닉 피크들의 패턴 분석을 이용한 물질 분석 방법 및 장치
Lee et al. Dynamic nuclear polarization in the hyperfine-field-dominant region
Qiu et al. Low-field NMR measurement procedure when SQUID detection is used
Abou-Hamad et al. NMR in pulsed magnetic field
Sinibaldi et al. NMR Detection at 8.9 mT with a GMR based sensor coupled to a superconducting Nb Flux Transformer
FI80797C (fi) Maetningsprocedur foer att fraon nmr-signaler avlaegsna stoerningar genom jaemfoerelse av i olika riktningar roterande rf-faeltkomponenter.
Demachi et al. T 1-Weighted image by ultra-low field SQUID-MRI
Katz et al. Earth field NMR with chemical shift spectral resolution: Theory and proof of concept
EP4279938A1 (en) Apparatus and method for investigating a sample by using nuclear magnetic resonance and a magnetoresistance sensor
WO2006052236A1 (en) Nmr and mri apparatus and method involving a squid magnetometer
Giovannetti et al. Efficiency evaluation of a 13C Magnetic Resonance birdcage coil: Theory and comparison of four methods
US20190082997A1 (en) Ultra-low field nuclear magnetic resonance device
Pannetier-Lecoeur Superconducting-magnetoresistive sensor: Reaching the femtotesla at 77 K

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
AMND Amendment
E601 Decision to refuse application
AMND Amendment
X701 Decision to grant (after re-examination)
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20160630

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20170626

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20180625

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20190604

Year of fee payment: 7