R

Question

我正在尝试在栅格顶部绘制shapefile。坐标系都相同（重新投影之后），但是shapefile在栅格顶部稍微移动了。参见下图。

这是我的代码：

CB <- readOGR("S:/New_CB_shapefile/cb_shapefile_islands.shp")
> CB
class       : SpatialPolygonsDataFrame 
features    : 1 
extent      : -77.46415, -75.39099, 36.71462, 39.60927  (xmin, xmax, ymin, ymax)
crs         : +proj=utm +zone=18 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
variables   : 1
names       : dissolve 
value       :        1 
> 
rst <- raster("201605.tif")
> rst
class      : RasterLayer 
band       : 1  (of  11  bands)
dimensions : 1040, 600, 624000  (nrow, ncol, ncell)
resolution : 0.003514633, 0.00268125  (x, y)
extent     : -77.37981, -75.27103, 36.87297, 39.66147  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
source     : E:/S3_Backup_Anna/S3_L3/c2rcc/201605.tif 
names      : X201605 
values     : 0, 0.07534551  (min, max)

> 
CB_reproject <- spTransform(CB, CRS("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"))
> CB_reproject
class       : SpatialPolygonsDataFrame 
features    : 1 
extent      : -79.48944, -79.48942, 0.0003311427, 0.0003572506  (xmin, xmax, ymin, ymax)
crs         : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
variables   : 1
names       : dissolve 
value       :        1 
> 

plot(raster("201605.tif"))
plot(CB, add=TRUE, lwd=.1)

在转换过程中，我的ymin和ymax好像发生了变化？这在ArcMap或ESA的SNAP中不会发生，因此我不确定为什么在R中会发生这种情况。有人有什么建议吗？