Subset requires a fully set raster to work properly. Enable also the possibility to create a subset based on a rotated bounding box. Have a look at the "scidb4gdal" repository on a merge command, where non existing values are replaced by zero (or NA value)