如何獲得輪廓 plot 的 Isocontour 的 xy 坐標？

Question

我正在嘗試從 3D plot 獲取等值面的 xy 坐標。這是我的嘗試；

import matplotlib.pyplot as plt

from numpy import pi, cos, sin, linspace, meshgrid

x = linspace(0,50,1000)
y = linspace(0,50,1000)

n = 5
L = 50
t = 0

def gyroid(x, y, n, L, t):
    tanım1 = (sin(2*pi*n*x/L) * cos(2*pi*n*y/L) + sin(2*pi*n*y/L) + cos(2*pi*n*x/L))
    return tanım1*tanım1 - t**2

XX, YY = meshgrid(x, y)
z = gyroid(XX, YY, n, L, t)

thickness = 0.1
contour = plt.contour(XX, YY, z,levels=[thickness])
# Attempt to get x-y coordinates
dat0= contour.allsegs[0][0]
plt.plot(dat0[:,0],dat0[:,1])

陀螺儀 function 通常看起來像； 3D plot

我得到 z = 0.1 平面的等值線； 無效 plot

我需要 xy 對這些空隙。 但是當我嘗試時，代碼只會得到左下角的坐標。 顯然 function 是強非線性的，但是有沒有辦法檢索這些坐標？

提前感謝您的回復。

Answer 1

您指定contour.allsegs[0][0]以便獲得第一條輪廓線的第一行。

for lines in contour.allsegs:
    for line in lines:
        X, Y = line[:,0], line[:,1]
        plt.plot(X, Y)