將指定的檔位的期權，指定階段剩余到期日的期權數據合并，用于回測

import pandas as pd
import numpy as np
import akshare as ak
pd.set_option("display.max_rows",None)
pd.set_option("display.max_columns",None)nh_price = ak.nh_price_index('RM')
nh_price.head()nh_return = ak.nh_return_index('RM')
nh_return.head()nh_volatility = ak.nh_volatility_index('RM')
nh_volatility.head()# rm_fu = ak.get_futures_daily(start_date='20180101',end_date='20210624',market='CZCE')
# rm_fu = rm_fu.loc[rm_fu['symbol'].str.contains('RM')]
import datetime
rm_fu = pd.read_csv('rm_fu.csv')
rm_fu['date'] = rm_fu['date'].apply(lambda x:datetime.datetime.strptime(str(x),'%Y-%m-%d'))
rm_fu = rm_fu.iloc[:,6:]
#rm_fu.to_csv('rm_fu.csv')
rm_fu.head()rm = pd.read_csv('RM.csv')
import datetime
rm['date'] = pd.to_datetime(rm['date'])
rm = rm.iloc[:,6:]
#rm.to_csv('RM.csv')
rm.head()#rm['date'] = rm['date'].apply(lambda x:datetime.datetime.strptime(x, '%Y-%m-%d').strftime('%Y-%m-%d %H:%M:%S')).tolist()rm.rename(columns={'合約系列':'symbol'},inplace=True)
last = pd.merge(rm,rm_fu,how = 'left',on=['date','symbol'])
last.head()#根據長度，制作合約剩余交易日
#尋找每一個合約
"""
def fun1(df):for j in range(len(df['date'].unique())):df.iloc[j,-1] = len(df['date'])-jreturn df
"""
from tqdm import tqdm
final = pd.DataFrame()
for i in tqdm(last['品種代碼   '].unique()):#print(last['品種代碼   '].head(1))#print(i)last_1 = last.loc[last['品種代碼   '] == i]#print(last_1)last_1 = last_1.sort_values(by='date')last_1['rt'] = 0#print(last_1.date.head())#尋找每一個日期for j in tqdm(range(len(last_1['date'].unique()))):#print('j',j)#print(range(len(last_1['date'].unique())))if len(last_1['date'])-j-1 == 0:print('特殊')print(len(last_1['date']))print(j)print(last_1.iloc[j,:])else:last_1.iloc[j,-1] = len(last_1['date'])-j-1final = final.append(last_1)
#print(final)
final.to_csv('final_1.csv')#數據長度小于20日的刪除
from tqdm import tqdm
final = pd.DataFrame()
for i in tqdm(last['品種代碼   '].unique()):#print(last['品種代碼   '].head(1))#print(i)last_1 = last.loc[last['品種代碼   '] == i]#print(last_1)last_1 = last_1.sort_values(by='date')last_1['rt'] = 0#print(last_1.date.head())#尋找每一個日期for j in tqdm(range(len(last_1['date'].unique()))):#print('j',j)#print(range(len(last_1['date'].unique())))if len(last_1['date'])< 20:passelse:last_1.iloc[j,-1] = len(last_1['date'])-j-1final = final.append(last_1)
#print(final)
final.to_csv('final_2.csv')#選擇某一天的合約，
print(last['品種代碼   '].head(1))#選取某一天的全部合約
final = pd.read_csv('final_1.csv')
date = '2020-01-16'
test = final.loc[final['date']==date]
test.head()#根據CP，距離期貨價格數據長短，判定合約（也可以根據希臘字母，價格等）
final['date'][0]==date
final.to_csv('E:/final.csv')