Nikos Alexandris wrote:
[chewed up message by the "System" :D -- apologies for re-posting!]
# 4) i.pansharpen -- all of them, three methods -- from my history
### Brovey
18 31.Jul.13 [ 15:54:36 ] i.pansharpen sharpen=brovey pan=Pan_Radiance
ms1=Blue_Radiance ms2=Green_Radiance ms3=Red_Radiance
output=sharp_brovey_Radiance
19 31.Jul.13 [ 15:54:46 ] i.pansharpen sharpen=brovey pan=Pan_ToAR
ms1=Blue_ToAR
ms2=Green_ToAR ms3=Red_ToAR output=sharp_brovey_ToAR
24 31.Jul.13 [ 15:59:16 ] i.pansharpen sharpen=brovey pan=Pan_DNs
ms1=Blue_DNs
ms2=Green_DNs ms3=Red_DNs output=sharp_brovey_DNs
### PCA
59 31.Jul.13 [ 16:09:53 ] i.pansharpen sharpen=pca pan=Pan_DNs
ms1=Blue_DNs
ms2=Green_DNs ms3=Red_DNs output=sharp_pca_DNs
60 31.Jul.13 [ 16:22:03 ] i.pansharpen sharpen=pca pan=Pan_Radiance
ms1=Blue_Radiance
ms2=Green_Radiance ms3=Red_Radiance output=sharp_pca_Radiance
61 31.Jul.13 [ 16:27:28 ] i.pansharpen sharpen=pca pan=Pan_ToAR
ms1=Blue_ToAR
ms2=Green_ToAR ms3=Red_ToAR output=sharp_pca_ToAR
### IHS
62 31.Jul.13 [ 16:35:12 ] i.pansharpen sharpen=ihs pan=Pan_ToAR
ms1=Blue_ToAR
ms2=Green_ToAR ms3=Red_ToAR output=sharp_ihs_ToAR
63 31.Jul.13 [ 16:37:21 ] i.pansharpen sharpen=ihs pan=Pan_Radiance
ms1=Blue_Radiance
ms2=Green_Radiance ms3=Red_Radiance output=sharp_ihs_Radiance
64 31.Jul.13 [ 16:39:32 ] i.pansharpen sharpen=ihs pan=Pan_DNs
ms1=Blue_DNs
ms2=Green_DNs ms3=Red_DNs output=sharp_ihs_DNs
# 5) Let's look at the numbers
for SHARP in `g.mlist rast pat=sharp*`; do r.info -mr ${SHARP}; done
min=36.2130177514793
max=153
title= (sharp_brovey_DNs_blue)
min=0
max=112.211538461538
title= (sharp_brovey_DNs_green)
min=0
max=107.968559837728
title= (sharp_brovey_DNs_red)
min=7.93864844473524
max=102.4232751775
title= (sharp_brovey_Radiance_blue)
min=0
max=87.5889521289036
title= (sharp_brovey_Radiance_green)
min=0
max=74.435557812593
title= (sharp_brovey_Radiance_red)
min=0
max=0.398724812575096
title= (sharp_brovey_ToAR_blue)
min=0
max=0.335289700449259
title= (sharp_brovey_ToAR_green)
min=0
max=0.476377949923749
title= (sharp_brovey_ToAR_red)
min=0
max=255
title=Image blue (sharp_ihs_DNs_blue)
min=0
max=255
title=Image green (sharp_ihs_DNs_green)
min=0
max=255
title=Image red (sharp_ihs_DNs_red)
min=0
max=255
title=Image blue (sharp_ihs_Radiance_blue)
min=0
max=254
title=Image green (sharp_ihs_Radiance_green)
min=0
max=254
title=Image red (sharp_ihs_Radiance_red)
min=0
max=2
title=Image blue (sharp_ihs_ToAR_blue)
min=0
max=2
title=Image green (sharp_ihs_ToAR_green)
min=0
max=2
title=Image red (sharp_ihs_ToAR_red)
min=177.072681268903
max=490.188522572706
title= (sharp_pca_DNs_blue)
min=174.855236596567
max=415.384343723085
title= (sharp_pca_DNs_green)
min=-71.4712772318356
max=390.53877906274
title= (sharp_pca_DNs_red)
min=38.0505369148222
max=189.658927814102
title= (sharp_pca_Radiance_blue)
min=24.8653625573439
max=215.672403116314
title= (sharp_pca_Radiance_green)
min=-13.7969488642255
max=113.557246501694
title= (sharp_pca_Radiance_red)
min=0.0506311311052859
max=0.660054986167671
title= (sharp_pca_ToAR_blue)
min=0.0146679903931977
max=0.832078269086436
title= (sharp_pca_ToAR_green)
min=-0.0646181045289202
max=0.451095521435642
title= (sharp_pca_ToAR_red)
At least for some of the above results, it is obvious that something isn't
right.
# 6) I use the script copy-pasted in ... the next message... in G64 (...I
wish I could use frames in G7 :-p) to get some overviews, without and with
"i.landsat.rgb -p"
Will upload screenshots uploaded in trac -- if possible.
Summary
In almost all -- actually, in all of them! -- there are NULL cells introduced
there where they did not exist before. I don't use a MASK. I have set up the
region before running i.pansharpen to match the extent of the Red band, for
example (resolution is taken care by the module).
Nikos
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
