[javascript] 지점이 Google지도에서 육지 또는 물인지 확인

이 쿼리를 로컬에서 수행하는 것이 더 흥미 롭다고 생각했기 때문에 좀 더 자립 할 수 있습니다. 한 번에 25000 개의 임의의 토지 좌표를 생성하고 비용이 많이 드는 외부 API에 대한 호출을 피하고 싶습니다. 다음은 TomSchober가 언급 한 파이썬 예제를 사용하여 파이썬 으로 찍은 사진입니다. 기본적으로 모든 육지 좌표를 포함 하는 미리 만들어진 350MB 파일 에서 좌표를 조회 하고 좌표가 거기에 있으면 인쇄합니다.

import ogr
from IPython import embed
import sys

drv = ogr.GetDriverByName('ESRI Shapefile') #We will load a shape file
ds_in = drv.Open("land_polygons.shp")    #Get the contents of the shape file
lyr_in = ds_in.GetLayer(0)    #Get the shape file's first layer

#Put the title of the field you are interested in here
idx_reg = lyr_in.GetLayerDefn().GetFieldIndex("P_Loc_Nm")

#If the latitude/longitude we're going to use is not in the projection
#of the shapefile, then we will get erroneous results.
#The following assumes that the latitude longitude is in WGS84
#This is identified by the number "4236", as in "EPSG:4326"
#We will create a transformation between this and the shapefile's
#project, whatever it may be
geo_ref = lyr_in.GetSpatialRef()
point_ref=ogr.osr.SpatialReference()
point_ref.ImportFromEPSG(4326)
ctran=ogr.osr.CoordinateTransformation(point_ref,geo_ref)

def check(lon, lat):
    #Transform incoming longitude/latitude to the shapefile's projection
    [lon,lat,z]=ctran.TransformPoint(lon,lat)

    #Create a point
    pt = ogr.Geometry(ogr.wkbPoint)
    pt.SetPoint_2D(0, lon, lat)

    #Set up a spatial filter such that the only features we see when we
    #loop through "lyr_in" are those which overlap the point defined above
    lyr_in.SetSpatialFilter(pt)

    #Loop through the overlapped features and display the field of interest
    for feat_in in lyr_in:
        # success!
        print lon, lat

check(-95,47)

나는 12 개의 코디네이터를 시도했는데 훌륭하게 작동합니다. “land_polygons.shp”파일은 여기 에서 다운로드 할 수 있으며 OpenStreetMaps를 보완합니다. (첫 번째 WGS84 다운로드 링크를 직접 사용했으며 두 번째도 작동 할 수 있습니다)

답변