geopandas-projections

GeoPandas Spatial Projections

Overview

GeoPandas is built on top of Shapely and Fiona, enabling geographic data manipulation with proper coordinate reference systems (CRS). Using correct projections is critical for accurate distance calculations and spatial operations.

Installation

pip install geopandas shapely fiona pyproj

Key Concepts

Coordinate Reference Systems (CRS)

• EPSG:4326: WGS84 (lat/lon), commonly used for geographic data but NOT suitable for distance calculations
• EPSG:3857: Web Mercator, used for web mapping
• Regional Projected CRS: For accurate local distance calculations (e.g., UTM zones)

Distance Calculations

Always project to a projected CRS before calculating distances. Geographic CRS (like EPSG:4326) measure in degrees, not kilometers.

Code Examples

Creating GeoDataFrames from Points

import geopandas as gpd
from shapely.geometry import Point
import pandas as pd

# From earthquake data
earthquakes_df = pd.read_json('/root/earthquakes_2024.json')
geometry = [Point(xy) for xy in zip(earthquakes_df['longitude'], earthquakes_df['latitude'])]
gdf = gpd.GeoDataFrame(earthquakes_df, geometry=geometry, crs='EPSG:4326')

Loading GeoJSON with Boundaries

import json

# Load GeoJSON and convert to GeoDataFrame
with open('/root/PB2002_boundaries.json', 'r') as f:
    geojson_data = json.load(f)
boundaries_gdf = gpd.GeoDataFrame.from_features(geojson_data['features'], crs='EPSG:4326')

Projecting to Projected CRS

# Project to a suitable CRS for distance calculations
# Example: project to Mercator for global analysis
gdf_projected = gdf.to_crs('EPSG:3857')
boundaries_projected = boundaries_gdf.to_crs('EPSG:3857')

# Or use a specific UTM zone for a region
# EPSG:32633 is UTM zone 33N

Spatial Filtering (Point in Polygon)

# Check if points fall within polygons
earthquakes_in_plate = gpd.sjoin(gdf, plate_polygons, how='inner', predicate='within')