复杂的阶梯函数覆盖在时间序列码的matplotlib图上

时间:2013-02-08 14:47:17

标签: python function datetime matplotlib

是否有人知道使用步骤函数pyplot.step(x,y)时的问题,其中x值由datetime实例给出? 我可以轻松绘制函数,除非我使用datetime系列作为x值。 我有一个类似于以下内容的矩阵:

In [3]: Minor_trend_matrix[0:20]
Out[3]: 
matrix([[ 0,  1],
        [ 1,  2],
        [ 1,  3],
        [ 1,  4],
        [-1,  5],
        [ 0,  6],
        [ 1,  7],
        [ 2,  8],
        [ 0,  9],
        [ 1, 10],
        [ 1, 11],
        [ 1, 12],
        [ 1, 13],
        [ 1, 14],
        [-1, 15],
        [-1, 16],
        [-1, 17],
        [ 1, 18],
        [ 1, 19],
        [-1, 21]])

注意第一列仅包含-2,-1,0,1,2 这些代表了阶梯函数所需的不同动作; -2是指数范围的高位,然后是指数范围的低位 -1是指数范围的低点 0是没有动作的 1是指数范围高 2指数范围低,然后指数范围高

索引引用datetime对象。我认为这就是step函数遇到麻烦的地方。任何建议或替代方案将不胜感激。 此外,关于整数表示的动作的细节,你认为,因为例如矩阵的运行可能如下所示

 [ 1,  2],
 [ 1,  3],
 [ 1,  4]

但我只希望在更改发生之前对最终项目进行操作,排除零以外我应该删除相关项目之前的所有前面条目(在不同条目之前,即在这种情况下-1,2, - 2)如果是这样我将如何去做? 此外,任何有关实现适当代码以实现此目标的帮助也将受到赞赏。

# Package Imports

import csv
import datetime
import matplotlib.pyplot as plt
import numpy as np
import dateutil.relativedelta as rd
import bisect
import scipy as sp

## File Settings

ipath = "/home/kane/Downloads/Charts/Wheat/W2/"
fname = "W2_0_L1N.CSV"

dateList = []
openList = []
highList = []
lowList = []
closeList = []
x = []
g = []
portfolio_list = []
result = []

portfolio = csv.DictReader(open(ipath + fname, "rb"))
portfolio_list.extend(portfolio)


for data in portfolio_list:
    dateList.append(data['Date'])
    openList.append(data[' Open'])
    highList.append(data[' High'])
    lowList.append(data[' Low'])
    closeList.append(data[' Close'])



for item in dateList:
    d = datetime.datetime.strptime(item, '%Y%m%d')

    t = d.toordinal()
    x.append(d)
    g.append(t)


dataArray = np.asarray([x, openList, highList, lowList, closeList])
dataArrayOrd = np.array([g, openList, highList, lowList, closeList], 
dtype = float)

代码的绘图部分如下:

# Plot settings
fig = plt.figure(figsize=(10, 5))
ax = fig.add_subplot(111)

line0 = plt.plot(x, simple_slope, 'k', label = "High Slope")
line1 = plt.plot(x, dataArray[1,:], 'g', label = "Open")  
line2 = plt.plot(x, dataArray[2,:], 'r', label = "High")
line3 = plt.plot(x, dataArray[3,:], 'y', label = "Low")
line4 = plt.plot(x, dataArray[4,:], 'b', label = "Close")
lns = line0 + line1 + line2 + line3 + line4
labs = [l.get_label() for l in lns]
ax.legend(lns, 
          labs, 
          loc='upper left', 
          frameon = False, 
          shadow=False, 
          fancybox=False, 
          labelspacing=0.1, 
          ncol=4, 
          prop={'size':9}
          )
ax.set_title(fname , fontsize = 12)



plt.xticks(fontsize = 8)
plt.yticks(fontsize = 8)
ax.set_ylabel( 'Cents', fontsize = 10 )
ax.grid()
plt.show()

此图中的x轴是时间和日期时间值,因此步骤函数应如下工作:

plt.step(x[Minor_trend_matrix[:,1]],dataArray[2,Minor_trend_matrix[:,1]])

当然,确切的问题有点复杂,需要访问

dataArray[3,Minor_trend_matrix[:,1]]
步进功能垂直线时的值应该向下移动到范围的低点,当需要同时访问高和低或低,然后高时,值更加复杂。这将是开发正确的条件语句的问题。但是,当前的问题在于,在我可以进展之前,让一个简单的版本能够正确使用日期时间值。

1 个答案:

答案 0 :(得分:1)

解决方案:

我用步进功能解决了这个问题。这是由于步进功能的工作方式。我必须通过for循环生成datetime对象。我现在可以得到一个我想要实现的简单版本。

# Package Imports

import csv
import datetime
import matplotlib.pyplot as plt
import numpy as np
import dateutil.relativedelta as rd
import bisect
import scipy as sp

## File Settings

ipath = "/home/kane/tripkane_python/"
fname = "datetime_test_data.csv"


portfolio_list = []
dateList = []
valueList = []
x = []

portfolio = csv.DictReader(open(ipath + fname, "rb"))
portfolio_list.extend(portfolio)


for data in portfolio_list:
    dateList.append(data['date'])
    valueList.append(data[' value'])

for item in dateList:
    d = datetime.datetime.strptime(item, '%Y%m%d')
    t = d.toordinal()
    x.append(d)


dataArray = np.asarray([x, valueList])

stepList = [0, 3, 5]

TIda = []
TIv = []

for n in range(0, len(stepList)):
    TIda.append(x[stepList[n]])
    TIv.append(valueList[stepList[n]])


# Plot 

fig = plt.figure(figsize=(10, 5))
ax = fig.add_subplot(111)

line1 = plt.plot(x, valueList, 'b')
lns = line1

labs = [l.get_label() for l in lns]
ax.legend(lns, 
          labs, 
          loc='upper left', 
          frameon = False, 
          shadow=False, 
          fancybox=False, 
          labelspacing=0.1, 
          ncol=4, 
          prop={'size':9}
          )
ax.set_title(fname , fontsize = 12)



plt.xticks(fontsize = 8)
plt.yticks(fontsize = 8)
ax.set_ylabel( 'Cents', fontsize = 10 )
ax.grid()

plt.step(TIda, TIv, 'r')

plt.show()

现在的问题是如何产生适当的条件语句来实现目标。