SGB-CPRM — Geologia, Recursos Minerais e Hidrogeologia

import requests
from PIL import Image
from io import BytesIO
 
# WMS do Mapa Geológico do Brasil via GeoServer
url = "https://geoservicos.sgb.gov.br/geoserver/geosgb/wms"
params = {
    "service": "WMS",
    "request": "GetMap",
    "layers": "geosgb:unidades_litoestratigraficas",
    "bbox": "-55,-15,-45,-10",
    "width": 800,
    "height": 600,
    "srs": "EPSG:4326",
    "format": "image/png",
}
 
response = requests.get(url, params=params)
if response.status_code == 200 and "image" in response.headers.get("content-type", ""):
    img = Image.open(BytesIO(response.content))
    img.save("mapa_geologico.png")
    print("Mapa geológico salvo em mapa_geologico.png")
else:
    print(f"Erro: {response.status_code}")

import requests
import pandas as pd
 
# Consultar ocorrências minerais via REST API
url = "https://geoportal.sgb.gov.br/server/rest/services/geologia/afloramentos/MapServer/0/query"
params = {
    "where": "UF = 'MG'",
    "outFields": "*",
    "returnGeometry": "true",
    "f": "geojson",
    "resultRecordCount": 100,
}
 
response = requests.get(url, params=params)
response.raise_for_status()
dados = response.json()
 
features = dados.get("features", [])
print(f"Afloramentos geológicos em MG: {len(features)}")
 
if features:
    props = [f["properties"] for f in features]
    df = pd.DataFrame(props)
    print(f"Colunas: {list(df.columns)}")
    print(df.head())

import requests
 
# SIAGAS — poços tubulares via REST API
url = (
    "https://geoportal.sgb.gov.br/server/rest/services"
    "/hidrologia/SIAGAS_MODDAD/MapServer/0/query"
)
params = {
    "where": "UF = 'SP'",
    "outFields": "POCO_ID,LATITUDE,LONGITUDE,PROFUNDIDADE,VAZAO,AQUIFERO,MUNICIPIO",
    "returnGeometry": "false",
    "f": "json",
    "resultRecordCount": 50,
}
 
response = requests.get(url, params=params)
response.raise_for_status()
dados = response.json()
 
features = dados.get("features", [])
print(f"Poços em SP: {len(features)}")
for f in features[:5]:
    attrs = f["attributes"]
    print(f"  {attrs.get('MUNICIPIO')} — Prof: {attrs.get('PROFUNDIDADE')}m, Vazão: {attrs.get('VAZAO')} m³/h")