pandas移动窗口计算脉动

import os
import pandas as pd
import numpy as np
#每个文件只取前三行后生成的maidong_files
root = 'D:/maidong_files'

all_U = pd.DataFrame()
all_V = pd.DataFrame()
all_chuizhi = pd.DataFrame()

all_time = []

#定义计算函数，用到数组的元素级别计算
def calUV(jiaodu,sudu):

    u = (-1) * np.sin( float(jiaodu) * np.pi / 180) * float(sudu)
    v = (-1) * np.cos( float(jiaodu) * np.pi / 180) * float(sudu)

    return '%.2f' % u,'%.2f' % v
    #return u,v
for every_file in os.listdir(root):

#每个文件取一列
    df_jiaodu=pd.read_csv(os.path.join(root,every_file),sep=',').iloc[:,2].replace('/',np.nan)
    array_jiaodu =df_jiaodu.values
    df_sudu=pd.read_csv(os.path.join(root,every_file),sep=',').iloc[:,3].replace('/',np.nan)
    array_sudu = df_sudu.values
    df_chuizhi=pd.read_csv(os.path.join(root,every_file),sep=',').iloc[:,4].replace('/',np.nan)

    calUV_ufunc = np.frompyfunc(calUV, 2, 2)
    [U, V] = calUV_ufunc(array_jiaodu,array_sudu)
        #print U.dtype注意UV的类型为object
        #转换UV类型为float
    U = U.astype(np.float)
    V = V.astype(np.float)

        #矩阵重新转为dataframe
    Udf = pd.DataFrame(U)
    Vdf = pd.DataFrame(V)

    all_U = all_U.append(Udf.T, ignore_index=True)
    all_V = all_V.append(Vdf.T, ignore_index=True)
    all_chuizhi = all_chuizhi.append(df_chuizhi.T, ignore_index=True)

#解析文件名中所包含的日期并转换为datatime格式    
    time_index = pd.to_datetime(every_file.split('_')[4])
    all_time.append(time_index)

#把time作为index
all_U.index = all_time
all_V.index = all_time
all_chuizhi.index = all_time

#生成日期范围
rng = pd.date_range(start = '20130628000000', end = '20130630010000', freq='5min')
#reindex从而实现日期范围内未有的部分全部插空值
all_time_U = all_U.reindex(rng)
all_time_V = all_V.reindex(rng)
all_time_chuizhi = all_chuizhi.reindex(rng)
all_time_U.to_csv('all_time_U.txt')
all_time_V.to_csv('all_time_V.txt')
all_time_chuizhi.to_csv('all_time_chuizhi.txt')

def cal_maidong(df):
    maidong = df - pd.rolling_mean(df,window=12,center=6,min_periods=1)
    return maidong

U_maidong = cal_maidong(all_time_U)
V_maidong = cal_maidong(all_time_V)
chuizhi_maidong = cal_maidong(all_time_chuizhi)
U_maidong['20130628000000':'20130630000000'].to_csv('c:/U_maidong.txt')
V_maidong['20130628000000':'20130630000000'].to_csv('c:/V_maidong.txt')
chuizhi_maidong['20130628000000':'20130630000000'].to_csv('c:/chuizhi_maidong.txt')