153GIS2 - 10. cvičení - GRASS GIS: Porovnání verzí
(prazdna stranka zalozena) |
mBez shrnutí editace |
||
Řádek 2: | Řádek 2: | ||
Stránka obsahuje ukázky skriptů v [[Python|Pythonu]] pro [[GRASS GIS]] z [[153GIS2 2009 - 10. cvičení|10. cvičení GIS2]]. | Stránka obsahuje ukázky skriptů v [[Python|Pythonu]] pro [[GRASS GIS]] z [[153GIS2 2009 - 10. cvičení|10. cvičení GIS2]]. | ||
<source lang="python"> | |||
import arcgisscripting, sys #import modulu | |||
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: | |||
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 lang="python"> | |||
import arcgisscripting, sys #import modulu | |||
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: | |||
input = "Liberec" #String - nazev obce | |||
rows = gp.SearchCursor("obce.shp") #vytvori cursor | |||
row = rows.next() #promena row obsahuje jeden zaznam ze shapefilu | |||
while row: | |||
if row.GetValue("NAZEV") == input: | |||
pocet = row.GetValue("OB01") #pocet obyvatel zadane obce | |||
row = rows.next() | |||
radky = gp.SearchCursor("obce.shp") | |||
rad = radky.next() | |||
while rad: | |||
if rad.GetValue("OB01") > pocet: | |||
print rad.GetValue("NAZEV"), rad.GetValue("OB01") | |||
rad = radky.next() | |||
except: | |||
print gp.GetMessages() | |||
</source> | |||
<source lang="python"> | |||
import arcgisscripting, sys #import modulu | |||
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: | |||
input = sys.argv[1] #String - nazev obce | |||
output = sys.argv[2] #cesta k textovemu souboru | |||
f = open(output, 'w') #otevru file | |||
rows = gp.SearchCursor("obce.shp") #vytvori cursor | |||
row = rows.next() #promena row obsahuje jeden zaznam ze shapefilu | |||
while row: | |||
if row.GetValue("NAZEV") == input: | |||
pocet = row.GetValue("OB01") #pocet obyvatel zadane obce | |||
row = rows.next() | |||
radky = gp.SearchCursor("obce.shp") | |||
rad = radky.next() | |||
while rad: | |||
if rad.GetValue("OB01") > pocet: | |||
f.write(rad.GetValue("NAZEV")+" "+str(rad.GetValue("OB01"))+"\n") #do souboru pisu obce ktere maji vic obyvatel nez zadana | |||
rad = radky.next() | |||
f.close() | |||
except: | |||
print gp.GetMessages() | |||
</source> | |||
{{GIS}} | {{GIS}} | ||
{{GFOSS}} | {{GFOSS}} |
Verze z 6. 5. 2009, 09:33
Stránka obsahuje ukázky skriptů v Pythonu pro GRASS GIS z 10. cvičení GIS2.
import arcgisscripting, sys #import modulu
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:
data = gp.Describe("obce.shp")
fields = data.Fields
field = fields.next()
while field:
print field.Name
field = fields.next()
except:
print gp.GetMessages()
import arcgisscripting, sys #import modulu
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:
input = "Liberec" #String - nazev obce
rows = gp.SearchCursor("obce.shp") #vytvori cursor
row = rows.next() #promena row obsahuje jeden zaznam ze shapefilu
while row:
if row.GetValue("NAZEV") == input:
pocet = row.GetValue("OB01") #pocet obyvatel zadane obce
row = rows.next()
radky = gp.SearchCursor("obce.shp")
rad = radky.next()
while rad:
if rad.GetValue("OB01") > pocet:
print rad.GetValue("NAZEV"), rad.GetValue("OB01")
rad = radky.next()
except:
print gp.GetMessages()
import arcgisscripting, sys #import modulu
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:
input = sys.argv[1] #String - nazev obce
output = sys.argv[2] #cesta k textovemu souboru
f = open(output, 'w') #otevru file
rows = gp.SearchCursor("obce.shp") #vytvori cursor
row = rows.next() #promena row obsahuje jeden zaznam ze shapefilu
while row:
if row.GetValue("NAZEV") == input:
pocet = row.GetValue("OB01") #pocet obyvatel zadane obce
row = rows.next()
radky = gp.SearchCursor("obce.shp")
rad = radky.next()
while rad:
if rad.GetValue("OB01") > pocet:
f.write(rad.GetValue("NAZEV")+" "+str(rad.GetValue("OB01"))+"\n") #do souboru pisu obce ktere maji vic obyvatel nez zadana
rad = radky.next()
f.close()
except:
print gp.GetMessages()