R: Точка поиска в полигоне часового пояса

Question

Я написал функцию, чтобы определить, попадает ли данная точка lat/lon в определенный полигон из файла tz_world.shp. Я перехожу в lon = -77, lat = 42, который находится в штате Нью-Йорк. Идентификатор из 155-й записи - «Америка/Нью-Йорк». Но использование point.inpoly возвращает объект длины с нулевой строкой. Вот мой код. Пожалуйста помоги. Благодарю. 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) 

Экран-вид на типично показывает координатам ограничивающего LON/лат точек для NYS, которые должны окружают мой комплект поставки LON/лат точку.

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>