로봇이 되고픈 부엉이

K-space data는 numpy라이브러리에서 제공하는 numpy array의 자료형으로 저장되어 있다. 이 K-space data인 숫자 행렬들로 부터 영상을 보기 위해 다음과 같은 작업을 한다. import numpy as np from numpy import fft import matplotlib.pyplot as plt path = './T2_7_C1M1_1.npy' k_space = np.load(path) unshifted_image = fft.fft2(k_space) image = fft.fftshift(unshifted_image) fig = plt.figure(figsize=(15, 6)) ax1 = plt.subplot(1,3,1); ax1.set_title('K space') plt..

엠알아이 3번째 단원(Magnetic Resonance Hardware) 요약 Magnet 자석 1. 초전도 자석(솔레노이드)는 L-He( 거리 비례 주파수 변화를 줌: 공간정보를 알 수 있음. 2. 거리 비례 주파수로 a. 선형적 라모 frequency를 따르게 한다. b. Z 방향만 c. 코일의 전류에 따라 자기장의 세기가 변한다. 3. 냉각은 수냉식과 공냉식이 있다. 4. 노이즈: 두드리는 소리가 난다.

엠알아이 5번째 단원(Spatial Localization) 요약 Gradient 그래디언트, 경사자장 한계 1. 증폭기의 파워, heating (온도상승은 noise를 만듦) 2. G를 어느정도 올리면 신체 말단부에 저릿한 자극 3. G를 어느정도 올리면 Eddy Currents 와전류 4. 비선형성과 concomitant terms Gradient의 역할 위치정보를 넣어준다. 단, Strength (진폭)를 변화시키는 것이지 방향을 바꾸진 않는다. Gradient의 종류 1. Slice selection (z축) 스핀의 주파수를 변형시킨다. 원하는 위치(center)를 중심으로 음수와 양수로 나뉜다. Gradient를 세게 걸수록 더욱 얇은 band width를 갖게 된다. 결과적으로 특정부분만 se..

엠알아이 6번째 단원(k-space) 요약 k-space를 compose 하면 Image 가 된다. k-space의 특징 1. 중심과 평행한 선. 2. 중앙에서 멀어질수록 간격이 줄어든다. Resolution 해상도 Low-resolution은 큰 pixel이고, High-resolution은 작은 pixel이다. 해상도가 높은 영상을 k공간에 mapping할 경우, 더 많은 시간이 걸린다. 또한 큰 readout gradient를 쏘아야 한다. FOV field of view 시야 Small FOV는 면적이 작은 영상이고, Large FOV는 면적이 큰 영상이다. 따라서 Large FOV의 K-space는 small보다 픽셀수가 많다. #작은 FOV는 aliasing 문제가 발생한다. 작은 FOV를 s..

엠알아이 4번째 단원(Tissue contrast) 요약 Excitation 여기 H의 precession frequency를 가진 RF를 인가했을 때. Relaxation 회복, 이완 RF가 끊어지고 equilibrium으로 돌아갈 때. a. Transverse relaxation 횡축이완(T2) T2 : 횡축자화가 37% 남았을 때 걸리는 시간. Dephasing : 인접한 스핀들과 상대적 위상을 변동시킨다. T2* : T2 간섭과 비균일 B, M에 의한 T2’ 신호를 고려하는 신호 감소. 스핀 간의 Dephasing이 빨라진다. 만약 B가 균일하다면, T2* 와 T2 는 같다. b. Longitudinal relaxation 종축이완(T1) T1 : 종축자화가 67% 회복되었을 때 걸리는 시간. E..

문제 위 식으로 부터, B0를 증가시키면 같은 신호에서 해상도를 크게할 수 있습니다. B0 가 1.5 T 일 때, 다음 문제를 해결하세요. 3차원 공간에서 해상도를 2배로 올린다면, B0는 같은 강도의 신호에서 어떻게 되어야만 하나요? 답 (여러분들이 원하는 직접적인 답) 해상도를 2배로 올리는 것은 Vs(체적)을 2차원에선 체적 Vs를 1/4 Vs 로 3차원에선 1/8 Vs 로 하는 것 입니다. 따라서 위와 같은 식이 성립되게 되고 문제에서 B0 = 1.5 T 라고 했기 때문에 B0* (이것이 우리가 구해야하는 해상도를 좋게 하기 위한 "자기장")은 다음과 같은 식이 되어 구할 수 있습니다. 따라서 4.24 T 가 되어야만 3차원에서 2배의 해상도를 올릴 수 있습니다. 답 (여러분들이 원하지 궁금해하지..

저번 시간까지 우리는 세차 운동에 대해 다루었다 세차운동에서 중요한 점은 전자는 팽이와 같은 회전을 하는 것이 었다. 이번 시간은 그 팽이를 툭하고 건들였을 때 일어나는 현상을 이해해보는 시간이다. MRI 의 신호를 얻어야 한다. 그 이유는 1. 신체 조직마다 전자의 밀도가 다르다. 2. 전자의 밀도가 다르다는 것을 이용하여 신호를 계측할 수 있다. 3. 그 신호를 가지고 다시금 신체 조직을 만들어 낼 수 있다. 이기 때문이다. 그렇기 때문에 동일하게 있는 전자로는 우리가 원하는 정보를 얻어내기 힘들다. 따라서 Radio Frequency 를 우리 몸에 쏜다. 이전시간에 공명에 대해서 우리가 배운 바가 있다. 공명이란 1. 물질마다 고유한 진동수(주파수)를 가지고 있다. 2. 그 진동수와 같은 진동수를 ..

오늘도 연구소에서 알 수 없는 코딩만 주구장창 해대고.... 나중에 누구나 돈벌며 하루하루 살아갈텐데, 그나마 내가 오랬동안 해왔던 걸 사용하여 돈을 번다면 그래도 새로 배우는 일보다 능률이 좋고.. 뭔가 의미있고 하나의 스토리가 되지 않을까 생각해본다. 는 무슨.. 편해서 그런거다. 새로운 일을 배우기 뭐해서.. 솔직히 4년동안 대학생활하고 다른 은행일이라도 해보자.은행에서 날 원한다면, 날 가져요 앙앙 좀 그렇지 않은가... 아무튼 이런 허튼 생각 접어둬버리구다시금 MRI의 기초개념을 정리하는 두번째 시간으로 가보자. MRI의 기초개념들 두번째 1. 전자의 양자수(올게 와버렸다) 양자수는 어려운 게 맞다. 이렇게 생각해보면 접근하기 쉽다. 원자기호를 생각해보자. 다음 그림을 봐바라. 그림 출처 : h..