ran.sam.plot=function(data,Gx,Gy,size,bin,Num)

{#data里面含有gx,gy这两列

#Gx,Gy分别是整个样地的两个边长

#size是随机抽取的小样方边长的大小，

#bin是设计抽样距离

#Num是随机取样的个数

X=seq(0,Gx-size,bin)

Y=seq(0,Gy-size,bin)

a=sample(X,Num,replace=T)

b=sample(Y,Num,replace=T)

result=list()

for(i in 1:Num)

{result[[i]]=data[data$gx<=(a[i]+size)&data$gx>=a[i]&data$gy<=(b[i]+size)&data$gy>=b[i],]

}

return(result)

}

ran.sam.plot(gts,600,400,10,5,10)