PyProj in Python: Complete Beginner's Setup Guide

6 hours ago
3 min read

Accurate coordinate transformations, map projections, and calculations in geodesy are critical to geospatial applications. Coordinate reference systems (CRS) are essential to all applications that use GIS, GPS, conduct spatial analysis, and retrieve data from mapping platforms.

PyProj is one of the most advanced Python libraries available for executing coordinate transformations and carrying out geospatial computations. It is built on top of the industry-standard PROJ library and provides the ability to accurately convert coordinates between different spatial reference systems to developers, GIS analysts, data scientists, and researchers.

What Is Pyproj?

Pyproj is a Python interface to the PROJ library for cartographic projection and coordinate transformation.

It provides functionality for:

Coordinate Transformation
CRS (Coordinate Reference System) Management
Geodetic Calculations
Datum Transformations
Projection Conversions
Distance/Area Calculations
Integration with GIS Libraries

Pyproj is commonly used with:

GeoPandas
Shapely
Rasterio
Fiona
GDAL
PostGIS
OpenStreetMap Workflows

Why Use Pyproj?

Modern geospatial datasets often use different coordinate systems.

For example:

System	Description
EPSG:4326	WGS84 Latitude/Longitude
EPSG:3857	Web Mercator
EPSG:32633	UTM Zone 33N
EPSG:5070	USA Contiguous Albers Equal Area

Without a coordinate transformation, spatial data from multiple sources may not align correctly.

Pyproj provides the solution for this type of issue by offering:

High Precision Transformations
Access to Thousands of CRS Definitions
Fast Performance
Industry Standard Compatibility
Cross-Platform Support

Installing PyProj

Install the latest version using pip:

pip install pyproj

Verify installation:

pip show pyproj

Example output:

Name: pyproj
Version: 3.7.1

Installing Specific Versions

For compatibility requirements:

pip install pyproj==3.6.1

Check installed version:

import pyproj

print(pyproj.__version__)

Output:

3.6.1

Understanding Coordinate Reference Systems (CRS)

A CRS defines how geographic coordinates relate to real-world locations.

Examples:

WGS84

EPSG:4326

Uses:

GPS devices
Google Maps APIs
Web mapping services

Format:

Latitude, Longitude

Example:

40.7128, -74.0060

(New York City)

Web Mercator

EPSG:3857

Uses:

Web maps
Tile servers
Mapping applications

Coordinates are represented in meters rather than degrees.

Importing PyProj

Basic import:

from pyproj import CRS, Transformer

Verify functionality:

from pyproj import CRS

crs = CRS.from_epsg(4326)

print(crs)

Output:

EPSG:4326

Creating CRS Objects

Create CRS using EPSG codes:

from pyproj import CRS

wgs84 = CRS.from_epsg(4326)

print(wgs84)

Output:

EPSG:4326

Create CRS from projection strings:

crs = CRS.from_string("+proj=longlat +datum=WGS84")

Coordinate Transformation Basics

Coordinate transformation converts coordinates from one CRS to another.

Example:

from pyproj import Transformer

transformer = Transformer.from_crs(
    "EPSG:4326",
    "EPSG:3857",
    always_xy=True
)

x, y = transformer.transform(-74.0060, 40.7128)

print(x, y)

Output:

-8238310.23 4970071.57

Transforming Coordinates Back

Reverse transformation:

reverse = Transformer.from_crs(
    "EPSG:3857",
    "EPSG:4326",
    always_xy=True
)

lon, lat = reverse.transform(
    -8238310.23,
    4970071.57
)

print(lat, lon)

Output:

40.7128 -74.0060

Working with UTM Coordinates

Universal Transverse Mercator (UTM) is widely used in surveying and engineering.

Convert latitude and longitude to UTM:

from pyproj import Transformer

transformer = Transformer.from_crs(
    "EPSG:4326",
    "EPSG:32633",
    always_xy=True
)

easting, northing = transformer.transform(
    12.4924,
    41.8902
)

print(easting, northing)

Output:

291962.58 4640626.11

Geodetic Distance Calculations

PyProj includes accurate geodesic calculations.

Example:

from pyproj import Geod

geod = Geod(ellps="WGS84")

az12, az21, distance = geod.inv(
    -74.0060,
    40.7128,
    -118.2437,
    34.0522
)

print(distance / 1000)

Output:

3944.42

Distance in kilometers between New York and Los Angeles.

Calculating Polygon Area

PyProj can calculate geodesic areas.

from pyproj import Geod

geod = Geod(ellps="WGS84")

lons = [-74, -73, -73, 74]
lats = [40, 40, 41, 41]

area, perimeter = geod.polygon_area_perimeter(
    lons,
    lats
)

print(abs(area))

Useful for GIS boundary analysis.

Listing CRS Information

Inspect CRS metadata:

from pyproj import CRS

crs = CRS.from_epsg(4326)

print(crs.name)
print(crs.axis_info)

Output:

WGS 84

Additional metadata includes:

Datum
Ellipsoid
Units
Coordinate axes

PyProj Integration with GeoPandas

Example:

import geopandas as gpd

gdf = gpd.read_file("data.shp")

gdf = gdf.to_crs("EPSG:3857")

GeoPandas internally uses PyProj for CRS transformations.

The PyProj library is an important component for any developer looking to build a geospatial application. It enables users to transform coordinates and perform geodetic calculations using Coordinate Reference System (CRS) management, projections, and transform methodologies.

The greatest challenge that new users of PyProj will face is grasping the idea of CRS, as well as installing PyProj correctly, and using the provided Transformer class and CRS class properly. However, as geospatial applications become more complex, PyProj can deliver the performance, accuracy, and reliability expected of a production-quality spatial analysis application.

Developers involved with GPS, GIS, mapping applications, or scientific research will find that PyProj is one of the major components of their geospatial developer toolkit.

For more information or any questions regarding PyProj, please don't hesitate to contact us at

Email: info@geowgs84.com

USA (HQ): (720) 702–4849

GeoWGS84AI

(A GeoWGS84 Corp Company)

https://www.geowgs84.ai

https://www.geowgs84.com/services/deep-learning-with-geospatial-data

Visit www.lizardtech.com - The Creators of the MrSID Format