No subject

[Joyce Anne Nagle]

| We have a LandsatTM subset (900x900 pixels), bands 2-5 only, 
| as 4 separate files, one for each band (ie BSQ).  Can anyone 
| offer advice how to import this into GRASS, as *imagery*.
| 
| Note - we do not have a 9-track (1/2") tape drive. 
| 
| __________________________________________________________________
|                               Dr Simon Cox
|        __  L
|     ,~'  L_|\                 Department of Earth Sciences       
|  ,-'         \                Monash University    
|  (            \               Clayton  Vic  3168  Australia
|  \    ___     /
|   L,~'   "\_x/                Phone +61 (3) 565 5762
|             u                 Fax   +61 (3) 565 5062
|                               simon at cerberus.earth.monash.edu.au
| __________________________________________________________________
| 
|     Bouquets gratefully accepted, brickbats returned to sender
| __________________________________________________________________


This is something I struggle with and then Michael Shapiro helped me out.

Put all 4 bands in the cell directory.  Name each as you wish, e.g. tm.2,
tm.3, tm.4, tm.5.   Then run r.support for each file.   The header file
looks like this:
-----------------------------------------------------------------
proj:       0
zone:       0
north:      -0.5
south:      -900.5
east:       900.5
west:       0.5
cols:       900
rows:       900
e-w resol:  1
n-s resol:  1
format:     0
compressed: 0
------------------------------------------------------------------
You can either use r.support to create the header, in which case you have
rows = 900, cols = 900 and bytes/cell = 1, or just copy this header
into files called tm.2, tm.3, tm.4, tm.5 in the cellhd directory.

Run r.support to create all the support files and then run r.colors tm.#
color=grey.eq or i.grey.scale to create a histogram stretched color table.


It worked great for me.  Hope it does the same for you.


~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Joyce A. Nagle
U.S. Army Cold Regions Research and Engineering Laboratory (CRREL)
Geophysical Sciences Branch          email:  nagle at hanover-crrel.army.mil
72 Lyme Road                         phone:  (603) 646-4161 / 4100
Hanover, NH  03755-1290              fax:    (603) 646-4278
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~