Understanding the Landscape · Lesson 4 of 5
Land Cover & Coordinate Systems
~15 min
The Concept
With terrain understood, your supervisor adds another layer to the map: land cover, forest, wetland, farmland, and urban area across the basin. But before you can combine it with the DEM, she warns you about a mistake that has ruined more GIS projects than any other: assuming two maps use the same coordinate system when they don't.
A coordinate reference system (CRS, meaning the rulebook a map uses to turn a spot on Earth into two numbers) is the agreed-upon way of turning a location on Earth's curved surface into numbers you can plot on a flat map. Different datasets are often recorded in different CRSs, one might use latitude and longitude, another a local flat projection measured in meters.
If you plot two datasets in different CRSs on the same map without converting them, they will not line up, sometimes by a few meters, sometimes by hundreds of kilometers, even though every individual number is 'correct' in its own system.
The Analogy
It is like two people describing the same seat at a stadium, one using row and seat number, the other using GPS coordinates. Both are correct. But if you tried to plot them on the same seating chart without translating one into the other, you would think they were talking about two completely different seats.
Code (optional)
import geopandas as gpd
landcover = gpd.read_file("bluewater_landcover.geojson")
print("Original CRS:", landcover.crs)
# Reproject to match the DEM's coordinate system before combining them
landcover_matched = landcover.to_crs("EPSG:32619")
print("Reprojected CRS:", landcover_matched.crs)Line by line
- geopandas reads spatial data along with its coordinate reference system metadata.
- Printing .crs shows exactly which system the data currently uses, always check this before combining datasets.
- to_crs() mathematically converts every coordinate to match another system, here UTM zone 19N (EPSG:32619), which Bluewater Basin's DEM uses.
Why Real Researchers Care
CRS mismatches are one of the most common real-world GIS errors, they've caused everything from misaligned emergency-response maps to incorrect area calculations in published research. Checking and matching coordinate systems before any spatial analysis is standard professional practice, not an optional detail.
Quick Check
Q1. What happens if you plot two datasets with different, unconverted coordinate systems on the same map?
Your Goal
Bluewater Basin's rain gauge network was recorded in latitude/longitude (EPSG:4326), but the DEM uses UTM zone 19N (EPSG:32619). Write the one line of code you'd add before combining them.
Hint: You need to reproject one dataset to match the other, look at the to_crs() call above.
Teach It Back
Explain to a teammate why 'the coordinates are correct' doesn't guarantee 'the map is correct'.