10

背景

此处的文档示例中,可以使用代码片段轻松生成以下等高线图。

import matplotlib
import numpy as np
import matplotlib.cm as cm
import matplotlib.mlab as mlab
import matplotlib.pyplot as plt

matplotlib.rcParams['xtick.direction'] = 'out'
matplotlib.rcParams['ytick.direction'] = 'out'

delta = 0.025
x = np.arange(-3.0, 3.0, delta)
y = np.arange(-2.0, 2.0, delta)
X, Y = np.meshgrid(x, y)
Z1 = mlab.bivariate_normal(X, Y, 1.0, 1.0, 0.0, 0.0)
Z2 = mlab.bivariate_normal(X, Y, 1.5, 0.5, 1, 1)
# difference of Gaussians
Z = 10.0 * (Z2 - Z1)

# Create a simple contour plot with labels using default colors.  The
# inline argument to clabel will control whether the labels are draw
# over the line segments of the contour, removing the lines beneath
# the label
plt.figure()
CS = plt.contour(X, Y, Z)
plt.clabel(CS, inline=1, fontsize=10)
plt.title('Simplest default with labels')

在此处输入图像描述

我的目标

我已经获得了我的等高线图,同时获得了matplotlib.contour.QuadContourSet实例CS。在示例代码段中,CS仅用于clabel(). 但是对于我的情况,我需要获得等高线的方程或坐标集以进行进一步计算。

如何从实例中提取轮廓线的坐标CS?或者我怎样才能以其他方式实现它?

我敢打赌,一定有办法做到这一点。否则,轮廓的东西就只是一个“可视化的花瓶”。

4

1 回答 1

20

您可以从CS.allsegs列表中获取轮廓的坐标。

尝试:

dat0= CS.allsegs[0][0]
plt.plot(dat0[:,0],dat0[:,1])

绘制第一个 (-1) 等高线级别。

于 2013-10-17T06:40:44.680 回答