检查Lng和Lat是否属于ESRI Shape文件中的多边形

Question

我在R中尝试这段代码来检查哪个子区域具有特定的纬度，经度所在。想法是循环所有子区域名称并签入它们，而不是特定的子区域（在新加坡）。 wa.map $ SUBZONE_N给出了新加坡分区的名称。我在R中遇到以下错误：

  

（函数（classes，fdef，mtable）中的错误：无法找到   函数'proj4string＆lt; - '的继承方法用于签名   '＆＃34; data.frame＆＃34;，＆＃34;字符＆＃34;'

请帮助我如何做到这一点。

wa.map <- readOGR("C:/Users/somna/Desktop/master-plan-2014-subzone-boundary-   no-sea-shp/MP14_SUBZONE_NO_SEA_PL.shp", layer="MP14_SUBZONE_NO_SEA_PL")
dat <- data.frame(Longitude = c(103.83838894560753,103.83636995885735),
              Latitude =c(1.27744837946643, 1.27734776750189),
              names = c("Blk 3", "Bef Hosp Dr"))

for (i in 1:323)
{
  sodo <- wa.map[wa.map$SUBZONE_N[i],]

  proj4string(dat) <- proj4string(sodo)

  over(sodo, dat)

}