gdalwarpとprojectRasterの違い

ラスターを投影しようとしています。RにはprojectRaster()これを行うための関数があります（完全に再現可能な例の下に）：

# example Raster
require(raster)
r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40)
r <- setValues(r, 1:ncell(r))
projection(r)
# project to
newproj <- "+init=epsg:4714"


# using raster package to reproject
pr1 <- projectRaster(r, crs = CRS(newproj), method = 'bilinear')

それはうまくいきます。しかし、それはかなり遅いです。

速度を上げるために、gdalwarp代わりに使用します（SSD を使用すると、ディスク/ Rへの読み書きのコストはそれほど高くありません）。

ただし、をprojectRaster()使用した結果を再現できませんgdalwarp。

# using gdalwarp to reproject
tf <- tempfile(fileext = '.tif')
tf2 <- tempfile(fileext = '.tif')
writeRaster(r, tf)
system(command = paste(paste0("gdalwarp -t_srs \'", newproj, "\' -r bilinear -overwrite"), 
                       tf,
                       tf2))
pr2 <- raster(tf2)

動作するようですが、結果は異なります。

# Info
system(command = paste("gdalinfo", 
                       tf))
system(command = paste("gdalinfo", 
                       tf2))

# plots
plot(r)
plot(pr1)
plot(pr2)

#extents
extent(r)
extent(pr1)
extent(pr2)

# PROJ4
proj4string(r)
proj4string(pr1)
proj4string(pr2)

# extract value
take <- SpatialPoints(matrix(c(-100, 50), byrow = T, ncol = 2), proj4string = CRS(newproj))
plot(take, add = TRUE)
extract(pr1, take)
extract(pr2, take)

何が欠けている/間違っていますか？

他の（より速い）代替手段はありprojectRaster()ますか？

再現性のあるすばらしい質問です。個人的には、違いの理由は双線形再投影の実装にあると思います。あなたは明らかに2つのアプローチのソースコードを調べることができますが、私はそれが非常に過剰であることを期待します。
Rの実装では、生のGDALバージョンよりも大きな「エラー」/「変更」が導入されているようです（私のバージョンとテストでは少なくとも-projectRasterは+ -0.01前後の変更を導入し、GDALは+ -0.002前後の値を提供します）。

最近傍の再投影を使用して両方のアプローチを比較すると、期待どおりに一致します。