X
PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 11 Jan 2012 08:46 PM by  anon
Glint in WoldView 2 Image
 3 Replies
Sort:
You are not authorized to post a reply.
Author Messages

anon



New Member


Posts:
New Member


--
11 Jan 2012 08:46 PM
    Hello all, I got my hands recently on a WorldView 2 Image, sadly since I'm interested in marine applications (Coral reef/sea grass mapping) the entire image got a sun glint problem... So my qusetion is, is there a way to remove the glint and make the image easier to analysis and classify? BTW, Im using Envi 4.7 Thanks

    Deleted User



    New Member


    Posts:
    New Member


    --
    12 Jan 2012 10:53 AM
    Hi ramisaeed, ENVI doesn't have any tools specifically designed to remove sun glint from imagery. There may be some studies in the research literature that would suggest methods you could use, but I'm afraid that I don't know about them. If you do find a good method, would you share it with us here? Thanks, Peg

    Deleted User



    New Member


    Posts:
    New Member


    --
    17 Jan 2012 09:37 PM
    Here is what i got from lpinner at the ESRI forums when i posted this same question... sadly i didn't know how do any of the ones that he mentioned... can anyone help translate this into a more dumped down version? oh, when i tried using the code he provided in IDL, it gave me an error in line 12, or so the log says... "Have a read of this paper: Hedley et al. (2005) Simple and robust removal of sun glint for mapping shallow-water benthos. International Journal of Remote Sensing Vol. 26, No. 10, 2107–2112. Step by step instructions to implement this method are available for ENVI (not ArcGIS sorry). I use some IDL code from ENVI to do this instead of the above manual method, though I do generate the samples with the ROI tool by hand): nirs=samples[*,*,7] nird=data[*,*,7] d=size(nirs,/dim) i=where(nirs gt 0) o=where(nird eq 0) m=min(nirs[i]) m=nirs[i] m = m[(sort(m))[5 * N_ELEMENTS(m) / 100] ] out=0 help, m foreach b, bindgen(6) do begin band=samples[*,*,b] r=linfit(nirs[i],band[i]) band=data[*,*,b] deglint=long(band - r[1] * (nird - m)+0.5) deglint[o]=0 if size(out, /n_dimensions) eq 0 then $ ;Image stack doesn't exist yet out=temporary(deglint) $ else $ ;Concatenate image to stack out=[[[temporary(out)]],[[temporary(deglint)]]] endforeach out=[[[temporary(out)]],[[data[*,*,6]]]] out=[[[temporary(out)]],[[data[*,*,7]]]] end Again note that these methods are for ENVI, not ArcGIS. You could probably reimplement the code in numpy though."

    Deleted User



    New Member


    Posts:2
    New Member


    --
    29 Oct 2012 07:11 PM
    I've submitted an updated script to the Code Library that runs from the ENVI (Classic) or 4.8 spectral menu - "Sun glint removal tool". I've also replied to your original thread on the ESRI forum and attached the script.
    You are not authorized to post a reply.