Exemple #1
0
     imageutils.subset(inND4, contextvect, outND4, 'GTiff',
                       rsgislib.TYPE_8INT)
 except:
     input("imageutils.subset not working")
     gdwarpcmd = "gdalwarp -cutline {c} -crop_to_cutline {i} {out}".format(
         c=contextvect,
         i=inND4.replace("TOSHIBA EXT", "TOSHIBA\ EXT"),
         out=outND4)
     print(gdwarpcmd)
     subprocess.call(gdwarpcmd, shell=True)
 if not (MGSmode):
     try:
         # if possible, work in rsgislib but fall back on gdalwarp if needed
         if not (os.path.isfile(DTMband)):
             imageutils.selectImageBands(
                 inDTM, DTMband, 'GTiff', rsgislib.TYPE_32FLOAT,
                 [2])
         imageutils.subset(DTMband, contextvect, outDTM,
                           'GTiff', rsgislib.TYPE_32FLOAT)
     except:
         input(
             "imageutils.selectImageBands or imageutils.subset not working"
         )
         gdwarpcmd = "gdalwarp -cutline {c} -crop_to_cutline {i} {out}".format(
             c=contextvect,
             i=DTMband.replace("TOSHIBA EXT", "TOSHIBA\ EXT"),
             out=outDTM)
         print(gdwarpcmd)
         subprocess.call(gdwarpcmd, shell=True)
 else:
     try:
Exemple #2
0
def selectImageBands(inImage, outImage, bandList):
    gdalFormat = 'KEA'
    dataType = rsgislib.TYPE_32FLOAT
    
    imageutils.selectImageBands(inImage, outImage, gdalFormat, dataType, bandList)