2012-04-25 89 views
3

在R,說你希望在醫學成像如何覆蓋R中的等級圖?

example produced with Mango

每個層是不同的矩陣,以呈現圖像背景 - 上的統計數據等。 這裏是一個無義例如:

  • volcano在灰色
  • B<- (volcano>160) * rnorm(prod(dim(volcano))在 熱顏色。

你會怎麼樣?

+0

你熟悉保羅·馬雷爾這個JSS出版物,[在R圖形光柵圖像(http://journal.r-project.org/archive/2011-1/RJournal_2011 -1_Murrell.pdf)? – 2012-04-25 17:38:28

回答

3

image函數可以做到這一點,只需設置第二組顏色的透明度並使用add=TRUE覆蓋它。任何缺失的值都將是完全透明的,您可以爲其他顏色設置alpha級別(當然,您的圖形設備必須支持部分透明才能使其工作)。下面是一個簡單的例子:

A <- t(volcano)[ncol(volcano):1,] 
B <- (A>160) * rnorm(prod(dim(A))) 

B[ A < 120 ] <- NA # show some complete transparency 

image(A, col=grey(seq(0,1,length.out=12))) # initial plot 
hc <- sub('FF$','77',heat.colors(12)) # convert heat colors to have an alpha 
image(B, add=TRUE, col=hc) # overlay the new plot