On 09/01/12 19:16, Moritz Lennert wrote:
On 09/01/12 00:50, Markus Neteler wrote:
On Fri, Jan 6, 2012 at 5:42 PM, Moritz Lennert
<[email protected]> wrote:
Hi everyone,
I have some troubles with image fusion through IHS transformation of a
Landsat image.
Hi Moritz,
could you please post the procedure you used?
The same procedure as the one you describe in the book, i.e. the
classical i.rgb.his & i.his.rgb approach with grey.eq color table for
the resulting three bands.
I'll try to write up a reproducible example, including relevant data if
the NC set does not show the same behaviour.
Ok, here we go using nc_spm_08:
- Using "real-color" channels I get weird colors
g.region rast=lsat7_1002_20@PERMANENT
r.colors lsat7_1002_30@PERMANENT color=grey.eq
r.colors lsat7_1002_20@PERMANENT color=grey.eq
r.colors lsat7_1002_10@PERMANENT color=grey.eq
g.region -p rast=lsat7_1002_80@PERMANENT
i.rgb.his red_input=lsat7_1002_30@PERMANENT
green_input=lsat7_1002_20@PERMANENT blue_input=lsat7_1002_10@PERMANENT
hue_output=h intensity_output=i saturation_output=s --o
i.his.rgb hue_input=h@fusion_test
intensity_input=lsat7_1002_80@PERMANENT saturation_input=s@fusion_test
red_output=fr green_output=fg blue_output=fb --o
r.colors map=fr@fusion_test color=grey.eq
r.colors map=fg@fusion_test color=grey.eq
r.colors map=fb@fusion_test color=grey.eq
Color composite before fusion (d.rgb -o red=lsat7_2002_30@PERMANENT
green=lsat7_2002_20@PERMANENT blue=lsat7_2002_10@PERMANENT):
http://geog-pc40.ulb.ac.be/grass/fusion_test/d_rgb_30_20_10_grey_eq.png
Color composite after fusion (d.rgb -o red=fr@fusion_test
green=fg@fusion_test blue=fb@fusion_test):
http://geog-pc40.ulb.ac.be/grass/fusion_test/d_rgb_fr_fg_fb_grey_eq.png
- Using near-infrared false color composition the colors are as expected:
g.region rast=lsat7_2002_30@PERMANENT
r.colors lsat7_2002_40@PERMANENT color=grey.eq
r.colors lsat7_2002_30@PERMANENT color=grey.eq
r.colors lsat7_2002_20@PERMANENT color=grey.eq
g.region -p rast=lsat7_2002_80@PERMANENT
i.rgb.his red_input=lsat7_2002_40@PERMANENT
green_input=lsat7_2002_30@PERMANENT blue_input=lsat7_2002_20@PERMANENT
hue_output=h intensity_output=i saturation_output=s --o
i.his.rgb hue_input=h@fusion_test
intensity_input=lsat7_2002_80@PERMANENT saturation_input=s@fusion_test
red_output=pir_fr green_output=pir_fg blue_output=pir_fb --o
r.colors map=pir_fr@fusion_test color=grey.eq
r.colors map=pir_fg@fusion_test color=grey.eq
r.colors map=pir_fb@fusion_test color=grey.eq
Color composite before fusion (d.rgb -o red=lsat7_2002_40@PERMANENT
green=lsat7_2002_30@PERMANENT blue=lsat7_2002_20@PERMANENT):
http://geog-pc40.ulb.ac.be/grass/fusion_test/d_rgb_40_30_20_grey_eq.png
Color composite after fusion (d.rgb -o red=pir_fr@fusion_test
green=pir_fg@fusion_test blue=pir_fb@fusion_test):
http://geog-pc40.ulb.ac.be/grass/fusion_test/d_rgb_pir_fr_fg_fb_grey_eq.png
Can anyone reproduce this ? Any explanation ?
Moritz
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
