Skip to content
Snippets Groups Projects
Commit 11693b97 authored by Neha Hunka's avatar Neha Hunka
Browse files

Classifying the world into IPCC Table 4.7

parent a6607035
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
country_summaries/IPCC_classes_DPS/.ipynb_checkpoints/IPCC_GEDI_Table4.7-checkpoint.py
+ 15
− 11
View file @ 11693b97
......@@ -20,7 +20,7 @@ import argparse
# wget https://storage.googleapis.com/earthenginepartners-hansen/GFC-2022-v1.10/Hansen_GFC-2022-v1.10_lossyear_10S_060W.tif -O /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/Hansen_GFC-2022-v1.10_lossyear_10S_060W.tif
##### RUN SCRIPT FROM COMMAND LINE #####
# python /projects/ADE_biomass_harmonization/biomass_harmonization/country_summaries/IPCC_classes_DPS/IPCC_GEDI_Table4.7.py --FOREST_HEIGHT_2000 /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/2000_10S_060W.tif --FOREST_HEIGHT_2020 /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/2020_10S_060W.tif --FOREST_LOSS /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/netloss_10S_060W.tif --ESA_FOREST_COVER /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/treecover2010_10S_060W.tif --FOREST_OLD_LOSS /projects/my-public-bucket/Data/Harris_et_al_PAPER/Old_loss/10S_060W_OLD_LOSS.tif --FOREST_COVER_LOSSYEAR /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/Hansen_GFC-2022-v1.10_lossyear_10S_060W.tif --BOREAL_AGE /projects/my-public-bucket/Data/Harris_et_al_PAPER/Forest_age_boreal/BOREAL_AGE/10S_060W_BOREAL_AGE.tif --PLANTATIONS /projects/my-public-bucket/Data/Harris_et_al_PAPER/WRI_Planted_Forest_Type/10S_060W_plantation_type_oilpalm_woodfiber_other_unmasked.tif --IFL_rasters /projects/my-public-bucket/Data/Harris_et_al_PAPER/Intact_Forest_Landscapes/IFL_rasters/10S_060W_IFL_2000.tif --FII_raster /projects/my-public-bucket/Data/Harris_et_al_PAPER/Forest_Integrity_Index/10S_060W_FII.tif --primary_forest_asia /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_Asia.tif --primary_forest_sa /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_SouthAmerica.tif --primary_forest_africa /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_Africa.tif --EcoCont /projects/my-public-bucket/Data/Harris_et_al_PAPER/EcoCont_v2/10S_060W_EcoCont.tif --tile 10S_060W --output_file /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/IPCC_10S_060W.tif
# python /projects/ADE_biomass_harmonization/biomass_harmonization/country_summaries/IPCC_classes_DPS/IPCC_GEDI_Table4.7.py --FOREST_HEIGHT_2000 /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/2000_10S_060W.tif --FOREST_HEIGHT_2020 /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/2020_10S_060W.tif --FOREST_LOSS /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/netloss_10S_060W.tif --ESA_FOREST_COVER /projects/my-public-bucket/Data/Harris_et_al_PAPER/ESA_WorldCover2021_v2/WC_10S_060W_C.tif --FOREST_COVER_LOSSYEAR /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/Hansen_GFC-2022-v1.10_lossyear_10S_060W.tif --BOREAL_AGE /projects/my-public-bucket/Data/Harris_et_al_PAPER/Forest_age_boreal/BOREAL_AGE/10S_060W_BOREAL_AGE.tif --PLANTATIONS /projects/my-public-bucket/Data/Harris_et_al_PAPER/WRI_Planted_Forest_Type_v2/10S_060W_plantation_type_oilpalm_woodfiber_other.tif --IFL_rasters /projects/my-public-bucket/Data/Harris_et_al_PAPER/Intact_Forest_Landscapes/IFL_rasters_2020/10S_060W_IFL_2020.tif --FII_raster /projects/my-public-bucket/Data/Harris_et_al_PAPER/Forest_Integrity_Index/10S_060W_FII.tif --primary_forest_asia /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_Asia.tif --primary_forest_sa /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_SouthAmerica.tif --primary_forest_africa /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_Africa.tif --EcoCont /projects/my-public-bucket/Data/Harris_et_al_PAPER/EcoCont_v2/10S_060W_EcoCont.tif --tile 10S_060W --output_mask /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/Ages_10S_060W.tif --output_file /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/Classes_10S_060W.tif
######################################################################################
parse = argparse.ArgumentParser(description="Classifies world into GEZ, Continent and Forest Age classes")
......@@ -60,7 +60,7 @@ PRIMARY = FII_raster
del(FII_raster)
PRIMARY[PRIMARY <= 9600] = 0
PRIMARY[PRIMARY > 9600] = 1
PRIMARY[INTACT_FORESTS > 0] = 1
PRIMARY[INTACT_FORESTS == 1] = 1
del(INTACT_FORESTS)
PRIMARY[PRIMARY_FOREST_ASIA > 0] = 1
del(PRIMARY_FOREST_ASIA)
......@@ -93,6 +93,8 @@ ESA_FOREST_COVER = rasterio.open(args.ESA_FOREST_COVER).read(1)
PRIMARY[ESA_FOREST_COVER == 0] = 0
del(ESA_FOREST_COVER)
PRIMARY[(PRIMARY > 1) | (PRIMARY < 1)] = 0
# ###################################################################################
# ###################################################################################
# ##### YOUNG FOREST CONDITIONS - GIVEN VALUE 2 #####################################
......@@ -116,12 +118,9 @@ del(FOREST_COVER_LOSSYEAR,BACKUP_FOREST_HEIGHT_2020,FOREST_HEIGHT_2020,FOREST_HE
##### EXCEPTIONS FOR THE BOREAL #####
BOREAL_AGE = rasterio.open(args.BOREAL_AGE).read(1)
PRIMARY[(BOREAL_AGE <= 20) & (BOREAL_AGE > 0)] = 2 # LAYER ACTIVE ONLY IN BOREAL ZONES
YOUNG[(BOREAL_AGE <= 20) & (BOREAL_AGE > 0)] = 2 # LAYER ACTIVE ONLY IN BOREAL ZONES
del(BOREAL_AGE)
##### PRIMARY AND INTACT FOREST CONDITIONS #####
YOUNG[PRIMARY > 0] = 0
##### PLANTATIONS CONDITIONS #####
PLANTATIONS = rasterio.open(args.PLANTATIONS).read(1)
YOUNG[PLANTATIONS > 0] = 0
......@@ -132,6 +131,10 @@ ESA_FOREST_COVER = rasterio.open(args.ESA_FOREST_COVER).read(1)
YOUNG[ESA_FOREST_COVER == 0] = 0
del(ESA_FOREST_COVER)
##### PRIMARY AND INTACT FOREST CONDITIONS #####
YOUNG[PRIMARY > 0] = 0
YOUNG[(YOUNG > 2) | (YOUNG < 2)] = 0
###################################################################################
###################################################################################
##### OLD FOREST CONDITIONS - GIVEN VALUE 3 #######################################
......@@ -145,7 +148,7 @@ OLD = FOREST_HEIGHT_2020
del(FOREST_HEIGHT_2020)
OLD[OLD < 5] = 0
OLD[(OLD >= 5) & (FOREST_HEIGHT_2000 < 5)] = 0
OLD[(OLD >= 5) & (FOREST_HEIGHT_2000 > 5)] = 3
OLD[(OLD >= 5) & (FOREST_HEIGHT_2000 >= 5)] = 3
# FOREST_LOSS = rasterio.open(args.FOREST_LOSS).read(1)
# OLD[FOREST_LOSS > 0] = 0
......@@ -161,10 +164,6 @@ FOREST_COVER_LOSSYEAR = rasterio.open(args.FOREST_COVER_LOSSYEAR).read(1)
OLD[FOREST_COVER_LOSSYEAR > 0] = 0
del(FOREST_COVER_LOSSYEAR)
##### PRIMARY AND INTACT FOREST CONDITIONS #####
OLD[PRIMARY > 0] = 0
OLD[YOUNG > 0] = 0
##### PLANTATIONS CONDITIONS #####
PLANTATIONS = rasterio.open(args.PLANTATIONS).read(1)
OLD[PLANTATIONS > 0] = 0
......@@ -175,6 +174,11 @@ ESA_FOREST_COVER = rasterio.open(args.ESA_FOREST_COVER).read(1)
OLD[ESA_FOREST_COVER == 0] = 0
del(ESA_FOREST_COVER)
##### PRIMARY AND INTACT FOREST CONDITIONS #####
OLD[PRIMARY > 0] = 0
OLD[YOUNG > 0] = 0
OLD[(OLD > 3) | (OLD < 3)] = 0
###################################################################################
###################################################################################
##### AGGREGATING GLOBAL FOREST CLASSES AND OUTPUT TO RASTER ######################
......
country_summaries/IPCC_classes_DPS/IPCC_GEDI_Table4.7.py
+ 15
− 11
View file @ 11693b97
......@@ -20,7 +20,7 @@ import argparse
# wget https://storage.googleapis.com/earthenginepartners-hansen/GFC-2022-v1.10/Hansen_GFC-2022-v1.10_lossyear_10S_060W.tif -O /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/Hansen_GFC-2022-v1.10_lossyear_10S_060W.tif
##### RUN SCRIPT FROM COMMAND LINE #####
# python /projects/ADE_biomass_harmonization/biomass_harmonization/country_summaries/IPCC_classes_DPS/IPCC_GEDI_Table4.7.py --FOREST_HEIGHT_2000 /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/2000_10S_060W.tif --FOREST_HEIGHT_2020 /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/2020_10S_060W.tif --FOREST_LOSS /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/netloss_10S_060W.tif --ESA_FOREST_COVER /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/treecover2010_10S_060W.tif --FOREST_OLD_LOSS /projects/my-public-bucket/Data/Harris_et_al_PAPER/Old_loss/10S_060W_OLD_LOSS.tif --FOREST_COVER_LOSSYEAR /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/Hansen_GFC-2022-v1.10_lossyear_10S_060W.tif --BOREAL_AGE /projects/my-public-bucket/Data/Harris_et_al_PAPER/Forest_age_boreal/BOREAL_AGE/10S_060W_BOREAL_AGE.tif --PLANTATIONS /projects/my-public-bucket/Data/Harris_et_al_PAPER/WRI_Planted_Forest_Type/10S_060W_plantation_type_oilpalm_woodfiber_other_unmasked.tif --IFL_rasters /projects/my-public-bucket/Data/Harris_et_al_PAPER/Intact_Forest_Landscapes/IFL_rasters/10S_060W_IFL_2000.tif --FII_raster /projects/my-public-bucket/Data/Harris_et_al_PAPER/Forest_Integrity_Index/10S_060W_FII.tif --primary_forest_asia /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_Asia.tif --primary_forest_sa /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_SouthAmerica.tif --primary_forest_africa /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_Africa.tif --EcoCont /projects/my-public-bucket/Data/Harris_et_al_PAPER/EcoCont_v2/10S_060W_EcoCont.tif --tile 10S_060W --output_file /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/IPCC_10S_060W.tif
# python /projects/ADE_biomass_harmonization/biomass_harmonization/country_summaries/IPCC_classes_DPS/IPCC_GEDI_Table4.7.py --FOREST_HEIGHT_2000 /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/2000_10S_060W.tif --FOREST_HEIGHT_2020 /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/2020_10S_060W.tif --FOREST_LOSS /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/netloss_10S_060W.tif --ESA_FOREST_COVER /projects/my-public-bucket/Data/Harris_et_al_PAPER/ESA_WorldCover2021_v2/WC_10S_060W_C.tif --FOREST_COVER_LOSSYEAR /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/Hansen_GFC-2022-v1.10_lossyear_10S_060W.tif --BOREAL_AGE /projects/my-public-bucket/Data/Harris_et_al_PAPER/Forest_age_boreal/BOREAL_AGE/10S_060W_BOREAL_AGE.tif --PLANTATIONS /projects/my-public-bucket/Data/Harris_et_al_PAPER/WRI_Planted_Forest_Type_v2/10S_060W_plantation_type_oilpalm_woodfiber_other.tif --IFL_rasters /projects/my-public-bucket/Data/Harris_et_al_PAPER/Intact_Forest_Landscapes/IFL_rasters_2020/10S_060W_IFL_2020.tif --FII_raster /projects/my-public-bucket/Data/Harris_et_al_PAPER/Forest_Integrity_Index/10S_060W_FII.tif --primary_forest_asia /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_Asia.tif --primary_forest_sa /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_SouthAmerica.tif --primary_forest_africa /projects/my-public-bucket/Data/Harris_et_al_PAPER/Turubanova_Humid_tropical_primary_forest/10S_060W_Africa.tif --EcoCont /projects/my-public-bucket/Data/Harris_et_al_PAPER/EcoCont_v2/10S_060W_EcoCont.tif --tile 10S_060W --output_mask /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/Ages_10S_060W.tif --output_file /projects/my-public-bucket/Data/Harris_et_al_PAPER/_TEST/Classes_10S_060W.tif
######################################################################################
parse = argparse.ArgumentParser(description="Classifies world into GEZ, Continent and Forest Age classes")
......@@ -60,7 +60,7 @@ PRIMARY = FII_raster
del(FII_raster)
PRIMARY[PRIMARY <= 9600] = 0
PRIMARY[PRIMARY > 9600] = 1
PRIMARY[INTACT_FORESTS > 0] = 1
PRIMARY[INTACT_FORESTS == 1] = 1
del(INTACT_FORESTS)
PRIMARY[PRIMARY_FOREST_ASIA > 0] = 1
del(PRIMARY_FOREST_ASIA)
......@@ -93,6 +93,8 @@ ESA_FOREST_COVER = rasterio.open(args.ESA_FOREST_COVER).read(1)
PRIMARY[ESA_FOREST_COVER == 0] = 0
del(ESA_FOREST_COVER)
PRIMARY[(PRIMARY > 1) | (PRIMARY < 1)] = 0
# ###################################################################################
# ###################################################################################
# ##### YOUNG FOREST CONDITIONS - GIVEN VALUE 2 #####################################
......@@ -116,12 +118,9 @@ del(FOREST_COVER_LOSSYEAR,BACKUP_FOREST_HEIGHT_2020,FOREST_HEIGHT_2020,FOREST_HE
##### EXCEPTIONS FOR THE BOREAL #####
BOREAL_AGE = rasterio.open(args.BOREAL_AGE).read(1)
PRIMARY[(BOREAL_AGE <= 20) & (BOREAL_AGE > 0)] = 2 # LAYER ACTIVE ONLY IN BOREAL ZONES
YOUNG[(BOREAL_AGE <= 20) & (BOREAL_AGE > 0)] = 2 # LAYER ACTIVE ONLY IN BOREAL ZONES
del(BOREAL_AGE)
##### PRIMARY AND INTACT FOREST CONDITIONS #####
YOUNG[PRIMARY > 0] = 0
##### PLANTATIONS CONDITIONS #####
PLANTATIONS = rasterio.open(args.PLANTATIONS).read(1)
YOUNG[PLANTATIONS > 0] = 0
......@@ -132,6 +131,10 @@ ESA_FOREST_COVER = rasterio.open(args.ESA_FOREST_COVER).read(1)
YOUNG[ESA_FOREST_COVER == 0] = 0
del(ESA_FOREST_COVER)
##### PRIMARY AND INTACT FOREST CONDITIONS #####
YOUNG[PRIMARY > 0] = 0
YOUNG[(YOUNG > 2) | (YOUNG < 2)] = 0
###################################################################################
###################################################################################
##### OLD FOREST CONDITIONS - GIVEN VALUE 3 #######################################
......@@ -145,7 +148,7 @@ OLD = FOREST_HEIGHT_2020
del(FOREST_HEIGHT_2020)
OLD[OLD < 5] = 0
OLD[(OLD >= 5) & (FOREST_HEIGHT_2000 < 5)] = 0
OLD[(OLD >= 5) & (FOREST_HEIGHT_2000 > 5)] = 3
OLD[(OLD >= 5) & (FOREST_HEIGHT_2000 >= 5)] = 3
# FOREST_LOSS = rasterio.open(args.FOREST_LOSS).read(1)
# OLD[FOREST_LOSS > 0] = 0
......@@ -161,10 +164,6 @@ FOREST_COVER_LOSSYEAR = rasterio.open(args.FOREST_COVER_LOSSYEAR).read(1)
OLD[FOREST_COVER_LOSSYEAR > 0] = 0
del(FOREST_COVER_LOSSYEAR)
##### PRIMARY AND INTACT FOREST CONDITIONS #####
OLD[PRIMARY > 0] = 0
OLD[YOUNG > 0] = 0
##### PLANTATIONS CONDITIONS #####
PLANTATIONS = rasterio.open(args.PLANTATIONS).read(1)
OLD[PLANTATIONS > 0] = 0
......@@ -175,6 +174,11 @@ ESA_FOREST_COVER = rasterio.open(args.ESA_FOREST_COVER).read(1)
OLD[ESA_FOREST_COVER == 0] = 0
del(ESA_FOREST_COVER)
##### PRIMARY AND INTACT FOREST CONDITIONS #####
OLD[PRIMARY > 0] = 0
OLD[YOUNG > 0] = 0
OLD[(OLD > 3) | (OLD < 3)] = 0
###################################################################################
###################################################################################
##### AGGREGATING GLOBAL FOREST CLASSES AND OUTPUT TO RASTER ######################
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment