Following a shell script written by Glenn Lawson, I have been able to tint
fields on orthographic image to represent crops or HEL. Now I am trying to
tint soils for capability class. This involves reclassing, resampling,
rescaling and then patching the component capability classes. When I try
this, the script seems to run correctly but I end up with the wrong colors.
Could this be caused by the categories not being continuous? My code
follows. I am running grass4.1 on a sparc 2 with sunos 4.1.3. The script
is in csh.
## lets see if we can tint soil classes
#use number of soil capability classes
set numcats=8
echo "The tint map, $tint, has $numcats soil categories."
#######do fields- loop through soil categories
set catno = 1
while ($catno <= $numcats)
r.reclass in=$tint out=MASK$catno << EOF
$catno = $catno
EOF
g.copy rast=MASK$catno,MASK >/dev/null
echo "resampling capability class $catno ....."
r.resample in=$ortho out=T_ORTHO$catno
g.remove rast=MASK,MASK$catno >/dev/null
#target range
if ($catno == 1) then
set rangeBegin = 1
else @ rangeBegin = $catno * 40 - 39
endif
@ rangeEnd = $catno * 40
r.rescale in=T_ORTHO$catno out=T_ORTHO$catno.1 to=$rangeBegin,$rangeEnd
@ catno += 1
end
#end loop through fields
####patching maps
set ct = 1
set pmaps = T_ORTHO1.1
@ catno = $numcats
while($ct < $catno)
@ ct += 1
echo "ct = $ct "
set pmaps=$pmaps,T_ORTHO$ct.1
end
echo "patching maps together......"
r.patch in=$pmaps out=soilsNE
set ct = 1
while ($ct <= $numcats + 1)
g.remove rast=T_ORTHO$ct,T_ORTHO$ct.1 >/dev/null
@ ct += 1
end
### rescale background and patch to tinted fields
if (-e $LOCATION/cats/rescaleOrtho) then
### DO NOTHING ###
else
r.rescale in=bridgene out=rescaleOrtho to=361,400
endif
r.patch in=soilsNE,rescaleOrtho, out=nresultSoilsNE
endif
#set proper colors for soil capability classes
cat ${colorrules} | r.colors map=nresultSoilsNE co=rules
#display tinted soils/home/grass4.1data
d.rast nresultSoilsNE
d.vect 16.nw color=black
d.vect streams.16 color=blue
d.vect soils.16 color=white
~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*
Rebecca Scheckler rebecca@vt.edu(internet)
Information Systems and Insect Studies rebeccas@vtvm1(bitnet)
Virginia Tech fax: 703-231-9131
202 Price Hall phone: 703-231-7216
Blacksburg, VA 24061-0319
~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*