R - 同じcrs、extent、dimension、resolutionのアライメントが合っていないラスタ

Question

メープルシロップの平均生産日数を求めています。私のメープル分布データはASCIIファイルにあります。 brick.Tmaxと呼ばれるラスタ（NetCDFファイルから作成）があります。私はbrick.Tmaxのスペックを私のメープル分布データにマッチさせたいと思っています。

## These are the specs I want to use for my maple distribution 
brick.Tmax 
class  : RasterBrick 
dimensions : 222, 462, 102564, 366 (nrow, ncol, ncell, nlayers) 
resolution : 0.125, 0.125 (x, y) 
extent  : -124.75, -67, 25.125, 52.875 (xmin, xmax, ymin, ymax) 
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
data source : E:\all_files\gridded_obs.daily.Tmax.1980.nc 
names  : X1980.01.01, X1980.01.02, X1980.01.03, X1980.01.04, X1980.01.05, X1980.01.06, X1980.01.07, X1980.01.08, X1980.01.09, X1980.01.10, X1980.01.11, X1980.01.12, X1980.01.13, X1980.01.14, X1980.01.15, ... 
Date  : 1980-01-01, 1980-12-31 (min, max) 
varname  : Tmax 

## reading in red maple data from ascii file into rasterLayer 
red_raster <- raster("E:/all_files/Maple_Data/redmaple.asc") 
red_raster 
class  : RasterLayer 
dimensions : 140, 150, 21000 (nrow, ncol, ncell) 
resolution : 20000, 20000 (x, y) 
extent  : -1793092, 1206908, -1650894, 1149106 (xmin, xmax, ymin, ymax) 
coord. ref. : NA 
data source : E:\all_files\Maple_Data\redmaple.asc 
names  : redmaple 
values  : -2147483648, 2147483647 (min, max) 

まだred_rasterの値を維持しながら、どのように私は、red_rasterへbrick.Tmaxからすべての仕様（寸法、CRS、resoluion、および範囲）を投影できますか？ 2つは互いに排他的であるようです。

注：私の質問を合理化しようとすると、元の投稿から質問を編集したので、下記のコメントが現在の状況で混乱している場合は謝ります。 （私は仲介者のように動作していたラスタprodavg_rastを削除しました）。

Answer 1

2つのラスタは、明らかに同じ程度にはありません。実際には、異なる世界にあります（座標参照系）。 brick.Tmaxは、角度（経度/緯度）座標が+proj=longlat +datum=WGS84 ですが、red_rasterは明らかにextent : -1793092, 1206908, -1650894, 1149106と指定されていません。したがって、これらのデータを一緒に使用するには、2つのうちの一方を変換する必要があります（もう一方の座標参照系に投影する）。問題は、red_rasterのcrsが何であるかわからないことです（esri asciiファイルはその情報を格納しません！）。したがって、データソースから何があるかを知る必要があります。また、対象となる領域や慣習を推測する必要があります。見つけたら、次のようなことをすることができます：

library(raster) 
tmax <- raster(nrow=222, ncol=462, xmn=-124.75, xmx=-67, ymn=25.125, ymx=52.875, crs="+proj=longlat +datum=WGS84") 

red <- raster(nrow=140, ncol=150, xmn=-1793092, xmx=1206908, ymn=-1650894, ymx=1149106, crs=NA) 
crs(red) <- " ??????  " 

redLL <- projectRaster(red, tmax) 

プロジェクションラスタに時間がかかります。あなたがcrsを見つけたかどうかをテストする良い方法は、物事が整列しているかどうかを示すことができるポリゴンを変換することです。