[GRASS-user] Pan-sharpening problem

Hamish hamish_b at yahoo.com
Tue Feb 5 21:53:49 EST 2008


Reno Bladergroen  wrote:
> Since my last post I learned a lot in the few past days, but keep
> stumbling on problems :-/. According to an excellent site 
> http://www.oryx.com/ams/landsat.html (tip from hamish_b... thanks!) i
> followed their instructions almost exactly and tried to pan-sharpen
> an area.
> I called my imported output files simply "Tm-2/4/5 and pan-80". Now
> when I start the i.fusion.brovey command the following error pops up:
> Calculating brov.red: ...(step 1/4)
> Invalid map <Tm>
> Parse error
> 
> An error occurred while running r.mapcalc
> 
> Then I thought: lets try to run r.mapcalc on the first three band
> files with the formula A+B+C. That worked fine!
> What am I doing wrong? All files are in the same area (and accepted
> that way by grass) and also imported in the same location.

r.mapcalc (GRASS's powerful map calculator module) is interpreting the
"-" in your map name as a minus sign and getting confused. In 6.3svn I
have now quoted the map name variables in the script to avoid those
problems, but it needs testing.

If you are using 6.3 you can download the new version and give it a
try.

http://trac.osgeo.org/grass/browser/grass/trunk/scripts/i.fusion.brovey

If you are using 6.2 for now you can just g.rename the maps with "-"
replaced with a less problematic character like "_" or "."s.


Hamish



      ____________________________________________________________________________________
Never miss a thing.  Make Yahoo your home page. 
http://www.yahoo.com/r/hs



More information about the grass-user mailing list