임시 블로그 이름

이런 코드는 진짜 짜지 말자 나는 코드를 잘 짜는 방법은 잘 모르겠다. 여기저기 검색 해 보면 좋은 조언들이 많이 나온다. 그 대신에 나는 정말 이런식으로 짜면 안되는 케이스들을 보여주고 싶다. 혼자서 공부하는 코드라면 모르겠지만, 회사에서 직장에서 코딩을 업으로 하면서, 아래 소개된 식으로 코드를 짜는건 정말 여러사람에게 민폐를 끼치는 행동이다. C++기준으로 써본다. 1. 이미 있는 타입에 대한 재정의에.... 대한 재정의...무슨 말이냐, 예를 들면 아래 코드 같이 짠거 말이다. typedef int TData; ... typedef TData DataPoint; ... DataPoint x = 0; 이건 정말 typedef를 남발한 코드다.지금이야 한번에 저 세 줄을 같이 볼 수 있으니 상관 없지..

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을 기준으로 해서 다른 특이점들과의 간격등..

MATLAB의 연산 결과를 자동으로 저장해 보자. 매우 오래 걸리는 연산을 하거나, 매번 수행하는 연산 결과를 알아서 자동으로 따로따로 저장하게 만들고 싶을때 유용하다. 날짜별로 폴더를 만들고, 그 밑에 파일 이름을 시간으로 해서 저장을 해보자. 예시:clear all;close all; s = 0; for j= 1: 1000 s = s + j; % 1에서 1000까지 합을 계산하는 연산end DIRNAME = datestr(now,'yyyymmdd');FILENAME = sprintf('./%s/SAVE_%s.mat',DIRNAME,datestr(now,'HHMMSS')); mkdir(DIRNAME); save(FILENAME); datestr 함수는 날짜와 시간을 문자열로 반환해 주는 함수고, now..

서론 이미지 워핑이라는것은 쉽게 말해, 영상을 이렇게 찌그러트리는 기술이다. [그림 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신호를 비롯한 다양한 생체신호를 수집하고, 그것을 전송하고, 의사의 진단을 받고, 그것을 다시 환자에게 피드백 해주는 이런 시스템을 생각해 볼 수 있다. 학회에 가보면 실제로 많은 기업들이 이..

(위의 그림은 http://sixflow.net/osx86/2646 여기서 가져옴) 맥북프로 2010중반 모델, 최근에 Lion업데이트하고 잘 써오고 있었는데..... 갑자기 앱스토어랑, iCloud계정에 문제가 생겼다 앱스토어에서는 아이디랑 비밀번호를 넣으면 장치 또는 컴퓨터를 확인할 수 없습니다 라는 메시지가 떴다..... 이건 분명 아이디랑 비밀번호는 맞았다는 이야긴데.... 게다가, 앱스토어에서 받은 몇가지 앱들도 계속 아이디랑 비번을 요구하면서 결국 안되고.... iCloud에서는 로그인을 하면 이 Apple ID 는 유효하지만 iCloud 계정이 아닙니다 라고 메시지가 나오면서 iCloud계정 로그인이 안된다..... 내 아이폰에서는 문제가 없었다. 결국 계정 문제라기 보다는 내 컴퓨터가 문제..

유입로그를 보아하니 ㅋㅋㅋㅋ 나도 빨리 끝났으면 좋겠네 ㅠㅠ

(이미지는 위키피디아 optical flow 설명 페이지에서 퍼왔다. http://en.wikipedia.org/wiki/File:Opticfloweg.png) 영상에서 Motion Estimation을 하는 방법 중 아마 가장 유명한것이 Block-Matching-Algorithm (BMA) 와 이 Optical Flow가 아닐까 한다. BMA는 뭐 그냥 블록을 나누고 두 영상에서 가장 비슷한 블록을 찾아내는 거니까 이해하는데 별 큰 어려움이 없고, 근데 이 Optical Flow는 좀 아리송 하다. 설명을 해야하는데, 귀찮기도 하고 위키피디아에 설명이 잘 되어있으니까 그걸 가지고 하겠다. 아래는 그냥 위키피디아에 있는걸 요점 정리만 해서 쓴거다. 원본은 여기 있다. http://en.wikipedia..

Matlab에서 주파수 분석을 할때 일일이 fft를 하고, 다시 fftshift를 하고, 거기다 abs를 취해서 plot을 하기가 귀찮아서 만들어 보았다. 사용법은 간단하다. >> freqshow(signal) >> freqshow(signal, Fs, 'option', 'scale') signal은 주파수 분석을 할 신호이고 Fs는 sampling frequency, 'option'은, 'amp'와 'phase', 'both' 중에 하나를 선택하면 되는데, 각각 amplitude만 보여줄건지, phase만 보여줄건지, 둘 다 보여줄건지 선택하는 것이다. 기본값은 'amp'이다. 'scale'은 'normal'과 'log'가 있다. amplitude spectrum의 스케일을 조정하는 것이다. 첫 줄 처..

블로그 유입 키워드를 보면 matlab lowpass filter가 꽤 있다. 아마도 컴퓨터나 전기전자 공학 학생들인것 같고, 처음에 이런 필터를 어떻게 설계해야 하는지 찾다 찾다 여기 까지 왔으리라... 사실 이런 필터를 설계하는 방법은 여러가지가 있고, 직접 필터 설계하는 방법 및 공식을 코드로 작성할 수 도 있다. 그렇지만 매번 이렇게 코드를 작성하는 것은 너무 귀찮은 일이고 magnitude response, phase response, pole-zero plot등을 보기 위해서는 또 코드를 작성해야 할 것이다., 매틀랩에서는 다행히도 좋은 툴을 제공하고 있다. 커맨드 라인에서 fdatool을 입력해 보자. (단 signal processing toolbox가 설치되어 있어야 한다.) 왼쪽 아래 ..