<div dir="ltr">Buongiorno a tutti,<div>sto lavorando con QGIS 3.32.0-Lima su Windows 11.</div><div>Sto cercando di fare il porting di uno script di R su QGIS creando un R script (<a href="https://docs.qgis.org/3.28/en/docs/user_manual/appendices/qgis_r_syntax.html">https://docs.qgis.org/3.28/en/docs/user_manual/appendices/qgis_r_syntax.html</a>)</div><div><br></div><div>Se lo faccio girare però ottengo un errore di enconding che non colgo:<br></div><div><br></div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote"><span style="color:rgb(255,0,0)">Traceback (most recent call last):<br></span><span style="color:rgb(255,0,0)">File "C:\Users/giano/AppData/Roaming/QGIS/QGIS3\profiles\SGI/python/plugins\processing_r\processing\algorithm.py", line 396, in processAlgorithm<br></span><span style="color:rgb(255,0,0)">output = RUtils.execute_r_algorithm(self, parameters, context, feedback)<br></span><span style="color:rgb(255,0,0)">File "C:\Users/giano/AppData/Roaming/QGIS/QGIS3\profiles\SGI/python/plugins\processing_r\processing\utils.py", line 281, in execute_r_algorithm<br></span><span style="color:rgb(255,0,0)">for line in iter(proc.stdout.readline, ''):<br></span><span style="color:rgb(255,0,0)">File "C:\OSGeo4W\apps\Python39\lib\codecs.py", line 322, in decode<br></span><span style="color:rgb(255,0,0)">(result, consumed) = self._buffer_decode(data, self.errors, final)<br></span><span style="color:rgb(255,0,0)">UnicodeDecodeError: 'utf-8' codec can't decode byte 0xe8 in position 21: invalid continuation byte<br></span><span style="color:rgb(255,0,0)">Execution failed after 37.93 seconds</span></blockquote></div></blockquote><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div><p style="margin:0px"><i><br></i></p></div></blockquote><div><p style="margin:0px">Il codice che sto usando è il seguente e ho due raster in ingresso, uno con valori Integer e l'altro è float.</p><p style="margin:0px">Di fatto l'obiettivo è di fare delle statistiche zonali tra raster calcolando la mediana dei valori del raster A per ciascuna classe del raster B (cosa che il tool "Raster Layer Zonal statistics" non fa) e riclassifica i valori del raster A in funzione del valore calcolato.</p></div><i><br>##My scripts=group</i><br><div><p style="margin:0px"><i><br>##LC=raster<br>##issue=raster<br>##output= output raster<br><br>library(raster)<br>library(sf)<br>library(sp)<br><br>x<-resample(issue,LC, method='ngb')<br>issue <-x<br><br>issue <- crop(issue, extent(LC))  <br>res(issue) <- res(LC) <br>extent(issue) <- extent(LC)<br> <br>s<- issue<br>r<- LC<br>result <- zonal(s, r, 'median')<br><br>r2 <- reclassify(r, result)<br><br></i></p><p style="margin:0px"><i>rc <- function(x) {<br>  <br>  x1 <- x[1]  ##valore pixel della prima banda<br>  x2 <- x[2] ##valore pixel 2 banda<br>  ifelse(x1 <=  x2, 0,<br>  ifelse(x1 >  x2, 1, NA)) }<br><br>s2 <- stack(s,r2)<br><br>beginCluster(6) ##ASSIGN CORE<br>z1 <- clusterR(s2, calc, args=list(fun=rc))<br>endCluster()</i></p><p style="margin:0px"><i><br>output<- z1</i><span style="color:rgb(255,0,0)"><br></span></p><p style="margin:0px"><i><br></i></p><p style="margin:0px">Qualcuno ha avuto problemi simili di encoding o di script R su QGIS? </p><p style="margin:0px">Grazie mille</p><p style="margin:0px"><br></p><p style="margin:0px">Federico</p></div></div>