Več

Kako uporabiti GRASS v.kernel za izračun gostote točk vzdolž omrežja?

Kako uporabiti GRASS v.kernel za izračun gostote točk vzdolž omrežja?


V cestnem omrežju moram izbrati odseke, kjer na razdalji 100 m se pojavi najmanj 5 pik (1 pika = 1 naslov) (priložena slika ponazarja moje vprašanje). Poskušam uporabiti GRASS v.kernel za to, z a enakomerna porazdelitev vendar ne vem, kako ga konfigurirati. Z uporabo teh parametrov:

v.kernel -v vhod = MY_LOCATIONS net = MY_NETWORK output = MY_OUTPUT stddeviation = 100 distmax = 1 mult = 200 vozlišče = split jedro = uniform

Posledično dobim segmente z vrednostmi, kjer se pojavi več pik, vendar ne vem, kako jih razlagati.

Na priloženi sliki: segmenti, kot jih označuje zeleni krog, so tisto, kar iščem.


Morda sem našel rešitev za svoje vprašanje, vendar nisem prepričan:

v.kernel -v vhod = MY_LOCATIONS neto = MY_NETWORK izhod = MY_OUTPUT stddeviation = 57.735 distmax = 1 več = 200 jedro = enotno

Krepki parametri in vrednosti so pomembni. Od kod prihajajo? Jedro, ki obdaja vsako piko, ima enakomerno porazdelitev z razponom 100 m na obeh straneh lokacije pik. Zato je „doseg“ (b - a) enotne porazdelitve 200 m. Standardni odmik enotne porazdelitve 1 / (b - a) je (b - a) / sqrt (12), zato je parameter standardnega odklona 200 / sqrt (12) = 57,735 (referenca)

Ker je višina enotne porazdelitve 1/200, rezultat pomnožim z 200 (torej parameter mult = 200).

Sčasoma dobim odseke vrstic z vrednostmi, ki označujejo število pik na dolžini 100 m.

Zdi se, da je rezultat pravilen, vendar nisem prepričan.

Uredi:

No, nisem prepričan, da je moja uporaba modula navsezadnje pravilna. Na splošno se ocena gostote jedra poveča na segmentih, kjer se pojavi več pik, kot bi moralo biti. Vendar pa obstajajo nekateri protislovni rezultati, ki jih ne morem razložiti. Spodaj je posnetek zaslona, ​​ki prikazuje omrežje (tanka črtkana črta), lokacije v omrežju (rdeče pike) in rezultat ocene gostote jedra, ki ima vrednost 5 ali več (debele vijolične črte). Ne razumem rezultata za dve regiji, označeni z oranžno. Manj kot 3 pike na razponu 100 m nikoli ne smejo imeti jedra več kot 3 (vendar je v mojem rezultatu 5). Tudi segment z okoli 8 pikami nad nekaj več kot 100 m bi moral imeti oceno gostote jedra vsaj 5, vendar je v mojem izhodu manj (torej brez debele vijolične črte). Mi lahko kdo to razloži? Ali pravilno uporabljam modul GRASS?

Hvala.

Uredi (20.01.2018):

Modul v.kernel se je spremenil v GRASS 7.2.2. Tu je posodobitev, kako bi zdaj uporabljal modul:

SET XMIN = 227511 SET XMAX = 240244 SET YMIN = 202932 SET YMAX = 213221 g. Regijska res = 1 n =% YMAX% s =% YMIN% e =% XMAX% w =% XMIN% v.in.ogr dsn =. /network.shp output = V_NETWORK --overwrite --quiet v.in.ogr dsn =./events.shp output = V_EVENTS -overwrite --quiet SET input = V_EVENTS SET net = V_NETWORK SET net_output = V_KERNEL_OUTPUT SET polmer = 100 SET dsize = 0 SET segmax = 10 SET distmax = 1 SET multiplikator = 1000 SET vozlišče = razdeljeno jedro SET = enakomerno v.kernel -overwrite --verbose input =% input% net =% net% net_output =% net_output% radius = % polmer% dsize =% dsize% segmax =% segmax% distmax =% distmax% multiplikator =% množitelj% vozlišče =% vozlišče% jedro =% jedro% v.out.ogr -sce vnos = V_KERNEL_OUTPUT dsn =. olayer = v_kernel format = ESRI_Shapefile type = line --overwrite --quiet g.remove -f type = ime vektorja = V_KERNEL_OUTPUT --quiet

Poglej si posnetek: SANS Institute