Как построить карту плотности в полярных координатах с помощью R?

Я хочу построить фигуру типа карты плотности в полярных координатах (мне не разрешено чтобы вставить цифры в мой пост, пожалуйста, перейдите по ссылке, чтобы увидеть это). Это карта плотности в полярных координатах. Я не знаком с R, поэтому даже я нашел этот пост, я до сих пор не знаю не знаю, как получить то, что я хочу.

Прямо сейчас у меня есть данные о разбросе в декартовых координатах. Буду признателен, если кто-нибудь сможет мне помочь. Большое спасибо.

======================== Обновление: мое решение ===================== =========

cart2pol <- function(x){
# x: (x,y)
y <- x[2]
x <- x[1]
r <- sqrt(x^2 + y^2)
t <- atan2(y,x)/pi*180
c(r,t)
}

angle <- apply(cockpit.data[c('x1','y1')],1,cart2pol)[2,]
r <- apply(cockpit.data[c('x1','y1')],1,cart2pol)[1,]
observations <-table(cut(angle,breaks=c(seq(-180,180,by=15))),cut(r,breaks=c(seq(0,sight,by=25))))

mm <- melt(observations,c('angle','r'))

labels <- seq(-172.5,172.5,length.out = 24) - 90
labels[labels<=0] <- labels[labels<=0] + 360
labels.y <- as.vector(rbind('', seq(0,sight,by=50)[-1]))
rosedensity <- ggplot(mm,aes(angle,r,fill=value))+geom_tile()+
    coord_polar(start=pi/2, direction = -1) + ggtitle('Rose Density') +
    scale_fill_gradientn(name="Frequency", colours = rev(rainbow(32)[1:23])) + #terrain.colors(100) , brewer.pal(10,'Paired')
    scale_x_discrete(labels = labels) + scale_y_discrete(labels = labels.y) +
    xlab('Angle') + ylab('R') +
    theme(
        plot.title = element_text(color="red", size=28, face="bold.italic"),
        axis.title.y = element_text(color="black", size=24, face="bold"),
        axis.title.x = element_text(color="black", size=24, face="bold"),
        axis.text=element_text(size=20),
        legend.justification=c(1,0), legend.position=c(1,0),
        legend.background = element_rect(fill="gray90", size=.1, linetype="dotted")
    )
ggsave(rosedensity, file=paste(dirOutput,'rosedensity.png',sep=''), width=8, height=8)

Вот мой выходной рисунок.

Я нашел решение из этот ответ.


r density-plot contour polar-coordinates
person Li Jiang    schedule 27.05.2016    source источник
comment
Добро пожаловать в StackOverflow! Чтобы мы могли вам помочь, я рекомендую вам сначала прочитать этот пост: stackoverflow.com/help/mcve. Затем разместите код, который вы написали, чтобы попытаться решить проблему.   -  person gariepy    schedule 27.05.2016

Ответы (1)


arrow_upward
2
arrow_downward

Вы можете попробовать это так:

library(ggplot2)
ggplot(faithful, aes(x = eruptions, y = waiting)) + 
  stat_density_2d(
    geom = "tile", 
    aes(fill = ..density..),
    n=c(40, 10), 
    contour = F
  ) + 
  scale_fill_gradientn(colours=rev(rainbow(32)[1:23])) + 
  coord_polar()

введите здесь описание изображения

person lukeA    schedule 27.05.2016