First it loops through all of the features in our stream gauges feature class (SGID93.WATER.StreamGaugesNHD).
For each feature, it uses a USGS id (SourceFeature_ID) to build a url to hit their Instantaneous Values web service.
# get json objectdata = json.loads(urllib2.urlopen(r'http://waterservices.usgs.gov/nwis/iv?format=json&site=' + id).read())
This is an example of one of the urls: http://waterservices.usgs.gov/nwis/iv?format=json&site=09413700. It then uses the json library to parse the data and get the values that we are interested in. These values are used to populate the appropriate fields in our feature class.
def getJsonValue(variableCode, data): for ts in data['value']['timeSeries']: if ts['variable']['valueType'] == variableCode: value = ts['values'][0]['value'][0]['value'] return value
The NOAA data is served up via an rss feed which means xml. The minidom object from the xml.dom library came in handy here for parsing the xml data.
# get noaa data gaugeID = row.getValue('GuageID') if gaugeID: ndata = minidom.parse(urllib2.urlopen('http://water.weather.gov/ahps2/rss/fcst/' + gaugeID.lower() + '.rss')) descriptionText = ndata.getElementsByTagName('description')[2].firstChild.nodeValue descriptionList = descriptionText.split('<br />') row.setValue('HIGHEST_FORECAST', descriptionList[5].split()[2].strip()) row.setValue('HIGHEST_FORECAST_DATE', getNOAADate(descriptionList[6].split('Time:')[1].strip())) row.setValue('LAST_FORECAST', descriptionList[8].split()[2].strip()) row.setValue('LAST_FORECAST_DATE', getNOAADate(descriptionList[9].split('Time:')[1].strip()))
Hi Scott,
ReplyDeleteCan you possibly re-host the full script? I am looking to do something similar and it would be helpful if I can look at your example.
Thanks!
Judson
Hey Judson,
DeleteThat script is long gone. Sorry! Did you have any specific questions?