仗剑走天涯

简要说明:推荐学习暗原色去雾的原理两篇：第一篇为pdf格式 900k可以下载：https://zkxb.jsu.edu.cn/CN/article/downloadArticleFile.do?attachType=PDF&id=2170第二篇为：https://blog.youkuaiyun.com/yangdashi888/article/details/54571669本文代码参考上面两篇文章，码出来的。效果不错。figure,f=imread("小树林.png");%f=imresiz

       图像的频率是表征图像中灰度变化剧烈程度的指标，是灰度在平面空间上的梯度。如：大面积的沙漠在图像中是一片灰度变化缓慢的区域，对应的频率值很低；而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域，对应的频率值较高。傅里叶变换在实际中有非常明显的物理意义，设f是一个能量有限的模拟信号，则其傅里叶变换就表示f的频谱。从纯粹的数学意义上看，...

原文：http://blog.sina.com.cn/s/blog_818aa1b90102vgor.html为什么matlab中的fft运算后想要得到信号的幅值需要对运算结果的幅值进行乘2除N对于一个时域信号x(t)，采样频率为Fs,信号频率F,采样点数为N，进行N采样后为y[n],则fft(y,N)计算所得的结果就是一个为N点的复数a+bi.每一个点就对应着一个频率点.这个点的模值,就是...

简介：搞不清surf函数和plot3函数当个日志1.画一个曲面使用函数 surfclc,clear,m=1:1:200;n=1:200;sigma = 200;f = fspecial('gaussian', [length(m),length(n)], sigma); %创建预定义的二维滤波器figure,surf(m,n,f); %m为横坐标，n为纵坐标。f点组成面如果上...

1.语法:imshow(I)imshow(I,[low high])imshow(I,[])imshow(RGB)imshow(BW)imshow(X,map)imshow(filename)imshow(___,Name,Value)himage = imshow(___)2.说明:imshow(I) 在图窗中显示灰度图像 I。imshow 优化图窗、坐标区和图像对象属性以...

fspecial创建预定义的二维滤波器语法h = fspecial(type)h = fspecial('average',hsize)h = fspecial('disk',radius)h = fspecial('gaussian',hsize,sigma)h = fspecial('laplacian',alpha)h = fspecial('log'...

网上关于傅里叶变换的解释特别多，但大部分都比较偏理论，导致我看来N多教程也还是懵懵懂懂。在某本书（信号完整性分析？？？）中看到一句震耳发聩的话：”每个工程师都应该亲自动手计算一遍傅里叶变换”。我知道很多工具可以直接给出傅里叶变换结果，但不清不楚一直是我心里过不去的坎，今天终于把坎踏平了，已经忘记了这是第几次冲锋。。。20200324本文主要对一个连续周期信号进行采样所得的有限离散周期...

注：将一个带有噪音的时域信号转化为频域信号clc,clear;figure,Hz=0.02;%周期；也可以说x轴的步长Fs=1/Hz; %频率%定义函数x=0:Hz:10-1/50; y=5*sin(2*pi*15*x) + 10*sin(2*pi*20*x)+ 17*cos(2*pi*5*x);%此处可以添加噪音y=y+10*randn(size(x));%上面的参数定义好，...

1.语法X = randnX = randn(n)X = randn(sz1,…,szN)X = randn(sz)X = randn(___,typename)X = randn(___,‘like’,p)2.说明X = randn 返回一个从标准正态分布中得到的随机标量。X = randn(n) 返回由正态分布的随机数组成的 n×n 矩阵。X = randn(sz1,…,s...

1.语法：y = linspace(x1,x2)y = linspace(x1,x2,n)2.说明：y = linspace(x1,x2) 返回包含 x1 和 x2 之间的 100 个等间距点的行向量。y = linspace(x1,x2,n) 生成 n 个点。这些点的间距为 (x2-x1)/(n-1)。linspace 类似于冒号运算符“:”，但可以直接控制点数并始终包括端点。“l...

1.语法：Y = fftn(X)Y = fftn(X,sz)2.说明：Y = fftn(X) 使用快速傅里叶变换算法返回 N 维数组的多维傅里叶变换。N 维变换等于沿 X 的每个维度计算一维变换。输出 Y 的大小与 X 相同。Y = fftn(X,sz) 将在进行变换之前根据向量 sz 的元素截断 X 或用尾随零填充 X。sz 的每个元素定义对应变换维度的长度。例如，如果 X 是一个 ...

1.语法：        X = ifftshift(Y)        X = ifftshift(Y,dim)2.说明：1.X = ifftshift(Y) 将进行过零频平移的傅里叶变换 Y 重新排列回原始变换输出的样子。换言之...

1.语法：        Y = fftshift(X)        Y = fftshift(X,dim)2.说明：Y = fftshift(X) 通过将零频分量移动到数组中心，重新排列傅里叶变换 X。  ...

1.语法：        Y = fft2(X)        Y = fft2(X,m,n)2.说明：        Y = fft2(X) 使用快速...

        冈萨雷斯版<图像处理>里面的解释非常形象：一个恰当的比喻是将傅里叶变换比作一个玻璃棱镜。棱镜是可以将光分解为不同颜色的物理仪器，每个成分的颜色由波长（或频率）来决定。 傅里叶变换可以看作是数学上的棱镜，将函数基于频率分解为不同的成分。当我们考虑光时,讨论它的光谱或频率谱。同样, 傅立叶变换使...

应用：二维离散傅里叶变换是将图像从空间域转至频域，在图像增强、图像去噪、图像边缘检测、图像特征提取、图像压缩等等应用中都起着极其重要的作用。**理论基础：**任意函数都可以表示成正弦函数的线性组合的形式。二维离散傅里叶变换公式如下：式中 f(x,y) 代表一幅大小为 M x N 的矩阵，其中 x = 0,1,2,···,M-1 和 y = 0,1,2,···,N-1 ，F(u,v) 表示 ...

1.三角函数和双曲函数2.指数函数注：        pow2(x)=x^2               pow2(x,y)=x^y3.复数函数4.圆整函数...

高斯滤波器是空间滤波器的一种，学习高斯滤波器，需要空间滤波的知识作为先验知识（各种数字图像处理的书中都有介绍，冈萨雷斯的经典教材说的还是比较清楚的，建议看英文版，个人感觉比翻译过来的更容易理解）。高斯滤波器是一种平滑线性滤波器，使用高斯滤波器对图像进行滤波，其效果是降低图像灰度的“尖锐”变化，也就是使图像“模糊”了。高斯滤波对于抑制服从正态分布的噪声效果非常好，其代价是使图像...

前言在野外数据采集中，虽然单个仪器采集的是一维信号，但是当把多台仪器数据汇总并生成做二维剖面的图像时，噪声可不只有一维的，更有x,y两个方差同时存在的"二维噪声"！我们已经知道一维噪声可以用一维傅里叶变换到频域滤波，同理二维噪声也可以用二维傅里叶变换到"频率滤波"。二维傅里叶正变换的原理笔者很讨厌一上来就看到一连串复杂的公式！因此当我看懂一个原理后，我就会用最好理解的方式来重述它，毕竟我更偏...