img_b5_11 = functions.create_band(band=b5_11) img_b7_11 = functions.create_band(band=b7_11) img_ndvi_11 = functions.create_band(band=ndvi_11) img_ndwi_11 = functions.create_band(band=ndwi_11) img_mndwi1_11 = functions.create_band(band=mndwi1_11) img_mndwi2_11 = functions.create_band(band=mndwi2_11) img_ndbi_11 = functions.create_band(band=ndbi_11) img_mndbi_11 = functions.create_band(band=mndbi_11) img_11 = functions.combine_bands_sf( b1=img_b1_11, b2=img_b2_11, b3=img_b3_11, b4=img_b4_11, b5=img_b5_11, b7=img_b7_11, ndvi=img_ndvi_11, ndwi=img_ndwi_11, mndwi1=img_mndwi1_11, mndwi2=img_mndwi2_11, ndbi=img_ndbi_11, mndbi=img_mndbi_11, multiband_array_file='L2 imagery/img_11_sf.npy') #----------------------------------------------------------------------------------------------- #train and predict with RF #400 trees result_11_300 = functions.train_rf( trees=300, maxfeatures=None, train_array=train_array_11, gt_array=train11,
#2009 b1_09_t = functions.extract_values(shp = file_test09_shp, raster = b1_09) b2_09_t = functions.extract_values(shp = file_test09_shp, raster = b2_09) b3_09_t = functions.extract_values(shp = file_test09_shp, raster = b3_09) b4_09_t = functions.extract_values(shp = file_test09_shp, raster = b4_09) b5_09_t = functions.extract_values(shp = file_test09_shp, raster = b5_09) b7_09_t = functions.extract_values(shp = file_test09_shp, raster = b7_09) ndvi_09_t = functions.extract_values(shp = file_test09_shp, raster = ndvi_09) ndwi_09_t = functions.extract_values(shp = file_test09_shp, raster = ndwi_09) mndwi1_09_t = functions.extract_values(shp = file_test09_shp, raster = mndwi1_09) mndwi2_09_t = functions.extract_values(shp = file_test09_shp, raster = mndwi2_09) ndbi_09_t = functions.extract_values(shp = file_test09_shp, raster = ndbi_09) mndbi_09_t = functions.extract_values(shp = file_test09_shp, raster = mndbi_09) test_09 = functions.combine_bands_sf(b1 = b1_09_t, b2 = b2_09_t, b3 = b3_09_t, b4 = b4_09_t, b5 = b5_09_t, b7 = b7_09_t, ndvi = ndvi_09_t, ndwi = ndwi_09_t, mndwi1 = mndwi1_09_t, mndwi2 = mndwi2_09_t, ndbi = ndbi_09_t, mndbi = mndbi_09_t, multiband_array_file = 'Samples/SLC off/C_test_array_sf_09.npy') #2011 b1_11_t = functions.extract_values(shp = file_test11_shp, raster = b1_11) b2_11_t = functions.extract_values(shp = file_test11_shp, raster = b2_11) b3_11_t = functions.extract_values(shp = file_test11_shp, raster = b3_11) b4_11_t = functions.extract_values(shp = file_test11_shp, raster = b4_11) b5_11_t = functions.extract_values(shp = file_test11_shp, raster = b5_11) b7_11_t = functions.extract_values(shp = file_test11_shp, raster = b7_11) ndvi_11_t = functions.extract_values(shp = file_test11_shp, raster = ndvi_11) ndwi_11_t = functions.extract_values(shp = file_test11_shp, raster = ndwi_11) mndwi1_11_t = functions.extract_values(shp = file_test11_shp, raster = mndwi1_11) mndwi2_11_t = functions.extract_values(shp = file_test11_shp, raster = mndwi2_11) ndbi_11_t = functions.extract_values(shp = file_test11_shp, raster = ndbi_11) mndbi_11_t = functions.extract_values(shp = file_test11_shp, raster = mndbi_11)
img_b5_09 = functions.create_band(band=b5_09) img_b7_09 = functions.create_band(band=b7_09) img_ndvi_09 = functions.create_band(band=ndvi_09) img_ndwi_09 = functions.create_band(band=ndwi_09) img_mndwi1_09 = functions.create_band(band=mndwi1_09) img_mndwi2_09 = functions.create_band(band=mndwi2_09) img_ndbi_09 = functions.create_band(band=ndbi_09) img_mndbi_09 = functions.create_band(band=mndbi_09) img_09 = functions.combine_bands_sf( b1=img_b1_09, b2=img_b2_09, b3=img_b3_09, b4=img_b4_09, b5=img_b5_09, b7=img_b7_09, ndvi=img_ndvi_09, ndwi=img_ndwi_09, mndwi1=img_mndwi1_09, mndwi2=img_mndwi2_09, ndbi=img_ndbi_09, mndbi=img_mndbi_09, multiband_array_file='L2 imagery/img_09_sf.npy') #----------------------------------------------------------------------------------------------- #train and predict with RF #400 trees result_09_400 = functions.train_rf( trees=400, maxfeatures=None, train_array=train_array_09, gt_array=train09,