QGIS 2. proccesing e output degli scripi di esempio di R

[Ummarino]

Sto sperimentando Processing e gli script di R.
In linux Ubuntu (LTS "Precise") fila tutto liscio.
Con Windows 7 ...dopo un rapido aggiornamnento dei pachetti di R (3.0.2) e il flag sulla opozione "64bit version"  dei setting di R-sextante (Proccesing) ..fila tutto liscio
In ufficio Windows XP (bleach..!) R (3.0.2) QGIS 2.0.1 (aggiornatissima via OSGEO) mi da problemi: finestra di output (R Output) comletamente bianca...
Qualcuno ha lo stesso problema

[Mattia De Agostino]

Ciao,
io purtroppo posso darti una risposta solo parziale.

Ho un sistema operativo WindowsXP, con QGIS 2.0.1-Dufour (aggiornata via OSGEO), ma con R 2.15.2 (non così aggiornata ma...mi sono affezionato!) e sistema operativo 32 bit, e funziona tutto.

Per esperienza, il pacchetto di Processing si interfaccia con R anche per lo scarico e l'aggiornamento dei pacchetti di R richiesti (RGDAL e simili), quindi verifica che proxy e impostazioni internet siano correttamente settate sia in QGIS che in R.
Mi ricordo di aver tribolato un po' per inserire il proxy in R. Avevo risolto seguendo queste istruzioni:
http://cran.r-project.org/bin/windows/rw-FAQ.html#The-Internet-download-functions-fail_002e

Ciao!
Mattia

[Ummarino]

Grazie..la soluzione del proxy lo adottata per prima, a causa del proxy "aziendali" il problema sussisteva già per il comando install.packages che poi ho risolto come da guida.
Ho installato il pachetto proj2 i grafici li visualizza (ad esmpio l'example script Field histogram) ma l'output dell' example script Field summary statistc no....
Se puoi mi fai un listato dei pachetti di R che hai istallato...così vedo se mi manca qualcosa

[Mattia De Agostino]

Ciao, questi i pacchetti installati sul mio R (come detto, ricorda che si tratta della versione 2.15).

ape, base, boot, chron, class, cluster, coda, codetools, colorspace, compiler, DAAG, datasets, deldir, dichromat, digest, evaluate, foreign, Formula, fts, ggplot2, graphics, grDevices, grid, gtable, hexbin, Hmisc, igraph, igraphdata, its, KernSmooth, labeling, lattice, latticeExtra, mapproj, maps, maptools, MASS, matlab, Matrix, methods, mgcv, , mondate, munsell, ncdf, nlme, nnet, nws, parallel, plyr, proto, proj4, quadprog, quantreg, raster, rasterVis, RColorBrewer, reshape2, rgdal, rgeos, rgl, rlecuyer, rpart, rsprng, sandwich, scales, snow, sp, SparseM, spatial, splines, stats, stats4, stringr, strucchange, survival, tcltk, testthat, timeDate, timeSeries, tis, tools, tripack, tseries, utils, XML, xts, zoo

Ho provato il comando, ed in effetti restituisce l'output HTML vuoto nella finestra di Processing.
Non sono sicuro però che debba riempire qualcosa (nelle versioni da te provate sugli altri SO si vedeva qualcosa?). Nel codice dello script non mi sembra che venga dato nulla in output a video.

Le statistiche vengono calcolate, e le trovi nella finestra di Log dello script:


Ciao!
Mattia

[Ummarino]

Ciao...stasera riprovo a casa su Ubuntu e W7, così ti faccio lo screenshot, purtroppo sono in ufficio e ho solo xp.
A casa cmq mi riportava i valori nella finestra di output (sia linux che Win7) mentre in ufficio (XP) mi succede lo stesso che a te ..log delllo script ok ma restituzione di html no.
Ne approfitto per chiederti una cosa...su R, ti ricordo che sono un newbyissimo: ho modificato un pachetto, ho contattato via mail il "Maintainer" per chiedergli un aiutino (una verifica) ma non mi ha risposto.
Se hai tempo, e se puoi farlo, mi daresti una occhiata al mio "pachetto".

[Ummarino]

Ecco l'output di linux Kubuntu ...come vedi quà funziona

[Ummarino]

Ti illustro il mio vero obbiettivo e "problema" di R + QGIS

Sto provando con Processing uno script di R ma mi da questo risultato:

writeOGR(prova,"/tmp/processing/067a68c63d9b450db70141f6dc25dd3a/prova.shp","prova", driver="ESRI Shapefile")
Errore in inherits(obj, "Spatial") : oggetto "prova" non trovato
Calls: writeOGR -> stopifnot -> inherits
Esecuzione interrotta

Lo script originale di R in realtà già genera uno shape file, quindi credo che devo modificarlo e tagliare le ultime istruzioni dello script e sostituirle con delle istruzioni che generino un oggetto che possa poi essere trasformato e importato dalla funzione writeOGR.
Come posso generare questo "oggetto"...? quali sono le librerie e le istruzioni ?
Ti faccio vedere qua lo script (nella versione R scripts di "Processing-QGis"):

1##[Crime Point Analisys] = group
2##points= vector
3##prova=output vector
4library(maptools)
5library(aspace)
6P1 <- (points)
7P1P <- as (P1, "SpatialPoints")
8P2P_df <- as.data.frame(P1P)
9calc_2sde(id=1, filename="/home/alessandro/RLab/Tempo/P2P_ouput.text", centre.xy=NULL, calccentre=TRUE, weights=NULL, points=P2P_df, verbose=FALSE)
10sde_shp <- convert.to.shapefile(sdeloc, sdeatt, "id", 5)
11write.shapefile(sde_shp, "/home/alessandro/RLab/Tempo/SED2_ouput", arcgis=T)

Le righe incriminate sono la 10 e la 11, ma a me interessa disegnare il file P2P_output.txt della riga 9  (4 colonne di cui le ultime due sono coordinate di punti da interpolare e formano un poligono ellissoidico).

Ciao ...