#meanrefs = find_files("/Users/phoetrymaster/Documents/School/Geography/Thesis/Data/MODIS_KANSAS_2007-2012/reprojected/Refs/2012", ".ref", recursive=False) #reffiles.append(meanrefs) rootout = "/Users/phoetrymaster/Documents/School/Geography/Thesis/Data/MODIS_KANSAS_2007-2012/reprojected/Classified/test1_envicurves/fullpxonly/geometricmean" #outfolders = [os.path.join(rootout, "clip1refs")]#, os.path.join(rootout, "clip2refs"), os.path.join(rootout, "clip3refs"), os.path.join(rootout, "clip4refs"), os.path.join(rootout, "clip5refs"), os.path.join(rootout, "clip6refs"), outfolders = [os.path.join(rootout, "clip1refs")] imagelist = [] clip1imgs = ["/Users/phoetrymaster/Documents/School/Geography/Thesis/Data/MODIS_KANSAS_2007-2012/reprojected/clips/KansasNDVI_2012_clip1.tif"]#, "/Users/phoetrymaster/Documents/School/Geography/Thesis/Data/MODIS_KANSAS_2007-2012/reprojected/clips/KansasEVI_2012_clip1.tif"] imagelist.append(clip1imgs) searchstrings = ["soy", "corn", "wwheat", "sorghum", "wwheatsoydbl"] fitmethods = ["SLSQP"]#, "TNC"] fullpixels = get_px_coords_from_points(clip1imgs[0], "/Users/phoetrymaster/Documents/School/Geography/Thesis/Data/MODIS_KANSAS_2007-2012/reprojected/clips/clip1_fullcells_points.shp") for method in fitmethods: for images in imagelist: for img in images: name = os.path.splitext(os.path.basename(img))[0] if "NDVI" in name: type = "NDVI" else: type = "EVI" for reffile, outfolder in zip(reffiles, outfolders): refs = {} for f in reffile: if type in f:
outroot = "/Users/phoetrymaster/Documents/School/Geography/Thesis/Data/MODIS_KANSAS_2007-2012/reprojected/Refs/2012clip1test2" samplepointsroot = "/Users/phoetrymaster/Documents/School/Geography/Thesis/Data/MODIS_KANSAS_2007-2012/SampleAreas/clip1/envipoints" clips = [] #clip1 #wwheat = [(39, 37), (45, 50), (11, 54), (4, 19), (75, 34)] #corn = [(48, 44), (10, 68), (8, 37), (67, 66), (83, 21)] #soy = [(42, 45), (31, 86), (7, 64), (55, 83), (73, 73)] #wwheatsoydbl = [(32, 62), (38, 55), (25, 56), (39, 30), (97, 3)] #sorghum = [(82, 52), (70, 38), (49, 33), (36, 28), (28, 27)] clip1imgs = ["/Users/phoetrymaster/Documents/School/Geography/Thesis/Data/MODIS_KANSAS_2007-2012/reprojected/clips/KansasNDVI_2012_clip1.tif", "/Users/phoetrymaster/Documents/School/Geography/Thesis/Data/MODIS_KANSAS_2007-2012/reprojected/clips/KansasEVI_2012_clip1.tif"] clip1out = os.path.join(outroot, "clip1") wwheat = get_px_coords_from_points(clip1imgs[0], os.path.join(samplepointsroot, "wwheat.shp")) corn = get_px_coords_from_points(clip1imgs[0], os.path.join(samplepointsroot, "corn.shp")) soy = get_px_coords_from_points(clip1imgs[0], os.path.join(samplepointsroot, "soy.shp")) wwheatsoydbl = get_px_coords_from_points(clip1imgs[0], os.path.join(samplepointsroot, "wwheatsoydbl.shp")) sorghum = get_px_coords_from_points(clip1imgs[0], os.path.join(samplepointsroot, "sorghum.shp")) clip1locs = {"wwheat": wwheat, "corn": corn, "soy": soy, "wwheatsoydbl": wwheatsoydbl, "sorghum": sorghum} clips.append((clip1locs, clip1imgs, clip1out)) ##clip2 #wwheat = [(2, 11), (6, 28), (63, 29), (23, 55), (6, 55), (3, 36)] #corn = [(3, 36), (79, 53), (50, 50), (66, 25), (70, 4), (65, 63)] #soy = [(6, 25), (6, 14), (97, 35), (85, 49), (60, 25)] #wwheatsoydbl = [(41, 53), (73, 61), (44, 32), (47, 64)] #sorghum = []