02ClippingSubscene
Created viernes 19 febrero 2016
R functions
S2Subsc_log.R
- /dataext3/dataalobo/Sentinel2/RSentinel2/S2Subsc_log.R
- Based in gdalutilities
- Can use davfs for ictjabox (but that can be slooow)
- Steps
- 1. Read and crop aoi (from shape) using gdal_translate
- Results in 3 images for 10, 20 and 60 m bands
- 2. Merge and Resample to 10m resolution using gdal_merge.py
- 3. Write multi-spectral image with correct band order
- reord <- c(11,1,2,3,5,6,7,4,8,12,13,9,10)
- 1. Read and crop aoi (from shape) using gdal_translate
- make spreasheet with bands
- CIR can be done with 8,4,2
S2Subsc_withR_log.R
Version trying to use R as much as possible (but R cannot read the format, thus the initial step must be the same as above
/dataext3/dataalobo/Sentinel2/RSentinel2/S2Subsc_withR_log.R
Improvements
- use pansharpening for interpolating to 10m
- push it to github
Example
/dataext3/dataalobo/Sentinel2/S2aoi1_10.tif