Hi,

I think you are talking about striping, which occurs quite often in remote
sensing images : http://landsat.usgs.gov/science_an_detectorstriping.php
It occurs especially in Landsat 7 images since its Scan Line Corrector
(slc) is off (it also occurs with other sensors).

I've failed to find an easy way to destripe images using open source
software like R, Grass and Sextante. You may have a look on other softwares
like : ossim (http://www.dpi.inpe.br/spring/english/), saga gis.

You also may find tips on this web page : http://landsat.usgs.gov/sci_an.php

HTH,

Alban


2013/3/19 Alessandro Samuel Rosa <alessandrosamuel@yahoo.com.br>

> Hi Robert,
> I made a wrong description of my problem. My Landsat images have saturated
> scan lines. When I make the RGB composition using bands 3, 2 and 1, I see
> strong green lines. There also are clusters of pixels with DN = 255 or DN =
> 0. I know these are spurious values because I have field experience in the
> study area. There are no NAs. Any idea?
>
>
>
> Alessandro Samuel-Rosa
> Postgraduate Course in Agronomy - Soil Science
> Federal Rural University of Rio de Janeiro
> Seropédica, Rio de Janeiro, Brazil
> soil-scientist.net
>
>
> ________________________________
>  De: Robert J. Hijmans <r.hijmans@gmail.com>
>
> Cc: R-sig-Geo <r-sig-geo@r-project.org>
> Enviadas: Sábado, 16 de Março de 2013 9:52
> Assunto: Re: [R-sig-Geo] Correcting missing scan lines in satelite images
>
>
> Alessandro,
>
> It is probably easier to do
>
> d <- focal(x, ..., NAonly=TRUE)
>
> # or
>
>
> a <- focal(x, ...)
> b <- cover(x, a)
>
> Robert
>
> On Fri, Mar 15, 2013 at 6:59 AM, Alessandro Samuel Rosa <
> alessandrosamuel@yahoo.com.br> wrote:
>
> Hi all,
> >
> >I am working with a few Landsat images that has missing scan lines. I
> want to apply a low-pass filter to the image to correct these missing scan
> lines. However, I do not want to apply the filter to the entire image.
> Thus, I'll draw polygons (shp or kml files) to delineate the areas were the
> filter shall be applied. I can set these polygons as masks using the
> function 'mask' of the 'raster' package. I know that the low-pass filter is
> applied using the function 'focal' of package 'raster'. However, I don't
> know how to apply the low-pass filter using the mask. I hope someone can
> give me some advice, please.
> >
> >Cheers,
> >
> >
> >Alessandro Samuel-Rosa
> >Postgraduate Course in Agronomy - Soil Science
> >Federal Rural University of Rio de Janeiro
> >Seropédica, Rio de Janeiro, Brazil
> >soil-scientist.net
> >
> >        [[alternative HTML version deleted]]
> >
> >
> >_______________________________________________
> >R-sig-Geo mailing list
> >R-sig-Geo@r-project.org
> >https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >
> >
>         [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>


-- 
____________________________
       Alban Thomas
Géomatique & Environnement
http://www.alban-thomas.azerttyu.net/

	[[alternative HTML version deleted]]

