高校生と高齢者の利用で
中央商店街を活性化しよう！！

import os
import numpy as np
import pandas as pd
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

plt.rcParams['font.family'] = 'Hiragino Sans'   # Windowsは 'Yu Gothic'
plt.rcParams['axes.unicode_minus'] = False
plt.rcParams['figure.dpi'] = 150

FIG_DIR = 'html/figures'
os.makedirs(FIG_DIR, exist_ok=True)

# 原論文 図2 の報告値をそのまま定義（%）
years = ['H10', 'H12', 'H13', 'H15', 'H17', 'H19', 'H21', 'H23', 'H25']
vac = pd.DataFrame({
    '市全体':     [6.7, 10.4, 11.1, 11.9, 10.6, 14.2, 15.5, 16.0, 16.2],
    '中心市街地': [5.2,  5.4,  5.7,  6.2,  7.6,  8.9, 11.2, 13.8, 14.0],
    '中央商店街': [1.6,  2.6,  2.4,  2.3,  2.7,  6.5,  9.8,  9.6, 11.3],
}, index=years)

print(vac.to_string())
ratio = vac['中央商店街'].iloc[-1] / vac['中央商店街'].iloc[0]
print(f"中央商店街: {ratio:.1f}倍 ≒ 約7倍")

fig, ax = plt.subplots(figsize=(10, 5.5))
colors_v = {'市全体': '#1F3B99', '中心市街地': '#D63FA0', '中央商店街': '#E8A33D'}
markers  = {'市全体': 'o', '中心市街地': 's', '中央商店街': '^'}
for col in vac.columns:
    ax.plot(years, vac[col], marker=markers[col], color=colors_v[col],
            linewidth=2, markersize=7, label=col)
    for x, y in zip(years, vac[col]):        # 各点に報告値を印字
        ax.annotate(f'{y}', (x, y), textcoords='offset points',
                    xytext=(0, 7), ha='center', fontsize=8)
ax.set_ylabel('空き店舗率（%）'); ax.set_ylim(0, 18); ax.legend()
fig.savefig(f'{FIG_DIR}/2020_H5_2_fig1.png', bbox_inches='tight')

# 原論文 図5・図6 の報告値（本文に明記された代表値のみ）
area = pd.DataFrame({
    '65歳以上人口の割合': [30, 25],
    '単身世帯の割合':     [65, 40],
}, index=['商店街周辺エリア', '松山市全体'])
print(area.to_string())

# 原論文 図9 の報告値（半径約1km圏の7校）
schools = pd.DataFrame({
    '高校': ['松山東', '松山南', '松山工業', '松山商業', '聖カタリナ学園', '済美', '松山東雲'],
    '全校生徒数': [1050, 1050, 950, 950, 950, 1800, 300],
})
print(schools.to_string(index=False))
print(f"合計: 約{schools['全校生徒数'].sum():,}人")

# 原論文 図10（利用頻度, n=1023）の報告値を円グラフに
freq_labels = ['毎日', '週4〜6回', '週1〜3回', 'それ以下', '通らない']
freq_vals   = [3, 3, 9, 57, 28]   # %（合計100）

fig, ax = plt.subplots(figsize=(7, 7))
ax.pie(freq_vals,
       labels=[f'{l}\n{v}%' for l, v in zip(freq_labels, freq_vals)],
       startangle=90, counterclock=False,
       wedgeprops=dict(edgecolor='white', linewidth=1.5))
fig.savefig(f'{FIG_DIR}/2020_H5_2_fig4.png', bbox_inches='tight')
print('週1回以下（それ以下+通らない）=', 57 + 28, '%')

# 図11（自習室の利用希望）と 図12（希望価格帯）の報告値
wish_vals  = [58, 42]              # はい / いいえ
price_vals = [2, 4, 19, 16, 59]   # 1000円以上 … 払いたくない
print('自習室に賛成:', wish_vals[0], '%')
print('料金を払いたくない:', price_vals[-1], '%（最多）')