[Gfoss] QGIS - georeferencer

[Pietro Blu Giandonato]

Saluti a tutti...

Spero di poter dare un contributo utile con quella che è la mia esperienza.
In effetti, la questione ruota tutta attorno al RMS "tollerabile", che poi è
in funzione della scala della carta scansionata e della sua risoluzione di
acquisizione. Giustamente, buona norma è mantenersi attorno al pixel. Dato
che il RMS viene espresso in unità di mappa (nel nostro caso in metri), è
necessario calcolare a quanti metri equivale un pixel dell'immagine da
rettificare, ovvero la sua risoluzione geometrica reale.

Quando si georiferiscono immagini raster di qualunque tipo, che si tratti di
foto aeree stampate o cartografia di qualunque genere scansionata, bisogna
assolutamente conoscerne la risoluzione di acquisizione e la scala di
rappresentazione.

Mando ad esempio, oltre alla risoluzione di acquisizione (72dpi) della sua
carta storica, non ha citato qual'è la scala di rappresentazione, ma
possiamo fare un esempio ipotizzando quest'ultima. Mettiamo essa sia 1:
50.000, per calcolare la risoluzione reale del pixel procediamo in questo
modo:

1) sapendo che in scala 1:50.000, 1 cm = 500 m e 1 inch = 2,54 cm
2) per una risoluzione di 72dpi, 1 pixel = 0,01 inch = 0,0254 cm
3) la risoluzione reale di un pixel sarà quindi Rr = (0,0254 cm * 500 m) / 1
cm = 12,7 m

E' consigliabile dunque in fase di georeferenziazione sceglie i GCP
ottenendo un RMS che stia attorno a questo valore, ma non sempre è cosa
facile. Bisogna infatti tenere conto delle deformazioni che la carta ha
subito nel tempo (soprattutto se è antica) o se è una riproduzione di un
originale (collage di più carte o fotocopia di fotocopia di fotocopia...), o
ancora della metodologia di acquisizione utilizzata (scanner piano di grande
formato (e già va bene) o peggio più scansioni mosaicate successivamente con
un SW di fotoritocco).

Mi pare quindi fondamentale avere ben chiara l'origine dell'immagine da
processare, e nel caso questa abbia subito molte "deformazioni", incaponirsi
nel cercare di ottenere un RMS che stia attorno al pixel è pura chimera. SI
possono effettuare tre o quattro procedure di georef, magari scegliendo
punti diversi, ma se il risultato rimane invariato è meglio rassegnarsi.

A presto,
Pietro Blu




Il 19/08/07, Marco Cerruti < cerruti.marco a gmail.com> ha scritto:
>
> Non ricampionare, semplicemente calcola a quanto corrisponde il singolo
> pixel (conosci le dimensioni dell'area ed il numero di pixel dell'immagine)
> ed individua un set di GCP che ti avvicina l'errore alle dimensioni del
> pixel.
>
> Tra l'altro se porti l'immagine da 72 a 400 dpi il pixel corrisponde a
> dimensioni reali più piccole, e lo stesso root mean square error nella
> georef diventa ancora più grande rispetto al pixel ... no? Comunque faccio
> un po' fatica a parlare di dpi, mi sembra un passaggio mentale in più, per
> me è meglio calcolare le dimensioni del pixel e parlare di metri e di
> ricampionamenti da x metri a y metri.
>
>
> M
>
> PS: per Steko, sì una pagina wiki sarebbe l'ideale, vedo cosa posso fare.
> Mi pare che anche il wiki di Open Aerial Map stesse creando pagine
> didattiche su argomenti del genere.
>
> Il 19/08/07, mando <mandoluca a gmail.com> ha scritto:
> >
> >
> >
> >
> > mando wrote:
> > >
> > > Dunque, per abbassare l'errore sotto il pixel:
> > >
> > > la cartografia raster 1:5000 che uso come riferimento per prendere le
> > > coordinare è un file a con risoluzione al 400 pxel/inch, mentre la
> > carta
> > > storica me l'hanno passata a 72 pixel/inch.
> > > E' meglio ricampionare l'immagine e portarla a 400?
> > >
> >
> > Ho provato ad aumentare la risoluzione:
> >
> > Helmert:
> > Se si aumenta la risoluzione del file e si riesegue l'inserimento dei
> > valori
> > dei GCP, la helmert riporta il file modificato alla risoluzione di 72
> > pixel/inch e il file viene caricato con dimensioni abnormi rispetto al
> > rapporto con la cartografia caricata.
> >
> > Lineare:
> > il file non è modificato nelle sue dimensioni ne ruotato (come dovrebbe,
> > se
> > ho ben capito)
> > --
> > View this message in context: http://www.nabble.com/QGIS---georeferencer-tf4248844.html#a12222158
> >
> > Sent from the Gfoss mailing list archive at Nabble.com.
> >
> >
> > _______________________________________________
> > Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
> >
> > Gfoss a faunalia.com
> > http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss
> >
>
>
> _______________________________________________
> Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
>
> Gfoss a faunalia.com
> http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss
>
>
-------------- parte successiva --------------
Un allegato HTML ? stato rimosso...
URL: http://www.faunalia.com/pipermail/gfoss/attachments/20070819/65705f7f/attachment-0001.htm