Dsp Tian

摘要： 这里的变换公式是模拟领域的公式，这里虽然是数字图像，不过也可以这样写，YCbCr采样另外的公式。clear all;close all;clc;img=imread('lena_color.jpg');img=mat2gray(img); %任意区间映射到[0,1];[m n dim]=size(img);imshow(img);%%图像的RGBR=img(:,:,1);G=img(:,:,2);B=img(:,:,3);%%RGB2YUVY=zeros(m,n); %亮度U=zeros(m,n); %彩度V=zeros(m,n); %浓度matrix=[0.299 0.58. 阅读全文

摘要： 色相饱和度亮度和三原色的相互转换function main img=imread('lena_color.jpg'); img=mat2gray(img); %任意区间映射到[0,1]; [m n dim]=size(img); imshow(img); %%图像的RGB R=img(:,:,1); G=img(:,:,2); B=img(:,:,3); %%图像RGB2HSL H=zeros(m,n); %色相角 S=zeros(m,n); %饱和度 L=zeros(m,n); %亮度 for i=... 阅读全文

摘要： 色相饱和度明度和三原色的相互转换clear all;close all;clc;img=imread('lena_color.jpg');img=mat2gray(img); %任意区间映射到[0,1];[m n dim]=size(img);imshow(img);%%图像的RGBR=img(:,:,1);G=img(:,:,2);B=img(:,:,3);%%图像的RGB2HSVH=zeros(m,n); %色相角S=zeros(m,n); %饱和度V=zeros(m,n); %明度for i=1:m for j=1:n r=R(i,j); g... 阅读全文

摘要： 老是搞灰度图像我都不好意思了，所以现在开始看看彩色图像方面的东西。首先就是彩色空间了。clear all;close all;clc;img=imread('lena_color.jpg');img=mat2gray(img); %任意区间映射到[0,1];[m n dim]=size(img);%%图像的RGBR=img(:,:,1);G=img(:,:,2);B=img(:,:,3);%%图像的RGB2CMYC=1-R;M=1-G;Y=1-B;%%图像的CMY2RGBR=1-C;G=1-M;B=1-Y;%%如果正反变换都没错的话，那么图像是不变的img(:,:,1)=R;i 阅读全文

摘要： 直方图匹配或叫做直方图规定化都可以，是把原图像的直方图按照给定的直方图加以映射，使新图像的直方图的分布类似于给定的函数。 总共有以下几步： 1.求给定的函数的累积直方图s。 2.求原图像的累积直方图G。 3.求s中每一个值在G中距离最小的位置index。 4.求原图像每个像素通过index映射到的新 阅读全文

摘要： 双边滤波模板主要有两个模板生成，第一个是高斯模板，第二个是以灰度级的差值作为函数系数生成的模板。然后这两个模板点乘就得到了最终的双边滤波模板。 第一个模板是全局模板，所以只需要生成一次。第二个模板需要对每个像素都计算一次，所以需要放到循环的里面来生成，这很像表面模糊啊。哦，表面模糊就是用了一个截尾滤 阅读全文

摘要： 程序运行时如果太慢了，那么就需要有一些标示来说明程序到底运行到哪里了，为了给一个直观的程序运行进度，就需要用到waitbar了。我也是看人家程序有这个东西，感觉挺有用的，就记下来了。h=waitbar(0,'wait...');for i=1:1000 waitbar(i/1000);endclose(h) 阅读全文

摘要： 过去也写过透视变换，当时算法真是弱爆了，我竟然会通过两次变换。不过那引用的三篇文章都是非常好的文章，直到今天我才看明白。所谓的倾斜校正，一定要有标定点，将一个倾斜的矩形变为不倾斜的。因此可以从原四边形四个点和新矩形四个点得到一个变换矩阵，根据这个矩阵再作用到全局图像就可以了。详细原理在这里，MIT的 阅读全文

摘要： SUSAN算子既可以检测角点也可以检测边缘，不过角点似乎比不过harris，边缘似乎比不过Canny。不过思想还是有点意思的。 主要思想就是：首先做一个和原图像等大的目标图像。然后用一个圆形的模板，用模板去遍历原图像每个像素，把模板内的每个像素都和模板中心像素比较，如果灰度小于一个阈值，那么就对目标 阅读全文

摘要： 一种评价压缩后图像质量的方法，和PSNR类似。具体论文看这里。 代码：function main img=imread('lena.jpg'); [m n]=size(img); imgn=imresize(img,[floor(m/8) floor(n/8)]); imgn=imresize(imgn,[m n]); img=double(img); imgn=double(imgn); tmp=[]; for i=9:8:m for j=9:8:n tmp=[tmp SSIM(img(i-8:i,j-8:j... 阅读全文