问题描述
我对使用空间数据非常陌生，所以我要说的大部分内容都是我尝试说一门外语。现在我正在尝试学习如何在 R 中完成这一切（我在 QGIS 中处理这些数据的能力稍强一些，但对于这个解决方案，我只在寻找 R）。
我的研究涉及宾夕法尼亚州 (PA) 的生态数据，因此我正在尝试将美国 NLCD 数据集裁剪为宾夕法尼亚州。我有 NLCD 的栅格图层和宾夕法尼亚州边界的 shapefile。我能够成功地将较大的美国栅格裁剪为 PA，如下所示：
library(raster)
library(rgdal)
pabound <- readOGR(dsn="...",layer="PAbound")
nlcdRast <- raster(".../NLCD_2016_Land_Cover_L48_20190424.img")
pabound <- spTransform(pabound,CRS(proj4string(nlcdRast)))
PAnlcd <- raster::crop(nlcdRast,pabound)
如果我为 nlcdRast 和 PAnlcd 运行简单的绘图命令（即 plot(nlcdRast)，它们保持相同的配色方案。但是当我通过 tmap 运行它时，它似乎以不同的方式查看裁剪后的数据，但我不确定如何解决这个问题。请参阅下面的图表：
library(tmap)
tm_shape(nlcdRast) +
tm_raster()
然后当我在 tmap 中绘制裁剪后的版本时：
tm_shape(PAnlcd) +
tm_raster()
如您所见，改变的不仅仅是调色板（我相信我能解决这个问题），但真正的问题是我丢失了图例中的重要信息。尽管完整图实际上显示了栅格 NLCD 的分类值，但裁剪后的版本现在似乎只显示了一些未知的数值范围。尽管目前看起来很糟糕，但我希望拥有与美国完整地图中看到的相同的图例/信息。
我很抱歉没有更可重复的示例，但我完全不了解这里发生的事情，所以我无法完全复制它。我想现在我只是在寻找在哪里寻找并尝试找出发生了什么变化。提前致谢。
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）