如何在matplotlib中绘制3d曲面上的直线

Question

我有一个3d图，想在图的表面上画几条线。我不清楚我应该如何组织这些线的数据，使其落在表面上。

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

函数plot_TemperatureEffectOnRuBisCOKinetics绘制3D绘图。现在我想看看在表面上的函数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

从你的问题中还不能完全清楚你正在尝试绘制哪些“线”，但我猜你想要计算Ha的一些固定值的TemperatureEffectOnRuBisCOKinetics (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。