Fontes de dados para PostGIS

IBGE oficial atualizada

... principal

OpenStreetMap/planet

OpenStreetMap/overpass

NaturalEarth/10m-cultural

See download page, used wget https://www.naturalearthdata.com/http//www.naturalearthdata.com/download/10m/cultural/ne_10m_admin_1_states_provinces.zip. Them unzip and shp2pgsql ne_10m_admin_1_states_provinces.shp | psql

SELECT iso_3166_2, name_pt, round(st_area(geom,true)::numeric/1000000.0) area_km2  
FROM  ne_10m_admin_1_states_provinces where iso_3166_2 like 'BR-%';

iso_3166_2	name_pt	area_km2
BR-RS	Rio Grande do Sul	272283
BR-RR	Roraima	224746
BR-PA	Pará	1230800
BR-AC	Acre	153863
BR-AP	Amapá	139181
BR-MS	Mato Grosso do Sul	356042
BR-PR	Paraná	198666
BR-SC	Santa Catarina	95097
BR-AM	Amazonas	1565685
BR-MT	Mato Grosso	902952
BR-RO	Rondônia	237942
BR-MA	Maranhão	326134
BR-PI	Piauí	252011
BR-CE	Ceará	150228
BR-RN	Rio Grande do Norte	52812
BR-PB	Paraíba	56516
BR-PE	Pernambuco	97630
BR-AL	Alagoas	27619
BR-RJ	Rio de Janeiro	43978
BR-SP	São Paulo	248638
BR-SE	Sergipe	21606
BR-BA	Bahia	560050
BR-ES	Espírito Santo	45806
BR-MG	Minas Gerais	586372
BR-TO	Tocantins	278973
BR-GO	Goiás	341247
BR-DF	Distrito Federal	5792

Usando correcao

SELECT iso_3166_2, name_pt, round(st_area(st_buffer(geom,0.003),true)::numeric/1000000.0) area_km2,
    get_utmzone_bydump(geom)  utm_urid, st_area(st_transform(st_setsrid(geom,4326),32724))/1000000  
FROM  ne_10m_admin_1_states_provinces where iso_3166_2='BR-SE';

A área esperada (Wikipedia) é maior, foi necessário buffer de 0.003 (geography) para equiparar.