Giáo trình xử lý ảnh y tế Tập 1b P3

Xử lý ảnh có quan hệ mật thiết với nhận thức về ảnh của con người. Nói một cách khác, “thị giác máy dựa trên phép xử lý ảnh bằng sự phân tích của máy. Theo một số tác giả, hai lĩnh vực: xử lý ảnh số và “thị giác máy được liên kết chặt chẽ với nhau.

Hình 2.11 Ví dụ 2.4 Hình 2.12 Ví dụ 2.4 Biểu thức biểu diễn thành h(n1 ,0)  R J ( Rn1 ) 2n1 J1(x) = hàm Bessel loại Vì : h(n1 , n )  RJ1 ( R n12  n22 ) 2 n12  n22 Không phải lúc dễ dàng rút biểu thức phân tích đáp ứng xung Ví dụ minh hoạ việc sử dụng phép tích phân số để thu h(n1 , n2 ) Ví dụ 2.5 Tính đáp ứng xung lọc thơng thấp Butterworth đối xứng vịng trịn có đáp ứng tần số cho H (1 ,  )  1 1 R (1 ,  ) Do2 16    1   ,       , R(1 ,  )  12   22 D0 = 3-dB (Giả thiết Do  0.3 ) Dùng biểu thức cơng thức 2.17 để tính đáp ứng xung với n1  5, ,5, n2  5, ,5 Giải Vì H (1 ,  ) thực đối xứng vịng trịn, cơng thức (2.17) viết đơn giản lại h(n1 , n )  4     H ( ,  ) cos(1 n1   n )d 1d   đơn giản ta viết    h(n1 , n )   cos( 1n1 )    H (1 ,  ) cos( n2 )d d1  0  (2.20) Lời giải toán đưa thơng qua chương trình viết ngơn ngữ C Ngơn ngữ C sử dụng rộng rãi tính linh hoạt C phù hợp cho ứng dụng khác Trước bạn viết chương trình, bạn cần nghiên cứu phương pháp tính tích phân kép Sau ta phát triển thêm qui tắc Simpson để tính tích phân kép Qui tắc tính tích phân Simpson viết sau: xm I  f ( x)dx x0   m 1 m x    f ( x0 )   f ( xi )   f ( xi )  f ( xm )  i 1, 3, i  , 4,   ch½n lỴ    x  ( xm  x ) / m; m số lẻ Với tích phân kép xm I yn   x0 f ( x, y )dxdy y0 Sử dụng qui tắc Simpson ta viết sau:  I  x yn  y0   m 1 m 2    f ( x0 , y )   f ( xi , y )   f ( xi , y )  f ( xm , y) dy i 1, 3, i  2, ,   lỴ ch½n   17 viết y  n I  x  y0  f ( x0 , y )     f ( x1 , y)  [14242 241]   f ( x2 , y )  dy       f ( xm , y )   Lại áp dụng qui tắc Simpson vào công thức ta có I x  y [14242 241]  3 1   4    f ( x0 , y0 ) f ( x0 , y1 ) f ( x0 , y2 ) f ( x0 , yn )  2  f ( x , y ) f ( x , y ) f ( x , y ) f ( x , y )    1 1 n    4     2                f ( xm , y0 ) f ( xm , y1 ) f ( xm , y2 ) f ( xm , yn ) 2     4 1    (2.21)  x  ( x m  x0 ) / m  y  ( y m  y0 ) / n m n số lẻ Nếu H (1 ,  ) thực có tính đối xứng vịng trịn, đáp ứng xung tính từ tích phân kép h( n1 , n2 )  2     H (1 ,  ) cos(1n1 ) cos( n2 ) d1d Dùng công thức viết chương trình C để thu đáp ứng xung từ đáp ứng tần số Dưới chương trình viết Turbo C 2.0 Chương trình 2.1 "FILTERD.C" Tính đáp ứng xung từ đáp ứng tần số /*Program to compute the impulse from the frequency response.*/ 18 /* Chuong trinh dung de tinh dap ung xung cua bo loc 2-D tuan hoan doi xung Dap ung tan so cua bo loc duoc cho boi nguoi dung chuong trinh tinh ham H(float w1, float w2) */ #include #include #include #include #include #define pi 3.1415926 int n1,n2; float f(float,float); float H(float,float); float Simpson2(float(*)(float,float),float,float,float,float,int,int); void main(void) { float f(float,float),xmin,xmax,ymin,ymax; int M,N; int NS,xt,yt; float h[11][11],temp; char file_name[40],ch; FILE *fptr; clrscr(); xmin=ymin=0.0; xmax=ymax=pi; M=N=20; /* So cac diem chia theo phuong x va y cac gia tri phai la chan.*/ gotoxy(1,2); printf (" This program calculates the impulse response from the"); printf("\n frequency response."); printf("\n The impulse response is calculated over a window"); printf("\n centered around the origin."); 19 printf ("\n Enter the number of samples in the impulse response."); printf("\n e.g 7*7 (max 11*11) (number have to be odd.) >"); scanf("%d*%d",&NS,&NS); NS=(NS-1)>>1; /* chieu rong va chieu dai cua mot cua so*/ xt=wherex(); yt=wherey(); gotoxy(70,25); textattr(WHITE+(GREEN