INCRA Estrutura Fundiária — Imóveis Rurais e Assentamentos

import pandas as pd
 
# Exportação via Painel de Assentamentos (https://painel.incra.gov.br/)
# Suporta download em PDF, XLS e ODS
# Alternativamente, baixar o CSV do dados.gov.br
 
df = pd.read_csv(
    "assentamentos_incra.csv",
    sep=";",
    encoding="utf-8",
    dtype=str,
)
 
print(f"Total de assentamentos: {len(df):,}")
print(f"Colunas: {list(df.columns)}")
 
# Assentamentos por UF
por_uf = df["UF"].value_counts()
print("\nAssentamentos por UF (top 10):")
print(por_uf.head(10))
 
# Total de famílias assentadas
df["FAMILIAS"] = pd.to_numeric(df["FAMILIAS"], errors="coerce")
print(f"\nTotal de famílias assentadas: {df['FAMILIAS'].sum():,.0f}")

import geopandas as gpd
import matplotlib.pyplot as plt
 
# Leitura de parcelas certificadas (download do Acervo Fundiário)
gdf = gpd.read_file("parcelas_sigef_sp.shp")
 
print(f"Parcelas certificadas em SP: {len(gdf):,}")
print(f"Área total (ha): {gdf['area_ha'].sum():,.0f}")
print(f"CRS: {gdf.crs}")
 
# Mapa de parcelas
fig, ax = plt.subplots(figsize=(10, 10))
gdf.plot(ax=ax, edgecolor="green", facecolor="lightgreen", linewidth=0.2)
ax.set_title("Parcelas certificadas — SIGEF/INCRA (SP)")
ax.set_axis_off()
plt.tight_layout()
plt.savefig("parcelas_sigef_sp.png", dpi=150)

import pandas as pd
 
# Dados do SNCR via Base dos Dados (basedosdados.org)
# pip install basedosdados
import basedosdados as bd
 
query = """
SELECT
    sigla_uf,
    COUNT(*) AS total_imoveis,
    SUM(area_total_hectares) AS area_total_ha,
    AVG(area_total_hectares) AS area_media_ha
FROM `basedosdados.br_incra_sncr.imoveis_rurais`
GROUP BY sigla_uf
ORDER BY area_total_ha DESC
"""
 
df = bd.read_sql(query, billing_project_id="SEU_PROJETO_BIGQUERY")
print("Estrutura fundiária por UF:")
print(df.to_string(index=False))