[GRASS-user] Reading 'ecw'-files

[Joop Goedbloed]

Hello Usersgroup

Have a mapset of vectorfiles with a square of 1000x1000 m in a location
There is also en set of 'ecw' files on dvd-rom with a square of 450x450m
 to read the ecw-files in grass and mapping them with the vevtos-files 
of 1000x1000m I wrote scripts
to do this job.

Here is an example of the script:

#!/usr/bin/sh
r.in.gdal -o 'input=/media/Luchtfoto'\''s Ven/207900_367650_75_2007.ecw' \
output=207900_367650_75_2007 --overwrite
r.in.gdal -o 'input=/media/Luchtfoto'\''s Ven/207900_368100_75_2007.ecw' \
output=207900_368100_75_2007 --overwrite
r.in.gdal -o 'input=/media/Luchtfoto'\''s Ven/207900_368550_75_2007.ecw' \
output=207900_368550_75_2007 --overwrite
r.in.gdal -o 'input=/media/Luchtfoto'\''s Ven/208350_368550_75_2007.ecw' \
output=208350_368550_75_2007 --overwrite
g.region vect=gbknH02 op venlo
r.patch input=207900_367650_75_2007.red op venlo,\
207900_368100_75_2007.red op venlo,\
207900_368550_75_2007.red op venlo,\
208350_368550_75_2007.red op venlo\
 output=gbknH02.red --overwrite
r.patch input=207900_367650_75_2007.green op venlo,\
207900_368100_75_2007.green op venlo,\
207900_368550_75_2007.green op venlo,\
208350_368550_75_2007.green op venlo\
 output=gbknH02.green --overwrite
r.patch input=207900_367650_75_2007.blue op venlo,\
207900_368100_75_2007.blue op venlo,\
207900_368550_75_2007.blue op venlo,\
208350_368550_75_2007.blue op venlo\
 output=gbknH02.blue --overwrite
g.remove rast=207900_367650_75_2007.red op venlo,\
207900_368100_75_2007.red op venlo,\
207900_368550_75_2007.red op venlo,\
208350_368550_75_2007.red op venlo
g.remove rast=207900_367650_75_2007.green op venlo,\
207900_368100_75_2007.green op venlo,\
207900_368550_75_2007.green op venlo,\
208350_368550_75_2007.green op venlo
g.remove rast=207900_367650_75_2007.blue op venlo,\
207900_368100_75_2007.blue op venlo,\
207900_368550_75_2007.blue op venlo,\
208350_368550_75_2007.blue op venlo
r.composite -d red=gbknH02.red op venlo \
green=gbknH02.green op venlo \
blue=gbknH02.blue op venlo \
levels=32 output=gbknH02 --overwrite
g.remove rast=gbknH02.red op venlo,gbknH02.green op venlo,gbknH02.blue op venlo


this scrip works fine but the resulting map is very poor resollusion
is ther any way making a better quality map (Higher resolution)


thanks

Joop Goedbloed