Given latitude and longitude values entered as decimal degrees (DD), this function returns a list (as an atomic vector) of station ID values, which can be used in get_GSOD to query for specific stations as an argument in the station parameter of that function.

nearest_stations(LAT, LON, distance)



Latitude expressed as decimal degrees (DD) (WGS84)


Longitude expressed as decimal degrees (DD) (WGS84)


Distance in kilometres from point for which stations are to be returned.


By default a class character vector object of station identification numbers. in order from nearest to farthest in increasing order. If return_full is TRUE, a data.table with full station metadata including the distance from the user specified coordinates is returned.


The GSOD data, which are downloaded and manipulated by GSODR stipulate that the following notice should be given. “The following data and products may have conditions placed on their international commercial use. They can be used within the U.S. or for non- commercial international activities without restriction. The non-U.S. data cannot be redistributed for commercial purposes. Re-distribution of these data by others must provide this same notification.”


Adam H. Sparks, [email protected]


if (FALSE) { # interactive() # Find stations within a 100km radius of Toowoomba, QLD, AUS n <- nearest_stations(LAT = -27.5598, LON = 151.9507, distance = 100) n }