我能正确的表面变形通过减去当地平均水平?
24日视图(30天)
显示旧的评论
嗨!
我试图重建一个振动表面从一些图片。然而,我的图片似乎扭曲了一个设置错误。现在我试图摆脱与下面的脚本设置错误,我做了一个简单的例子:
10×10厘米,表面变形在以下假设:
由于额外的表面高度失真大数量级。例如
% %例子;扭曲的图像= vib +失真
x = 0:0.1:10;
y = 0:0.1:10;
distortedImage = [];
振动= [];
为i = 1:长度(x)
disIm_x = [];
Vib_x = [];
为j = 1:长度(y)
Vib_xy = 0.09 * sin (4 * x (i)) + 0.15 *罪(y (j) + 0.5 * x (i));
disIm_xy y = x(我)* (j) + Vib_xy;
Vib_x = [Vib_x;Vib_xy];
disIm_x = [disIm_x;disIm_xy];
结束
distortedImage = [distortedImage, disIm_x];
振动=(振动、Vib_x);
结束
图(1)
次要情节(1、2、1)
冲浪(x, y, distortedImage);在;
标题(“扭曲的形象”)
次要情节(1、2、2)
冲浪(x, y,振动)
zlim (5 [5])
标题(“隐藏信号”)
知道我隐藏的将不完美的正弦信号,在这个例子中,我知道我不能使用傅里叶变换。
我决定从扭曲的图像减去当地平均使用卷积:
半径= 1;
sqr = 1 + 2 *半径;%的平方大小平均将3 x3
= 1 (sqr) / (sqr * sqr);
S = conv2 (distortedImage,“相同”);% S (x, y) =平均distIm (x, y)和邻国8输入
子= (distortedImage-S);
子=子(1 +半径:长度(x)半径,半径1 +:长度(x)半径);%修剪边缘附近的曲线玲珑
图(2)
次要情节(1、3、1)
显示亮度图像(年代);在;
colorbar
标题(的平均表面)
次要情节(1、3、2)
显示亮度图像(子);在;
colorbar
标题(“distored形象——平均”)
次要情节(1、3、3)
显示亮度图像(振动)
colorbar
标题(“隐蔽的信号”)
图(3)
次要情节(1、2、1)
冲浪(子);在;
zlim ([-0.5 - 0.5])
标题(“distored形象——平均”)
次要情节(1、2、2)
冲浪(振动)
zlim (5 [5])
标题(“隐蔽的信号”)
你知道为什么encovered形象不符合隐藏信号?频率和相对振幅encovered映像中是可见的,但绝对振幅要小得多。有更好的方法吗?注意,扭曲和隐藏信号并不在我的数据。我只有一个“扭曲的形象”。
非常感谢阅读我的问题。我希望它是清晰和简洁。
