使用 imshow 时出现 xticks 和 yticks 问题

Question

我在修复imshow图的xticks和yticks时遇到问题。 当我将单个图用于imshow时，范围会更好（见下图）。 但是，当我使用在imshow图上插入局部最大值的坐标时，刻度发生了变化（第二个图）。 任何人都可以在绘制局部最大值后帮助我修复xticks和yticks 。

import numpy as np
import matplotlib.pyplot as plt
from skimage.feature import peak_local_max

# Input of a 2D array/ 
data = np.real(data) # data link : https://www.file.io/LvDn/download/WqdRCoW2dUPV

# Find the coordinates of local maxima in the data
coordinates = peak_local_max(data, min_distance=8,  num_peaks=2) #, threshold_abs = 0.0004) #
# Display results

kf=2*2*np.pi/10
resolution = 40
N=9
kx = np.linspace(-kf, kf, N * resolution)  # space vector
ky = np.linspace(-kf, kf, N * resolution)  # space vector


fig, axes = plt.subplots(1, 2, figsize=(13, 8), sharex=True, sharey=True)
ax = axes.ravel()
ax[0].imshow(data, cmap='seismic', vmin=data.min(), vmax=data.max(),
                       extent=[kx.min(), kx.max(), ky.min(), ky.max()],
                       interpolation='nearest', origin='lower')
ax[0].axis('on')
ax[0].set_title('Original')

上面的代码片段生成了这个图：

在这里， imshow 范围可以正常工作，这对我来说是理想的。 当我使用脚本绘制局部最大值时，刻度会更改，如下所示。

ax[1].imshow(data, cmap='seismic')
ax[1].autoscale(False)
ax[1].plot( coordinates[:, 1], coordinates[:, 0], 'k.', marker='o', markersize=12)
ax[1].axis('on')
ax[1].set_title('Peak_local_max')
fig.tight_layout()
plt.show()

[

如何修复“峰值局部最大值”图中的xticks和yticks ？

如果我将范围用于第二个绘图的脚本，则该图如下所示：

Answer 1

很难确定您要的是什么，并且您提供的示例代码不起作用。 例如，像spec这样的变量是未定义的。

但我猜主要问题是您正在混合像素（来自peak_local_max ）和数据坐标。 您可以使用以下方法修复：

extent = [kx.min(), kx.max(), ky.min(), ky.max()]

fig, ax = plt.subplots(1, 2, figsize=(13, 8), sharex=True, sharey=True)

ax[0].set_title('Original')
ax[0].imshow(data, cmap='seismic', extent=extent)

ax[1].set_title('Peak_local_max')
ax[1].imshow(data, cmap='seismic', extent=extent)
ax[1].plot(kx[coordinates[:, 1]], ky[coordinates[:, 0]], 'go', markersize=12)

Matplotlib 使用的范围也应该是外部范围，我不确定您的坐标数组kx和ky实际上是这样的。 这可能会导致分辨率一半的不匹配/偏移。