arcpy 生成等时圈isochron

#coding:utf-8
import json
from datetime import datetime
import time
import os
import urllib
import urllib2
import requests
import arcpy

time.sleep(0.5)

#119.422428,32.397402
x='119.422428'
y='32.397402'
mapbox_url='https://api.mapbox.com/isochrone/v1/mapbox/cycling/xxx%2Cyyy?contours_minutes=10&polygons=true&access_token='
token ='pk.eyJ1IjoieWFuc2MyMDIyIiwiYSI6ImNramF6c283ajJlNWMyeW5xenxxxx6qMhSItFGTw'
req_url=mapbox_url.replace('xxx',x).replace('yyy',y)+token

request=urllib2.Request(req_url)
response=urllib2.urlopen(request)
bodydata=response.read()
#print bodydata
result = json.loads(bodydata)  # 将字符串转换为json


spatialref=arcpy.SpatialReference(4326)
arcpy.env.workspace = r'E:	est	est2.gdb' #shp 执
fc=r'aaa'
cursor = arcpy.da.InsertCursor(fc, ['SHAPE@']) 
arr = arcpy.Array()
for i in result['features'][0]['geometry']['coordinates'][0]:
    arr.add(arcpy.Point(i[0],i[1]))
arr.remove(0)
polygon = arcpy.Polygon(arr, spatialref)
cursor.insertRow([polygon])

arr.removeAll()
del cursor

#https://account.mapbox.com/access-tokens

原文地址:https://www.cnblogs.com/yansc/p/14212669.html