joshobrien / gdalutilities Goto Github PK

m <- matrix(rnorm(n=10^4,m=0,sd=2),nrow=100)
m[1:10,] <- NA
imatest <- raster(m)
cellStats(imatest,range)
imatest <- brick(imatest,imatest*2,imatest+10)
extent(imatest) <- c(0,100,0,100)
projection(imatest) <- CRS("+init=epsg:25831")
imatest
class       : RasterBrick 
dimensions  : 100, 100, 10000, 3  (nrow, ncol, ncell, nlayers)
resolution  : 1, 1  (x, y)
extent      : 0, 100, 0, 100  (xmin, xmax, ymin, ymax)
coord. ref. : +init=epsg:25831 +proj=utm +zone=31 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs 
data source : in memory
names       :    layer.1,    layer.2,    layer.3 
min values  :  -7.844195, -15.688390,   2.155805 
max values  :   7.681265,  15.362530,  17.681265 

writeRaster(imatest,filename="imatest",format="GTiff",NAflag=-9999,overwrite=TRUE)
qt <- quantile(imatest,probs=c(0.02,0.98),na.rm=TRUE)
qt
scal <- 1
for(i in 1:nlayers(imatest)){
  scal <- c(scal,c(qt[i,],c(1,255)))
}
scal <- scal[-1]

> scal
        2%        98%                               2%        98%                       
 -4.133383   4.107563   1.000000 255.000000  -8.266766   8.215126   1.000000 255.000000 
        2%        98%                       
  5.866617  14.107563   1.000000 255.000000 
> gdal_translate(src_dataset="imatest.tif",
+                dst_dataset="imateststretchgdal2.tif",
+                ot="Byte",scale=scal,dryrun=TRUE)

> gdal_translate(src_dataset="imatest.tif",
+                dst_dataset="imateststretchgdal2.tif",
+                ot="Byte",scale=scal)
Warning message:
In CPL_gdaltranslate(source, destination, options) :
  GDAL Error 6: Unknown option name '-8.26676644515004'

> mscal <- matrix(scal,ncol=3, byrow=FALSE)
> mscal
           [,1]       [,2]       [,3]
[1,]  -4.133383  -8.266766   5.866617
[2,]   4.107563   8.215126  14.107563
[3,]   1.000000   1.000000   1.000000
[4,] 255.000000 255.000000 255.000000
> gdal_translate(src_dataset="imatest.tif",
+                dst_dataset="imateststretchgdal2.tif",
+                ot="Byte",scale=mscal,dryrun=TRUE)
[1] "gdal_translate imatest.tif imateststretchgdal2.tif -ot Byte -scale -4.13338322257502 4.10756297485811 1 255 -8.26676644515004 8.21512594971622 1 255 5.86661677742498 14.1075629748581 1 255

--->  Building R-gdalUtilities
--->  Testing R-gdalUtilities
* using log directory ‘/opt/local/var/macports/build/_opt_PPCSnowLeopardPorts_R_R-gdalUtilities/R-gdalUtilities/work/gdalUtilities/gdalUtilities.Rcheck’
* using R version 4.3.1 (2023-06-16)
* using platform: powerpc-apple-darwin10.0.0d2 (32-bit)
* R was compiled by
    gcc-mp-12 (MacPorts gcc12 12.3.0_0) 12.3.0
    GNU Fortran (MacPorts gcc12 12.3.0_0) 12.3.0
* running under: OS X Snow Leopard 10.6
* using session charset: UTF-8
* checking for file ‘gdalUtilities/DESCRIPTION’ ... OK
* checking extension type ... Package
* this is package ‘gdalUtilities’ version ‘1.2.5’
* package encoding: UTF-8
* checking package namespace information ... OK
* checking package dependencies ... OK
* checking if this is a source package ... OK
* checking if there is a namespace ... OK
* checking for executable files ... OK
* checking for hidden files and directories ... OK
* checking for portable file names ... OK
* checking for sufficient/correct file permissions ... OK
* checking whether package ‘gdalUtilities’ can be installed ... OK
* checking installed package size ... OK
* checking package directory ... OK
* checking DESCRIPTION meta-information ... OK
* checking top-level files ... OK
* checking for left-over files ... OK
* checking index information ... OK
* checking package subdirectories ... OK
* checking R files for non-ASCII characters ... OK
* checking R files for syntax errors ... OK
* checking whether the package can be loaded ... OK
* checking whether the package can be loaded with stated dependencies ... OK
* checking whether the package can be unloaded cleanly ... OK
* checking whether the namespace can be loaded with stated dependencies ... OK
* checking whether the namespace can be unloaded cleanly ... OK
* checking loading without being on the library search path ... OK
* checking dependencies in R code ... OK
* checking S3 generic/method consistency ... OK
* checking replacement functions ... OK
* checking foreign function calls ... OK
* checking R code for possible problems ... OK
* checking Rd files ... OK
* checking Rd metadata ... OK
* checking Rd cross-references ... OK
* checking for missing documentation entries ... OK
* checking for code/documentation mismatches ... OK
* checking Rd \usage sections ... OK
* checking Rd contents ... OK
* checking for unstated dependencies in examples ... OK
* checking examples ... OK
* checking PDF version of manual ... OK
* DONE

Status: OK

 *** caught segfault ***
address 0xc, cause 'invalid permissions'

Traceback:
 1: CPL_gdalwarp(source, destination, options, oo, doo, config_options,     quiet, "-overwrite" %in% options)
 2: gdal_utils("warp", source = srcfile, destination = dstfile, options = opts,     config_options = config_options)
 3: gdalwarp(srcfile = f[X], dstfile = outPath, t_srs = "EPSG:32617",     te = c(extCommon[1], extCommon[3], extCommon[2], extCommon[4]),     tr = c(resMin, resMin), tap = TRUE, r = "cubic", overwrite = TRUE,     dryrun = FALSE)
 4: FUN(X[[i]], ...)
 5: lapply(X = X, FUN = FUN, ...)
 6: sapply(1:length(orth$f), standardizeOrthos, f = orth$f, orthos = orth$orthos,     pathOriginal, pathStandard, extCommon, resMin)

Possible actions:
1: abort (with core dump, if enabled)
2: normal R exit
3: exit R without saving workspace
4: exit R saving workspace
Selection: 1
R is aborting now ...
zsh: segmentation fault  R