zScoreThresh,shadowSumThresh,\
contractPixels,dilatePixels,\
correctIllumination,correctScale,\
exportComposites,outputName,exportPathRoot,crs,transform,scale,resampleMethod,
preComputedCloudScoreOffset,preComputedTDOMIRMean,preComputedTDOMIRStdDev,
landsatCollectionVersion = landsatCollectionVersion)
# Separate into scenes and composites for subsequent analysis
processedScenes = lsAndTs['processedScenes']
processedComposites = lsAndTs['processedComposites']
# Map.addLayer(processedComposites.select(['NDVI','NBR']),{'addToLegend':'false'},'Time Series (NBR and NDVI)',False)
for year in range(startYear + timebuffer, endYear + 1 - timebuffer):
t = processedComposites.filter(ee.Filter.calendarRange(year, year,
'year')).first()
Map.addLayer(t.float(), getImagesLib.vizParamsFalse, str(year), False)
####################################################################################################
# Load the study region
Map.addLayer(studyArea, {'strokeColor': '0000FF'}, "Study Area", True)
Map.centerObject(studyArea)
####################################################################################################
####################################################################################################
# View map
Map.turnOnInspector()
Map.view()
####################################################################################################
####################################################################################################
# If exporting composites, track the exports
if exportComposites:
taskManagerLib.trackTasks2()
transform = [1000,0,-2361915.0,0,-1000,3177735.0]
# Specify scale if transform is null
scale = None
####################################################################################################
#End user parameters
####################################################################################################
####################################################################################################
####################################################################################################
#Start function calls
####################################################################################################
####################################################################################################
#Call on master wrapper function to get climate composites
climateComposites = getImagesLib.getClimateWrapper(collectionName,studyArea,startYear,endYear,startJulian,endJulian,timebuffer,weights,compositingReducer,exportComposites,exportPathRoot,crs,transform,scale,exportBands)
Map.addTimeLapse(climateComposites.select(exportBands),{},'Climate Composite Time Lapse')
####################################################################################################
#Load the study region
Map.addLayer(studyArea, {'strokeColor': '0000FF'}, "Study Area", True)
Map.centerObject(studyArea)
####################################################################################################
####################################################################################################
# View map
Map.turnOnInspector()
Map.view()
####################################################################################################
####################################################################################################
# If exporting composites, track the exports
if exportComposites:taskManagerLib.trackTasks2()
#Set properties for asset
# ccdc = ccdc.copyProperties(processedScenes)
# ccdc = ccdc.setMulti(ccdcParams)
# ccdc = ee.Image(ccdc)
Map.addLayer(ccdc,{},'CCDC Output',False);
#Export output
if exportCCDC:
print('Exporting CCDC output')
exportName = outputName + '_'+str(startYear) + '_' + str(endYear) + '_' + str(startJulian) + '_' + str(endJulian)
exportPath = exportPathRoot + '/'+ exportName
#Export output
getImagesLib.exportToAssetWrapper(ccdc,exportName,exportPath,{'.default':'sample'},studyArea,scale,crs,transform)
####################################################################################################
#Load the study region
Map.addLayer(studyArea, {'strokeColor': '0000FF'}, "Study Area", False)
Map.centerObject(studyArea)
####################################################################################################
# View map
Map.turnOnInspector()
Map.view()
####################################################################################################
####################################################################################################
# Track the export
if exportCCDC:taskManagerLib.trackTasks2()
# }
#Export image
if exportCoefficients:
if not detrend:
coeffsOut = coeffs.multiply(1000).int16()
else:coeffsOut = coeffs.float()
coeffsOut = coeffsOut.copyProperties(coeffs).copyProperties(coeffs,['system:time_start'])
outName = outputName + str(startYearT) + '_'+ str(endYearT)
outPath = exportPathRoot + '/' + outName
getImagesLib.exportToAssetWrapper(coeffs,outName,outPath,'mean',studyArea,scale,crs,transform)
coeffCollection.append(coeffs)
####################################################################################################
#Load the study region
Map.addLayer(studyArea, {'strokeColor': '0000FF'}, "Study Area", False)
Map.centerObject(studyArea)
####################################################################################################
####################################################################################################
# View map
Map.turnOnInspector()
Map.view()
####################################################################################################
####################################################################################################
# If exporting composites, track the exports
if exportCoefficients:taskManagerLib.trackTasks2()
#Set up proper resampling for each band
#Be sure to change if the band names for the exported image change
pyrObj = {'_yr_':'mode','_dur_':'mode','_mag_':'mean','_slope_':'mean'}
possible = ['loss','gain']
how_many_list = ee.List.sequence(1,howManyToPull).getInfo()
outObj = {}
for p in possible:
for key in pyrObj.keys():
for i in how_many_list:
i = int(i)
kt = indexName + '_LT_'+p + key+str(i)
outObj[kt]= pyrObj[key]
#Export output
getImagesLib.exportToAssetWrapper(ltOutStack,exportName,exportPath,outObj,studyArea,scale,crs,transform)
####################################################################################################
#Load the study region
Map.addLayer(studyArea, {'strokeColor': '0000FF'}, "Study Area", False)
Map.centerObject(studyArea)
####################################################################################################
####################################################################################################
# View map
Map.turnOnInspector()
Map.view()
####################################################################################################
####################################################################################################
# If exporting LT stack, track the exports
if exportLTStack:taskManagerLib.trackTasks2()