153GIS2 - 10. cvičení - GRASS GIS: Porovnání verzí

Z GeoWikiCZ
VizuálníWikitext
mBez shrnutí editace
 
(Není zobrazeno 12 mezilehlých verzí od stejného uživatele.)
Řádek 1: Řádek 1:
{{Upravit}}
Stránka obsahuje ukázky skriptů v [[Python|Pythonu]] pro [[GRASS GIS]] z [[153GIS2 2009 - 10. cvičení|10. cvičení GIS2]]. Viz také [[Skriptování v GRASS GIS]].
 
Stránka obsahuje ukázky skriptů v [[Python|Pythonu]] pro [[GRASS GIS]] z [[153GIS2 2009 - 10. cvičení|10. cvičení GIS2]].


=== 1 ===
=== 1 ===


<source lang="python">
<source lang="python">
import arcgisscripting, sys    #import modulu
from grass.script import vector as grass
gp = arcgisscripting.create()  #vytvoreni objektu geoprocessingu
 
gp.Overwriteoutput = 1          #nastaveno prepisovani
for nazev in grass.vector_columns(map = 'obce').keys():
gp.workspace = "d:/jirka/skola/vyuka/gis/model/"    #nastaveni workspace - pracovni adresar nebo geodatabaze
     print nazev
try:
    data = gp.Describe("obce.shp")  
    fields = data.Fields
    field = fields.next()
    while field:
      print field.Name
      field = fields.next()
except:
     print gp.GetMessages()
</source>
</source>


Řádek 25: Řádek 13:


<source lang="python">
<source lang="python">
import arcgisscripting, sys    #import modulu
import grass.script as grass
gp = arcgisscripting.create()  #vytvoreni objektu geoprocessingu
gp.Overwriteoutput = 1          #nastaveno prepisovani
gp.workspace = "d:/jirka/skola/vyuka/gis/model/"    #nastaveni workspace - pracovni adresar nebo geodatabaze
   
   
try:
select = grass.vector_db_select(map = "obce",
  input = "Liberec"    #String - nazev obce
                                where = "NAZEV = 'Liberec'")
  rows = gp.SearchCursor("obce.shp")  #vytvori cursor
 
  row = rows.next()      #promena row obsahuje jeden zaznam ze shapefilu
idx = select['columns'].index('OB01')
  while row:
ob01 = 0
      if row.GetValue("NAZEV") == input:
for record in select['values'].itervalues():
          pocet = row.GetValue("OB01")     #pocet obyvatel zadane obce
    ob01 += int(record[idx])
      row = rows.next()
 
  radky = gp.SearchCursor("obce.shp")
nazev_idx = select['columns'].index('NAZEV')
  rad = radky.next()
for record in grass.vector_db_select(map = "obce",
  while rad:
                                    where = "OB01 > %d" % ob01)['values'].itervalues():
      if rad.GetValue("OB01") > pocet:
    print unicode(record[nazev_idx], 'latin2').ljust(15), record[idx]
          print rad.GetValue("NAZEV"), rad.GetValue("OB01")
      rad = radky.next()
except:
  print gp.GetMessages()  
</source>
</source>


Řádek 51: Řádek 32:


<source lang="python">
<source lang="python">
import arcgisscripting, sys     #import modulu
import sys
gp = arcgisscripting.create()  #vytvoreni objektu geoprocessingu
import grass.script as grass
gp.Overwriteoutput = 1          #nastaveno prepisovani
 
gp.workspace = "d:/jirka/skola/vyuka/gis/model/"    #nastaveni workspace - pracovni adresar nebo geodatabaze
input = sys.argv[1]    # String - nazev obce
try:
output = sys.argv[2]    # cesta k textovemu souboru
  input = sys.argv[1]    #String - nazev obce
 
  output = sys.argv[2]    #cesta k textovemu souboru
