Có hai điều sẽ tác động đến độ mượt của cốt truyện, băng thông được sử dụng cho ước tính mật độ hạt nhân của bạn và các khoảng nghỉ bạn gán màu cho trong ô.
Theo kinh nghiệm của tôi, để phân tích khám phá, tôi chỉ điều chỉnh băng thông cho đến khi tôi có được một âm mưu hữu ích. Trình diễn dưới đây.
library(spatstat)
set.seed(3)
X <- rpoispp(10)
par(mfrow = c(2,2))
plot(density(X, 1))
plot(density(X, 0.1))
plot(density(X, 0.05))
plot(density(X, 0.01))
Chỉ cần thay đổi lược đồ màu mặc định sẽ không giúp ích gì cả, cũng sẽ không thay đổi độ phân giải của pixel (nếu bất cứ điều gì độ phân giải mặc định quá chính xác và bạn nên giảm độ phân giải và làm cho pixel lớn hơn). Mặc dù bạn có thể muốn thay đổi bảng màu mặc định cho mục đích thẩm mỹ, nhưng nó có ý định phân biệt đối xử cao.
Những điều bạn có thể làm để giúp màu sắc thay đổi mức độ tỷ lệ thành logarit (sẽ thực sự chỉ có ích nếu bạn có một quy trình rất không đồng nhất), thay đổi bảng màu để thay đổi nhiều hơn ở cấp thấp hơn (sai lệch về đặc điểm kỹ thuật của dải màu trong R), hoặc điều chỉnh chú giải để có các thùng rời rạc thay vì liên tục.
Các ví dụ về sự thiên vị trong truyền thuyết được điều chỉnh từ đây và tôi có một bài đăng khác trên trang web GIS giải thích việc tô màu các thùng rời rạc trong một ví dụ khá đơn giản ở đây . Mặc dù vậy, những điều này sẽ không giúp ích nếu mô hình kết thúc hoặc được làm mịn mặc dù bắt đầu bằng.
Z <- density(X, 0.1)
logZ <- eval.im(log(Z))
bias_palette <- colorRampPalette(c("blue", "magenta", "red", "yellow", "white"), bias=2, space="Lab")
norm_palette <- colorRampPalette(c("white","red"))
par(mfrow = c(2,2))
plot(Z)
plot(logZ)
plot(Z, col=bias_palette(256))
plot(Z, col=norm_palette(5))
Để làm cho màu trong suốt trong hình ảnh cuối cùng (trong đó thùng màu đầu tiên là màu trắng), người ta có thể chỉ cần tạo đoạn màu và sau đó thay thế thông số RGB bằng màu trong suốt. Ví dụ dưới đây sử dụng cùng một dữ liệu như trên.
library(spatstat)
set.seed(3)
X <- rpoispp(10)
Z <- density(X, 0.1)
A <- rpoispp(100) #points other places than density
norm_palette <- colorRampPalette(c("white","red"))
pal_opaque <- norm_palette(5)
pal_trans <- norm_palette(5)
pal_trans[1] <- "#FFFFFF00" #was originally "#FFFFFF"
par(mfrow = c(1,3))
plot(A, Main = "Opaque Density")
plot(Z, add=T, col = pal_opaque)
plot(A, Main = "Transparent Density")
plot(Z, add=T, col = pal_trans)
pal_trans2 <- paste(pal_opaque,"50",sep = "")
plot(A, Main = "All slightly transparent")
plot(Z, add=T, col = pal_trans2)