Population by Region

import matplotlib.pyplot as plt
import numpy as np
from mpl_toolkits.mplot3d import Axes3D

fig = plt.figure(figsize=(12, 8), facecolor='white')
ax = fig.add_subplot(111, projection='3d')
ax.set_facecolor('white')

# Data: 5 regions, 4 age groups
regions = 5
age_groups = 4
xpos = np.arange(regions)
ypos = np.arange(age_groups)
xpos, ypos = np.meshgrid(xpos, ypos)
xpos = xpos.flatten()
ypos = ypos.flatten()
zpos = np.zeros_like(xpos)

dx = dy = 0.65
dz = np.array([120, 85, 95, 45, 180, 110, 75, 60,
               150, 95, 88, 52, 200, 130, 92, 68, 90, 65, 55, 38])

# Green gradient
colors = plt.cm.colors.LinearSegmentedColormap.from_list('', ['#27F5B0', '#6CF527'])
bar_colors = [colors(v/max(dz)) for v in dz]

ax.bar3d(xpos, ypos, zpos, dx, dy, dz, color=bar_colors, alpha=0.9, edgecolor='#000000', linewidth=0.3)

ax.set_xlabel('Region', fontsize=11, color='#1f2937', labelpad=10)
ax.set_ylabel('Age Group', fontsize=11, color='#1f2937', labelpad=10)
ax.set_zlabel('Population (K)', fontsize=11, color='#1f2937', labelpad=10)
ax.set_title('Population Distribution by Region', fontsize=14, color='#1f2937', fontweight='bold', pad=20)

ax.set_xticks(range(5))
ax.set_xticklabels(['North', 'South', 'East', 'West', 'Central'])
ax.set_yticks(range(4))
ax.set_yticklabels(['0-18', '19-35', '36-55', '55+'])
ax.tick_params(colors='#000000', labelsize=9)

ax.xaxis.pane.fill = False
ax.yaxis.pane.fill = False
ax.zaxis.pane.fill = False
ax.xaxis.pane.set_edgecolor('#000000')
ax.yaxis.pane.set_edgecolor('#000000')
ax.zaxis.pane.set_edgecolor('#000000')
ax.grid(True, alpha=0.5, linewidth=0.5)
ax.xaxis._axinfo['grid']['color'] = '#000000'
ax.yaxis._axinfo['grid']['color'] = '#000000'
ax.zaxis._axinfo['grid']['color'] = '#000000'
ax.xaxis._axinfo['tick']['color'] = '#000000'
ax.yaxis._axinfo['tick']['color'] = '#000000'
ax.zaxis._axinfo['tick']['color'] = '#000000'


ax.xaxis.line.set_color('#000000')
ax.yaxis.line.set_color('#000000')
ax.zaxis.line.set_color('#000000')

plt.tight_layout()
plt.show()