- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我尝试使用 python 自学一些 FFT 基础知识。目前,我正在尝试使用 python 和相关库重现 matlab FFT2 衍射结果。
我尝试模拟的示例可以在这里找到: http://www.mathworks.ch/ch/help/matlab/math/fast-fourier-transform-fft.html#brfb2vw-3
在页面底部您将看到以下 matlab 代码:
D1 = fft2(M);
D2 = fftshift(D1);
imagesc(abs(D2))
axis image
colormap(hot)
title('{\bf Diffraction Pattern}')
这会将圆形孔径图像(如果有区别的话,他们实际上将其构建为阵列)转换为衍射图案。
[圆形光圈]www.mathworks.ch/help/releases/R2013b/matlab/math/aperture.gif
【Matlab计算结果】www.mathworks.ch/help/releases/R2013b/matlab/math/diffraction1.gif
问题是,当我拍摄他们的圆形孔径图像(作为 gif 图形)并应用等效的 scipy FFT2 指令时,我得到了不同的结果。为什么会出现这种差异以及创建衍射图案的正确方法是什么。从衍射图案返回到生成它的图像是否存在类似的问题?只有将光圈构建为像这样的数组时,这才有效吗?
我的python代码如下:
import sys
import numpy as np
import pylab as py
from scipy import misc, fftpack
image = misc.imread(sys.argv[1])
D1 = fftpack.fft2(image)
D2 = fftpack.fftshift(D1)
abs_image = np.abs(D2)
py.imshow(abs_image)
py.show()
我下载并使用了他们的圆形光圈图像作为输入(在裁剪尺度等之后),但结果我在蓝色背景上得到了一个彩色十字线。十字线中间有一个看起来像小点的东西。这个结果是由于我使用下载的带有方形边缘的圆形孔径图像还是由于代码所致?
提前非常感谢
最佳答案
我的 scipy.misc.imread
版本将为您提供一个 W*H*4 RGBA 颜色数组,而不是强度值,并且其上的 FFT 可能会给您不同的结果。
我已将用于生成 M
的 Matlab 示例代码移植到 numpy
中:
import numpy as np
import pylab as py
from scipy import misc, fftpack
n = 2**10
I = np.arange(1, n)
x = I - n / 2
y = n / 2 - I
R = 10
X = x[:, np.newaxis]
Y = y[np.newaxis, :]
M = X**2 + Y**2 < R**2
D1 = fftpack.fft2(M)
D2 = fftpack.fftshift(D1)
abs_image = np.abs(D2)
py.imshow(abs_image)
py.show()
这会产生更好的结果:
关于python - 将 Matlab fft2 衍射示例转换为 Python,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/20517309/
我正在尝试将一个字符串逐个字符地复制到另一个字符串中。目的不是复制整个字符串,而是复制其中的一部分(我稍后会为此做一些条件......) 但我不知道如何使用迭代器。 你能帮帮我吗? std::stri
我想将 void 指针转换为结构引用。 结构的最小示例: #include "Interface.h" class Foo { public: Foo() : mAddress((uint
这有点烦人:我有一个 div,它从窗口的左上角开始过渡,即使它位于文档的其他任何位置。我试过 usign -webkit-transform-origin 但没有成功,也许我用错了。有人可以帮助我吗?
假设,如果将 CSS3 转换/转换/动画分配给 DOM 元素,我是否可以检测到该过程的状态? 我想这样做的原因是因为我正在寻找类似过渡链的东西,例如,在前一个过渡之后运行一个过渡。 最佳答案 我在 h
最近我遇到了“不稳定”屏幕,这很可能是由 CSS 转换引起的。事实上,它只发生在 Chrome 浏览器 上(可能还有 Safari,因为一些人也报告了它)。知道如何让它看起来光滑吗?此外,您可能会注意
我正在开发一个简单的 slider ,它使用 CSS 过渡来为幻灯片设置动画。我用一些基本样式和一些 javascript 创建了一支笔 here .注意:由于 Codepen 使用 Prefixfr
我正在使用以下代码返回 IList: public IList FindCodesByCountry(string country) { var query =
如何设计像这样的操作: 计算 转化 翻译 例如:从“EUR”转换为“CNY”金额“100”。 这是 /convert?from=EUR&to=CNY&amount=100 RESTful 吗? 最佳答
我使用 jquery 组合了一个图像滚动器,如下所示 function rotateImages(whichHolder, start) { var images = $('#' +which
如何使用 CSS (-moz-transform) 更改一个如下所示的 div: 最佳答案 你可以看看Mozilla Developer Center .甚至还有例子。 但是,在我看来,您的具体示例不
我需要帮助我正在尝试在选中和未选中的汉堡菜单上实现动画。我能够为菜单设置动画,但我不知道如何在转换为 0 时为左菜单动画设置动画 &__menu { transform: translateX(
我正在为字典格式之间的转换而苦苦挣扎:我正在尝试将下面的项目数组转换为下面的结果数组。本质上是通过在项目第一个元素中查找重复项,然后仅在第一个参数不同时才将文件添加到结果集中。 var items:[
如果我有两个定义相同的结构,那么在它们之间进行转换的最佳方式是什么? struct A { int i; float f; }; struct B { int i; float f; }; void
我编写了一个 javascript 代码,可以将视口(viewport)从一个链接滑动到另一个链接。基本上一切正常,你怎么能在那里看到http://jsfiddle.net/DruwJ/8/ 我现在的
我需要将文件上传到 meteor ,对其进行一些图像处理(必要时进行图像转换,从图像生成缩略图),然后将其存储在外部图像存储服务器(s3)中。这应该尽可能快。 您对 nodejs 图像处理库有什么建议
刚开始接触KDB+,有一些问题很难从Q for Mortals中得到。 说,这里 http://code.kx.com/wiki/JB:QforMortals2/casting_and_enumera
我在这里的一个项目中使用 JSF 1.2 和 IceFaces 1.8。 我有一个页面,它基本上是一大堆浮点数字段的大编辑网格。这是通过 inputText 实现的页面上的字段指向具有原始值的值对象
ScnMatrix4 是一个 4x4 矩阵。我的问题是什么矩阵行对应于位置(ScnVector3),旋转(ScnVector4),比例(ScnVector3)。第 4 行是空的吗? 编辑: 我玩弄了
恐怕我是 Scala 新手: 我正在尝试根据一些简单的逻辑将 Map 转换为新 Map: val postVals = Map("test" -> "testing1", "test2" -> "te
输入: This is sample 1 This is sample 2 输出: ~COLOR~[Green]This is sample 1~COLOR~[Red]This is sam
我是一名优秀的程序员,十分优秀!