Buongiorno Luca e grazie della risposta. Ho letto e riletto quel post della Graser, e provato e riprovato ad adattarlo al mio caso, che poi è il seguente.
In un cvs ho il capo ORIG con la sezione di censimento da cui si parte, DEST con la sezione di censimento a cui si arriva e poi il campo TOTALE con l'ampiezza delle varie linee che dovrebbero collegare partenza e arrivo.
Poi ho un shp con le coordinate dei punti di ciascuna sezione di censimento di mio interesse. Il campo che unisce i due file è SEZ nello shp, che dovrei unire con ORIG e DEST (i cui valori sono inclusi in quelli di SEZ).
La Graser dice di creare un file virtuale e poi scrivere queste linee di comando:
SELECT StartID, DestID, Weight,
make_line(a.geometry, b.geometry)
FROM edges
JOIN nodes a ON edges.StartID = a.ID
JOIN nodes b ON edges.DestID = b.ID
WHERE a.ID != b.ID
Nel mio caso sarebbe:
SELECT ORIG, DEST, TOTALE make_line(a.geometry, b.geometry)
-- il mio primo dubbio è su quello che dovrei scrivere nella parentesi di make_line. Con geometry si indicano in generale i centroidi, credo. ma 'a' e 'b' occorrono o devo adattarli al mio caso? --
FROM Bari2011_i -- mio file csv
JOIN SezIta_2011_Point_WGS84 a ON Bari2011_i.ORIG = a.SEZ
JOIN SezIta_2011_Point_WGS84 b ON Bari2011_i.DEST = b.SEZ
WHERE a.SEZ != b.SEZ
Nelle linee di JOIN c'è qualcosa che non va, ovviamente. Python mi dà errore e le mie conoscenze sono estremamente scarse.
Se qualcuno potesse aiutarmi, gliene sarei estremamente grato. Mi interessa, per ricerca, capire come si spostano le persone verso la mia area di interesse nel comune di Bari. Posso anche allegare file.
Buona domenica
qgis-it-user-request@lists.osgeo.org ha scritto:
Invia le richieste di iscrizione alla lista QGIS-it-user all'indirizzo
qgis-it-user@lists.osgeo.orgPer iscriverti o cancellarti attraverso il web, visita
https://lists.osgeo.org/mailman/listinfo/qgis-it-user
oppure, via email, manda un messaggio con oggetto `help' all'indirizzo
qgis-it-user-request@lists.osgeo.orgPuoi contattare la persona che gestisce la lista all'indirizzo
qgis-it-user-owner@lists.osgeo.orgSe rispondi a questo messaggio, per favore edita la linea dell'oggetto
in modo che sia più utile di un semplice "Re: Contenuti del digest
della lista QGIS-it-user..."Argomenti del Giorno:
1. Re: Due cose: mostrare flussi Istat pendolarismo e comprare
dominio per webGIS (luca barbadoro)
2. PYQGIS riferirsi a una banda in in raster multibanda
(Giacomo Fontanelli)
3. Re: PYQGIS riferirsi a una banda in in raster multibanda
(Luigi Pirelli)----------------------------------------------------------------------
Message: 1
Date: Fri, 4 Oct 2019 07:54:24 +0200
From: luca barbadoro <luca.barbadoro@gmail.com>
To: qgis-it-user@lists.osgeo.org
Subject: Re: [QGIS-it-user] Due cose: mostrare flussi Istat
pendolarismo e comprare dominio per webGIS
Message-ID: <85a618bf-0750-2b37-2ae1-87d3ea52198b@gmail.com>
Content-Type: text/plain; charset=utf-8; format=flowedmi inserisco alla discussione!
per collegare visivamente i centroidi (o comunque gli elementi puntuali
a cui devono essere associati i dati del pendolarismo, siano comuni, ma
anche con zone censuarie è possibile), si può realizzare una FLOW map
direttamente in qgis senza plug in come spiegato da Anita Graser in
questo suo post:https://anitagraser.com/2019/05/04/flow-maps-in-qgis-no-plugins-needed/
Luca B
Il 29/09/19 22:12, Iacopo ha scritto:
Ciao Gabriele,
anni fa ho fatto così:
- ho creato un grafo unendo i centridi dei comuni interessati e creando
un campo con la combinazione dei codici dei 2 comuni collegati dalla
linea (in postgres ST_MakeLine (a.the_geom, b.the_geom) AS the_geom,
a.codcom || b.codcom AS codice, a.nome AS origine, b.nome AS
destinazione)
- dalla tabella ISTAT ho creato un campo concatenando codice origine e
codice destinazione
- poi basta una jointNon so se è il sistema miglire, ma ha funzionato.
Iacopo
Il giorno dom, 29/09/2019 alle 11.44 +0200, gabriele.dipalma@uniba.it
ha scritto:Cari tutti,
sono per favore avrei bisogno di una mano per due cose: la prima è
mostrare i flussi dei pendolari presi dai censimenti Istat e la
seconda sono le caratteristiche di uno spazio web dove vorrei
esportare il mio GIS.Per la prima, conosco Flow Mapper ma, per quantità di dati come
quelli
che vorrei mostrare (più o meno una matrice 200x200), faccio prima
a
disegnarla a mano. Ci sarebbero modi più semplici di
visualizzazione?
Ho visto che può essere usato Python, ma non conosco il suo
linguaggio. Non so, se qualcuno potesse darmi un suggerimento
gliene
sarei grato.La seconda richiesta è questa: sto creando un GIS con
caratteristiche
socio-demografiche di un territorio italiano e ricostruzione di un
programma di rigenerazione urbana di una amministrazione comunale.
Il
GIS è frutto di dati di enti pubblici (ISTAT, Agenzia delle
Entrate),
e frutto di un'analisi qualitativa e quantitativa che sto
conducendo
sul campo. Vorrei comprare uno spazio sul web dove tra pochi mesi
esportare il mio lavoro, creando un Public GIS, potreste suggerirmi
caratteristiche del servizio da acquistare ed eventualmente proprio
qualche provider di questo tipo di servizi (spazio web compatibile
con
QGIS). Essendo economista e sociologo non ne capisco granché.Ad ogni modo grazie e una buona domenica
Gabriele Di Palma
PHD Candidate in Human Relations Sciences
University of Bari Aldo Moro_______________________________________________
QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/qgis-it-user_______________________________________________
QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/qgis-it-user------------------------------
Message: 2
Date: Fri, 4 Oct 2019 12:20:41 +0200
From: Giacomo Fontanelli <giacomofontanelli76@gmail.com>
To: Utenti QGIS it <qgis-it-user@lists.osgeo.org>
Subject: [QGIS-it-user] PYQGIS riferirsi a una banda in in raster
multibanda
Message-ID:
<CAO5DHxFPUNaN4A0jKrcxXsntOm9o3F+TJh+Zb6PGS5b_FJ=00w@mail.gmail.com>
Content-Type: text/plain; charset="utf-8"salve a tutti
Ho caricato un raster multibanda in qgis tramite
rasterLayer = QgsRasterLayer("/path/file.tif")
ho controllato il numero di bande tramite
nBand = rasterLayer.bandCount()
adesso vorrei potermi riferire a una particolare banda del rasterLayer,
tramite qualcosa come ad esempiorasterLayer.band(i)
ma 'QgsRasterLayer' non ha l'attributo 'band'
Sapete coe fare?
Grazie
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
priva di virus. www.avast.com
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
-------------- parte successiva --------------
Un allegato HTML è stato rimosso...
URL: <http://lists.osgeo.org/pipermail/qgis-it-user/attachments/20191004/9d39c6fa/attachment-0001.html>------------------------------
Message: 3
Date: Fri, 4 Oct 2019 12:27:30 +0200
From: Luigi Pirelli <luipir@gmail.com>
To: Giacomo Fontanelli <giacomofontanelli76@gmail.com>
Cc: Utenti QGIS it <qgis-it-user@lists.osgeo.org>
Subject: Re: [QGIS-it-user] PYQGIS riferirsi a una banda in in raster
multibanda
Message-ID:
<CAFO80_q2xD+aZuqyQPe5Fw4Lj6k93uq7WycXuXGA3S+p-7apWg@mail.gmail.com>
Content-Type: text/plain; charset="utf-8"https://docs.qgis.org/testing/en/docs/pyqgis_developer_cookbook/raster.html
Luigi Pirelli
**************************************************************************************************
* LinkedIn: https://www.linkedin.com/in/luigipirelli
* Stackexchange: http://gis.stackexchange.com/users/19667/luigi-pirelli
* GitHub: https://github.com/luipir
* Book: Mastering QGIS3 - 3rd Edition
<https://www.packtpub.com/eu/application-development/mastering-geospatial-development-qgis-3x-third-edition>
* Hire a team: http://www.qcooperative.net
**************************************************************************************************On Fri, 4 Oct 2019 at 12:20, Giacomo Fontanelli <
giacomofontanelli76@gmail.com> wrote:salve a tutti
Ho caricato un raster multibanda in qgis tramite
rasterLayer = QgsRasterLayer("/path/file.tif")
ho controllato il numero di bande tramite
nBand = rasterLayer.bandCount()
adesso vorrei potermi riferire a una particolare banda del rasterLayer,
tramite qualcosa come ad esempiorasterLayer.band(i)
ma 'QgsRasterLayer' non ha l'attributo 'band'
Sapete coe fare?
Grazie
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> Mail
priva di virus. www.avast.com
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
<#m_3388865954445822882_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
_______________________________________________
QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/qgis-it-user-------------- parte successiva --------------
Un allegato HTML è stato rimosso...
URL: <http://lists.osgeo.org/pipermail/qgis-it-user/attachments/20191004/bedcaef6/attachment-0001.html>------------------------------
Subject: Chiusura del digest
_______________________________________________
QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/qgis-it-user------------------------------
Fine di Digest di QGIS-it-user, Volume 48, Numero 4
***************************************************
Gabriele Di Palma
PHD Candidate in Human Relations Sciences
University of Bari Aldo Moro