csv读取后生成这样一个数据框
第一行，第一列相当于我的坐标值XY，矩阵相当于Z

tpx = unlist(data[1,-1])
tpy = data[,1][-1]
tpz = unlist(data[-1,-1])
xy = expand.grid(tpx, tpy)
xyz = data.frame(xy, tpz)


或者楼主的函数改成
matrix.grid<-function(matrix1)
{
  data.<-numeric()
  x<-numeric()
  y<-numeric()
  
  for(i in 2:dim(matrix1)[1])
  {
    for(j in 2:dim(matrix1)[2])
    {
      x<-c(x,matrix1[i, 1])
      y<-c(y,matrix1[1, j])
      data.<-c(data.,matrix1[i,j])
    }
  }
  data.frame(x,y,z=data.)
}

这个问题很简单，就是把矩阵转换成一一对应的data.frame。有reshape包里的melt函数即可解决，代码如下：

也可以这样
dt=read.csv("2line14dayrelable.csv",head=F)
dtrow=dt[-1,1]
dtcol=t(as.matrix(dt[1,-1]))
dtcore=stack(dt[-1,-1])[,1]
dtpos=merge(dtrow,dtcol)
outdata=cbind(merge(dtrow,dtcol),dtcore)
names(outdata)=c("x","y","z")

上面代码中 dtpos=merge(dtrow,dtcol) 这一句不需要，忘了删除