Come si usa v.kernel di GRASS?

Sono confuso su come usare il kernel v di GRASS.

Ho uno strato vettoriale di circa 2,5 milioni di punti. Voglio creare una mappa di calore usando v.kernel per mostrare le concentrazioni, dal momento che ho istanze variabili con punti sovrapposti, a volte enormi sovrapposizioni.

Ho già ottenuto questo livello vettoriale in GRASS e viene visualizzato bene.

Ho provato ad usare il comando v.kernel di GRASS sulla base di ciò che ho visto qui e su altri forum, e non riesco a farlo fare altro oltre all'output di un raster che è solo un quadrato rosa.

Ecco il comando che sto usando:

v.kernel --verbose input=master_grass7 output=master_grass7a_heatmap stddeviation=.0001

Ho variato lo stddeviation a tutti i tipi di valori da 1000000 a .000001 e non ha avuto alcun effetto.

Ho letto ripetutamente la documentazione di v.kernel e non capisco davvero a cosa stia arrivando. Almeno, le istruzioni sono su concetti esoterici, niente di pratico. Ho anche controllato il codice sorgente e non lo capisco nemmeno davvero. Sì, posso leggere C. Il problema è che dipende da molte cose definite altrove in GRASS GIS.

Ho anche fatto molte ricerche su Google e non riesco a trovare una guida completa. Tutto quello che sto ricevendo sono copie sparse della pagina doc / man di v.kernel o persone che apparentemente sono riuscite a farlo funzionare senza problemi.

Ho anche controllato il concetto di stima della densità del kernel (KDE), e anche allora non vedo come usare il comando v.kernel. Quel comando sembra essere un'interpretazione specifica di KDE; i suoi interruttori non sembrano corrispondere bene ai concetti generici di KDE.

Quindi torniamo alla domanda principale qui: come può qualcuno che non è intimo con lo sviluppo del prodotto GRASS usare il comando v.kernel? È disponibile una traduzione in lingua normale?

L' v.kernelalgoritmo calcola la densità dei punti vettoriali per ogni cella di una mappa raster. Se finora hai utilizzato solo oggetti vettoriali, è probabile che tu non abbia impostato regioncorrettamente il tuo (che non solo definisce le estensioni, ma anche la risoluzione raster): probabilmente hai il tuo regionset su una sola riga e una colonna, che significa che l' v.kernelalgoritmo calcolerà la densità del kernel solo come un singolo valore sull'intera mappa. Un'impostazione della regione come questa va bene se non si utilizzano mappe raster, perché le mappe vettoriali non si preoccupano delle impostazioni di risoluzione. Controllare le impostazioni della propria regione utilizzando g.region -pe, se rowse colsè impostato su 1, aumentare la risoluzione utilizzando

g.region rows=<number of rows> cols=<number of columns>

o

g.region res=<length>

dove <length>è la lunghezza di una cella di risoluzione in unità della mappa.