for rast in os.listdir(r'/mnt/local_drive/britton/PRISM_data/PRISM_daily_tmax'):
if rast[-4: ] == '.tif':
tmax = rasterio.open(r'/mnt/local_drive/britton/PRISM_data/PRISM_daily_tmax' + '//' + rast)
tmax_array = tmax.read(1)
affine = tmax.transform
tract_average = zonal_stats(tract_polygon[(n-1):n],
tmax_array,
affine = affine,
stats = ['mean'],
all_touched = True,
geojason_out = False)
tract_average = tract_average[0]['mean']
if n == 1:
data_map.loc[i]['Date'] = rast[11:-4]
data_map.iloc[(i-1):i, n] = tract_average * 1.8 + 32 # optional conversion to F
i = i + 1
for rast in os.listdir(r'/mnt/local_drive/britton/PRISM_data/PRISM_daily_tmax'):
if rast[-4: ] == '.tif':
tmax = rasterio.open(r'/mnt/local_drive/britton/PRISM_data/PRISM_daily_tmax' + '//' + rast)
tmax_array = tmax.read(1)
affine = tmax.transform
tract_average = zonal_stats(tract_polygon[(value-1):value],
tmax_array,
affine = affine,
stats = ['mean'],
all_touched = True,
geojason_out = False)
tract_average = tract_average[0]['mean']
if value == 1:
data_map.loc[i]['Date'] = rast[11:-4]
data_map.iloc[(i-1):i, value] = tract_average * 1.8 + 32 # optional conversion to F
i = i + 1
pool_obj = multiprocessing.Pool()
process = pool_obj.map(spatialaverage, range(1, len(tract_id_list)))
data_map
