gibbs采样 合成图像

clear all; close all; clc %-------------------------------------------------------------------------- % 设置参数 %-------------------------------------------------------------------------- height = 256; % 图像的尺寸 width = 256; num_label = 3; % 标号的个数，标号为{1,2,..., num_label} gray = [55, 115, 225]'; % 灰度级 %gray = [24, 83, 130, 180, 230]'; % 灰度级 intensity = zeros(height, width); % 灰度矩阵 theta = 0.23; % 高斯分布的标准差 gauss_list = zeros(num_label, num_label); guass_list = calc_gauss_list(gray, theta); neighbor = 2; % 1阶邻域 ngb_len = 2 * neighbor; num_iter = 30; % 迭代次数 prior = zeros(num_label, 1); % 先验概率 %-------------------------------------------------------------------------- % Gibbs 采样 %-------------------------------------------------------------------------- tic % 初始化图像 intensity = initial_intensity(height, width, gray); ratio = calc_ratio(intensity, gray) pad_intensity = zeros(height + ngb_len, width + ngb_len); i = 0; j = 0; ngb_label = zeros(ngb_len + 1, ngb_len + 1); for iter = 1:num_iter pad_intensity = padarray(intensity, [neighbor, neighbor], 'symmetric'); for i = 1:height for j = 1:width ngb_label = pad_intensity(i:i+ngb_len, j:j+ngb_len); prior = calc_prior(ngb_label, gray); intensity(i, j) = gibbs_sample(gray, prior, guass_list); end end toc ratio = calc_ratio(intensity, gray) end intensity = uint8(intensity); imshow(intensity);