从DEM创建累积的hypsograhic曲线

Question

我想从两个数字高程模型创建累积的hypsometric / hypsographic曲线，最好使用matplotlib。 我已经将DEM的值读入两个独立的数组，得到了一个简单的累积直方图，其中y轴上的计数/频率和x上的高程。现在，我需要y轴上的高程值，然后是x轴上的百分比，如下所示：https://de.wikipedia.org/wiki/Hypsografische_Kurve#/media/File:Hypsografische_Kurve_der_Schweiz.png

有人有建议如何实现这个目标吗？

编辑： 到目前为止，我有这个：

f,ax = pylab.subplots(figsize=(10,10))
ax.spines["top"].set_visible(False)
ax.spines["left"].set_visible(False)
ax.spines["right"].set_visible(False)
bins = [500.0, 1000.0, 1500.0, 2000.0, 2500.0, 3000.0, 3500.0, 4000.0, 4500.0, 5000.0, 5500.0, 6000.0, 6500.0]
pylab.xlabel("Elevation in [m]")
pylab.ylabel("Count")
pylab.xticks(bins)
ax.hist([height1,height], bins, rwidth=0.9, color=["r","b"], alpha=0.5, label=['dem_c','dem_t'],histtype="bar", cumulative = True)
pylab.legend(loc='upper right')

height和height1是包含通过gdal从DEM读取的高程值的列表。到目前为止，这是我的输出。