如何区分gplot中的背景颜色和NA值

Question

我正在生成湖泊叶绿素图。我想用蓝色填充湖泊，其中叶绿素浓度非常低，NA值为淡蓝色。我使用的是下面的代码

gplot(Chlorophyll_map_5) + geom_tile(aes(fill=value)) + scale_fill_gradient(low = 'blue', high = 'red', na.value='blue',name="Chl-a (ug/l)",limits=c(0,1000)) + coord_equal()+theme_bw() 

这让我对na.value='blue'这样的情节一样给出：

na.value='blue'

当我使用na.value='transparent'我得到这个图片：

na.value='transparent'

如果我改变na.value的颜色，它也会改变背景。有没有办法在不改变背景的情况下用颜色填充湖面？

我的数据的输出：`我的数据的输出：

Formal class 'RasterLayer' [package "raster"] with 12 slots 
..@ file :Formal class '.RasterFile' [package "raster"] with 13 slots 
.. .. ..@ name  : chr "/private/var/folders/68/hm_5ts9x7psb6j3wnb91_bfr0000gn/T/RtmpZ3BLZD/raster/r_tmp_2017-07-18_133827_28365_34843.grd"  
.. .. ..@ datanotation: chr "FLT8S" 
.. .. ..@ byteorder : Named chr "little" 
.. .. .. ..- attr(*, "names")= chr "value" 
.. .. ..@ nodatavalue : num -1.7e+308 
.. .. ..@ NAchanged : logi FALSE 
.. .. ..@ nbands  : int 1 
.. .. ..@ bandorder : Named chr "BIL" 
.. .. .. ..- attr(*, "names")= chr "value" 
.. .. ..@ offset  : int 0 
.. .. ..@ toptobottom : logi TRUE 
.. .. ..@ blockrows : int 0 
.. .. ..@ blockcols : int 0 
.. .. ..@ driver  : chr "raster" 
.. .. ..@ open  : logi FALSE 
..@ data :Formal class '.SingleLayerData' [package "raster"] with 13 slots 
.. .. ..@ values : logi(0) 
.. .. ..@ offset : num 0 
.. .. ..@ gain  : num 1 
.. .. ..@ inmemory : logi FALSE 
.. .. ..@ fromdisk : logi TRUE 
.. .. ..@ isfactor : logi FALSE 
.. .. ..@ attributes: list() 
.. .. ..@ haveminmax: logi TRUE 
.. .. ..@ min  : num 0.00335 
.. .. ..@ max  : num 3870657 
.. .. ..@ band  : int 1 
.. .. ..@ unit  : chr "" 
.. .. ..@ names  : chr "layer" 
..@ legend :Formal class '.RasterLegend' [package "raster"] with 5 slots 
.. .. ..@ type  : chr(0) 
.. .. ..@ values : logi(0) 
.. .. ..@ color  : logi(0) 
.. .. ..@ names  : logi(0) 
.. .. ..@ colortable: logi(0) 
..@ title : chr(0) 
..@ extent :Formal class 'Extent' [package "raster"] with 4 slots 
.. .. ..@ xmin: num 35.8 
.. .. ..@ xmax: num 36.7 
.. .. ..@ ymin: num 2.4 
.. .. ..@ ymax: num 4.65 
..@ rotated : logi FALSE 
..@ rotation:Formal class '.Rotation' [package "raster"] with 2 slots 
.. .. ..@ geotrans: num(0) 
.. .. ..@ transfun:function() 
..@ ncols : int 3240 
..@ nrows : int 8321 
..@ crs  :Formal class 'CRS' [package "sp"] with 1 slot 
.. .. ..@ projargs: chr "+proj=longlat +ellps=WGS84 +no_defs" 
..@ history : list() 
..@ z  : list() 

Answer 1

x <- trim(Chlorophyll_map_5) 
ggplot(Chlorophyll_map_5) + 
geom_tile(aes(fill=value)) + 
scale_fill_gradient(low = 'blue', high = 'red', na.value='blue',name="Chl-a (ug/l)",limits=c(0,1000)) + 
coord_equal()+theme_bw() 

每文档中，trim功能“通过去除外围行和列 只包含NA值作物RasterLayer “