Jak wykreślić za pomocą png jako tło?

Question

Zrobiłem fabułę z 3 milionami punktów i zapisałem ją jako PNG. Zajęło mi to kilka godzin i chciałbym uniknąć ponownego rysowania wszystkich punktów.

Jak mogę wygenerować nowy spisek, który ma tę PNG jako tło?

Answer 1

Spróbuj tego:

library(png) 

#Replace the directory and file information with your info 
ima <- readPNG("C:\\Documents and Settings\\Bill\\Data\\R\\Data\\Images\\sun.png") 

#Set up the plot area 
plot(1:2, type='n', main="Plotting Over an Image", xlab="x", ylab="y") 

#Get the plot information so the image will fill the plot box, and draw it 
lim <- par() 
rasterImage(ima, lim$usr[1], lim$usr[3], lim$usr[2], lim$usr[4]) 
grid() 
lines(c(1, 1.2, 1.4, 1.6, 1.8, 2.0), c(1, 1.3, 1.7, 1.6, 1.7, 1.0), type="b", lwd=5, col="white") 

Poniżej znajduje się działka.

Answer 2

Podczas użytkownika @ bill_080 odpowiedź bezpośrednio odpowiedzi na pytanie, czy to jest naprawdę to, co chcesz? Jeśli chcesz na to spisać, musisz dokładnie wyrównać swoje układy współrzędnych. Zobacz np. Houston Crime Map jak to zrobić z ggplot2.

Dla twojego problemu wydaje mi się, że może być łatwiejsze rozwiązanie: binning, czyli kończenie 2d histogramów.

> df <- data.frame (x = rnorm (1e6), y = rnorm (1e6)) 
> system.time (plot (df)) 
     User  System verstrichen 
    54.468  0.044  54.658 
> library (hexbin) 
> system.time (binned <- hexbin (df, xbins=200)) 
     User  System verstrichen 
     0.252  0.012  0.266 
> system.time (plot (binned)) 
     User  System verstrichen 
     0.704  0.040  0.784 

hexbin współpracuje bezpośrednio z siatki i ggplot2, ale centrum współrzędne pojemniki są w binned@xcm i binned@ycm, więc można też wykreślić wynik w grafice bazowych. Z dużej liczby pojemników, można uzyskać szybki wersji oryginalnej fabuły:

> system.time (plot (binned@xcm, binned@ycm, pch = 20, cex=0.4)) 
     User  System verstrichen 
     0.780  0.004  0.786 

ale łatwo można mieć kolory kodujące gęstość:

> plot (binned@xcm, binned@ycm, pch = 20, cex=0.4, col = as.character (col)) 

> col <- cut (binned@count, 20) 
> levels (col) <- grey.colors (20, start=0.9, end = 0) 
> plot (binned@xcm, binned@ycm, pch = 20, cex=0.4, col = as.character (col))