def createAOIName(AOI_Name,
AOI_Num,
srcImageryListOrig,
srcVectorAOIFile,
srcVectorFileList,
outputDirectory,
clipImageryToAOI=True,
windowSizeMeters=200,
clipOverlap=0.0,
minpartialPerc=0.0,
vrtMosaic=True,
createPix=False,
createSummaryCSVChallenge=True,
csvLabel='All',
featureName='Buildings'):
srcImageryList = []
if clipImageryToAOI:
for srcImagery in srcImageryListOrig:
print(srcImagery)
AOI_HighResMosaicName = os.path.join(
outputDirectory,
'AOI_{}_{}_{}.vrt'.format(AOI_Num, AOI_Name, srcImagery[1]))
if vrtMosaic:
AOI_HighResMosaicClipName = AOI_HighResMosaicName.replace(
'.vrt', 'clipped.vrt')
else:
AOI_HighResMosaicClipName = AOI_HighResMosaicName.replace(
'.vrt', 'clipped.TIF')
if os.path.isfile(AOI_HighResMosaicClipName):
os.remove(AOI_HighResMosaicClipName)
if vrtMosaic:
command = 'gdalwarp -of VRT ' + \
'-cutline ' + srcVectorAOIFile + ' ' + \
srcImagery[0] + ' ' + \
AOI_HighResMosaicClipName
else:
command = 'gdalwarp ' + \
'-cutline ' + srcVectorAOIFile + ' ' + \
srcImagery[0] + ' ' + \
AOI_HighResMosaicClipName
print(command)
os.system(command)
srcImageryList.append([AOI_HighResMosaicClipName, srcImagery[1]])
else:
srcImageryList = srcImageryListOrig
# rasterFileList = [['rasterLocation', 'rasterDescription']]
# i.e rasterFileList = [['/path/to/3band_AOI_1.tif, '3band'],
# ['/path/to/8band_AOI_1.tif, '8band']
# ]
chipSummaryList = gT.cutChipFromMosaic(srcImageryList,
srcVectorFileList,
outlineSrc=srcVectorAOIFile,
outputDirectory=outputDirectory,
outputPrefix='',
clipSizeMX=windowSizeMeters,
clipSizeMY=windowSizeMeters,
clipOverlap=clipOverlap,
minpartialPerc=minpartialPerc,
createPix=createPix,
baseName='AOI_{}_{}'.format(
AOI_Num, AOI_Name),
imgIdStart=1)
outputCSVSummaryName = 'AOI_{}_{}_{}_{}_solutions.csv'.format(
AOI_Num, AOI_Name, csvLabel, featureName)
createCSVSummaryFile(chipSummaryList,
outputCSVSummaryName,
rasterChipDirectory='',
replaceImageID='',
createProposalsFile=False,
pixPrecision=2)
def createAOIName(AOI_Name, AOI_Num,
srcImageryListOrig,
srcVectorAOIFile,
srcVectorFileList,
outputDirectory,
clipImageryToAOI=True,
windowSizeMeters=200,
clipOverlap=0.0,
minpartialPerc=0.0,
vrtMosaic=True,
createPix=False,
createSummaryCSVChallenge=True,
csvLabel='All',
featureName='Buildings'):
srcImageryList = []
if clipImageryToAOI:
for srcImagery in srcImageryListOrig:
print(srcImagery)
AOI_HighResMosaicName = os.path.join(outputDirectory, 'AOI_{}_{}_{}.vrt'.format(AOI_Num, AOI_Name, srcImagery[1]))
if vrtMosaic:
AOI_HighResMosaicClipName = AOI_HighResMosaicName.replace('.vrt', 'clipped.vrt')
else:
AOI_HighResMosaicClipName = AOI_HighResMosaicName.replace('.vrt', 'clipped.TIF')
if os.path.isfile(AOI_HighResMosaicClipName):
os.remove(AOI_HighResMosaicClipName)
if vrtMosaic:
command = 'gdalwarp -of VRT ' + \
'-cutline ' + srcVectorAOIFile + ' ' + \
srcImagery[0] + ' ' + \
AOI_HighResMosaicClipName
else:
command = 'gdalwarp ' + \
'-cutline ' + srcVectorAOIFile + ' ' + \
srcImagery[0] + ' ' + \
AOI_HighResMosaicClipName
print(command)
os.system(command)
srcImageryList.append([AOI_HighResMosaicClipName, srcImagery[1]])
else:
srcImageryList = srcImageryListOrig
# rasterFileList = [['rasterLocation', 'rasterDescription']]
# i.e rasterFileList = [['/path/to/3band_AOI_1.tif, '3band'],
# ['/path/to/8band_AOI_1.tif, '8band']
# ]
chipSummaryList = gT.cutChipFromMosaic(srcImageryList, srcVectorFileList, outlineSrc=srcVectorAOIFile,
outputDirectory=outputDirectory, outputPrefix='',
clipSizeMX=windowSizeMeters, clipSizeMY=windowSizeMeters, clipOverlap=clipOverlap,
minpartialPerc=minpartialPerc, createPix=createPix,
baseName='AOI_{}_{}'.format(AOI_Num, AOI_Name),
imgIdStart=1)
outputCSVSummaryName = 'AOI_{}_{}_{}_{}_solutions.csv'.format(AOI_Num, AOI_Name, csvLabel,featureName)
createCSVSummaryFile(chipSummaryList, outputCSVSummaryName, rasterChipDirectory='', replaceImageID='',
createProposalsFile=False,
pixPrecision=2)