如何使用传单过滤光栅图层并仅打印R中某个值以上的单元格？

我有一个光栅图层，它包含了全美24小时的积雪。可以从这里提取数据：

我只想在传单地图上绘制光栅中值大于或等于4（英寸）的单元。这是我当前地图的外观：

我将所有小于4的值更改为NA，认为光栅单元不会显示在地图上。我想删除地图上所有变灰的单元格。函数subset（）和filter（）不适用于光栅图层。有什么想法吗？我的代码如下，以供参考：

library(dplyr)
library(rgdal)
library(raster)
library(ncdf4)
library(leaflet)
library(leaflet.extras)

download.file(obsvSnow_Link, destfile = file.path(folderpath, 'observedSnow.tif'))
obsvSnow <- raster(file.path(folderpath, 'observedSnow.tif'))
names(obsvSnow) <- 'snowfall'

obsvSnow[obsvSnow < 4,] <- NA 
          
colores <- c("transparent","#99CCFF","#3399FF","#0000FF","#FFE066", "#FF9900", "#E06666","#CC0000","#990033")
at <- c(4,8,seq(12,42,6),100)
cb <- colorBin(palette = colores, bins = at, domain = at)
          
mp <- leaflet(width = "100%",options = leafletOptions(zoomControl = FALSE)) %>%
          addTiles() %>%
          addRasterImage(x=obsvSnow$snowfall, 
                           colors = cb, 
                           opacity = 0.6) %>%
          addLegend(title = 'Inches',
                    position='bottomright',
                    pal = cb, values = at) %>%
          leaflet.extras::addSearchUSCensusBureau(options = searchOptions(autoCollapse=TRUE, minLength=10)) %>%
          addScaleBar(position='bottomleft') %>%
          addFullscreenControl()
          
mp

库（dplyr）
图书馆（rgdal）
图书馆（光栅）
图书馆（ncdf4）
图书馆（单张）
图书馆（单张、附加资料）
下载.file（obsvSnow_链接，destfile=file.path（folderpath，'observedSnow.tif'））
obsvSnow请提供一个最小的、自包含的、可复制的示例。也就是说，使用R附带的一些数据，或者使用代码创建的数据。