In fact this is not a problem with the implementation but maybe a problem of comprehension.
I dont understand why the posterior probability isnt equal to zero whereas there is not landslide in this area.
import numpy as np
import rasterio
#import math
import csv
import os
from itertools import product
def raster_to_numpy(fileraster):
with rasterio.open(fileraster) as src:
if src.meta['nodata'] == None:
print("Warnings : Attention, valeur nodata non renseignée dans le fichier entrée")
array = src.read(1, masked=True)
nodata = np.array(array.fill_value, dtype=array.dtype)
array = array.filled()
else:
array = src.read(1)
meta = src.meta
meta = meta.values()
nodata = meta[-1]
return array,nodata
def variable(raster,nodata):
variable = list(np.unique(raster))
if (nodata in variable):
variable = np.delete(variable,variable.index(nodata))
return variable
FileDossierRaster = r"C:\Users\user\Desktop\Boulot\Dossier charge"
fileglissement = r"C:\Users\user\Desktop\Boulot\landslides_tests\ASCII\scarp_rb.asc"
fileout = r"C:\Users\user\Desktop\Boulot\Fichiers_csv\sdmuc1.csv"
def attributes(FileDossierRaster,fileglissement,fileout):
liste_raster = []
liste_nodata = []
liste_var = []
for i in os.listdir(FileDossierRaster):
if not os.path.isdir(i):
new_path = os.path.join(FileDossierRaster,i)
raster_var,nodata = raster_to_numpy(new_path)
liste_raster.append(raster_var)
liste_nodata.append(nodata)
v_ras = variable(raster_var,nodata)
liste_var.append(v_ras)
raster_gli,nodata_gli = raster_to_numpy(fileglissement)
taille = [len(i) for i in liste_var]
n = np.prod(taille)
valeurs = np.zeros((n,11+len(liste_raster)))
combi = list(product(*liste_var)) ## ici je fais les combinaisons
for i in range(n):
valeurs[i,0] = i+1
valeurs[i,1:len(combi[0])+1] = combi ## ici je les mets dans mon tableau
t = 1
for k in range(len(liste_raster)):
t *= np.where(liste_raster == valeurs[i,k+1],1,0) ## possibilé erreur au dernier
c = np.count_nonzero(t == 1)
g = np.where(raster_gli == 1,1,0)
s = np.count_nonzero(t*g == 1)
valeurs[i,len(combi[0])+1] = c
valeurs[i,len(combi[0])+2] = s
###########writing csv file#################
entete = ['numero']
for i in range(len(liste_var)):
v = 'variable' + str(i)
entete.append(v)
entete = entete + ['intersection','Cellules_affectees','Proba_posteriori','c','d','e','f','g']
with open(fileout,'w') as f:
w = csv.writer(f, delimiter = ';')
w.writerow(entete)
for line in zip(valeurs[:,0],valeurs[:,1],valeurs[:,2],valeurs[:,3],valeurs[:,4],valeurs[:,5],valeurs[:,6],valeurs[:,7]):
w.writerow(line)
#######################################
attributes(FileDossierRaster,fileglissement,fileout)
