How can I parse GeoJSON with Python

Question

I have have geojson data from a query which I now want to parse and print on screen. My current code is:

import urllib
import geojson

while True:

    url = 'https://earthquake.usgs.gov/fdsnws/event/1/query?format=geojson&starttime=2012-01-01&endtime=2017-03-01&minmagnitude=4.0&maxmagnitude=9.0&minlongitude=5.95&maxlongitude=10.50&minlatitude=45.81&maxlatitude=47.81'
    uh = urllib.urlopen(url)
    data = uh.read()
    print data
    break

It seems that data is a simple string. However, I thought it could be parsed like a json parameter. How do I have to handle geojson data in order to print a single point, e.g. to extract the coordinates of the first point only?

I guess the `while True` and the `break` can be removed completely? — Martin Thoma, Jul 25 '18 at 06:30

score 51 · Answer 1 · answered Dec 13 '17 at 11:54

51

import geojson
with open(path_to_file) as f:
    gj = geojson.load(f)
features = gj['features'][0]

answered Dec 13 '17 at 11:54

4

How does this differ from loading the file with the base package `json`? – duhaime Mar 05 '21 at 00:17
Works! Should be marked as the accepted answer! – hzitoun Apr 01 '21 at 22:00
2

The `geojson` library didn't have any releases since 2019, so I recommend to use more popular and recent `geopandas` instead – lucidyan May 29 '21 at 16:06

score 31 · Accepted Answer · answered Mar 12 '17 at 22:41

31

You can read it like any json:

import json
data = json.loads(datastring)
data['features'][0]['geometry'] #Your first point

answered Mar 12 '17 at 22:41

Jesse Bakker

2,403
13
25

2

if loading from a file: use data = json.load(fileobject) – onxx Apr 15 '21 at 18:04
This is what I am searching for. When I use your solution, I get a dict which looks like this `0 MultiLineString [[[8.5551415, 49.8598025], [8.5563399, 49.8598...` with 3 columns. As you can see the third column has the coordinates I want. Do you know by any chances how I can transfer the coordinates into a dataframe with latitude and longitude as seperate colums since the values represent this? – Gobrel Jan 04 '22 at 23:09

score 22 · Answer 3 · answered Apr 24 '19 at 06:17

22

You can also use geopandas:

import geopandas as gpd
earthquake = gpd.read_file('earthquake.geojson')
print(earthquake.head())

answered Apr 24 '19 at 06:17

sagarr

1,152
1
9
16

8

Caveat: this method is not reading data as geojson. It reads the data as geodataframe. – Yanni Cao Jun 06 '20 at 03:08
2

This is often the quickest and easiest way to read GeoJSON and make sense of what it contains. – Dahn Nov 16 '20 at 19:05

score 5 · Answer 4 · answered Jun 28 '19 at 20:21

5

You can read it with json import, and file open:

import json
with open(path) as f:
    data = json.load(f)
for feature in data['features']:
    print(feature['properties'])

answered Jun 28 '19 at 20:21

poetyi

236
4
13

I had to use utf-8 encoding: open(path, encoding='utf-8') – Stefan May 20 '22 at 10:45

score 4 · Answer 5 · answered Sep 11 '19 at 05:30

4

You can use pandas library directly

import pandas as pd
data = pd.read_json('File.geojson')

The important this is to understand the structure of this json file and manipulate the dictionaries in there

answered Sep 11 '19 at 05:30

Himanshu Aggarwal

163
1
1
5

How can I parse GeoJSON with Python

5 Answers5

Linked