Khôi phục ảnh part 2

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

Thêm vào BST

Báo xấu

68
lượt xem 8
download

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

Đ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.

Chủ đề:

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

Lưu

Nội dung Text: Khôi phục ảnh part 2

/* This program calculates the inverse magnitude response of the extracted blurring function.Theprogram requires the file names of the blurred imageand the edge strength image for the blurred image.To obtain the edge strength image you can use Kirsh.exe or Sobel.exe. */ #include #include #include #include #include #include #define pi 3.1414 void 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;i
for(j=1;j
for(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;i
gi¶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íi ph¬ng ph¸p FFT h¬n, so víi c¸ch tiÕp cËn nµy. C¸c hÖ sè thu ®îc cña bé 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óc x¹ thay ®æi theo ¸p xuÊt khÝ quyÓn, vµ còng cã thÓ lµ do më èng kÝnh trong 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, nhng nã 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Ø ra mé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ïng thñ tôc vµ dïng bé läc FIR 7  7 thiÕt kÕ sö dông cöa sæ Blackmann chó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 216
H×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".` 217
H×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Ët thÓ trªn bÒ mÆt ®Òu chÞu mét t¸c ®éng b»ng nhau cña c¸c vÕt mê. §iÒu nµ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