Dúvidas sobre cálculo de população e gráfico bloxpot em python

Bom dia a todos.

Estou tentando obter o cálculo de uma população e o respectivo gráfico bloxpot deste conforme o programa á baixo, depois de obter as latitudes e longitudes eu não estou conseguindo obter o cálculo de da população total de cada região e gerar um gráfico boxplot com a população dos municípios tendo, no eixo x cada uma das regiões e no eixo y a população.

GERAR UM GRÁFICO BOXPLOT COM A
POPULAÇÃO DOS MUNICÍPIOS TENDO, NO EIXO X,CADA UMA DAS REGIÕES DA CATEGORIA ALVO E, NO EIXO Y
A POPULAÇÃO

import pandas as pd

tabela1 = pd.read_csv('C:\\scripts\\pasta 1\\pasta 2\\pasta 3\\cadmun_simple.csv')
tabela2 = pd.read_csv('C:\\scripts\\pasta 1\\pasta 2\\pasta 3\\BP-2010.csv')

import matplotlib as mt
import matplotlib.pyplot as plt
import matplotlib_venn
from matplotlib_venn import venn2

tabela=tabela1.merge(tabela2, how='outer', left_on='MUNCODDV', right_on='geocode')
tabela.to_csv('C:\\scripts\\pasta 1\\pasta 2\\pasta 3/tabela.csv')
b = tabela.loc[138:152, :]
b.to_csv('C:\\scripts\\pasta 1\\pasta 2\\pasta 3/tabela_final.csv')

latitude = b.LATITUDE.median()
longitude = b.LONGITUDE.median()
quantis = b.population.quantile([.25,.50,.75])
pop_t = b.population.sum()

b['LATITUDE'] == latitude
b['LONGITUDE']== longitude
b.index[b.LATITUDE == latitude, b.LONGITUDE ==longitude] # Esta parte ocorre um erro