我有一個循環遍歷 50 個數據幀併計算p90_annual
每個調用的變量。在這個循環中,我想p90_annual
在 pdf 中繪製一個頁面,其中每個頁面都對應p90_annual
於每個數據框。(我想要一個 50 頁的 pdf,其中每一頁都是p90_annual
每個數據框的圖)我目前正在使用:
with PdfPages('90thPercentile.pdf') as pdf:
plt.figure
plt.plot(p90_annual)
plt.title(j)
plt.ylabel('Days Above 90th Percentile')
pdf.savefig()
plt.close()
當我這樣做時,我只會得到一頁,其中包含最後一個p90_annual
繪製的實例。我如何修改它以便p90_annual
它在循環時為每個實例添加一個新頁面?
對於上下文......下面是我試圖讓它在其中工作的更大的循環
# Huge Loop
for j in TempDict:
#Make Baseline
df=TempDict[j]
df=pd.to_numeric(df.tmax, errors='coerce')
mask = (df.index >= '1900-01-01') & (df.index <= '1940-12-31')
Baseline=df.loc[mask]
Tmax=Baseline.astype(np.float)
Index=Baseline.index
DailyBase=pd.DataFrame(data={'date':Index,'tmax':Tmax})
#pivot dataframe
DailyBase['year']=DailyBase.date.dt.year
DailyBase['day']=DailyBase.date.dt.strftime('%m-%d')
BaseResult=DailyBase[DailyBase.day!='02-29'].pivot(index='year',columns='day',values='tmax')
#Calculate Percentiles
BaseResult.index=list(range(1,42))
BaseResult.insert(0,'12-31_',BaseResult['12-31'])
BaseResult.insert(0,'12-30_',BaseResult['12-30'])
BaseResult['01-01_'] = BaseResult['01-01']
BaseResult['01-02_'] = BaseResult['01-02']
p90_todict = {}
for i in range(len(BaseResult.columns)-4):
index = i+2
p90_todict[BaseResult.columns[index]] = np.quantile(BaseResult.iloc[:,index-2:index+3].dropna(),.9)
np.quantile(BaseResult.iloc[:,index-2:index+3].dropna(),.98)
#Make POR dataframe
#pull tmax and dates from original ACIS data
FullTmax=df.astype(np.float)
FullIndex=df.index
#create and rotate data frame
DailyPOR=pd.DataFrame(data={'date':FullIndex,'tmax':FullTmax})
DailyPOR['year']=DailyPOR.date.dt.year
DailyPOR['day']=DailyPOR.date.dt.strftime('%m-%d')
PORResult=DailyPOR[DailyPOR.day!='02-29'].pivot(index='year',columns='day',values='tmax')
#Compare POR and baseline
import copy
#eliminate leap years from POR daily data
noleap_DailyPOR = copy.copy(DailyPOR[DailyPOR.day != '02-29'])
noleap_DailyPOR.index = noleap_DailyPOR.date
#Use only winter months
only_winter = noleap_DailyPOR[(noleap_DailyPOR.index.month >= 12) | (noleap_DailyPOR.index.month <= 2)]
#set results to 0 for counts
p90results = pd.DataFrame(index = only_winter.date)
p90results['above90'] = 0
#Compare POR and percentiles
for index, row in only_winter.iterrows():
if row.tmax > p90_todict[row.day]:
p90results.loc[row.date,'above90'] = 1
#Sum annual counts above percentiles
p90_annual=p90results.groupby(p90results.index.year).sum()
with PdfPages('90thPercentile.pdf') as pdf:
plt.rcParams['text.usetex'] = False
plt.figure
plt.plot(p90_annual)
plt.title(j)
plt.ylabel('Days Above 90th Percentile')
pdf.savefig()
plt.close()
將 for 循環放在裡面with PDFPages() as pdf:
會有所幫助:
import numpy as np
from matplotlib.backends.backend_pdf import PdfPages
import matplotlib.pyplot as plt
import pandas as pd
p90_annual = pd.DataFrame({"data": [19,20,13]})
with PdfPages('90thPercentile.pdf') as pdf:
for i in range(0,50):
plt.figure
plt.plot(p90_annual)
plt.title(i)
plt.ylabel('Days Above 90th Percentile')
pdf.savefig()
plt.close()
更新:我想我沒有說得足夠清楚,重點是移到with PDFPages() as pdf:
你的循環之外,然後一切都會被歸檔。
with PdfPages('90thPercentile.pdf') as pdf:
#your loop here
for j in TempDict:
....
plt.figure
plt.plot(p90_annual)
plt.title(i)
plt.ylabel('Days Above 90th Percentile')
pdf.savefig()
plt.close()
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句