%作业4：选择长江三角洲区域的遥感图像，通过阈值分割法提取地物信息
R=imread('右洞庭.jpg'); %读入图像
I=rgb2gray(R); %由于阈值分割法主要用于单色图像，所以先将图像转换成灰度图像
hgram=0:255; %制定灰度变化范围
J=histeq(I,hgram); %实现图像直方图规定化
G=roicolor(J,150,255); %根据图像灰度直方图选择150和255这两个阈值作为分割门限
figure,imshow(I) %显示图像灰度图
figure,imshow(J); %显示均衡化后的图像
figure,imhist(I,64) %对原始图像进行灰度直方统计并显示
figure,imhist(J,64) %对均衡化后的图像进行直方图统计并显示
figure,imshow(G); %显示经阈值分割后得到的图像
colorbar %显示颜色条
