R : 표준 시간대 폴리곤의 찾아 내기

Question

주어진 lat/lon 점이 tz_world.shp 파일의 특정 다각형 내에 있는지를 결정하는 함수를 작성했습니다. 나는 뉴욕 주에있는 lon = -77, lat = 42로 전달합니다. 155 번째 항목의 ID는 'America/New_York'입니다. 그러나 point.inpoly를 사용하면 0 행 길이의 객체가 반환됩니다. 여기 내 코드가있다. 도와주세요. 감사. BSL.

i = 155 

timeZoneList = split(timeZonesShpFile, timeZonesShpFile$TZID) 

ps = lapply(timeZoneList, Polygon) 

p1 = lapply(seq_along(ps), function(i) Polygons(list(ps[[i]]), ID = names(timeZoneList)[i])) 

my_spatial_polys = SpatialPolygons(p1, proj4string = CRS("+proj=longlat +datum=WGS84")) 

polyNames = sapply(slot(my_spatial_spdf, 'polygons'), function(i) slot(i, 'ID')) 

pt = SpatialPointsDataFrame(cbind(lon,lat), data.frame(row=1), proj4string=CRS("+proj=longlat +datum=WGS84")) 

thisPoly = my_spatial_polys[ i ] 

thisList = getSpPPolygonsIDSlots(thisPoly) 

thisDf = data.frame(row = 1, row.names = thisList) 

thisSpdf = SpatialPolygonsDataFrame(thisPoly, thisDf) 

pIp = point.in.poly(pt, thisSpdf) 

> pIp 

[1] coordinates row   row.1  
<0 rows> (or 0-length row.names) 

coords의 화면보기는 제공되는 lon/lat 포인트를 둘러싸는 NYS의 일반적인 경계 l/lat 점을 보여줍니다.

Answer 1

여기 다른 접근법 : 위에

library(maptools) 
library(rgdal) 
library(rgeos) 
download.file("http://efele.net/maps/tz/world/tz_world.zip", tf <- tempfile(fileext = ".zip")) 
unzip(tf, exdir = tempdir()) 
shp <- readShapeSpatial(file.path(tempdir(), "world", "tz_world.shp")) 
over(
    SpatialPoints(data.frame(lon = c(13, -77, 51), lat = c(52, 42, 0))), # Berlin, New York, London 
    shp[shp$TZID %in% c("America/New_York", "Europe/Berlin"), ] 
) 
#    TZID 
# 1 Europe/Berlin 
# 2 America/New_York 
# 3    <NA>