请教各位大神：
distance <- c()
for (i in 1:nrow(cz_bts_3g))
{
  for (j in 1:nrow(cz_yyt_wz))
  {
    dist <- round(distm(cz_bts_3g[i,9:10],cz_yyt_wz[j,4:5],fun = distVincentyEllipsoid),0)
    b <- data.frame(cz_bts_3g[i,1:10],cz_yyt_wz[j,2],dist)
    distance <- rbind(distance,b)
  }
}
write.csv(distance,"C:/Users/hhx/Desktop/jz_4g_alter/cz_3g_yyt_dist.csv")
现在有cz_bts_3g，这是基站的表，包含经纬度信息，cz_yyt_wz这个是营业厅的表，也包含经纬度信息，我想操作的就是求出每一个基站和每一个营业厅之间的距离，然后再求出距离这个基站最近的营业厅是哪个，以及他们的距离。我写了一个for循环，但是数据量大的话（30000*500，也就是1500万条数据），根本跑不动，所有想请教各位大神，有没有快速的方法。

跪求各位大神。。。。。