光栅:如何使用python获得lat/long的高程

Question

我也在SO的GIS部分发布了这个问题。由于我不确定这是否是一个"纯粹"的python问题，我也在这里再次提问。

我想知道是否有人在不使用ArcGIS的情况下从光栅获取高程数据，而是将信息作为python list或dict获取？

我将XY数据作为元组列表获取。

我想循环遍历列表，或者将其传递给函数或类方法，以获得xy对的相应提升。

我对这个主题做了一些研究，gdalAPI听起来很有前途。有人能建议我如何处理事情、陷阱和示例代码吗？其他选择？

感谢您的努力，LarsVegas

Answer 1

我建议查看Google Elevation API

使用起来非常简单：

http://maps.googleapis.com/maps/api/elevation/json?locations=39.7391536，-104.9847034&传感器=真_假

{
   "results" : [
      {
         "elevation" : 1608.637939453125,
         "location" : {
            "lat" : 39.73915360,
            "lng" : -104.98470340
         },
         "resolution" : 4.771975994110107
      }
   ],
   "status" : "OK"
}

请注意，免费版本限制为每天2500次请求。

Answer 2

我们使用此代码来获取给定纬度/经度的高程（注意：我们只要求打印高程以及舍入的纬度和经度值）。

import urllib.request
import json
lati = input("Enter the latitude:")
lngi = input("Enter the longitude:")
# url_params completes the base url with the given latitude and longitude values
ELEVATION_BASE_URL = 'http://maps.googleapis.com/maps/api/elevation/json?'
URL_PARAMS = "locations=%s,%s&sensor=%s" % (lati, lngi, "false")
url=ELEVATION_BASE_URL + URL_PARAMS
with urllib.request.urlopen(url) as f:
response = json.loads(f.read().decode())    
status = response["status"]
result = response["results"][0]
print(float(result["elevation"]))
print(float(result["location"]["lat"]))
print(float(result["location"]["lng"]))

Answer 3

看看高度计，它是Google Elevation API 的包装器

Answer 4

下面是我构建的另一个很好的API：https://algorithmia.com/algorithms/Gaploid/Elevation

import Algorithmia
input = {
    "lat": "50.2111",
    "lon": "18.1233"
}
client = Algorithmia.client('YOUR_API_KEY')
algo = client.algo('Gaploid/Elevation/0.3.0')
print algo.pipe(input)