图像纹理合成

关于纹理合成最经典的论文应该就数EfrosTexture Synthesis by Non-parametric Sampling这篇论文了,引用量近2000。

这里的合成是基于样例的,就是先有一个小的纹理图像,然后合成一个大的。

我个人的理解是,先产生一个大的随机图像,然后对随机图像中每个像素及其邻域像素和小的纹理图像做对比,找出随机图像当前邻域和小的纹理图像最接近的邻域,然后将小的纹理图像最接近的像素赋值给大的随机图像。这样每处理一个随机图像中的像素,就需要遍历小的纹理图像的所有像素,所以运行时间是很慢的。

先看下效果:

小的纹理图像:

合成的纹理图像:

matlab代码如下:

clear all;
close all;
clc;


mask=mat2gray(imread('wen.jpg'));   %小的纹理图像
[m n]=size(mask);


imgn=mat2gray(rand(256,256));   %最终需要的大的纹理图像,现在是随机图像


w=2;            %L邻域长的半径
I=[1 1 1 1 1;
1 1 1 1;
1 0 0 0];  %L邻域,这里是5*3的,当然也可以是7*4或9*5的
for i=1+w:256
    for j=1+w:256-w
        
        mi=inf;
        for p=1+w:m
            for q=1+w:n-w
                tmp=mean2(abs(imgn(i-w:i,j-w:j+w).*I-mask(p-w:p,q-w:q+w).*I));             
                if tmp<mi       %取差别最小的像素
                   mi=tmp;
                   ii=p;
                   jj=q;
                end
            end
        end        
        imgn(i,j)=mask(ii,jj);
        
    end
end


figure;
imshow(imgn,[]);

纹理合成Texture Synthesis)是指基于已有的纹理图像,生成一个与该图像类似但又不完全一样的新纹理图像的过程。在Matlab中,可以使用Image Processing Toolbox中的相关函数来实现纹理合成。 以下是一个简单的纹理合成Matlab例程: ```matlab % 读取纹理图像 texture = imread('texture.jpg'); % 设置合成后的图像大小 outSize = [512 512]; % 设置纹理块的大小 blockSize = [64 64]; % 计算纹理块的数量 numBlocks = ceil(outSize./blockSize); % 初始化合成后的图像 out = zeros(outSize, 'uint8'); % 随机选择一个纹理块 randBlock = randi([1 size(texture,1)-blockSize(1)], [1 2]); % 将选择的纹理块复制到合成后的图像中 out(1:blockSize(1), 1:blockSize(2)) = texture(randBlock(1):(randBlock(1)+blockSize(1)-1), randBlock(2):(randBlock(2)+blockSize(2)-1)); % 逐块合成 for i = 1:numBlocks(1) for j = 1:numBlocks(2) % 如果不是第一个纹理块,则计算此块与上一块的匹配度 if i>1 || j>1 % 计算上一块的特征向量 lastBlockFeatures = getFeatures(out((1:blockSize(1))+(i-2)*blockSize(1), (1:blockSize(2))+(j-2)*blockSize(2))); % 在纹理图像中寻找最匹配的块 bestBlock = findBestBlock(texture, lastBlockFeatures, blockSize); % 将最匹配的块复制到合成后的图像中 out((1:blockSize(1))+(i-1)*blockSize(1), (1:blockSize(2))+(j-1)*blockSize(2)) = texture(bestBlock(1):(bestBlock(1)+blockSize(1)-1), bestBlock(2):(bestBlock(2)+blockSize(2)-1)); end end end % 显示合成后的图像 figure; imshow(out); % 特征提取函数 function features = getFeatures(img) % 将RGB图像转换为灰度图像 img = rgb2gray(img); % 计算梯度幅值和方向 [Gmag, Gdir] = imgradient(img); % 将方向量化为8个方向 Gdir = floor((Gdir+180)/45)+1; Gdir(Gdir==9) = 1; % 计算每个方向的平均梯度幅值 features = zeros(1,8); for i = 1:8 features(i) = mean(Gmag(Gdir==i)); end end % 查找最匹配的块函数 function bestBlock = findBestBlock(img, features, blockSize) % 初始化最匹配块的位置和匹配度 bestBlock = [1 1]; bestMatch = inf; % 逐个块计算匹配度 for i = 1:(size(img,1)-blockSize(1)) for j = 1:(size(img,2)-blockSize(2)) % 计算当前块的特征向量 curFeatures = getFeatures(img(i:(i+blockSize(1)-1), j:(j+blockSize(2)-1), :)); % 计算当前块与目标块的匹配度 curMatch = sum((features-curFeatures).^2); % 如果匹配度更好,则更新最匹配块的位置和匹配度 if curMatch<bestMatch bestMatch = curMatch; bestBlock = [i j]; end end end end ``` 这个例程中使用了两个函数`getFeatures`和`findBestBlock`,它们分别用于提取纹理块的特征向量和查找最匹配的块。这些函数的实现比较简单,可以根据需要进行修改和优化。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值