这里有一个例子,我思考你在做什么。与加权平均不确定自己在做什么,但我想我做了一件,你可以适应吗?不太清楚你的意思什么阴谋这层的新矩阵但给你。
% % Matlab答案测试脚本
% % noiseonf matlab函数。米从< //www.tatmou.com/matlabcentral/answers/217417-how-fourier-synthesis-is-performed-using-matlab >检索
%使用noiseonf常规的物质
%彼得Kovesi生成虚拟云图片:
% (c) 1996 - 2005年彼得Kovesi版权
%学院计算机科学与软件工程
西澳大利亚大学的%
% http://www.csse.uwa.edu.au/
%下载链接似乎被打破,所以访问matlab上面的答案
%这里是一些虚拟的图像数据生成的,因为我们只是为了演示
%现在可以运行代码!
行= 258;
列= 2570;
系数= 1.7;
我= randn(行、列);
imfft = fft2 (im);%的fft的形象。
imfft = fftshift (imfft);% 0频率转移到中间。
杂志= abs (imfft);%得到大小
阶段= imfft. /杂志;%和阶段
x = 1(行,1)*(列/ 2:(列/ 2 - 1));
y =(行/ 2:(行/ 2 - 1))的*的(1、列);
半径=√x。^ 2 + y ^ 2);%矩阵值包含来自中心的半径。
半径(行/ 2 + 1,列/ 2 + 1)= 1;% . .避免除零。
过滤器= 1. /(半径。^因素);%构造滤波器。
newfft =过滤器。*阶段;
grayImage =实际(ifft2 (fftshift (newfft)));%转化获得最终的噪声图像
grayImage = grayImage-min (grayImage (:));
grayImageN = grayImage. / max (grayImage (:));
%现在虚拟图像设置为值在0和1之间
%可以展示如何改变颜色映射到暗= 1 = 0
图(1)中,clf,
colormap灰色的%不确定这里的需要但是网络实现需要灰度
斧子(1)=次要情节(211);%保存次要情节轴柄colormap分配颜色映射
显示亮度图像(grayImageN) colormap (ax (1),“灰色”),%的原始colormap 0黑暗:1点燃
标题(“原始地图”)
colorbar
提出= colormap;
斧子(2)=情节(212);;
显示亮度图像(grayImageN)、标题(“翻转彩色地图”),colorbar
colormap (ax (2), flipud城市规划机构(cmap))%翻转colormap desination颠倒
%层不需要循环检查每个像素
zPlaneVals = [。1。3。5。7。9);%的层使用
Splane = 0(行、列元素个数(zPlaneVals));%为层预先分配内存堆栈
图(2)中,
为林德= 1:元素个数(zPlaneVals)%层指数
TempBinImg = 0(大小(grayImageN));%临时零形象
%,因为在你的例子只有当原始图像的值设置为1
否则我们已经% <阈值设置为0。
TempBinImg (grayImageN < zPlaneVals(林德))= 1;%替代临时图像索引值<阈值
Splane(:,:,林德)= TempBinImg;%这个存储到堆栈
Sax(林德)=次要情节(元素个数(zPlaneVals), 1,林德);%创建次要情节
显示亮度图像(Splane(:,:,林德)),colormap (Sax(林德)flipud城市规划机构(cmap));%显示图像和翻转colorbar
标题([“< Zplane踹谷”num2str (zPlaneVals(林德))))
WeightSplanes(:,:,林德)= zPlaneVals(林德)* Splane(:,:,林德);%保存一些假体重应用层
结束
WeightedPlane = (WeightSplanes 3) /金额总和(zPlaneVals);%应用图像分层的重量之和除以重量的总和
图,显示亮度图像(WeightedPlane), colormap灰色的colormap (flipud城市规划机构(cmap)), colorbar
