加快坐标转换的循环代码

Speed up looping code for coordinate conversion

我有一个树干及其位置的数据集,最初是在绘图坐标中。我需要将这些坐标转换为 utm 坐标。我已经编写了以下函数来将我们的绘图坐标转换为 utm 坐标,并将以下脚本转换为整个数据集上的 运行。我遇到的问题是我总共有 261403 个词干,而脚本已经花费了难以置信的长时间 运行。我知道循环通常很慢,应该避免。关于如何在不使用循环的情况下执行此操作的任何建议?

定义函数:

    newcoords_x=function(gx,gy)
    {  
      rot_x <- gx*cos(-0.031989084) - gy*sin(-0.031989084)
      rot_y<- gx*sin(-0.031989084) + gy*cos(-0.031989084) 
      utm_x<-rot_x+625774
      utm_y<-rot_y+1011776  
      return(utm_x)
    }

    newcoords_y=function(gx,gy)
    { 
      rot_x <- gx*cos(-0.031989084) - gy*sin(-0.031989084)
      rot_y<- gx*sin(-0.031989084) + gy*cos(-0.031989084)  
      utm_x<-rot_x+625774
      utm_y<-rot_y+1011776
      return(utm_y)
    }

循环:

    for (i in 1:length(x$tag))
    {
      x$utm_x[i]=newcoords_x(x$gx[i], x$gy[i])
      x$utm_y[i]<-newcoords_y(x$gx[i], x$gy[i])
    }

这是一些示例数据,x,

    tag    gx    gy 
    2  994.1 488.3 
    4  990.5 488.9   
    6  993.5 498.3  
    7  992.7 469.3 
    8  981.9 473.5 
    13 983.0 452.6 

非常感谢!

最好的办法是向量化操作而不是使用循环。你可以这样做:

transform(x, utm_x=gx*cos(-0.031989084) - gy*sin(-0.031989084) + 625774
           , utm_y=gx*sin(-0.031989084) + gy*cos(-0.031989084) + 1011776)

#  tag    gx    gy    utm_x   utm_y
#1   2 994.1 488.3 626783.2 1012232
#2   4 990.5 488.9 626779.6 1012233
#3   6 993.5 498.3 626782.9 1012242
#4   7 992.7 469.3 626781.2 1012213
#5   8 981.9 473.5 626770.5 1012218
#6  13 983.0 452.6 626771.0 1012197