raster::sampleRandom采样超出掩膜区域（含NA）问题求助

我在构建物种分布模型时，使用WorldClim和SoilGrids的30秒分辨率环境预测因子，通过raster包的sampleRandom函数采样背景点，遇到了异常情况：多数采样点落在海洋区域（即掩膜外的NA值区域）。

我的采样代码如下：

bg<-as.data.frame(sampleRandom(env.pred[[1]], size=10000, na.rm=F, xy=TRUE))

其中env.pred是已经完成掩膜处理的RasterStack，其对应对象的信息如下：

> pred.mask
class      : RasterBrick 
dimensions : 21600, 43200, 933120000, 24  (nrow, ncol, ncell, nlayers)
resolution : 0.008333333, 0.008333333  (x, y)
extent     : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : r_tmp_2024-05-03_142200_7224_49931.grd 
names      :       bio1,      bio10,      bio11,      bio12,      bio13,      bio14,      bio15,      bio16,      bio17,      bio18,      bio19,       bio2,       bio3,       bio4,       bio5, ... 
min values :   9.283334,  11.466667,   6.466667, 502.000000,  58.000000,   0.000000,   5.105296, 157.000000,   1.000000,  79.000000,   1.000000,   5.566667,  36.323528,  86.614578,  16.600000, ... 
max values :   26.15417,   27.68333,   24.81667, 3249.00000,  403.00000,  163.00000,  103.06243, 1157.00000,  549.00000, 1157.00000,  841.00000,   14.12500,   78.41881,  465.92282,   33.30000, ... 

当我设置na.rm=T时，仅能得到61个采样点，但我原本期望所有采样点都位于掩膜区域内。有没有人遇到过这种情况？

内容的提问来源于stack exchange，提问作者AmandaKM