def execute(self, parameters, messages):
""" Run tool sourcecode by calling metric.py """
metric_workspace = str(parameters[15].value)
metrci_workspace = metric_workspace.replace("\\", "/")
landsat_files = [
str(parameters[1].value),
str(parameters[2].value),
str(parameters[3].value),
str(parameters[4].value),
str(parameters[5].value),
str(parameters[6].value),
str(parameters[7].value),
str(parameters[8].value)
]
landsat_meta = str(parameters[0].value)
dem_path = str(parameters[9].value)
hot_shp_path = str(parameters[10].value)
cold_shp_path = str(parameters[11].value)
weather_path = str(parameters[12].value)
crop = str(parameters[13].value)
timezone = float(parameters[14].value)
try:
wx_elev = float(parameters[16].value)
except:
wx_elev = None
try:
wx_zom = float(parameters[17].value)
except:
wx_zom = None
try:
LE_ref = float(parameters[18].value)
except:
LE_ref = None
try:
mounts = bool(parameters[19].value)
except:
mounts = None
testflag = str(parameters[20].value)
recalc = bool(parameters[21].value)
# print parameters to screen
for param in parameters:
message = "{0} : {1}".format(
str(param.name).ljust(12, " "), param.value)
print(message)
arcpy.AddMessage(message)
metric_py.run(metric_workspace, landsat_files, landsat_meta, dem_path,
hot_shp_path, cold_shp_path, weather_path, testflag,
recalc, crop, timezone, wx_elev, wx_zom, LE_ref, mounts)
def execute(self, parameters, messages):
""" Run tool sourcecode by calling metric.py """
metric_workspace = str(parameters[15].value)
metrci_workspace = metric_workspace.replace("\\","/")
landsat_files = [str(parameters[1].value),
str(parameters[2].value),
str(parameters[3].value),
str(parameters[4].value),
str(parameters[5].value),
str(parameters[6].value),
str(parameters[7].value),
str(parameters[8].value)]
landsat_meta = str(parameters[0].value)
dem_path = str(parameters[9].value)
hot_shp_path = str(parameters[10].value)
cold_shp_path = str(parameters[11].value)
weather_path = str(parameters[12].value)
crop = str(parameters[13].value)
timezone = float(parameters[14].value)
try: wx_elev = float(parameters[16].value)
except: wx_elev = None
try: wx_zom = float(parameters[17].value)
except: wx_zom = None
try: LE_ref = float(parameters[18].value)
except: LE_ref = None
try: mounts = bool(parameters[19].value)
except: mounts = None
testflag = str(parameters[20].value)
recalc = bool(parameters[21].value)
# print parameters to screen
for param in parameters:
message = "{0} : {1}".format(str(param.name).ljust(12," "), param.value)
print(message)
arcpy.AddMessage(message)
metric_py.run(metric_workspace, landsat_files, landsat_meta, dem_path,
hot_shp_path, cold_shp_path, weather_path, testflag,
recalc, crop, timezone, wx_elev, wx_zom, LE_ref, mounts)
landsat_files =[r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B2.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B3.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B4.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B5.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B6.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B7.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B10.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B11.TIF.tif"]
landsat_meta = r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_MTL.txt"
dem_path = r"C:\Users\Jeff\Desktop\metric_testing\SRTM_mosaic.tif"
hot_shp_path = r"C:\Users\jwely\Desktop\metric_trub\test_inputs_NC\input_ref_pixels\nonirrigated_hot.shp"
cold_shp_path = r"C:\Users\jwely\Desktop\metric_trub\test_inputs_NC\input_ref_pixels\irrigated_cool.shp"''
weather_path = r"C:\Users\jwely\Desktop\metric_trub\test_inputs_NC\input_weather\6564006671082dat.txt"
testflag = "ALL" # valid inputs are "ALL", "LIMITED", "ET-ONLY" (ET-ONLY 25% less time than ALL)
recalc = True # recalc = True will prevent all
crop = "alfalfa" # reference crop. possible values are "grass" , "alfalfa"
timezone = -5.0 # -5 or -4 for Eeastern time depending on daylight savings
wx_elev = 1.000 # elevation of the weather station
wx_zom = 0.010 # estimated zom at weather station (temporary)
LE_cold_cal_fac = 1.050 # LE calibration factor, should probably always be 1.05
mountains = True # set true for mountainous regions
metric_py.run(workspace, landsat_files, landsat_meta, dem_path,
hot_shp_path, cold_shp_path, weather_path, testflag,
recalc, crop, timezone, wx_elev, wx_zom, LE_cold_cal_fac, mountains)
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B2.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B3.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B4.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B5.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B6.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B7.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B10.TIF",
r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_B11.TIF.tif"
]
landsat_meta = r"C:\Users\Jeff\Desktop\metric_testing\LC80440282015144LGN00\LC80440282015144LGN00_MTL.txt"
dem_path = r"C:\Users\Jeff\Desktop\metric_testing\SRTM_mosaic.tif"
hot_shp_path = r"C:\Users\jwely\Desktop\metric_trub\test_inputs_NC\input_ref_pixels\nonirrigated_hot.shp"
cold_shp_path = r"C:\Users\jwely\Desktop\metric_trub\test_inputs_NC\input_ref_pixels\irrigated_cool.shp" ''
weather_path = r"C:\Users\jwely\Desktop\metric_trub\test_inputs_NC\input_weather\6564006671082dat.txt"
testflag = "ALL" # valid inputs are "ALL", "LIMITED", "ET-ONLY" (ET-ONLY 25% less time than ALL)
recalc = True # recalc = True will prevent all
crop = "alfalfa" # reference crop. possible values are "grass" , "alfalfa"
timezone = -5.0 # -5 or -4 for Eeastern time depending on daylight savings
wx_elev = 1.000 # elevation of the weather station
wx_zom = 0.010 # estimated zom at weather station (temporary)
LE_cold_cal_fac = 1.050 # LE calibration factor, should probably always be 1.05
mountains = True # set true for mountainous regions
metric_py.run(workspace, landsat_files, landsat_meta, dem_path,
hot_shp_path, cold_shp_path, weather_path, testflag, recalc,
crop, timezone, wx_elev, wx_zom, LE_cold_cal_fac, mountains)
