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

Chia sẻ: Cinny Cinny | Ngày: | Loại File: PDF | Số trang:10

Thêm vào BST

Báo xấu

643
lượt xem 9
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Qua chiến lược chia để trị ta đạt được hiệu suất tính toán trên máy tính của giải thuật 1-D FFT. Thuật toán FFT vector 2-D sau đâylà cùng một chiến lược. Giải thuật DFT 2-D được xen kẽ với những giải thuật DFT 2-D nhỏ hơn, cuối cùng chỉ DFT 2-D của phần tử đơn được tính. Chúng ta sẽ kiểm tra vector FFT.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Giáo trình xử lý ảnh y tế Tập 1a P9

printf("\nPress any key to exit."); getch() ; exit(1); } nsq=filelength(fileno(fptr)); printf("Is this a square image ?"); printf("\n i.e. Is image_length=image width (y or n)? -> "); while(((ch=tolower(getch()))!='y')&&(ch!='n')); putch(ch); switch(ch) { case 'y': image_length=image_width=sqrt(nsq); printf("\n Image size = %d x %d",image_length, image_width); break; case 'n': printf("\nEnter image_width-->"); scanf("%d",&image_width); image_length=nsq/image_width; printf("image length is %d", image_length); break; } printf("\nEnter file name for filtered image ->"); scanf("%s",file_name); ind=access(file_name,0); while(!ind) { gotoxy(1,6); printf("File exists. Wish to overwrite? (y or n)-->"); while(((ch=tolower(getch()))!='y')&&(ch!='n')); putch(ch); switch(ch) { case 'Y': ind=1; break; case 'n': gotoxy(1,6); 70
printf(" "); gotoxy(1,5); printf(" "); gotoxy(1,5); printf("Enter file name -->"); scanf("%s",file_name); ind=access(file_name,0); } } fptr1=fopen(file_name,"wb"); gotoxy(70,25); textattr(WHITE+(GREEN
{ y1=y2=(int)0; for(k1=-1; k1
K irsh Compass operator. */ #include #include #include #include #include #include #include /* Kirsh operators. */ int T[8][3][3]={ { {5,5,5}, {-3,0,-3}, {-3,-3,-3},}, {{5,5,-3}, {5,0,-3}, {-3,-3,-3}, }, {{5,-3,-3}, {5,0,-3}, {5,-3,-3},}, {{-3,-3,-3}, {5,0,-3}, {5,5,-3},}, {{-3,-3,-3}, {-3,0,-3}, {5,5,5},}, {{-3,-3,-3}, {-3,0,5}, {-3,5,5},}, {{-3,5,5}, {-3,0,5}, {-3,-3,5},}, {{-3,5,5}, 73
{-3,0,5}, {-3,-3 ,-3},} }; void main() { int i,j,n1,n2,image_width, image_length,k1,k2,ind; char file_name[14],ch; unsigned char **w; unsigned char *temp,tmp; int y[8] ; float nsq; unsigned int zn2; FILE *fptr, *fptr1; clrscr (); printf("Enter file name for input image ->"); scanf("%s",file_name); if((fptr=fopen(file_name,"rb"))==NULL) { p rintf("%s does not exist.", file_name); p rintf("\nPress any key to exit."); getch (); exit(1); } nsq=filelength(fileno(fptr)); printf("Is this a square image ?"); printf ("\n i.e. Is image_length=image_width (y or n)? ->"); while(((ch=tolower(getch()))!='y')&&(ch!='n')); putch(ch); switch(ch) { case 'y': image_length=image_width=sqrt(nsq); printf("\n Image size = %d x %d",image_length, image_width); break; case 'n': printf("\nEnter image_width-->"); scanf("%d",&image_width); 74
image_length=nsq/image_width; printf("image length is %d", image_length); break; } printf("\nEnter file name for filtered image ->"); scanf("%s",file_name); ind=access(file_name,0); while(!ind) { gotoxy(1,6); printf("File exists. Wish to overwrite? (y or n)-->"); while(((ch=tolower(getch()))!='Y')&&(ch!='n')); putch(ch); switch(ch) { case 'y': ind=1; break; case 'n': gotoxy(1,6); printf (" "); gotoxy(1,5); printf(" "); gotoxy(1,5); printf("Enter file name -->"); scanf("%s",file_name); ind=access(file_name,0); } } fptr1=fopen(file_name,"wb"); gotoxy(70,25); textattr(WHITE+(GREEN
for(i=0;i
for(j=0; j
Hình 5.5 (a) Kết quả thu được khi dùng xử lý Sobel trên hình 4.19. (b) K ết quả của xử lý Kirsh trên hình 4.19. 78
CHƯƠNG 6 BIẾN ĐỔI FOURIER RỜI RẠC 6.1 Chỉ dẫn Trong chương 2,chúng ta đã chứng minh rằng đáp ứng tần số của hệ thống của hệ thống tuyến tính bất biến (LSI ) 2-D được cho bởi:    j (1k1  2 k 2 )   h ( k , k )e (6.1) H (1 ,  2 )  1 2 k1   k 2   Nếu h(k1,k2) chỉ có chỉ tồn tại với k1  0, k  0 và tổng quát được xác định 2 trong miền hữu hạn có kích thước N  N thì N 1 N 1   h(k1, k 2 )e  j ( k  k ) (6.2) H (1 ,  2 )  11 22 k1 0k 2 0 Công thức này chứng tỏ rằng H ( 1 ,  2 ) là tuần hoàn, chu kỳ tuần hoàn là 2. Nếu chúng ta lấy mẫu dưới dạng  1, 2, và miền xác định là (0   1  2) và (0  2  2), N  N mẫu, chúng ta có thể viết: 2 2 và  2  (6.3) 1  n1 n2 N N N 1 N 1  j 2  n1k1  n2 k 2  H (n1 , n 2 )    h(k 1 , k 2 )e vì thế (6.4) N k1  0 k 2  0 Biểu thức (6.4) được gọi là biến đổi Fourier rời rạc 2-D hay còn gọi là DFT. Công thức này được áp dụng vào nhiều ứng dụng như lọc, nén ảnh, phóng đại ảnh. Trong chương này chúng ta sẽ nghiên cứu 2-D DFT và các kỹ thuật tính toán. Đầu tiên, chúng ta sẽ xem xét 1-D DFT, sau đó mở rộng ra cho 2-D. 6.2 Biến đổi Fourier 1-D Biến đổi Fourier 1-D cho tín hiệu thời gian rời rạc f(kT) tính theo công thức : N 1  j 2 nk F (n)   f (kT )e (6.5) N k 0 Công thức này có thể viết lại dưới dạng 75