matlab
文章平均质量分 73
神海之龙
这个作者很懒,什么都没留下…
展开
-
傅里叶光学(二)一维光栅与二维光栅的夫琅禾费衍射
%% 一维情况下 《傅里叶光学》 吕乃光 p104func1=@(x) convn(rect(x),comb_o(x/2),'same').*rect(x/4); x=linspace(-4,4,801);y=func1(x);[out,xx]=fourier1(func1);figure,plot(x,y);%% 二维情况下 《傅里叶光学》 吕乃光 P105-106L=1;原创 2014-06-07 15:20:23 · 15983 阅读 · 0 评论 -
傅里叶光学(三) 狄拉克与高斯函数
%GAUS exp(-pi*x.*x) function.[x,y]=meshgrid(linspace(-4,4,250));r=sqrt(x.^2+y.^2);z1=gaus(r);% gaus functionz2=gaus((r-3)/0.2);%ring delta functionz3=gaus((r-3)/3); subplot(1,3,1),mysurf(x,y,原创 2014-06-07 15:42:33 · 3456 阅读 · 0 评论 -
稀疏矩阵 sparse 函数 matlab
a=fix(rand(10)*10)+1b=fix(rand(10)*10)+1c=full(sparse(a(:),b(:),1,11,11))sum(c(:))原创 2014-05-26 15:34:21 · 1179 阅读 · 0 评论 -
The Haar Transform 哈尔变化
Probably the simplest useful energy compression process is the Haar transform. In 1-dimension, this transforms a 2-element vector (x(1)x(2))T into (y(1)y(2))Tusing:(y(1)y(2))=T(x(1)x(2))(1翻译 2014-05-27 21:07:46 · 2502 阅读 · 0 评论 -
傅里叶加密
理论基础:原创 2014-05-28 11:47:19 · 1741 阅读 · 0 评论 -
真彩色图像小波分解
clc;clear;a=imread('mm.jpg');figure,imshow(a);map=colormap;nbcol=size(map,1);R=double(a(:,:,1)); % 真彩色图像的三个分量转换成double 格式G=double(a(:,:,2));B=double(a(:,:,3));tic;wname='db1'; % 选取的小波t1=[];原创 2014-05-29 16:03:30 · 4654 阅读 · 0 评论 -
用小波对图像分解,和特征分析
path = 'C:\Program Files\MATLAB\R2013a\bin\Original_Images\DIP3E_Original_Images_CH10\';fileExt = '*.tif';files = dir(fullfile(path,fileExt));len = size(files,1);n=0;for i=1:len fileName = st原创 2014-05-29 16:13:16 · 6769 阅读 · 1 评论 -
gabor function
function [gb_R,gb_I]=gabor_fn(lambda,theta,gamma,psi,b) %[gb_R,gb_I]=gabor_fn(lambda,theta,gamma,psi,b)% lambda Wavelength in pixels of Gabor filter to construct% psi phase % theta原创 2014-05-31 21:29:47 · 1240 阅读 · 0 评论 -
两种mask-边缘提取-类拉普拉斯
clc;clear;f=imread('mm.jpg');f=rgb2gray(f);ff=double(f);W1=[1,2,1;2,-12,2;1,2,1];W2=[2,1,2;1,-12,1;2,1,1];f1=conv2(ff,W1,'same');f2=conv2(ff,W2,'same');subplot(2,2,1),imshow(f);title(' ori img原创 2014-05-31 22:24:08 · 2224 阅读 · 0 评论 -
gabor filter 已修改
function [Aim,Eim,Oim]=gabor_filter(img,gamma,theta,lambda,psi,b) % [Aim,Eim,Oim]=gabor_filter3(img,gamma,theta,lambda,psi,b)% lambda Wavelength in pixels of Gabor filter to construct% psi原创 2014-05-31 22:00:09 · 1151 阅读 · 0 评论 -
gabor wavelet 初步代码 已修改
function [gb_R,gb_I]=gabor_wt(lambda,theta,m,gamma,a,psi,b,delta)%[gb_R,gb_I]=gabor_wt(kx,ky,theta,lambda,m,a,psi,delta)%lambda 为函数波长%theta 为gabor 核函数的方向%psi 为相位偏移%sigma 为高斯函数的标准差%gamma Th原创 2014-05-31 22:02:38 · 1482 阅读 · 0 评论 -
傅里叶光学(四)衍射与干涉
optics5Contentscylindershifted cylinderphase-shifted cylinderinterference of two aperturesmoving apertures closer togethereffect of aperture shapeinterference of three aperturesinterferenc翻译 2014-06-07 15:54:48 · 6169 阅读 · 1 评论 -
傅里叶光学 (一)
Functions used in opticsContentscyl function (circular aperture)somb (sombrero) functionring delta functioncomposite function exampleanamorphic functions2D triangle function (separable)2D tr转载 2014-06-01 21:59:08 · 3927 阅读 · 1 评论 -
Convolution
ConvolutionConvolution is the most important and fundamental concept in signal processing and analysis. By using convolution, we can construct the output of system for any arbitrary input signal,转载 2014-05-07 16:37:41 · 1070 阅读 · 0 评论 -
Proof of Separable Convolution 2D 证明 二维卷积的可分性
原文地址:http://www.songho.ca/dsp/convolution/convolution2d_separable.htmlBy the definition of Convolution 2D; Since convolution is commutative (x[n] * y[n] = y[n] * x[n]), swap the order of c翻译 2014-05-07 16:31:26 · 1984 阅读 · 0 评论 -
空间滤波(三)
clc;clear;ori_img=imread('mm.jpg');ori_img=mat2gray(ori_img);r=ori_img(:,:,1);g=ori_img(:,:,2);b=ori_img(:,:,3);w1=double([0,-1,0;-1,5,-1;0,-1,0]); %sharpen 锐化 算子w2=double([0,0,0;-1,1,0;0,0,0原创 2014-05-07 22:37:59 · 661 阅读 · 0 评论 -
matlab 中各种算子集合与应用(一)
clearori_mag=imread('photo3.bmp');r=double(ori_mag(:,:,1));g=double(ori_mag(:,:,2));b=double(ori_mag(:,:,3));figure,imshow(ori_mag,[]);title('original image')cc=cat(3,r,g,b);% w=double([0,-1,0;原创 2014-05-07 20:47:57 · 1916 阅读 · 0 评论 -
空间滤波(二)
%de_noi_img 空间滤波ori_img=imread('photo3.bmp');% figure,imshow(ori_img,[]);title('original image')noi_img=imnoise(ori_img,'salt & pepper',0.002);r=double(noi_img(:,:,1));g=double(noi_img(:,:,2));原创 2014-05-07 21:47:05 · 884 阅读 · 0 评论 -
A Comparison Between Haar Wavelet Transform and Fast Fourier Transform in Analyzing Financial Time S
Research Journal of Applied SciencesYear: 2010 | Volume: 5 | Issue: 5 | Page No.: 352-360DOI: 10.3923/rjasci.2010.352.360 A Comparison Between Haar Wavelet Transform a转载 2014-05-07 23:04:12 · 1611 阅读 · 0 评论 -
空间滤波(四) sobel 算子
ori_img=imread('Fig0342(a)(contact_lens_original).tif');ori_img=mat2gray(ori_img); % 归一化 转化为double sobel_v=double([1,2,1;0,0,0;-1,-2,-1]); %Sobel operator 3x3 convolution kernels% 垂直 Sobel_v kron原创 2014-05-08 11:14:18 · 3583 阅读 · 2 评论 -
频域滤波(一) 傅里叶谱与相位
ori_rec=imread('C:\Program Files\MATLAB\R2013a\bin\Original_Images\DIP3E_Original_Images_CH04\Fig0424(a)(rectangle).tif');ori_rec=mat2gray(ori_rec(257:768,257:768)); % 归一化 转化为double fft_rec=fftshift原创 2014-05-09 15:25:02 · 3366 阅读 · 0 评论 -
空间滤波(五) 综合运用
ori_img=imread('C:\Program Files\MATLAB\R2013a\bin\Original_Images\DIP3E_Original_Images_CH03\Fig0343(a)(skeleton_orig).tif');stretch_img=intrans(ori_img,'stretch',mean2(im2double(ori_img)),0.9);w_原创 2014-05-08 12:10:07 · 1212 阅读 · 0 评论 -
频率滤波(2)
function [f_filt]=freq_filter(f,h)f=mat2gray(f); % 归一化 转化为double h=double(h);pad_h=zeros(size(f)+size(h)-1); % 在频率域滤波 空域的图像和滤波器都要填充,如果不填充会产生缠绕误差,填充后的大小pad_f=pad_h; % 填充后滤波器的大小与填充后图像的大小是相同的%但是填充零原创 2014-05-09 18:40:14 · 927 阅读 · 0 评论 -
haar_kron 用克罗内克积求haar 矩阵
function y=haar_kron(N)% 2^n=N ,n=1,2,3·······y=1;a=[1,1];b=[1,-1];k=1;while(k<N)yy=diag(ones(1,k)); y=[kron(y,a);kron(yy,b)];k=k*2;endfor i=1:Ny(i,1:end)=y(i,1:end)/sqrt(sum(abs(y(i,1:en原创 2014-05-11 12:16:47 · 1430 阅读 · 1 评论 -
相似图片搜索的原理(matlab)
根据Neal Krawetz博士的解释,原理非常简单易懂。我们可以用一个快速算法,就达到基本的效果。这里的关键技术叫做"感知哈希算法"(Perceptual hash algorithm),它的作用是对每张图片生成一个"指纹"(fingerprint)字符串,然后比较不同图片的指纹。结果越接近,就说明图片越相似。下面是一个最简单的实现:第一步,缩小尺寸。将图片缩小到原创 2014-06-07 15:02:24 · 4933 阅读 · 3 评论