メキシコ北部、カナダの一部、および米国本土の気候データを含むラスター レイヤーがあります。気候データを米国本土のゾーンのみに制限しようとしています。これを行うには、米国本土の地図をインポートして、ラスター データでオーバーレイするのは簡単だと思いました (もちろん、これは思ったよりもずっと難しいことがわかりました)。これまでのところ、rgdalライブラリを使用して、米国本土のマップと州の区分を含むシェープファイルをインポートすることができました。これをラスター レイヤーに変換して、最終的に最初の気候ラスター レイヤーにオーバーレイできるようにします。これは私が使用しているコードです:
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
私は何を間違っていますか?? どんなアドバイスでも大歓迎です!!!!