[Utenti] [R script] output sempre vuoto
Totò Fiandaca
pigrecoinfinito a gmail.com
Mar 1 Dic 2015 19:29:12 CET
SO - windows 10 64bit.
R 3.2.2 da cui:
> .libPaths()
[1] "C:/Users/nomeutente/Documents/R/win-library/3.2"
[2] "C:/Program Files/R/R-3.2.2/library"
alcune osservazioni (prove effettuate con geoalgoritmi: Histogram, Summary
statics):
1. layer 'postgis_view' con tipo di dato 'qlonglong' dà sempre output vuoto;
2. opzioni di processing R utilizzate:
a. Cartella R : C:\Program Files\R\R-3.2.2;
b. Cartella della libreria di R
: C:\Users\nomeutente\Documents\R\win-library\3.2;
c. Cartella script
R: C:\Users\nomeutente\.qgis2\processing\rscripts.
3. nei 'log' è sempre presente:
per il geoalgoritmo *Summary_statistics:*
2015-12-01T18:44:32 0 Comandi di esecuzione R
options("repos"="http://cran.at.r-project.org/")
.libPaths("C:/Users/Salvatore/Documents/R/win-library/3.2")
tryCatch(find.package("rgdal"), error=function(e) install.packages("rgdal",
dependencies=TRUE))
tryCatch(find.package("raster"), error=function(e)
install.packages("raster", dependencies=TRUE))
library("raster")
library("rgdal")
Layer =
readOGR("C:/Users/SALVAT~1/AppData/Local/Temp/processinga6c0ea2d8af34a348d818540a1ed682e/6d8cc8a91b8b4dcd8538863e7c77c7c6",layer="
*ganagrafiche*")
Field="nri_anagr"
Summary_statistics<-data.frame(rbind(sum(Layer[[Field]]),
length(Layer[[Field]]),
length(unique(Layer[[Field]])),
min(Layer[[Field]]),
max(Layer[[Field]]),
max(Layer[[Field]])-min(Layer[[Field]]),
mean(Layer[[Field]]),
median(Layer[[Field]]),
sd(Layer[[Field]])),row.names=c("Sum:","Count:","Unique values:","Minimum
value:","Maximum value:","Range:","Mean value:","Median value:","Standard
deviation:"))
colnames(Summary_statistics)<-c(Field)
*Summary_statistics*
per il geoalgoritmo *Histogram:*
2015-12-01T18:40:31 0 Comandi di esecuzione R
options("repos"="http://cran.at.r-project.org/")
.libPaths("C:/Users/Salvatore/Documents/R/win-library/3.2")
tryCatch(find.package("rgdal"), error=function(e) install.packages("rgdal",
dependencies=TRUE))
tryCatch(find.package("raster"), error=function(e)
install.packages("raster", dependencies=TRUE))
library("raster")
library("rgdal")
Layer =
readOGR("C:/Users/SALVAT~1/AppData/Local/Temp/processinga6c0ea2d8af34a348d818540a1ed682e/e3bde974e7414e1e8503f803196e8dd3",layer="
*vnroconcessionitomba2*")
Field="nro"
png("C:/Users/SALVAT~1/AppData/Local/Temp/processinga6c0ea2d8af34a348d818540a1ed682e/bd2ac4e27919468fb6bec4ff328031bd/RPLOTS.html.png")
hist(Layer[[Field]],main=paste("Histogram of",Field),xlab=paste(Field))
dev.off()
la cosa che non capisco negli script di sopra: perchè scrive in modo errato
i nomi dei layer nel readOGR? nel primo caso il layer è: 'g_anagrafiche' e
non 'ganagrafiche' ; nel secondo è: 'v_nro_concessioni_tomba2' e non '
*vnroconcessionitomba2'*; praticamente toglie '_'.
concludendo:
se utilizzo layer postigis che non sono viste tutto procede bene anche se
ottengo i log di sopra;
per il resto sembra funzionare tutto bene!!!
saluti.
Il giorno 1 dicembre 2015 18:12, Totò Fiandaca <pigrecoinfinito a gmail.com>
ha scritto:
>
> Il giorno 1 dicembre 2015 17:22, matteo <matteo.ghetta a gmail.com> ha
> scritto:
>
>> ti funziona con questi suggerimenti?
>>
>> cosi in caso si possono aggiungere al manuale di QGIS e fare un favore a
>> tutti dopo la tua faticata.. :)
>>
>
> la risposta è articolata:
> credo di aver capito i percorsi da mettere per configurare R in QGIS, ma
> resta il dubbio: perchè per alcuni layer funziona bene e per altri no?
> perchè per alcuni layer dà sempre pagina bianca e nel log segnala errori su
> 'rgdal' e 'raster'?
> la cosa strana è che per alcuni layer dove dà l'output corretto compare
> sempre, nel 'log' ,il messaggio di errore di 'rgdal' e 'raster' .
>
> ci sto lavorando ancora, spero di trovare presto la soluzione.
>
>
>
> --
> *Salvatore Fiandaca*
> *mobile*.:+39 327.493.8955
> *m*: *pigrecoinfinito a gmail.com <pigrecoinfinito a gmail.com>*
> 43°51'0.54"N 10°34'27.62"E - EPSG:4326
>
--
*Salvatore Fiandaca*
*mobile*.:+39 327.493.8955
*m*: *pigrecoinfinito a gmail.com <pigrecoinfinito a gmail.com>*
43°51'0.54"N 10°34'27.62"E - EPSG:4326
-------------- parte successiva --------------
Un allegato HTML è stato rimosso...
URL: <http://lists.qgis.it/pipermail/utenti/attachments/20151201/2df4cbb9/attachment-0001.html>
Maggiori informazioni sulla lista
Utenti