导入投影的GeoTIFF并转换为R中的纬度和经度坐标

Question

虽然我不是R的新手，但我不熟悉其中的空间数据，而不是ArcGIS。

我正在尝试从已经在Mercator中投影的公共访问GIS网站导入GeoTIFF图像（图像可以找到here）。使用raster中的R包，我可以成功将GeoTIFF导入R并绘制它（见下文）。

> library(raster)
> mb <- raster("/Maps/mb_srelief10m/mb_srelief10m.tif")
> mb
class       : RasterLayer 
dimensions  : 8060, 7220, 58193200  (nrow, ncol, ncell)
resolution  : 10, 10  (x, y)
extent      : -48598.54, 23601.46, 3850014, 3930614  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=merc +lon_0=-70.31666666666668 +lat_ts=41.65 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0 
data source : C:\Users\Connor\Documents\ArcGIS\Maps\mb_srelief10m\mb_srelief10m.tif 
names       : mb_srelief10m 
values      : 0, 255  (min, max)
attributes  :
       ID    COUNT BinValues Value
 from:  0  1460797         0     0
 to  : 51 14795988       255   255

如您所见，x轴和y轴标签表示光栅文件中的NAD83值（图像中的红色框）。然而，当我尝试使用crop函数仅使用此光栅文件的一部分时，我无法选择适当的区域，因为我的用户定义的extent对象是纬度和经度坐标。我知道我可以使用drawExtent函数手动选择区域，但我想自动执行此过程。

我的问题是：

如何变换光栅文件以使我的轴在纬度和经度上，我可以使用纬度和经度坐标裁剪文件？

如果这篇文章在其他地方更合适，请告诉我，我会删除该文件。谢谢！