Re: [FreeGeoCZ] skript - vyhladanie - region

[Martin Landa <landa.martin AT gmail.com>]

Dobry den,

2009/5/4 Dusan Senko <dusan AT xlan.sk>:
> obraciam sa na Vas s otazkou. Mame vektor s atrib. tabulkou, kde v
> jednom stlpci je nazov geomorf. celku. Zaujimalo by nas, ake prikazy
> (len strucne) by mal obsahovat skript, kt by vyhladal urcity nazov
> (napr. Male Karpaty) a podla neho urcil g.region? Dalej pokracuje

v GRASS >= 6.4.0 pouzijte prepinac -r modulu v.db.select

v.db.select map=geology where="GEO_NAME='Qp'" -r

v nizsich verzich

v.extract
v.info -g

> vypocet, kt uz mame skriptovo zastreseny. Kriterium na vyber bude bud 1)
> postupne kazdy celok v tejto vektorovej vrstve alebo 2) vyber bude
> cerpat z textaku, kde bude pod sebou uvedene nazvy.

priklad skriptu v Pythonu (pouze GRASS >= 6.4.0)

#!/usr/bin/env python

import sys

import grass

def main():
    map = 'geology'
    where = "GEO_NAME = 'Qp'"

    ret = grass.read_command('v.db.select',
                             flags = 'r',
                             map = map,
                             where = where)

    if not ret:
        grass.fatal('Dotaz skoncil chybou')

    reg = grass.parse_key_val(ret)
    grass.run_command('g.region',
                      n = reg['n'],
                      s = reg['s'],
                      e = reg['e'],
                      w = reg['w'])

    return 0

if __name__ == "__main__":
    sys.exit(main())

Martin

-- 
Martin Landa <landa.martin gmail.com> * http://gama.fsv.cvut.cz/~landa