Visualización Geospacial Avanzada con GeoPandas

1. Tienes los estados de México en EPSG:4326 y necesitas su área en km². ¿Qué haces?

• Proyectar con to_crs(6372) y dividir .area entre 1e6
• Calcular .area directamente y multiplicar por 111
• Usar set_crs(6372) y después .area
• Calcular .length en lugar de .area

EPSG:4326 está en grados, así que .area devuelve grados cuadrados sin sentido. Primero se proyecta a un CRS en metros, como EPSG:6372 para México, y luego se divide entre 1e6.

2. ¿Qué hace gpd.sjoin(ciudades, estados, predicate="within")?

• Agrega a cada ciudad los datos del estado que la contiene
• Une las dos tablas por la columna name
• Recorta los estados a la forma de las ciudades
• Calcula la distancia de cada ciudad al estado más cercano

sjoin es una unión espacial: a cada ciudad le agrega las columnas del estado cuyo polígono la contiene, usando la geometría en lugar de una columna en común.

3. ¿Cuál es la diferencia entre set_crs y to_crs?

• set_crs etiqueta sin cambiar coordenadas; to_crs las transforma
• Son equivalentes; to_crs es el nombre moderno
• set_crs transforma y to_crs solo etiqueta
• set_crs solo funciona con Shapefile

set_crs solo etiqueta los datos con un CRS sin tocar las coordenadas; to_crs transforma las coordenadas al nuevo sistema. Etiquetar mal con set_crs desplaza el mapa sin dar error.

4. Quieres unir los 32 estados en tres regiones y sumar su área. ¿Qué método usas?

• dissolve(by="region", aggfunc={"area_km2": "sum"})
• groupby("region").sum(), que también une las geometrías
• sjoin(regiones, predicate="within")
• buffer(0) sobre cada región

dissolve(by=...) agrupa por una columna, une las geometrías de cada grupo en una sola y agrega los atributos con aggfunc, como un groupby espacial.

5. Haces un buffer de 100 km alrededor de cinco ciudades cercanas entre sí y sumas sus áreas. ¿Qué problema tiene el resultado?

• Las zonas traslapadas se cuentan más de una vez; hay que unirlos con union_all() antes
• Ninguno, la suma de áreas siempre es correcta
• El buffer solo funciona con polígonos, no con puntos
• El área de un buffer siempre es cero

Los buffers se traslapan y esas zonas se cuentan varias veces. Primero hay que unirlos con union_all() y después medir el área.