Dsp Tian

摘要： cl;w=0.6; %放大或缩小的宽度h=1.4; %放大或缩小的高度img=imread('Corner.png');imshow(img);[m n]=size(img);imgn=zeros(h*m,w*n);rot=[h 0 0;0 w 0;0 0 1]; %变换矩阵x=h*u,y=w*vinv_rot=inv(rot);for x=1:h*m for y=1:w*n pix=[x y 1]*inv_rot; ... 阅读全文

摘要： cl;w=0.7; %放大或缩小的宽度h=1.5; %放大或缩小的高度img=imread('Corner.png');imshow(img);[m n]=size(img);imgn=zeros(h*m,w*n);rot=[h 0 0;0 w 0;0 0 1]; %变换矩阵x=h*u,y=w*vfor u=1:m for v=1:n pix=[u v 1]*rot; if (pix(1)>=0.5 && pix(2)>=0.5) ... 阅读全文

摘要： cl;bw=[eye(3) eye(3);eye(3) eye(3)]; %要处理的矩阵f=eye(3); %模板矩阵,其实是需要使用到的逻辑lut=makelut(@(nbood) isequal(nbood,f),3); %生成查表向量，将f矩阵这样的逻辑应用到3*3的nbood上bw1=applylut(bw,lut); %应用查表向量。方法：遍历bw所有pixel,如果哪个pixel和f完全相同的,则bw1中相应位置的pixel置1，否则置0bw,f,bw1 ... 阅读全文

摘要： cl;img1=imread('gaosi.jpg');[m,n]=size(img1);w=fspecial('gaussian',[3 3]);img2=imresize(imfilter(img1,w),[m/2 n/2]);img3=imresize(imfilter(img2,w),[m/4 n/4]);img4=imresize(imfilter(img3,w),[m/8 n/8]);img5=imresize(imfilter(img4,w),[m/16 n/16]);imshow(img1);figure,imshow(img2);figure, 阅读全文

摘要： #include <iostream>#include <stdlib.h>using namespace std;typedef struct tree{ int data; struct tree *left; struct tree *right; struct tree *prior;}node;void preorder_tree(node *head){ if (head==NULL) return; cout<<head->data<<""; preorder_tree(head->left); pr 阅读全文

摘要： cl;img=imread('mask.jpg');imshow(img);[x,y]=size(img);img_man=zeros(x,y);img_com=zeros(x,y);%% 直方图均衡化算法Max=max(max(img));Min=min(min(img));Hist=zeros(1,256);for i=1:x for j=1:y Hist(img(i,j)+1)=Hist(img(i,j)+1)+1; endendfigure,plot(Hist);p=zeros(1,256);for i=1:256 p(i)=Hist(i)/(x*y);... 阅读全文

摘要： cl;img=imread('fupeng.jpg');z=3;[x y]=size(img);imshow(img);img=double(img);img_cai=zeros(x,y,z);Max=max(max(img));Min=min(min(img));img=(255/(Max-Min))*img-(255*Min)/(Max-Min);figure,imshow(uint8(img));r=1;g=2;b=3;for i=1:x for j=1:y temp=(2*pi/(Max-Min))*img(i,j)-(2*pi*Min)/(Max-Min); ... 阅读全文

摘要： cl;img=imread('TestDatabase/1.bmp');imshow(img);img=double(img);figure,imshow(mat2gray(img))ma=double(max(max(img)));mi=double(min(min(img)));img=(255/(ma-mi))*img-(255*mi)/(ma-mi);img=uint8(img);figure,imshow(img) 阅读全文

摘要： cl;m=31;n=31;img=zeros(m+1,n+1);img=double(img);pi=3.1415926;sigma=10;for i=-(m/2):m/2 for j=-(n/2):n/2 img(i+m/2+1,j+n/2+1)=(1/(2*pi*sigma*sigma))*exp(-(i*i+j*j)/(2*sigma*sigma)); endendimg=mat2gray(img);imshow(img);imwrite(img,'pic.bmp');%一个更好的实现，2012.12.18修改h=5;w=2;[x y]=meshgrid(-w... 阅读全文

摘要： cl;raw=zeros(200,256,30);for i=1:30 filename=strcat('F:\算法实验\data\seq3\',int2str(i),'.bmp'); raw(:,:,i)=imread(filename);end方法二：cl;img=cell(1,30);for i=1:30 filename=strcat('F:\算法实验\data\seq3\',int2str(i),'.bmp'); img{i}=imread(filename);end 阅读全文