[QGIS-it-user] mostrare flussi Istat pendolarismo
gabriele.dipalma a uniba.it
gabriele.dipalma a uniba.it
Dom 6 Ott 2019 01:25:46 PDT
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 a lists.osgeo.org ha scritto:
> Invia le richieste di iscrizione alla lista QGIS-it-user all'indirizzo
> qgis-it-user a lists.osgeo.org
>
> Per 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 a lists.osgeo.org
>
> Puoi contattare la persona che gestisce la lista all'indirizzo
> qgis-it-user-owner a lists.osgeo.org
>
> Se 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 a gmail.com>
> To: qgis-it-user a 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 a gmail.com>
> Content-Type: text/plain; charset=utf-8; format=flowed
>
> mi 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 joint
>>
>> Non so se è il sistema miglire, ma ha funzionato.
>>
>> Iacopo
>>
>> Il giorno dom, 29/09/2019 alle 11.44 +0200, gabriele.dipalma a 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 a lists.osgeo.org
>>> https://lists.osgeo.org/mailman/listinfo/qgis-it-user
>> _______________________________________________
>> QGIS-it-user mailing list
>> QGIS-it-user a 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 a gmail.com>
> To: Utenti QGIS it <qgis-it-user a lists.osgeo.org>
> Subject: [QGIS-it-user] PYQGIS riferirsi a una banda in in raster
> multibanda
> Message-ID:
> <CAO5DHxFPUNaN4A0jKrcxXsntOm9o3F+TJh+Zb6PGS5b_FJ=00w a 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 esempio
>
> rasterLayer.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>
> <#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 a gmail.com>
> To: Giacomo Fontanelli <giacomofontanelli76 a gmail.com>
> Cc: Utenti QGIS it <qgis-it-user a lists.osgeo.org>
> Subject: Re: [QGIS-it-user] PYQGIS riferirsi a una banda in in raster
> multibanda
> Message-ID:
> <CAFO80_q2xD+aZuqyQPe5Fw4Lj6k93uq7WycXuXGA3S+p-7apWg a 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 a 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 esempio
>>
>> rasterLayer.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 a 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 a 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
Maggiori informazioni sulla lista
QGIS-it-user