Python：從緯度和經度值獲取高程

Question

首先，這個問題有代碼。 但是，它似乎只適用於美國。 正如他們的網站 ( https://nationalmap.gov/epqs/ ) 上所述，當無法找到值時，它將返回 -1000000，這就是我的情況。

我試圖在瑞典獲得海拔。

使用此示例數據：

lat = [57.728905, 57.728874, 57.728916, 57.728836, 57.728848]
lon = [11.949309, 11.949407, 11.949470, 11.949342, 11.949178]

# create df
df = pd.DataFrame({'lat': lat, 'lon': lon})

如何使下面的代碼適用於全球范圍而不僅限於美國？ 或者至少，有什么方法可以讓我把重點放在瑞典上？ 美國地質調查局不應該有全球數據嗎？

def make_remote_request(url: str, params: dict):
   """
   Makes the remote request
   Continues making attempts until it succeeds
   """

   count = 1
   while True:
       try:
           response = requests.get((url + urllib.parse.urlencode(params)))
       except (OSError, urllib3.exceptions.ProtocolError) as error:
           print('\n')
           print('*' * 20, 'Error Occured', '*' * 20)
           print(f'Number of tries: {count}')
           print(f'URL: {url}')
           print(error)
           print('\n')
           count += 1
           continue
       break

   return response


def elevation_function(x):
   url = 'https://nationalmap.gov/epqs/pqs.php?'
   params = {'x': x[1],
             'y': x[0],
             'units': 'Meters',
             'output': 'json'}
   result = make_remote_request(url, params)
   return result.json()['USGS_Elevation_Point_Query_Service']['Elevation_Query']['Elevation']

要運行該函數：

df['elevation'] = df.apply(elevation_function, axis=1)


# Print output
df

源代碼： 使用簡單的python腳本從緯度經度坐標獲取高程

更新：

使用已接受答案中的建議並在高程函數中添加 time.sleep(1) 允許我獲得所有觀測的高程數據。 請注意，此 API 僅允許每個請求 100 個位置。 因此，對於更多的行，必須將其分解為不同的請求。

Answer 1

我一直在使用opentopodata.org的 api 來獲取高程值。 您可以使用涵蓋整個瑞典的EU-DEM數據集。 API 請求非常簡單，如下所示：

https://api.opentopodata.org/v1/eudem25m?locations=57.728905,11.949309

另一個高程 API 可以在https://open-elevation.com/找到。 一個請求看起來非常相似：

https://api.open-elevation.com/api/v1/lookup?locations=57.728905,11.949309

因此調整您的elevation_function ：

def elevation_function(x):
   url = 'https://api.opentopodata.org/v1/eudem25m?'
   # url = 'https://api.open-elevation.com/api/v1/lookup?'
   params = {'locations': f"{x[0]},{x[1]}"}
   result = make_remote_request(url, params)
   return result.json()['results'][0]['elevation']