On 26/02/15 15:13, Paulo van Breugel wrote:
Is there an (easy) way to upload raster category labels to a point file,
equivalent to uploading category values with v.what.rast?
That would be a nice addition to v.what.rast. Probably worth an
enhancement ticket in trac.
One workaround would be something like this:
g.copy vect=schools_wake,myschools
v.db.addcolumn myschools col="geolcat int"
v.db.addcolumn myschools col="geollabel varchar(4)"
v.what.rast myschools raster=geology_30m column=geolcat
and then a small Python script :
catlabels=grass.read_command('r.category', map='geology_30m').splitlines()
for catlabel in catlabels:
cat=catlabel.split('\t')[0]
label=catlabel.split('\t')[1]
where='geolcat='+cat
grass.run_command('v.db.update', map='myschools', column='geol_code',
value=label, where=where
Moritz
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
