La relación entre la investigación ionosférica y los sistemas de navegación por satélite global (GNSS) se puede analizar a través de dos enfoques. El enfoque directo utiliza modelos ionosféricos para mitigar su influencia, mientras que el enfoque indirecto aprovecha los datos de GNSS para modelar parámetros ionosféricos. Este estudio presenta un enfoque indirecto en el que el contenido total de electrones (TEC), un parámetro fundamental para las condiciones ionosféricas, se modela como una función armónica utilizando la expansión armónica esférica (SH). Se desarrollan modelos ionosféricos específicos de estación (locales) y regionales al descomponer la influencia ionosférica en componentes deterministas y estocásticos. Se utilizaron datos de GNSS de siete estaciones distribuidas uniformemente en Serbia para estimar los coeficientes de TEC. Los modelos locales se proporcionaron en formato ION como coeficientes SH, lo que permite la determinación de TEC en cualquier época, mientras que los modelos regionales tenían una resolución espacial y temporal de 2 horas. Los valores de raíz cuadrada media (RMS) de TEC variaron de 0.2 a 0.5 TECU (unidad de contenido total de electrones), con una media de 0.3 TECU. La validación contra mapas ionosféricos globales mostró un acuerdo dentro de 5.0 TECU. También se analizó el impacto del grado y orden de la expansión SH en los valores de TEC. Estos resultados refinan el modelado ionosférico regional, mejorando la precisión de posicionamiento de GNSS en Serbia y más allá.