writeRasterのtiffが開かない

Question

私はいくつかのラスタ計算をしていますが、結果を新しいtiffに書きたいと思います。しかし、RのwriteRasterからの出力.tifは、一般的に使用されているイメージ化ソフトウェアのいずれでも開かない。それを修正するために私のコードに追加すべきことについて何か提案してください。事前に おかげで、

setwd() 

library(sp) 
library(raster) 
library(rgdal) 

band530.2 <- raster("TTC001270001.tif", band=2) 
band530.3 <- raster("TTC001270001.tif", band=3) 
band530 <- 4 * band530.2 + (band530.2 - band530.3) 

lf<- function(x,a,b) { 
    v <- getValues(x) 
    v <- a * v + b 
    x <- setValues(x, v) 
    return(x) 
} 

calib.band530 <- lf(band530, 0.0009842088, -0.04702277) 

writeRaster(calib.band530, filename="cTTC001270001.tif", format="GTiff", overwrite=TRUE) 

Answer 1

あなたのコードでは、ArcGISでうまく開かれ、完全に使用可能な地理TIFFになります。あなたは、通常の画像ビューアで出力を表示したい場合は

、あなたは、通常のラスタに変換することができます。これは、地理空間情報を排除する

tiff("test.tif",dim(calib.band530)[2],dim(calib.band530)[1]) 
grid.raster(as.raster(calib.band530),interpolate=FALSE) 
dev.off() 

。

視認可能なジオティフを作成するには、人が見ることができるようにするにはデータを再スケーリングしてから、もう1つのパラメータを色空間に追加する必要があります。これは私の作品：

ただ、色空間パラメータは画像ビューアで開くようになりますが、おそらく、黒になります追加

calib.band_rezero <- calib.band530 - cellStats(calib.band530,min) 
calib.band_stretch <- calib.band_rezero/cellStats(calib.band_rezero,max) *254 
writeRaster(calib.band_stretch, filename="cTTC001270001_viewable.tif", format="GTiff",datatype='INT1U', overwrite=TRUE) 

：

writeRaster(calib.band530, filename="cTTC001270001.tif", format="GTiff",datatype='INT1U', overwrite=TRUE)