0

我有一个栅格图层，其中包含来自墨西哥北部、加拿大部分地区和美国大陆的气候数据。我试图将气候数据限制在美国大陆地区。为此，我认为导入美国大陆地图并用我的栅格数据覆盖它会很容易（当然这比我想象的要困难得多）。到目前为止，使用该rgdal库，我已经成功导入了一个 shapefile，其中包括美国大陆地图和各州分部。现在我想将其转换为栅格图层，以便最终将其与我的初始气候栅格图层叠加。这是我正在使用的代码：

setwd ("C:/Climate_data/USA map")
ogrInfo(".", "USA_mainland_states")
usa_state = readOGR(dsn=".", layer="USA_mainland_states")

##Convert to Raster 
r_usa_state <- raster()
extent(r_usa_state) <- extent(usa_state)
rasterize(usa_state,r_usa_state, fun='last')
overlay (r_usa_state, sms_av, fun='mask')

但是，现在我收到以下错误：

Error in .readCells(x, cells, 1) : no data on disk or in memory

sms_av是气候栅格图层（103936 个元素，823.3 kb）。

另外，当我这样做时：

hasValues(r_usa_state)

我得到：

FALSE

我究竟做错了什么？？任何建议将不胜感激！！！！

4

0 回答 0