d.legend GLOBCOVER22_CA_lab use=30,40,50,60,70,90,100,110,120,130,140,150,160,170,180,190,200,210,200 -m r.category GLOBCOVER22_CA_lab rules=-<15%) broadleaved evergreen forest (>5m) ? 50:Closed (>40%) broadleaved deciduous forest (>5m) ? 60:Open (15-40%) broadleaved deciduous forest/woodland (>5m) ? 70:Closed (>40%) needleleaved evergreen forest (>5m) ? 90:Open (15-40%) needleleaved deciduous or evergreen forest (>5m) ? 100:Closed to open (>15%) mixed forest (>5m) ? 110:Mosaic forest or shrubland / grassland ? 120:Mosaic grassland / forest or shrubland ? 130:Closed to open (>15%) shrubland (<5m) ? 140:Closed to open (>15%) herbaceous vegetation ? 150:Sparse (<15%) vegetation ? 160:Closed to open (>15%) broadleaved forest regularly flooded ? 170:Closed (>40%) broadleaved forest or shrubland permanently flooded ? 180:Closed to open (>15%) grassland or woody vegetation on regularly flooded ? 190:Artificial surfaces and associated areas (Urban areas >50%) ? 200:Bare areas ? 210:Water bodies ? 220:Permanent snow and ice ? EOF d.vect bstreamsacc50t co=aqua nsfall d.rast srtm90ctr_shade d.vect bstreamsacc50t d.vect srtmstream90med3t co=blue width=2 d.vect nsfpts siz=18 icon=basic/marker co=red fcol=red d.vect pacorabasbigBW type=boundary co=yellow width=2 d.vect chagbasinBW type=boundary co=green width=2 d.vect pacorabasbigBW type=boundary co=yellow nsf_basins only v.overlay ain=bstreamsacc50t atype=line bin=pacorabasbigBW out=pacora_bstreamsacc50 op=and v.overlay ain=srtmstream90med3t atype=line bin=pacorabasbigBW out=pacora_stream90med op=and v.overlay ain=srtmstream90med3t atype=line bin=basin_chagres out=chag_stream90med op=and v.overlay ainput=bstreamsacc50t atype=line binput=basin_chagres out=chag_bstreamsacc50 operator=and --o d.vect chagbasinBW fcol=220:220:220 col=green d.vect pacorabasbigBW fcol=230:230:230 col=yellow d.vect chag_bstreamsacc50 co=aqua d.vect pacora_bstreamsacc50 co=aqua d.vect pacora_stream90med co=blue width=2 d.vect chag_stream90med co=blue width=2 d.vect nsfpts siz=18 icon=basic/marker co=red fcol=red d.rgb b=LSp012r054_20020528z17_30 g=LSp012r054_20020528z17_40 r=LSp012r054_20020528z17_50 awk 'BEGIN {OFS=":"} {print $1,$1,-1*$1,-1*$1}' Transpanam_map2009basrec.txt>Transpanam_map2009basrec1.txt r.watershed elev="srtm90panamalin" threshold=5000 accum="accum200" drain="drain200" bas="basin200 r.watershed elev="panama_srtm3_90m" threshold=5000 accum="accumv3_200" drain="drainiv3_200" bas="basiniv3_200 r.mapcalc "streams_srtm3_200m=if(abs(accumv3_200)>50,1,null())" r.mapcalc "streams_srtm3_200med=if(abs(accumv3_200)>500,1,null())" r.mapcalc "streams_srtm3_200maj=if(abs(accumv3_200)>5000,1,null())" r.watershed el=panama_srtm3_90m th=100000 accum=testaccumsrtm3_90m basin=testbasinsrtm3_90m r.mapcalc "teststreams_srtm3_90=if(abs(testaccumsrtm3_90m)>100,1,null())" r.watershed el=panama_srtm3_90m th=100000 accum=testaccumsrtm3_90ms basin=testbasinsrtm3_90ms r.mapcalc "teststreams_srtm3_90s=if(abs(testaccumsrtm3_90ms)>100,1,null())" r.watershed el=ifs_n10w080d th=100000 accum=testaccumifs_10m basin=testbasinifs_10m r.mapcalc "teststreams_ifs_10m=if(abs(testaccumifs_10m)>500,1,null())" r.thin teststreams_ifs_10m out=teststreams_ifs_10mt r.to.vect teststreams_ifs_10mt out=teststreams_ifs_10mt r.watershed el=ifs_n10w080d th=100000 accum=testaccumifs_30mnn basin=testbasinifs_30mnn r.mapcalc "teststreams_ifs_30mnn=if(abs(testaccumifs_30mnn)>100,1,null())" r.watershed el=panama_srtm3_90m th=100000 accum=testaccumsrtm3_30mnn basin=testbasinsrtm3_30mnn r.mapcalc "teststreams_srtm3_30mnn=if(abs(testaccumsrtm3_30mnn)>100,1,null())" r.terraflow -s ifsall_10m_patch filled=ifsall_10m_fill dir=ifsall_10m_dir swater=ifsall_10m_swatersh accum=ifsall_10m_accum tci=ifsall_10m_tci r.terraflow -s srtm90.8378 filled=srtm90.fill dir=srtm90.dir swater=srtm90.swatersh accum=srtm90.accum tci=srtm90.tci r.terraflow [-sq] elev=string filled=string direction=string swatershed=string accumulation=string tci=string [d8cut=value] [memory=value] [STREAM_DIR=string] [stats=string] r.mapcalc "srtmstream90maj=if(abs(srtm90.accum)>5000 && srtm90.8378 > 0,1,null())" r.mapcalc "srtmstream90med=if(abs(srtm90.accum)>1000 && srtm90.8378 > 0,1,null())" r.mapcalc "srtmstream90med2=if(abs(srtm90.accum)>2000 && srtm90.8378 > 0,1,null())" r.mapcalc "srtmstream90=if(abs(srtm90.accum)>100 && srtm90.8378 > 0,1,null())" r.mapcalc "srtmstream90med3=if(abs(accum90ctr) > 1000 && srtm90ctr > 0, 1, null())" r.mapcalc="if(bbasin.50m30K == 34 || bbasin.50m30K == 36 || bbasin.50m30K == 38 || bbasin.50m30K == 32, 1, null())" r.mapcalc "pacorabasin_90m=if(basinsctr90mask==122 || basinsctr90mask==120 || basinsctr90mask==114 || basinsctr90mask==118, 1, null())" r.mapcalc "pacorabasino_90m=if(basinsctr90mask==122 || basinsctr90mask==120 || basinsctr90mask==114 || basinsctr90mask==118 || basinsctr90mask==116 || basinsctr90mask==110, 1, null())" r.mapcalc "chagresbasin=if(basinsctr90mask==102 || basinsctr90mask==100 || basinsctr90mask==96 || basinsctr90mask==98, 1, null())" r.mapcalc "basinsmask=if(basinsctr90mask==102 || basinsctr90mask==100 || basinsctr90mask==96 || basinsctr90mask==98 ||basinsctr90mask==122 || basinsctr90mask==120 || basinsctr90mask==114 || basinsctr90mask==118 ||basinsctr90mask==116 || basinsctr90mask==110 || basinsctr90mask==92 || basinsctr90mask==34, 1, null())" r.mapcalc "basinsmask_inv=if(isnull(basinsmask), 1,null())" r.mapcalc "pequinibasin90m=if(basinsctr90mask==92, 1, null())" r.mapcalc "cuangobasin90m=if(basinsctr90mask==34, 1, null())" r.mapcalc "basinscanal=if(basin200mask==62 || basin200mask==60 || basin200mask==50 || basin200mask==44 ||basin200mask==46 || basin200mask==48 || basin200mask==52 || basin200mask==56 ||basin200mask==54 || basin200mask==366, 1, null())" r.mapcalc "boundaryctr=if(panama_srtm3_90m,1,null())" r.to.vect boundaryctr out=boundaryctr feature=area --o d.rast -o pacorabasin_90m d.rast -o cuangobasin90m d.rast -o pequinibasin90m d.rast -o chagresbasin d.rast -o pacorabasino_90m d.vect srtmstream90med3t d.vect PanamaMaster07_Feb05_GIS icon=basic/circle size=8 fcolor=black d.rast -o basinsmask_inv d.vect boundaryctr type=boundary width=2 d.barscale -mt d.vect PanamaMaster07_Feb05_GIS icon=basic/circle size=8 d.vect PanamaMaster07_Feb05_GIS icon=basic/marker size=16 width=2 d.vect sstreamctr90t d.vect srtmstream90medt d.vect sstreamctr90majt co=black width=2 d.rast pan_srtm3_300msh d.rast lake1_90m -o d.rast lake2 -o d.vect canal width=3 d.vect basinscanal type=boundary width=2 d.vect basin_cuango type=boundary d.vect basin_pacorabig type=boundary d.vect stream90m_canal d.vect stream90m_cuango d.vect stream90m_pacora v.overlay ain=srtmstream90majt bin=basinscanal out=stream90m_canal atype=line operator=and v.overlay ain=srtmstream90majt bin=basin_cuango out=stream90m_cuango atype=line operator=and v.overlay ain=srtmstream90majt bin=basin_pacorabig out=stream90m_pacora atype=line operator=and v.in.ascii feb05pts.txt r.mapcalc "lsfac=exp(fccum90ctr*90.,1.4)*exp(sin(slopectr_srtm3_90m),1.3)" r.mapcalc "lsfac2=exp(accum90ctr,1.4)*exp(sin(slope.30m),1.3)" r.average base=basinsctr90mask cover=slope30m_percc out=basins_slopeperc r.average base=basinsctr90mask cover=lsfac out=basins_lsfac r.average base=basinsctr90mask cover=lsfac2 out=basins_lsfac2 d.his i=srtm90ctr_shade h=basins_lsfac r.mapcalc "basins5ctrBW=if(isnull(pacorabasbigBW),0,1)+if(isnull(chagbasinBW),0,2) \ + if(isnull(pequinibasBW),0,3)+if(isnull(cuangobasBW),0,4)+if(isnull(diosbasinBW),0,5)" 0 white 1 170 170 170 2 110 110 110 3 210 210 210 4 150 150 150 5 130 130 130 551 20:53 d.rast -o pacorabasbigBW 552 20:53 d.rast -o pequinibasBW 553 20:54 d.rast -o chagbasinBW 554 20:54 d.rast -o cuangobasBW 555 20:55 history 556 20:55 r.colors diosbasinBW co=rules MAPS d.his i=pan_srtm3_300msh h=panama_srtm3_90m bright=20 d.vect Transpanama_map2009_pts attrcol=int_1 disp=attr lcol=blue lsiz=12 d.vect Transpanama_map2009_pts co=red siz=15 icon=basic/marker fcol=red width=2 d.vect srtmstream90majt co=blue d.vect sstreamsouth90maj co=blue d.grid r.colors basinsv3rec2_200m co=random d.rast -o basinsv3rec2_200m d.vect streams_srtm3_200medt d.vect Transpanama_map2009_pts co=black siz=20 icon=basic/marker fcol=red 1464 23:29 v.out.ascii Transpanama_map2009_pts fs=space | r.what basinsv3_200m 1464 23:29 v.out.ascii Transpanama_pub2009xy fs=space | r.what basinsv3_200m 1465 23:30 vi Transpanam_map2009_coordbasins.txt 1466 23:31 cp Transpanam_map2009_coordbasins.txt Transpanam_map2009basrec.txt 1467 23:31 vi Transpanam_map2009basrec.txt 1468 23:38 cp Transpanam_map2009basrec.txt Transpanam_map2009basrec.txtbak 1469 23:43 awk 'BEGIN {OFS=":"} {print $1,$1,-1*$1,-1*$1}' Transpanam_map2009basrec.txt > Transpanam_map2009basrec1.txt 1469 23:43 awk 'BEGIN {OFS=":"} {print $1,$1,-1*$1,-1*$1}' Transpanam_map2009basrec.txt > Transpanam_map2009basrec1.txt awk 'BEGIN {FS="|"; OFS=":"} {print $4,$4,-1*$4,-1*$4}' Transpanama_pub2009_coordbasins.txt >Transpanama_pub2009basrec.txt 1470 23:43 vi Transpanam_map2009basrec1.txt 1471 23:43 r.recode help 1472 23:44 r.recode basinsv3_200m out=basinsv3_200m rules=Transpanam_map2009basrec1.txt 1473 23:44 r.recode basinsv3_200m out=basinsv3rec_200m rules=Transpanam_map2009basrec1.txt 1473 23:44 r.recode basinsv3_200m2 out=basinsv3rec_200m rules=Transpanam_map2009basrec2.txt 1474 23:45 d.rast basinsv3rec_200m 1475 23:47 d.his h=basinsv3rec_200m i=pan_srtm3_300msh 1476 23:48 d.vect Transpanama_map2009_pts co=black siz=15 icon=basic/marker fcol=red r.water.outlet drain=drainiv3w2_90 basin=basinsmwest e=293515 n=913393. r.water.outlet drain=drainiv3n_90 basin=RioNobreDios e=666131 n=1057687 v.out.ascii Transpanama_pub2009xy fs=space | r.what basinsv3_200m vi Transpanama_pub2009_coordbasins.txt awk 'BEGIN {FS="|"; OFS=":"} {print $4,$4,-1*$4,-1*$4}' Transpanama_pub2009_coordbasins.txt > Transpanama_pub2009basrec.txt r.recode basinsv3_200m out=basinsv3recpub_200m rules=Transpanama_pub2009basrec.txt r.recode basinsv3_200m out=basinsv3recpub_200m3 rules=Transpanama_pub2009basrec3.txt r.recode basinsv3_200m out=basinsv3recpub_200m4 rules=Transpanama_pub2009basrec4.txt Data range is -458 to -34 > -458 20 20 20 > -362 80 80 80 > -200 120 120 120 > -100 150 150 150 > -34 200 200 200 > end -458 160 160 160 -440 20 20 20 -406 60 60 60 -394 80 80 80 -362 100 100 100 -243 120 120 120 -238 140 140 140 -196 150 150 150 -62 170 170 170 -34 190 190 190 MAP publ d.rast pan_srtm3_300msh d.rast basinsv3recpub_200m4 -o d.rast -o basins_smallw3 d.rast -o cuangobasin90m d.rast -o RioNobreDios d.vect streams_srtm3_200medt d.vect Transpanama_pub2009xy siz=7 icon=basic/circle fco=black d.out.file transpan_basinsall_temp4.png