如何在matplotlib中的3d表面上绘制一条线

时间:2014-02-03 13:57:27

标签: python numpy matplotlib

我有一个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()

1 个答案:

答案 0 :(得分:3)

您的问题并不完全清楚您要绘制哪些“线条”,但我猜您想要TemperatureEffectOnRuBisCOKinetics评估某些固定值HacHaOfGammaStarcHaOfK0等)和相同的温度值范围(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)

输出:

enter image description here

如果你想做一些更高级的注释,而不仅仅是使用图例,你应该看一下HYRY's answer here