zhifangtujunhenghua.zip_图像拉伸_非线性拉伸

clc; clear all; close all; img=imread('1.jpg'); img=rgb2gray(img); subplot(1,2,1); imshow(img,[]);title('Original Picture'); [M,N]=size(img); %Histogram Caculation img=double(img); h=zeros(256,1); for i=1:M, for j=1:N, f=img(i,j); h(f+1)=h(f+1)+1; end end subplot(1,2,2); bar(h);grid;title('Original Histogram'); %Histogram Equalizing hh=h/(M*N); for k=1:256 c(k)=uint8(sum(hh(1:k))*255); end for i=1:M for j=1:N f=img(i,j)+1; histeqimg(i,j)=c(f); end end figure; subplot(1,2,1); imshow(histeqimg,[]);title('Picture after Histogram Equalizing'); %Histogram Caculation of Equalized Picture histeqimg=double(histeqimg); h=zeros(256,1); for i=1:M, for j=1:N, f=histeqimg(i,j); h(f+1)=h(f+1)+1; end end subplot(1,2,2); bar(h);grid;title('Histogram of Equalized Picture'); %Histogram Equalizing by Function Histeq img=imread('1.jpg'); img=rgb2gray(img); img=histeq(img); figure; subplot(1,2,1); imshow(img,[]);title('Comparison:Picture by function histeq'); img=double(img); h=zeros(256,1); for i=1:M, for j=1:N, f=img(i,j); h(f+1)=h(f+1)+1; end end subplot(1,2,2); bar(h);grid;title('Comparison:Histogram by function histeq');