-
Notifications
You must be signed in to change notification settings - Fork 0
/
extract_values_to_points_grasspy.py
68 lines (45 loc) · 2.69 KB
/
extract_values_to_points_grasspy.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
import grass.script as grass
from grass.script import raster as grassR
import os
import string
import glob
import re
import fnmatch
lista_rasters_maps_old=grass.mlist_grouped ('rast', pattern='*Bin*') ['PERMANENT']
lista_vect_maps_old=grass.mlist_grouped ('vect', pattern='*id_2*') ['PERMANENT']
for i in lista_vect_maps_old:
check="id2_"+i[40:41]
for a in lista_rasters_maps_old:
if "FRAG60m_mata_clump_AreaHA" in a:
if check in a:
grass.run_command('g.region',rast=a)
grass.run_command('v.what.rast', vector=i,raster=a, column="FG060")
if "FRAG120m_mata_clump_AreaHA" in a:
if check in a:
grass.run_command('g.region',rast=a)
grass.run_command('v.what.rast', vector=i,raster=a, column="FG120")
if "patch_clump_mata_limpa_AreaHA" in a:
if check in a:
grass.run_command('g.region',rast=a)
grass.run_command('v.what.rast', vector=i,raster=a, column="PATCH")
lista_rasters_maps_new=grass.mlist_grouped ('rast', pattern='*segment*') ['PERMANENT']
lista_vect_maps_new=grass.mlist_grouped ('vect', pattern='*10SP*') ['PERMANENT']
for i in lista_vect_maps_new:
format=i[36:38]
format=format.replace("_","")
format="00"+format
format=format[-2:]
check="FID_"+format
for a in lista_rasters_maps_new:
if check in a:
if "FRAG60m_mata_clump_AreaHA" in a:
grass.run_command('g.region',rast=a)
grass.run_command('v.what.rast', vector=i,raster=a, column="FG060")
if "FRAG120m_mata_clump_AreaHA" in a:
if check in a:
grass.run_command('g.region',rast=a)
grass.run_command('v.what.rast', vector=i,raster=a, column="FG120")
if "patch_clump_mata_limpa_AreaHA" in a:
if check in a:
grass.run_command('g.region',rast=a)
grass.run_command('v.what.rast', vector=i,raster=a, column="PACTH")