Hi,

If only one scan line is wrong, there are two easy ways to deal with it,
depending your purpose :
1- if your image is for illustration use only, you can use 'bilinear
interpolation' (maybe
http://www.image.ucar.edu/GSP/Software/Fields/Help/interp.surface.html -
but I never tried it)
2- if you want to classify image, you should set the wrong scan line to NA,
classify image and use a majority filter

You mask should be at least 3 rows height and your {255,0,0} pixels should
be set to NA. Can you send your code ?

Cheers,

Alban



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

> Thank you for the references. But I don't think the problem I have falls
> into the striping concept. The RGB Landsat 5 image has only one scan line
> showing the feature I described. Because the scan line appears as a green
> line in the RGB image, I believe there is something wrong with band 2.
>
>
> 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:* David Stephens (Cefas) <David.Stephens@cefas.co.uk>
> *Para:* alban thomas <alban.thomas.fr@gmail.com>; Alessandro Samuel Rosa <
> alessandrosamuel@yahoo.com.br>
> *Cc:* R-sig-Geo <r-sig-geo@r-project.org>
> *Enviadas:* Quarta-feira, 20 de Março de 2013 11:04
> *Assunto:* RE: [R-sig-Geo] Correcting missing scan lines in satelite
> images
>
> If you are referring to stripping as described below I have had this
> problem using remotely sensed acoustic imagery. The solution I d is to use
> 2-fourier filtering as outlined in this paper:
>
> Wilken, D., Feldens, P., Wunderlich, T. And Heinrich, C. (2012).
> Application of 2D Fourier filtering for elimination of Strip noise in
> side-scan sonar mosaics. Geo-Mar Lett.
>
> I have had some quite good results (
> https://dl.dropbox.com/u/4683316/FF.png) using Fast fourier transform
> (fft) function in R.
>
> Cheers, Dave.
>
>
>
>
> -----Original Message-----
> From: r-sig-geo-bounces@r-project.org [mailto:
> r-sig-geo-bounces@r-project.org] On Behalf Of alban thomas
> Sent: 20 March 2013 10:52
> To: Alessandro Samuel Rosa
> Cc: R-sig-Geo
> Subject: Re: [R-sig-Geo] Correcting missing scan lines in satelite images
>
> 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]]
>
>
>
>
> This email and any attachments are intended for the na...{{dropped:25}}

