Using GeoAlchemy2 and PostGIS for Geospatial Data

Jan 25, 2025 | geospatial, postgis, geoalchemy2

In this blog, we will explore how GeoAlchemy2 and PostGIS work together to handle geospatial data in PostgreSQL. Geospatial data is crucial for applications that require location-based information, and these tools provide powerful capabilities for storing and querying such data.

Understanding GeoAlchemy2 and PostGIS

GeoAlchemy2 is an extension of SQLAlchemy that allows for the use of geospatial data types and functions in Python applications. PostGIS is a spatial database extender for PostgreSQL that adds support for geographic objects, enabling location queries to be run in SQL.

Example: Storing Latitude and Longitude in PostgreSQL

To demonstrate how to convert latitude and longitude into a PostgreSQL database using GeoAlchemy2 and PostGIS, follow these steps:

1. Install Required Packages: Make sure you have the necessary packages installed:

   pip install geoalchemy2 psycopg2
   

2. Define Your Model: Create a model that includes a geometry column for storing geospatial data:

   from sqlalchemy import create_engine, Column, Integer
   from geoalchemy2 import Geometry
   from sqlalchemy.ext.declarative import declarative_base
   
   Base = declarative_base()
   
   class Location(Base):
       __tablename__ = 'locations'
       id = Column(Integer, primary_key=True)
       geom = Column(Geometry(geometry_type='POINT', srid=4326))
   
   engine = create_engine('postgresql://user:password@localhost/mydatabase')
   Base.metadata.create_all(engine)
   

3. Insert Data: To insert latitude and longitude into the database:

   from sqlalchemy.orm import sessionmaker
   from geoalchemy2.shape import from_shape
   from shapely.geometry import Point
   
   Session = sessionmaker(bind=engine)
   session = Session()
   
   # Create a new location
   lat, lon = 40.7128, -74.0060  # Example coordinates for New York City
   point = Point(lon, lat)  # Note: Shapely uses (lon, lat) order
   location = Location(geom=from_shape(point, srid=4326))
   session.add(location)
   session.commit()
   

Conclusion

GeoAlchemy2 and PostGIS provide a robust framework for managing geospatial data in PostgreSQL. By following the steps outlined in this post, you can easily store and query geospatial information in your applications.

For more information, check out the GeoAlchemy2 documentation and the PostGIS documentation.