Форумы GIS-Lab.info

Denis Rykov писал(а):403

nickleb писал(а):исправил ссылку

Максим Дубинин писал(а):

nickleb писал(а):исправил ссылку

нет

gamm писал(а):

Код: Выделить всё

my.ellipse<-function(x0,y0,a,b,ell.angle,ell.segm=48,print=TRUE) {
  # horisonl ellipse
  xy<-matrix(ncol=2,nrow=ell.segm+1)
  for(i in 1:ell.segm) {
    ang<-2*pi/ell.segm
    x<-cos((i-1)*ang)*a
    y<-sin((i-1)*ang)*b
    xy[i,]<-c(x,y)
  }
  xy[ell.segm+1,]<-xy[1,]
  
  # turn ellipse
  ell.sin<-sin(ell.angle/180*pi)
  ell.cos<-cos(ell.angle/180*pi)
  mtr<-matrix(c(ell.cos,ell.sin,-ell.sin,ell.cos),ncol=2,nrow=2)
  xy1<-t((mtr %*% t(xy)) + c(x0,y0))
  if(print) points(xy1,type="l",lwd=3)
  invisible(xy1)
}

pos<-expand.grid(1:10*10,1:5*10)
plot(0,0,xlim=c(0,105),ylim=c(0,55),asp=1,type="n")
for(i in 1:nrow(pos)) {
  x0<-pos[i,1]        # center position 
  y0<-pos[i,2]        # center position 
  a<-runif(1,0.1,5)   # "main" semiaxis
  b<-runif(1,0.1,5)   # "secondary" semiaxis
  ang<-runif(1,0,360) # angle
  my.ellipse(x0,y0,a,b,ang,ell.segm=48,print=TRUE)
}

главная полуось горизонтальна, поворот - против часовой от горизонтали.