select = grass.vector_db_select(map = "obce",
  f = open(output, 'w')   #otevru file
                                where = "NAZEV = 'Liberec'")
  rows = gp.SearchCursor("obce.shp")   #vytvori cursor
 
  row = rows.next()       #promena row obsahuje jeden zaznam ze shapefilu
idx = select['columns'].index('OB01')
  while row:
ob01 = 0
      if row.GetValue("NAZEV") == input:
for record in select['values'].itervalues():
          pocet = row.GetValue("OB01")     #pocet obyvatel zadane obce
    ob01 += int(record[idx])
      row = rows.next()
 
          radky = gp.SearchCursor("obce.shp")
nazev_idx = select['columns'].index('NAZEV')
  rad = radky.next()
f = open(output, 'w')
  while rad:
try:
      if rad.GetValue("OB01") > pocet:
    for record in grass.vector_db_select(map = "obce",
          f.write(rad.GetValue("NAZEV")+" "+str(rad.GetValue("OB01"))+"\n")    #do souboru pisu obce ktere maji vic obyvatel nez zadana
                                        where = "OB01 > %d" % ob01)['values'].itervalues():
      rad = radky.next()
        f.write("%s%s\n" % (record[nazev_idx].ljust(15), record[idx]))
  f.close()
finally:
except:
    f.close()
  print gp.GetMessages()  
</source>
 
=== 4 ===
 
<source lang="python">
import grass.script as grass
 
for record in grass.vector_db_select(map = 'kraje',
                                    columns = "NAZEV_ENG")['values'].itervalues():
    val = unicode(record[0], 'ascii')
    print val
    grass.run_command('v.extract',
                      quiet = True,
                      overwrite = True,
                      input = 'kraje',
                      output = 'kraje_' + val,
                      where = "NAZEV_ENG = '%s'" % val)
</source>
</source>


{{GIS}}
{{GIS}}
{{GFOSS}}
{{GFOSS}}
{{GRASS}}
{{Python}}

Aktuální verze z 25. 2. 2010, 13:36

Stránka obsahuje ukázky skriptů v Pythonu pro GRASS GIS z 10. cvičení GIS2. Viz také Skriptování v GRASS GIS.

1

from grass.script import vector as grass

for nazev in grass.vector_columns(map = 'obce').keys():
    print nazev

2

import grass.script as grass
 
select = grass.vector_db_select(map = "obce",
                                where = "NAZEV = 'Liberec'")

idx = select['columns'].index('OB01')
ob01 = 0
for record in select['values'].itervalues():
    ob01 += int(record[idx])

nazev_idx = select['columns'].index('NAZEV')
for record in grass.vector_db_select(map = "obce",
                                     where = "OB01 > %d" % ob01)['values'].itervalues():
    print unicode(record[nazev_idx], 'latin2').ljust(15), record[idx]

3

import sys
import grass.script as grass

input = sys.argv[1]     # String - nazev obce
output = sys.argv[2]    # cesta k textovemu souboru

select = grass.vector_db_select(map = "obce",
                                where = "NAZEV = 'Liberec'")

idx = select['columns'].index('OB01')
ob01 = 0
for record in select['values'].itervalues():
    ob01 += int(record[idx])

nazev_idx = select['columns'].index('NAZEV')
f = open(output, 'w')
try:
    for record in grass.vector_db_select(map = "obce",
                                         where = "OB01 > %d" % ob01)['values'].itervalues():
        f.write("%s%s\n" % (record[nazev_idx].ljust(15), record[idx]))
finally:
    f.close()

4

import grass.script as grass

for record in grass.vector_db_select(map = 'kraje',
                                     columns = "NAZEV_ENG")['values'].itervalues():
    val = unicode(record[0], 'ascii')
    print val
    grass.run_command('v.extract',
                      quiet = True,
                      overwrite = True,
                      input = 'kraje',
                      output = 'kraje_' + val,
                      where = "NAZEV_ENG = '%s'" % val)
Citováno z „https://geo.fsv.cvut.cz/index.php?title=153GIS2_-_10._cvičení_-_GRASS_GIS&oldid=16151“