如何在matplotlib中在三维曲面上绘制直线

Question

我有一个三维绘图，想在绘图表面画几条线。我不清楚该如何组织这些线的数据，使之落到表面。在

以下代码的一些解释：我对描述Rubisco酶（光合作用中至关重要）活性的温度敏感参数进行了敏感性分析。活化能Ha是这个方程中唯一的参数。在

函数plot_TemperatureEffectOnRuBisCOKinetics绘制三维绘图。现在，我想在表面上看到函数TemperatureEffectOnRuBisCOKinetics中描述的4个参数中的每一个的行，每一行都被很好地标记。在

对于如何构造这些行的数据，我们将不胜感激！在

import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np

cRefTmp_C           = 25.                # [C]
cRefTmp_K           = cRefTmp_C + 273.15 # [K]
MolarGasConstant    = 8.314472           # [J mol-1 K-1]

def TemperatureEffectOnRuBisCOKinetics(Ha, LeafTemperature_C):
    """
        multiplier for temperature effects on Kc, K0, Ri and GammaStar [ - ]
        formula thesis Manfred Forstreuter p 66 (eq 2.41) 

        Parameter            ParameterValue
        cHaOfGammaStar   29000  
        cHaOfK0      35900
        cHaOfKc      59500
        cHaOfRi      46390   

        refs for equation:
        Harley P.C., Thomas R.B., Reynolds J.F., Strain B.R., 1992. 
            Modelling photosynthesis of cotton grown in elevated CO2. Plant, Cell Environ. 15: 271-282.
        Farquhar G.D., von Caemmerer S. & Berry J.A., 1980. 
            A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta 149: 78-90.

    """
    LeafTemperature_K = LeafTemperature_C + 273.15               # from Celsius to Kelvin
    return exp(Ha * (LeafTemperature_K - cRefTmp_K) / (MolarGasConstant * LeafTemperature_K * cRefTmp_K))


def plot_TemperatureEffectOnRuBisCOKinetics():
    Ha     = np.arange(25000., 60000., 1000.)
    T      = np.arange(0., 30., 1)
    Ha,T   = np.meshgrid(Ha,T)
    TEff   = TemperatureEffectOnRuBisCOKinetics(Ha, T)

    fig = plt.figure()
    fig = plt.figure(facecolor='White')
    ax = fig.gca(projection='3d')

    surf = ax.plot_surface(Ha,T,TEff, rstride=1, cstride=1, cmap=cm.coolwarm,
           linewidth=0, antialiased=False)

    ax.set_zlim(TEff.min() ,TEff.max())
    ax.zaxis.set_major_locator(LinearLocator(10))
    ax.zaxis.set_major_formatter(FormatStrFormatter('%.02f'))
    fig.colorbar(surf, shrink=0.5, aspect=5)

    ax.set_title('Effect of temperature on Michaelis Menten-parameters \n at different Ha values')
    ax.set_xlabel('Activation energy, Ha (J mol-1)')
    ax.set_ylabel('Leaf surface temperature (C)')
    ax.set_zlabel('T-multiplier to reference value')


    plt.show()

plot_TemperatureEffectOnRuBisCOKinetics()

Answer 1

从你的问题中还不完全清楚你想要画出哪些“线”，但我猜你想对TemperatureEffectOnRuBisCOKinetics的一些固定值Ha（cHaOfGammaStar，cHaOfK0等）和相同范围的温度值（0到30步，每1步）。在

例如，要绘制cHaOfGammaStar，可以执行以下操作：

cHaOfGammaStar = np.array([29000])
z = TemperatureEffectOnRuBisCOKinetics(cHaOfGammaStar, T)

# we need to hold the axes to plot on top of the surface
ax.hold(True)

# we multiply cHaOfGammaStar by a vector of ones to make it the same length
# as T and z
l, = ax.plot(cHaOfGammaStar * np.ones(T.size), T, z, ' k')

# create a figure legend
ax.figure.legend((l,), ('cHaOfGammaStar',), loc=4, fancybox=True)

输出：

如果你想做一些更华丽的注释而不是仅仅使用一个图形图例，你应该看看HYRY's answer here。在