diff --git a/clc_reclass_to_fragmented_adding_abandoment.py b/clc_reclass_to_fragmented_adding_abandoment.py deleted file mode 100644 index ccdb633300d8889714ace8e350ff8c0a487bef30..0000000000000000000000000000000000000000 --- a/clc_reclass_to_fragmented_adding_abandoment.py +++ /dev/null @@ -1,84 +0,0 @@ -# -*- coding: utf-8 -*- -""" -Created on Wed Feb 06 16:58:17 2019 - -@author: fw56moba -""" - -import numpy as np - -import MacPyver as mp - - -clc_path = r"I:\biocon\GIS_data\Corine_LandCover\g250_clc12_V18_5.tif" -aban_path = r"I:\biocon\Projects\RewildingMAP\DataSources\LandAbandonement\LandAbandonment_Scenarios\A1AbandAr2040_250m_3035.tif" - -#read in rasters -clc = mp.raster.tiff.read_tif(clc_path) -aba = mp.raster.tiff.read_tif(aban_path) - - -# values are the unique values form the column code_12 from the -#corine landcover shapefile -#using: shp_values = set(row[0] for row in arcpy.da.SearchCursor("clc12_Version_18_5", "code_12")) -shp_values = [u'111', u'112', u'121', u'122', u'123', u'124', u'131', u'132', u'133', u'141', - u'142', u'211', u'212', u'213', u'221', u'222', u'223', u'231', u'241', u'242', - u'243', u'244', u'311', u'312', u'313', u'321', u'322', u'323', u'324', u'331', - u'332', u'333', u'334', u'335', u'411', u'412', u'421', u'422', u'423', u'511', - u'512', u'521', u'522', u'523', u'999'] - -#get unique values from clc raster -uniq = np.unique(clc) - -#create translation from raster value to code_12 value -clc_dic = {} -for x in range(len(shp_values)): - clc_dic[uniq[x]] = int(shp_values[x]) - -"""result from the dictionary function above: -#dictionary contains the translation from raster value to shape code_12 value -''' -{1: 111, 2: 112, 3: 121, 4: 122, 5: 123, 6: 124, 7: 131, 8: 132, 9: 133, - 10: 141, 11: 142, 12: 211, 13: 212, 14: 213, 15: 221, 16: 222, 17: 223, - 18: 231, 19: 241, 20: 242, 21: 243, 22: 244, 23: 311, 24: 312, 25: 313, - 26: 321, 27: 322, 28: 323, 29: 324, 30: 331, 31: 332, 32: 333, 33: 334, - 34: 335, 35: 411, 36: 412, 37: 421, 38: 422, 39: 423, 40: 511, 41: 512, - 42: 521, 43: 522, 44: 523, 48: 999} -""" - -#copy raster -clc_code_12 = clc[:] -#translate raster values based on the dictionary -for x in clc_dic: - clc_code_12 = np.where(clc_code_12 == x, clc_dic[x], clc_code_12) - - -#reclass raster pastures as intensive agriculture -dicti = {'a':[111,113,999], 'b':[211,242,999], 'c':[512, 523, 999], 'd':[995,995,999]} - -inraster = clc_code_12[:] -outraster = inraster[:] -for x in dicti: - start,end,new_value = dicti[x] - outraster = np.where((inraster>=start)&(inraster<=end),new_value, outraster ) - -#translate all values to 1 which are not 999 or 255 -outraster = np.where((outraster != 999) & (outraster != 255), 1, outraster) - -#add abandomant area to landcover -outraster = np.where(aba == 1, 1, outraster) - -#save to tif file -mp.raster.tiff.write_tif(clc_path, r'O:\GIS_data\Corine_LandCover\g250_clc12_V18_5_reclass_binar_02.tif', outraster, 0, nodata = 255, option = "COMPRESS=DEFLATE") - - -outraster = np.where(outraster==999, 9, outraster) -outraster = np.where(clc == 44, 0, outraster) -outraster = np.where(outraster==255,0,outraster) -mp.raster.tiff.write_tif(clc_path, r'I:\biocon\Emmanuel_Oceguera\projects\2023_03_NaturaConnect\dataSources\clc2018_v2020_geoPackage_Raster\g250_clc12_V18_5_reclass_binar_excl_abandom_8bit.tif',outraster, 6, nodata=0, option="COMPRESS=DEFLATE") - - -""" -Continue in arcmap -convert raster to polygon -""" \ No newline at end of file