benedikt/rki-impfparser
1
0
Fork 0
Code Releases Activity

Compare commits

base: benedikt:55147eeed57ddb0496d7b0b4f21adcf6037cf8a1
Branches Tags
benedikt:main
..
compare: benedikt:78d4a48643caba6ca5d99662548d047534f78964
Branches Tags
benedikt:main

2 commits
55147eeed5 ... 78d4a48643

Author SHA1 Message Date
Benedikt Bastin
78d4a48643 feat: Added stats to title 
Added total vaccinations, vaccination percentage and mean daily vaccination rate to title
 2021-01-15 19:59:50 +01:00
Benedikt Bastin
386666d5b3 feat: Embed mean vaccination rate into plot 2021-01-15 19:47:30 +01:00
1 changed files with 15 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

19
plot.py
View file

@ -46,9 +46,14 @@ dates = impfungen['Datum']
daily = impfungen['Gesamtzahl Impfungen']
cumulative = np.cumsum(impfungen['Gesamtzahl Impfungen'])
mean_vaccinations_daily = np.mean(daily)
total_vaccinations = int(np.sum(daily))
total_vaccinations_percentage = float(total_vaccinations) / einwohner_deutschland
print(total_vaccinations_percentage)
to_be_vaccinated = einwohner_deutschland - np.sum(daily)
mean_vaccinations_daily = np.mean(daily)
mean_vaccinations_daily_int = int(np.round(mean_vaccinations_daily))
to_be_vaccinated = einwohner_deutschland - total_vaccinations
days_extrapolated = int(np.ceil(to_be_vaccinated / mean_vaccinations_daily))
extrapolated_dates = np.array([dates[0] + datetime.timedelta(days=i) for i in range(days_extrapolated)])
@ -56,6 +61,7 @@ extrapolated_dates = np.array([dates[0] + datetime.timedelta(days=i) for i in ra
extrapolated_vaccinations = mean_vaccinations_daily * range(days_extrapolated)
# Stand aus Daten auslesen
#stand = dates.iloc[-1]
#print_stand = stand.isoformat()
@ -78,7 +84,12 @@ def plot_extrapolation_portion(percentage):
plt.title(
'Tägliche Impfquote, kumulierte Impfungen und lineare Extrapolation bis {} % der Bevölkerung Deutschlands\n'
'Erstellung: {}, Datenquelle: RKI, Stand: {}'.format(print_percentage, print_today, print_stand)
'Erstellung: {}, Datenquelle: RKI, Stand: {}\n'
'Impfungen gesamt: {} ({} %), Durchschnittliche Impfrate: {} Impfungen/Tag'.format(
print_percentage,
print_today, print_stand,
total_vaccinations, np.round(total_vaccinations_percentage * 100, 2), mean_vaccinations_daily_int
)
)
ax2 = ax.twinx()
@ -89,7 +100,7 @@ def plot_extrapolation_portion(percentage):
ax2.set_xlim(xmax=dates[0] + datetime.timedelta(days=percentage * days_extrapolated))
ax2.grid(True)
ax2.plot(dates, cumulative, color='red', label='Kumulierte Impfungen')
ax2.plot(extrapolated_dates, extrapolated_vaccinations, color='orange', label='Kumulierte Impfungen (linear extrap.)')
ax2.plot(extrapolated_dates, extrapolated_vaccinations, color='orange', label='Extrap. kumulierte Impfungen\n({} Impfungen/Tag)'.format(mean_vaccinations_daily_int))
#ax2.plot()
ax.legend(loc='upper left')