3 Jul
2020
3 Jul
'20
10:10 a.m.
Message posté par : ntouati
----------------------------------------
Bonjour,
Je rencontre un problème que j'arrive pas à résoudre. J'aurai besoin d'un peu
d'aide.
Mon objectif est d'addition des couches raster pour ne former qu'une seule couche
qui sera mobiliser pour un ré-échantillonnage (SAGA algorithme)
Voici l'extrait de mon script qui ne fonctionne pas :
-----------------
Code :
#Use raster calculator to create city fraction
entries = [] # initialize input list
rast1 = QgsRasterCalculatorEntry() # create first raster
rast1.ref = 'rast@1' # define first raster name
objBuild = QgsRasterLayer(build['OUTPUT']) # get raster as QGIS object
rast1.raster = objBuild # define raster
rast1.bandNumber = 1 # define raster band to use
entries.append(rast1) # add raster to inputs
rast2 = QgsRasterCalculatorEntry() # same thing for second raster
objRoad = QgsRasterLayer(road['OUTPUT'])
rast2.ref = 'rast@2'
rast2.raster =objRoad
rast2.bandNumber = 1
entries.append(rast2)
rast3 = QgsRasterCalculatorEntry() # same thing for third raster
objEnv = QgsRasterLayer(env['OUTPUT'])
rast3.ref = 'rast@3'
rast3.raster = objEnv
rast3.bandNumber = 1
entries.append(rast3)
# Get a path to save outputs as temporary files
getPathSave = objEmp.dataProvider().dataSourceUri() # get access path from any temporary
file
split1 = getPathSave.split('|')
split2 = split1[0]
split3 = split2[:-10] # delete file name at the end of the string
# define path to save result as temporary file
outputCalc = split3 + 'addition.tif'
outputCalc2 = split3 + 'additionResampl.tif'
calc = QgsRasterCalculator('rast@1+rast@2+rast@3', outputCalc, 'GTiff',
objBuild.extent(),
objBuild.width(), objBuild.height(), entries) #
configure raster calculator
calc.processCalculation() # launch calculation
# Set raster dimensions
resamplFracVille = processing.run("saga:resampling", {'INPUT':
outputCalc,'KEEP_TYPE':True,'SCALE_UP':0,'SCALE_DOWN':0,'TARGET_USER_XMIN
TARGET_USER_XMAX TARGET_USER_YMIN TARGET_USER_YMAX':extent,
'TARGET_USER_SIZE':resolution,'TARGET_USER_FITS':1,'TARGET_TEMPLATE':None,'OUTPUT':outputCalc2})
addFraction = iface.addRasterLayer(outputCalc2, "FRACTION_DE_VILLE")
-----------------
Le script me renvoie ce message :
-----------------
Code :
Incapable d’exécuter l’algorithme
Impossible de charger la couche source pour INPUT:
memory?geometry=Point&crs=EPSG:2154&field=Jeu:string(254,0)&field=ID:int8(10,0)&field=Nom:string(254,0)&field=Lat:double(23,15)&field=Lon:double(23,15)&field=Alt:int8(10,0)&field=Name:string(15,0)&field=_13_00_00:string(254,0)&field=_14_00_00:string(254,0)&field=_15_00_00:string(254,0)&field=_16_00_00:string(254,0)&field=_17_00_00:string(254,0)&field=_18_00_00:string(254,0)&field=_19_00_00:string(254,0)&field=_20_00_00:string(254,0)&field=_21_00_00:string(254,0)&field=_22_00_00:string(254,0)&field=_23_00_00:string(254,0)&field=_00_00_00:string(254,0)&field=_01_00_00:string(254,0)&field=_02_00_00:string(254,0)&field=_03_00_00:string(254,0)&field=_04_00_00:string(254,0)&field=_05_00_00:string(254,0)&field=_06_00_00:string(254,0)&field=_07_00_00:string(254,0)&field=_08_00_00:string(254,0)&field=_09_00_00:string(254,0)&field=_10_00_00:string(254,0)&field=_11_00_00:string(254,0)&field=_12_00_00:string(254,0)&field=_13_00_00_:string(254,0)
&field=_14_00_00_:string(254,0)&field=_15_00_00_:string(254,0)&field=_16_00_00_:string(254,0)&field=_17_00_00_:string(254,0)&field=_18_00_00_:string(254,0)&field=_19_00_00_:string(254,0)&field=_20_00_00_:string(254,0)&field=_21_00_00_:string(254,0)&field=_22_00_00_:string(254,0)&field=_23_00_00_:string(254,0)&field=_00_00_00_:string(254,0)&field=_01_00_00_:string(254,0)&field=_02_00_00_:string(254,0)&field=_03_00_00_:string(254,0)&field=_04_00_00_:string(254,0)&field=_05_00_00_:string(254,0)&field=_06_00_00_:string(254,0)&field=_07_00_00_:string(254,0)&field=_08_00_00_:string(254,0)&field=_09_00_00_:string(254,0)&field=_10_00_00_:string(254,0)&field=_11_00_00_:string(254,0)&field=_12_00_00_:straddition.tif
introuvable
-----------------
Le problème vient-il du "split" ? ou du stockage temporaire
des"outputs" ?
En vous remerciant d'avance,
N.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=334109#p334109
Pour y répondre : qgis_fr(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net