히스토그램 임계 값에서 플로팅

Question

% OTSU THRESHOLDING 
% X=imread('cameraman.tif') 
% a=ones(256,1)*[1:256]; 
% c2=double(X).*(a/2+50)+(1-double(X)).*a/2; 
% c3=uint8(255*mat2gray(c2)); 
% IDX = otsu(c3,2); 
%---------------------------------------------------------- 
% 
% function [IDX,sep] = otsu(I,n) 

I = single(I); 

    %% Convert to 256 levels 
    I = I-min(I(:)); 
    I = round(I/max(I(:))*255); 

    %% Probability distribution 
    unI = sort((I)); 
    nbins = min(length(unI),256); 

    [histo,pixval] = hist(I(:),256); 

    P = histo/sum(histo); 


    %% Zeroth- and first-order cumulative moments 
    w = cumsum(P); 
    mu = cumsum((1:nbins).*P); 

    %% Maximal sigmaB^2 and Segmented image 
    if n==2 
     sigma2B =... 
      (mu(end)*w(1:end-1)-mu(1:end-1)).^2./w(1:end-1)./(1-w(1:end-1)); 

     [maxsig,k] = max(sigma2B); 
     k 

이것은 오츠 임계 값 (2 클래스 만)을 찾는 데 사용하는 코드입니다. 그래서 최적의 임계 값 즉 k를 찾은 후 임계 값 이미지를 계산하는 방법 ?? 각 이미지 픽셀 값을 변경하는 방법 < k = 0 및> = k = 1 그러면 이진 이미지가 생성됩니까?

Answer 1

임계 값 k이 주어지면 이진 연산을 사용하여 임계화된 이미지를 계산할 수 있습니다.

단순히 I이 이미지 인 경우, 역치 영상과 같이 계산 될 수 넣어 :

outThreshold = I >= k; 

또한 사용할 수있는 동일한 작업 할 것입니다 outThreshold = im2bw(I,k); (감사 @Daniel를!). 기본적으로 위에서 설명한 내용을 본질적으로 수행합니다. 어떤 경우

,이 반환 이진지도 곳보다 크거나 미만이다 k 것과 false/0가 설정되는 동안 true/1 설정 얻을 k 같다 강도가 픽셀. 다음은 graythresh을 사용한 예입니다. MATLAB에는 Otsu를 사용하여 graythresh에 의해 내장 된 임계 값 찾기 알고리즘이 있습니다.

I = imread('cameraman.tif'); %// Built-in to MATLAB - A grayscale image 
k = 255*graythresh(im); %// Output is scaled between 0 - 1. 
         %//Multiply by 255 to get 8-bit intensity 
subplot(1,2,1); 
imshow(I); 
outThresh = I >= k; 
subplot(1,2,2); 
imshow(outThresh); 

는, 당신은이 그림을 얻을 : 검은 색 위치가 그렇지 않은 사람들을 동안 임계 값을 능가하는 픽셀 흰색 나타낸다입니다

위치를.

의견에서

---

, 당신은 다중 임계 값을 수행하는 방법을 추구하고자하는 (당신은 정말 그냥 이 게시물을 업데이트해야한다). 이것은 여전히 ​​매우 간단합니다.

thresholds = [32 96 128 192]; 

이것이 말하는 것은 다섯 지역은 우리가 고려할 필요가 있다는 것입니다 : 0 사이의 강도 값

• 픽셀을 31
당신을 가정 할 것은 이렇게 같은 임계 값의 배열을했다
• 명암 값이 32에서 95 사이 인 픽셀
• 96에서 127 사이의 명도 값을 갖는 픽셀
• 픽셀 명도 va 매독 강도 값
191-128 사이의 픽셀이 예를 들어, 상부 단부가 배타적 것을

255-192 사이에 유의해야한다. 따라서 종료 간격은 포함하지 않지만 다음 간격에이 숫자를 포함시키고이를 간격의 하단으로 지정합니다.

for 루프를 사용하고 각 범위 내의 모든 픽셀을 찾아 ID 번호로 설정하는 지도 변수를 만들면됩니다.이 모든 ID를 순서대로 증가 시키십시오. 그러면 첫 번째 영역에 ID가 1, 다음에 ID가 할당되고 ID가 2 등으로 할당됩니다. 이 그것을 할 것입니다 코드입니다 :

map = zeros(size(I)); 

thresholdArray = [0 thresholds 255]; 

for k = 1 : numel(thresholdArray)-1 
    lower = thresholdArray(k); 
    upper = thresholdArray(k+1); 

    map(I >= lower & I < upper) = k; 
end 

%// Rescale to 0 - 255 
outMap = uint8(255*(map - min(map(:)))/(max(map(:)) - min(map(:)))); 

figure; 
subplot(2,2,1); 
imshow(I); 
subplot(2,2,2); 
imhist(I); 
subplot(2,2,3); 
imshow(outMap); 
subplot(2,2,4); 
imhist(outMap); 

이 출력 모습입니다 같은 :

코드는 다중 임계 값,뿐만 아니라 막대 그래프와 원본 이미지를 표시하는 수행 히스토그램이있는 임계 값 이미지뿐입니다. 출력 맵의 크기를 0에서 255로 변경했습니다. 따라서 클래스에는 (256/numThresholds)의 값이 할당됩니다. 이 경우, 64

Answer 2

im2bw(I,p) 위 내 의견을보고 몇 가지 방법을 제안 임계 P