how to get full active sensors id table

[TonyFly3000]

May I ask how to get full active sensors id table include senior id and city name?So that we could search certain city measurement data.

I try to use locations API to find sensors, but its require city coordinates and the active sensors may be outside city center 25000 radius range.

[russbiggs]

There is no construct of cities within the API, because this is highly subjective geographically. Instead of using the radius and coordinate method I recommend the bounding box parameter: https://docs.openaq.org/using-the-api/geospatial#bounding-box
This way you can define a bounding box that will fit around arbitrary geographic shapes, such as a city boundary. in practice this would look like:

locations = client.locations.list(limit = 1000, bbox=[-119.115567,33.514103,-117.138028,34.344752])

This bounding box covers most of Los Angeles.

Given a GeoJSON or shapefile of a city boundary theres an abundant number of libraries and tools that can help define the bounding box.

[TonyFly3000]

Thanks. I found there is a new R package [openaq]. I am trying to use it instead.

I want get daily PM2.5 data for a city(eg Los Angeles). There is how I do it.

But I think its may not be the best way. Do you have some recommendation:

load pacakge

library(jsonlite)
library(curl)
library(tidyverse)
library(janitor)
library(ggplot2)
library(plotly)
library(openaq)
library(tidygeocoder)

Step 1 get city box coordinate

find Los Angeles latitude and longitude

lat_longs <-tibble::tribble(~addr,"Los Angeles") %>%geocode(addr, method = 'osm', lat = latitude , long = longitude)

make a box

box_range=list(lat_longs$longitude-0.3,lat_longs$latitude-0.3,lat_longs$longitude+0.3,lat_longs$latitude+0.3)

Step 2 get location which have PM2.5 data at past 10 day

location001=list_locations(bbox=box_range,parameters_id=2)

location002=location001 %>% filter(is_monitor==TRUE, datetime_last %>% as.POSIXlt() %>% as.Date()>=Sys.Date()-10)

Step 3 get sensor with first location

using first location within box coordinate

sensor001=(location002 %>% head(1))$id %>% get_location_sensors() 

filter pm2.5 sensor

sensor002= sensor001 %>% filter(name=='pm25 µg/m³')

Step 4 get daily PM2.5 measurement data from sensor since 2024-01-01

measurement001=list_sensor_measurements(sensors_id=sensor002$id,datetime_from="2024-01-01"%>% as.POSIXct()
                                        ,datetime_to=Sys.Date()%>% as.POSIXct()
                                    ,limit=1000,rollup='daily'
                                        )