NUST COLLEGE OF ELECTRICAL AND MECHANICAL ENGINEERING
NUST COLLEGE OF ELECTRICAL AND MECHANICAL ENGINEERING

DIGITAL IMAGE PROCESSING
LAB REPORT 4 Submitted By:
NS Fatima Hassan
DE-34-CE-B-209
Submitted to:
Sir Saqib Ejaz

TASK1:
Code:
img=imread('fig01.jpg'); figure(1) imshow(img,[]); img=rgb2gray(img); rotatedim=imrotate(img,45,'bilinear'); threshold=graythresh(rotatedim);
BW = im2bw(rotatedim,threshold); figure(2) imshow(BW,[]) Output:

Figure [ 1 ]: Original gray scale image

Figure [ 2 ]: Rotated image after thresholding

TASK2:
Code:
img=imread('fig02.jpg'); figure(1) imshow(img,[]); img=rgb2gray(img) doubimg=im2double(img); %converting image to double figure(2) imshow(doubimg,[]); logtransform=log(doubimg); %log transform figure(3) imshow(logtransform,[]);

Output:

Figure 3: Original image

Figure 4: Image after log transformation

Figure 5: Image after Log Transform
TASK3:
Code: img=imread('fig03.jpg'); figure(1) imshow(img,[]); img=rgb2gray(img); doubimg=im2double(img); %power transformtaions powimg1=1*(doubimg.^5); figure(2) imshow(powimg1,[]); powimg2=1*(doubimg.^10); figure(3) imshow(powimg2,[]); powimg3=1*(doubimg.^2); figure(4) imshow(powimg3,[]); powimg4=1*(doubimg.^0.5); figure(5) imshow(powimg4,[]); powimg5=1*(doubimg.^0.1); figure(6) imshow(powimg5,[]); %log transformation logtransform=log10(doubimg); figure(7) imshow(logtransform,[]); %inverse log transformation inv=(10).^logtransform; figure(8) imshow(inv,[]);

Output:

Figure 5: Original image

Power Law Transformation

Figure 6a

Figure 7b

Figure 7c

Figure 7d

Figure 7e

Figure 7: Log Transformation

Figure 8: Inverse Log Transform
TASK4:
Code: img=imread('fig04b.jpg'); figure(1) imshow(img,[]); img=rgb2gray(img); doubimg=double(img); s=zeros(size(img,1),size(img,2)); for i=1:size(doubimg,1) for j=1:size(doubimg,2) if ((doubimg(i,j)<96 )||(doubimg(i,j)==96)) s(i,j)=doubimg(i,j)*(32/96); elseif (doubimg(i,j)>96 && (doubimg(i,j)<160 || doubimg(i,j)==160)) s(i,j)=((doubimg(i,j)-96)*((224-32)/(160-96)))+32; else s(i,j)=((doubimg(i,j)-160)*((225-254)/(255-160)))+224; end end end figure(2) imshow(s,[]); Output:

Figure [ 9 ]: Original Image

Figure [ 10 ]: High Contrast image
TASK5:
Code: img=imread('fig05.jpg'); figure(1) imshow(img,[]); img=rgb2gray(img); doubimg=double(img); for x=1:8 bp=bitget(img,x); figure(x+1) imshow(bp,[]); end Output:

Figure [ 11 ]: Bit planes 1 to 6

Figure [ 12 ]:Bit planes 7,8
TASK6:
Code: img=imread('fig03.jpg'); img1=imread('fig01.jpg'); figure(1) imshow(img,[]); img=rgb2gray(img); doubimg=double(img); g1=imadjust(img,[0.3 0.8],[0 1]); figure(2) imshow(g,[]);

g2=imadjust(img,[0.3 0.8],[0 1],10); figure(3) imshow(g2,[]); g3=imadjust(img1,[0 1],[1 0]); %inverting figure(4) imshow(g3,[]); g4=imadjust(img1,[],[],10); %fig4 figure(5) imshow(g4,[]); g5=imadjust(img,[],[],10); %fig3 figure(6) imshow(g5,[]);

Output:

Figure [ 13 ]a :Outputs Task 6

Figure 13b :Outputs Task 6