With GeoDjango we can find places in proximity to other places – this is very useful for things like a store locator. Let’s use a store locater as an example. Our store locator needs to be able to read in messy user input (zip, address, city, some combination). Then, locate any stores we have nearby.
General concept and theory
We have two problems to solve. One is to turn messy address input into a point on the globe. Then we need a way to query this point against other known points and determine which locations are close.
Set up known locations
Before we can really begin we need to set up GeoDjango. You can read the docs or use docker-compose.
It’s still a good idea to read the tutorial even if you use docker.
Let’s add a location. Something like:
class Location(models.Model): name = models.CharField(max_length=70) point = models.PointField() objects = models.GeoManager()
A PointField stores a point on the map. Because Earth is not flat we can’t use simple X, Y coordinates. Luckily you can almost think of Latitude and Longitude as X, Y. GeoDjango defaults to this. It’s also easy to get Latitude and Longitude from places like Google Maps. So if we want – we can ignore the complexities of mapping coordinates on Earth. Or you can read up on SRID if you want to learn more.
At this point we can start creating locations with points – but for ease of use add GeoModelAdmin to Django Admin to use Open Street Maps to set points.
from django.contrib import admin from django.contrib.gis.admin import GeoModelAdmin from .models import Location @admin.register(Location) class LocationAdmin(GeoModelAdmin): pass
Wow! We’re doing GIS!
Add a few locations. If you want to get their coordinates just type
location.point.x (or y).
Querying for distance.
Django has some docs for this. Basically make a new point. Then query distance. Like this:
from django.contrib.gis.geos import fromstr from django.contrib.gis.measure import D from .models import Location geom = fromstr('POINT(-73 40)') Location.objects.filter(point__distance_lte=(geom, D(m=10000)))
m is meters – you can pass all sorts of things though. The result should be a queryset of Locations that are near our “geom” location.
Already we can find locations near other locations or arbitrary points! In Part II I’ll explain how to use Open Street Maps to turn a fuzzy query like “New York” into a point. And from there we can make a store locator!