[GRASS-user] Reading 'ecw'-files

GRASS GRASS user
,
1

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@venlo
r.patch input=207900_367650_75_2007.red@venlo,\
207900_368100_75_2007.red@venlo,\
207900_368550_75_2007.red@venlo,\
208350_368550_75_2007.red@venlo\
output=gbknH02.red --overwrite
r.patch input=207900_367650_75_2007.green@venlo,\
207900_368100_75_2007.green@venlo,\
207900_368550_75_2007.green@venlo,\
208350_368550_75_2007.green@venlo\
output=gbknH02.green --overwrite
r.patch input=207900_367650_75_2007.blue@venlo,\
207900_368100_75_2007.blue@venlo,\
207900_368550_75_2007.blue@venlo,\
208350_368550_75_2007.blue@venlo\
output=gbknH02.blue --overwrite
g.remove rast=207900_367650_75_2007.red@venlo,\
207900_368100_75_2007.red@venlo,\
207900_368550_75_2007.red@venlo,\
208350_368550_75_2007.red@venlo
g.remove rast=207900_367650_75_2007.green@venlo,\
207900_368100_75_2007.green@venlo,\
207900_368550_75_2007.green@venlo,\
208350_368550_75_2007.green@venlo
g.remove rast=207900_367650_75_2007.blue@venlo,\
207900_368100_75_2007.blue@venlo,\
207900_368550_75_2007.blue@venlo,\
208350_368550_75_2007.blue@venlo
r.composite -d red=gbknH02.red@venlo \
green=gbknH02.green@venlo \
blue=gbknH02.blue@venlo \
levels=32 output=gbknH02 --overwrite
g.remove rast=gbknH02.red@venlo,gbknH02.green@venlo,gbknH02.blue@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

2

On Tue, Oct 14, 2008 at 6:05 PM, Joop Goedbloed <jlgoedbloed@hetnet.nl> wrote:

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@venlo

here you are missing
g.region vect=gbknH02@venlo align=207900_367650_75_2007.red

...

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

You need to set to the resolution of the raster map as well. Otherwise
it takes the previously set resolution.

Markus