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 [EMAIL PROTECTED]
r.patch [EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED]
output=gbknH02.red --overwrite
r.patch [EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED]
output=gbknH02.green --overwrite
r.patch [EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED]
output=gbknH02.blue --overwrite
g.remove [EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED]
g.remove [EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED]
g.remove [EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED],\
[EMAIL PROTECTED]
r.composite -d [EMAIL PROTECTED] \
[EMAIL PROTECTED] \
[EMAIL PROTECTED] \
levels=32 output=gbknH02 --overwrite
g.remove [EMAIL PROTECTED],[EMAIL PROTECTED],[EMAIL PROTECTED]
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
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
