일 | 월 | 화 | 수 | 목 | 금 | 토 |
---|---|---|---|---|---|---|
1 | 2 | 3 | 4 | |||
5 | 6 | 7 | 8 | 9 | 10 | 11 |
12 | 13 | 14 | 15 | 16 | 17 | 18 |
19 | 20 | 21 | 22 | 23 | 24 | 25 |
26 | 27 | 28 | 29 | 30 | 31 |
- DSP
- 미국
- Window Function
- 영상처리
- 이민
- Matlab
- Bessel Funciton
- lowpass filter
- Kaiser Window
- Signal Processing
- 신호처리
- 해외생활
- 됬으면 좋겠네
- 매틀랩
- Issaquah
- 차사기
- covid19
- 킹카운티
- 미국 차 구매
- cudnn
- IT
- Visual Studio
- 이사쿠아
- biosignal
- image processing
- User Control
- King County
- ecg
- c++
- Knob Control
- Today
- Total
목록엔지니어링 (24)
임시 블로그 이름
신호 분석, 필터 설계 등을 하다보면 다양한 Window Function등을 만나게 된다. Cosine, Raised Cosine, Hamming, Hanning, Blackman, Triangular, Gaussian 등등등 어떤 Window Function이 있는지는 위키피디아에 잘 정리되어 있으니 참고하자. (http://en.wikipedia.org/wiki/Window_function) 그러면 도데체 이렇게 많은 윈도우를 왜 쓰는 것일까? 간단하게 1차원 데이터에 대해서 생각해보자. 우리가 가지고있는 데이터를 $*$ s(t) $*$라고 하자.그리고 Window Function을 $*$ w(t) $*$라고 하자. Window Function을 쓰는 방법은 크게 2가지가 있는것 같다. (같다 라고 ..
FIR Filter Design에는 여러가지가 있는데, 그 중에서 가장 직관적이고 간단한 방법이 Frequency Sampling 방식이다. 방법은 이렇다. 1. 원하는 형태의 주파수 도메인의 Magnitude Response $*$H(k)$*$를 만든다.2. 만들어진 주파수 도메인 데이터를 Inverse DFT해서 시간축 필터 $*$h(n)$*$을 생성한다. 말했듯히 매우 간단하고 직관적이다. 간단하게 Low-Pass Filter를 하나 만들어보자. Filter의 조건은 다음과 같다.- Filter의 Order는 $*$N$*$는 30으로 하자. 그러면 Filter의 Tap수 $*$L$*$은 30 + 1 = 31이다.- Cutoff 주파수 $*${ \omega }_{ c }$*$는 0.25 라고 하자. ..
http://en.wikipedia.org/wiki/Kaiser_window Window Function은 종류도 많고, 다양한 경우에 쓰인다. 특히 Apodization등에 써서 주파수 축에서 Side Lobe를 줄이는데 많이 사용 된다.가장 유명한 Window Function은 아마 Hamming, Hanning 등 일텐데, 이러한 Window Function들은 따로 파라미터가 없이 고정되어있는 형태이다.대신에 Kaiser Window는 파라미터로 윈도우의 형태를 조절할 수 있다. 특히 Window-Method Filter Design에서 유용하게 쓰이는데, 이 Kaiser Window를 적용하면, Transition Band나, Ripple등을 조절할 수 있기 때문이다. Kaiser Window는..
http://mathworld.wolfram.com/ModifiedBesselFunctionoftheFirstKind.html Modified Bessel Function of the First Kind는 다음과 같은 수식으로 정의 된다. $*$ { I }_{ n }(z)={ i }^{ -n }{ J }_{ n }(ix) $*$ 여기에서 $*${n}$*$이 실수 $*${v}$*$인 경우엔 아래와 같이 계산할 수 있다. $*${ I }_{ v }(z)={ \left( \cfrac { 1 }{ 2 } z \right) }^{ v }\sum _{ k=0 }^{ \infty }{ \frac { { \left( \cfrac { 1 }{ 2 } z \right) }^{ 2k } }{ k!\Gamma \left( v..
Multimodal Non-Rigid Motion Artifact Correction with Concurrent Ultrasound - Presented at ISMRM2012, Melbourne, Australia. n37jan@gmail.com
Multimodal Non-Rigid Motion Artifact Correction with Concurrent Ultrasound - Presented at ISMRM2012, Melbourne, Australia. n37jan@gmail.com
R-peak 검출아마도 ECG 신호처리라는 키워드로 이 블로그 까지 온 사람들은, 대부분 'R-peak 검출'과 '노이즈 제거'라는 큰 키워드에 대해 궁금해하다가 여기까지 오지 않았을까 한다. 이 포스트에서는 R-peak 검출에 대해서 먼저 이야기 해 보자. 왜?일단 R-peak검출을 하기 전에 왜 R-peak 검출이 필요한 것인가에 대해서 이야기 해 보자. 이미지 출처: http://www.davita-shop.co.uk/ecg-instruments.html 내가 의사가 아니라 아주 자세하게 이야기는 못한다만 (ㅋㅋ)... 위의 그림을 보면 RR-interval, PR, QRS, QT, ST segment등의 용어가 나오는 것을 알 수 있다. 대부분 R-peak을 기준으로 해서 다른 특이점들과의 간격등..
서론 이미지 워핑이라는것은 쉽게 말해, 영상을 이렇게 찌그러트리는 기술이다. [그림 1] (좌) 원본 이미지 256x256 격자 이미지, (우) 워핑을 통해서 찌그러진 이미지 이쪽으로 유명한 책은 Randy Crane의 책이다. 워낙 유명한 책이다보니 한국어 번역판도 있다. C로 짜여진 소스코드를 공개 해서, 어떻게 실제로 구현해야 하는지 잘 설명하고 있다.영문판: A Simplified Approach to Image Processing: Classical and Modern Techniques in C한글판: 영상처리 이론과 실제 또, 이미지 워핑/몰핑의 대가인 Wolberg의 Digital Image Warping이란 책도 매우 유명하다. 다양한 예시와 함께, 실제로 영화나 공학쪽에서 어떠한 상황에..
심전도(ECG) 신호처리 사실 ECG신호에 대한 관심은 Einthoven할아버지가 ECG라는걸 집대성(?)한 이후 꾸준히 있어왔으나, 70년대 ~ 90년대가 절정이 아니었나 싶다. 논문들을 겁색해 보면 그 시절 논문들이 상당히 많다. 최근에는 U-Health, E-Health이런것 때문에 ECG신호에 대한 관심이 다시 매우 높아지고 있다. 사람의 생존상태, 스트레스 상태들을 직접/간접적으로 알려주고, 측정 장비 또한 매우 저가로 또 소형으로 만들 수 있기 때문이다. 한 가지 예를 들면, 모바일 또는 원격지에서 ECG신호를 비롯한 다양한 생체신호를 수집하고, 그것을 전송하고, 의사의 진단을 받고, 그것을 다시 환자에게 피드백 해주는 이런 시스템을 생각해 볼 수 있다. 학회에 가보면 실제로 많은 기업들이 이..
(이미지는 위키피디아 optical flow 설명 페이지에서 퍼왔다. http://en.wikipedia.org/wiki/File:Opticfloweg.png) 영상에서 Motion Estimation을 하는 방법 중 아마 가장 유명한것이 Block-Matching-Algorithm (BMA) 와 이 Optical Flow가 아닐까 한다. BMA는 뭐 그냥 블록을 나누고 두 영상에서 가장 비슷한 블록을 찾아내는 거니까 이해하는데 별 큰 어려움이 없고, 근데 이 Optical Flow는 좀 아리송 하다. 설명을 해야하는데, 귀찮기도 하고 위키피디아에 설명이 잘 되어있으니까 그걸 가지고 하겠다. 아래는 그냥 위키피디아에 있는걸 요점 정리만 해서 쓴거다. 원본은 여기 있다. http://en.wikipedia..