从功率谱中获取横截面

时间:2016-07-11 15:12:38

标签: matlab fft

我目前正在尝试使用光谱方法来分析地形景观。

当我对景观进行FFT并绘制功率谱时。从功率谱中可以找到景观中结构的方向。

2D功率谱: - 2D power-spectrum

在这个功率谱中,我想做一个横截面。 当峰值幅度方向沿x轴或y轴时,这很容易。 但是对于这个领域(以及其他领域),情况并非如此。

来自另一个区域的横截面 - 沿y轴方向: - Cross-section from another area - orientated along the y-axis

我的问题是我想在1中的山峰上做一个横截面,而我似乎无法弄清楚如何。

如果有人能指出我的解决方案。被困在这里几天了。

修改1  我希望横截面是沿着峰值方向的一条线。

编辑2 改进了第一张图片以显示我想要的横截面

2 个答案:

答案 0 :(得分:0)

我的解决方案是,正如GameOfThrows建议的那样:

  • 在我想要的方向上选择2个(或更多)点
  • 在点上使用最小二乘法来创建线
  • 为插值设置meshgrid。
  • 在新行上使用interp2功能。
  • 为部分
  • 定义适当的轴

在我的最后一个横截面中,我最终有多条线,这样我肯定会达到最大幅度。

我有点回答我的问题,但我一直很忙:)

答案 1 :(得分:0)

您可以使用ginput内置matlab函数存储功率谱的2(x,y)坐标,然后使用此值来分隔要插值的轮廓。