如何将100×100的meshgrid镜像扩展至(-200,200)范围？

针对你现有100×100、范围(0,200)的MeshGrid，要扩展至(-200,200)并获取(-200,0)区间的z值，分两种场景处理：

一、镜像法（适用于函数对称的场景） #

如果你的z对应的函数满足f(-x,y)=f(x,y)（X轴对称）或f(x,-y)=f(x,y)（Y轴对称），可以直接通过镜像现有z数组快速生成新区域的值，无需重新计算，效率极高。

代码实现 #

import numpy as np

# 原有网格与z数组（假设原Xmin=0，Xmax=200）
xmesh_original = np.linspace(0, 200, 100)
ymesh_original = np.linspace(0, 200, 100)
X_original, Y_original = np.meshgrid(xmesh_original, ymesh_original)
z_original = ...  # 你的100×100原z数组

# 构造扩展后的网格：合并(-200,0)与(0,200)区间，去除重复的0点
xmesh_neg = np.linspace(-200, 0, 100)
xmesh_extended = np.concatenate([xmesh_neg[:-1], xmesh_original])
ymesh_neg = np.linspace(-200, 0, 100)
ymesh_extended = np.concatenate([ymesh_neg[:-1], ymesh_original])
X_extended, Y_extended = np.meshgrid(xmesh_extended, ymesh_extended)

# 生成扩展后的z数组
# 1. 镜像X<0区域：沿列翻转原z数组，对应x从-200到0的对称值
z_x_neg = z_original[:, ::-1]
z_x_extended = np.concatenate([z_x_neg[:, :-1], z_original], axis=1)
# 2. 若需同时镜像Y<0区域：沿行翻转X扩展后的z数组
z_y_neg = z_x_extended[::-1, :]
z_extended = np.concatenate([z_y_neg[:-1, :], z_x_extended], axis=0)

关键说明 #

• [:,::-1]：沿列翻转数组，将原x从0→200的z值映射为x从-200→0的对称值
• [:-1]：去除重复的0点，避免网格中出现两个x=0或y=0的点

二、重新计算法（通用场景，不限函数对称性） #

如果函数无对称性，或需要精确计算新区域的z值，直接基于扩展后的网格重新计算即可。

代码实现 #

import numpy as np

# 保持与原网格一致的步长，生成(-200,200)范围的网格
original_step = (200 - 0) / (100 - 1)  # 原网格步长
xmesh_extended = np.arange(-200, 200 + original_step, original_step)
ymesh_extended = np.arange(-200, 200 + original_step, original_step)

# 或用linspace直接指定点数（与步长一致的话，总点数为199）
# xmesh_extended = np.linspace(-200, 200, 199)
# ymesh_extended = np.linspace(-200, 200, 199)

X_extended, Y_extended = np.meshgrid(xmesh_extended, ymesh_extended)

# 重新计算z值，替换为你的实际函数逻辑
def calculate_z(x, y):
    # 示例：替换成你生成原z的函数
    return np.sin(x/50) + np.cos(y/50)

z_extended = calculate_z(X_extended, Y_extended)

关键说明 #

• 保持步长一致：确保扩展后的网格密度与原网格相同，避免插值误差
• 直接计算：结果最准确，但如果原z的计算逻辑复杂，耗时会比镜像法高

内容的提问来源于stack exchange，提问作者tajurba experiment