首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >R光栅:裁剪后的光栅颜色(亮度)是否与原始光栅不同?

R光栅:裁剪后的光栅颜色(亮度)是否与原始光栅不同?
EN

Stack Overflow用户
提问于 2018-11-20 16:21:36
回答 1查看 40关注 0票数 0

我想要crop一个多波段光栅(4波段)由空间多边形(在SpatialPolygonsDataFrame)。当我在QGIS中显示原始栅格和裁剪后的栅格时,我发现裁剪后的栅格与原始栅格的颜色不同。下面是我的代码:

代码语言:javascript
复制
library(raster) 

mosaic_shp <- shapefile("mo_clipper.shp")
mosaic <- brick('orthomosaic.tif')

mosaic_sub <- crop(mosaic, extent(mosaic_shp))

writeRaster(mosaic_sub, 'mosaic_sub.tif', format = "GTiff", overwrite = TRUE)

QGIS中部分裁剪的栅格和原始栅格中的对应部分:

我不知道如何处理这个问题,任何帮助都将不胜感激。

EN

回答 1

Stack Overflow用户

发布于 2018-11-20 17:01:05

在QGIS中仔细比较了这两个栅格之后,我找到了答案。该问题与writeRaster函数中的dataType参数有关。所以我们只需要修改代码,如下所示:

代码语言:javascript
复制
library(raster) 

mosaic_shp <- shapefile("mo_clipper.shp")
mosaic <- brick('orthomosaic.tif')

mosaic_sub <- crop(mosaic, extent(mosaic_shp))
data_type <- unique(dataType(mosaic)) # get data type from original raster;

writeRaster(mosaic_sub, 'mosaic_sub.tif', format = "GTiff", overwrite = TRUE, 
            datatype = data_type) # set datatype;
票数 0
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/53388831

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档