To obtain earthquake events surrounding San Francisco (within 200 km) since 2013, minimum magnitude of 2.5, in geojson format:
import code.usgs
usgs.APIquery(starttime = "2013-01-01", endtime = "",
minmagnitude = "2.5",
latitude = "37.77", longitude = "-122.44",
minradiuskm = "0", maxradiuskm = "200",
format = "geojson")
- See exploration for an example of USGS data exploration with R.
- Results from the exploration are discussed here on abshinn.github.io.
Also within this repository is a shell script which uses curl to pull down the latest earthquakes from USGS's earthquake feed.
The simplest way to call the USGS feed is with the following curl command, which asks for earthquake events equal to or greater than 2.5 magnitude within the last 24 hours:
$ curl -s http://earthquake.usgs.gov/earthquakes/feed/v1.0/summary/2.5_day.csv
The script code/quakes.sh
uses the above curl, but allows a few options such as -s
for sorting by magnitude, and -H
for showing earthquake events within the last hour. For help, use the -h
option.
The script can be aliased within your respective shell profile such as ~/.bash_profile
, like so:
$ echo alias quakes="~/{wherever you cloned this repo}/usgs/code/quakes.sh" >> ~/.bash_profile
Then, whenever you would like to be up-to-date on earthquakes:
$ quakes -H
time latitude longitude depth mag place
2014-09-01T16:47:32.000Z 50.5447 -174.485 25.6 2.7 "184km S of Atka - Alaska"
Or, to see earthquake events that have occurred in California in the past 24 hours:
$ quakes | grep California
2014-09-01T16:10:37.400Z 36.9793 -121.464 5.4 3.5 "9km ESE of Gilroy - California"