Khôi phục ảnh part 2

Điều này dẫn chúng ta quay lại với các giả thiết ban đầu của chúng ta (coi PSF là bất biến khoảng cách), dùng một OTF duy nhất cho tất cả các trường hợp có thể không chấp nhận được trong một số trường hợp. Để khắc phục vấn đề này chúng ta sẽ xem xét giải thuật sau đây : 1. Chia ảnh thành các miền chữ nhật hoặc là vuông không chồng lên nhau.
Nội dung trích xuất từ tài liệu:
Khôi phục ảnh part 2/* This program calculates the inverse magnitude response of theextracted blurring function.Theprogram requires the file names of theblurred imageand the edge strength image for the blurred image.Toobtain the edge strength image you can use Kirsh.exe or Sobel.exe. */#include #include #include #include #include #include #define pi 3.1414void main() { FILE *fptr1,*fptr2,*fptr; char file_name1[12],file_name2[12]; unsigned char *buff; int M,M1,z,i,j,*wv,*wh,x0,y0,sign_max; float sum,A,*cv,*ch,mu1,mu2,pi2; float sigmax,sigmay,sx2,sy2,u,v,du,dv; float H[32][32],G,x; double nsq; char file_name[14],ch1;clrscr();printf (Enter file name for blurred image before obtaining edges-->);scanf(%s,file_name1);fptr1=fopen(file_name1,rb);if(fptr1=NULL) { printf(File does not exist.); exit(1); }printf(Enter file name for edge strength image-->);scanf(%s,file_name2);fptr2=fopen(file_name2,rb);if(fptr2=NULL) { printf(File does not exist.); 210 exit(1); }nsq=filelength(fileno(fptr1));M=sqrt(nsq);M1=M/4;buff=(char *)malloc(M);wv=(int *)malloc(11*sizeof(int));wh=(int *)malloc(11*sizeof(int));ch=(float *)malloc(11*sizeof(float));cv-(float *)malloc(11*sizeof(float));fseek(fptr2,(long)(M*M1),0);for(i-M/4;ifor(j=1;jfor(i=0;i H×nh 10.3 Hµm kh«i phôc ¶nh mê.for(i=0;i v+=dv; }fprintf(fptr,\n);u=-pi;for(i=0;igi¶m mê. §Çu tiªn bé läc ®îc thiÕt kÕ víi ph¬ng ph¸p Shanks vµ sau®ã lµ ph¬ng ph¸p lÆp cña ch¬ng 9. KÕt qu¶ ¸p dông bé läc IIR 2-D víi¶nh ®· nãi ë trªn cã sù c¶i thiÖn ®¸ng chó ý so víi ¶nh gèc. Tuy nhiªn,ph¬ng ph¸p läc th«ng thÊp m« t¶ trong ch¬ng 9 cho kÕt qu¶ gÇn víiph¬ng ph¸p FFT h¬n, so víi c¸ch tiÕp cËn nµy. C¸c hÖ sè thu ®îc cñabé läc IIR trong vÝ dô nµy ®îc ghi trong b¶ng 10.1. B©y giê chóng ta t¨ng cêng ¶nh sao ho¶ thu ®îc qua kÝnh thiªn v¨n®îc cho trªn h×nh 10.5. Nguyªn nh©n cña vÕt mê nµy lµ do hÖ sè khócx¹ thay ®æi theo ¸p xuÊt khÝ quyÓn, vµ còng cã thÓ lµ do më èng kÝnhtrong kho¶ng thêi gian dµi. ¸p dông thñ tôc ®êng biªn trong phÇn nµy,chóng ta thu ®îc ¶nh chØ ra trªn h×nh 10.6. MÆt dï ¶nh nÐt h¬n, nhngnã kh«ng cung cÊp cho ta nhiÒu chi tiÕt míi. Chóng ta cÇn chØ ra r»ng®©y lµ mét ¶nh cña toµn bé sao ho¶ víi ®é ph©n gi¶i chØ cã 256  256®iÓm ¶nh. §©y lµ mét ¶nh díi mÉu, vµ kh«ng cã th«ng tin nµo vÒ ¶nh®îc rót ra ë ®©y mµ l¹i kh«ng rót ra ®îc tõ ¶nh gèc. H×nh 10.7 chØ ramét ¶nh mê lµ kÕt qu¶ cña ¶nh cña vËt n»m ngoµi tiªu cù. ¸p dông cïngthñ tôc vµ dïng bé läc FIR 7  7 thiÕt kÕ sö dông cöa sæ Blackmannchóng ta thu ®îc ¶nh chØ trªn h×nh 10.8. B¶ng 10.1 c¸c hÖ sè cña bé läc kh«i phôc -0.034139 1.000000 -0.014686 - 0.650127 -0.012220 0.007701 -0.014683 - 0.650224 0.991926 0.474917 -0.655115 0.005353 -0.012214 0.007699 -0.655211 0.005351 0.468586 0.005273 216H×nh 10.4 §¸p øng biªn ®é cña bé läc 2-D IIR ®îc xÊp xØ bëi hµm kh«i phôc ¶nh mê. Bµi tËp 10.2 1. MARS.IMG vµ YOSSRA.IMG lµ c¸c tªn file t¬ng øng cña c¸c¶nh trong h×nh 10.5 vµ 10.7. Thùc hiÖn qu¸ tr×nh kh«i phôc. §Ó lµm bµi tËp nµy chóng ta cÇn ch¹y c¸c ch¬ng tr×nh sau theo thø tù: a. KIRSH ®Ó cã ®îc ®êng biªn râ nÐt cña ¶nh. b. DEBLUR ®Ó tÝnh hµm chuyÓn ®æi gi¶m mê. c. IMPULSE ®Ó thu ®îc ®¸p øng kh«ng gian cña hµm gi¶m mê. d. IIRD ®Ó thiÕt kÕ bé läc IIR tõ ®¸p øng kh«ng gian. e. IIRDITER ®Ó t¨ng cêng bé läc ®îc thiÕt kÕ (dïng mét hoÆc hai bíc lÆp). B¹n cÇn ph¶i tù viÕt ch¬ng tr×nh cho b¹n. 2. LÆp l¹i phÇn 1 dïng c¸c bé läc kiÓu FIR. 3. ¸p dông läc trung b×nh ®Ó kh«i phôc ¶nh YOSSRA.IMG. 4. ¸p dông bé läc thèng kª Wallis víi ¶nh YOSSRA.IMG.` 217H×nh 10.5 ¶nh sao ho¶ bÞ mê do ¶nh hëng cña khÝ quyÓn. H×nh 10.6 Kh«i phôc ¶nh h×nh 10.5. 218 H×nh 10.7 ¶nh mê do ngoµi tiªu cù. H×nh 10.8 Kh«i phôc ¶nh h×nh 10.7.10.6 Kh«i phôc l¹i ¶nh qua phÐp xö lý vïng C¸c phÐp gÇn ®óng ë phÇn trªn dùa trªn c¬ së coi r»ng tÊt c¶ c¸c vËtthÓ trªn bÒ mÆt ®Òu chÞu mét t¸c ®éng b»ng nhau cña c¸c vÕt mê. §iÒunµy sÏ ®óng nÕu chØ cã mét ®é s©u nhá trªn ¶nh hoÆc tÊt c¶ c¸c vËt thÓ 219