'Re: [R-sig-Geo] Correcting missing scan lines in satelite images'

[prev in list] [next in list] [prev in thread] [next in thread] 

List:       r-sig-geo
Subject:    Re: [R-sig-Geo] Correcting missing scan lines in satelite images
From:       Alessandro Samuel Rosa <alessandrosamuel () yahoo ! com ! br>
Date:       2013-03-25 12:34:52
Message-ID: 1364214892.94112.YahooMailNeo () web141405 ! mail ! bf1 ! yahoo ! com
[Download RAW message or body]

Hi, I try the second solution and send the code when I have it done.

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

________________________________
 De: alban thomas <alban.thomas.fr@gmail.com>

Cc: David Stephens (Cefas) <David.Stephens@cefas.co.uk>; R-sig-Geo \
                <r-sig-geo@r-project.org> 
Enviadas: Segunda-feira, 25 de Março de 2013 8:11
Assunto: Re: [R-sig-Geo] Correcting missing scan lines in satelite images

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

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 <al
> 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
> 
> 

> 
> > 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 < 

> > 
> > 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 named recipient only. Its \
> unauthorised use, distribution, disclosure, storage or copying is not permitted. If \
> you have received it in error, please destroy all copies and notify the sender. In \
> messages of a non-business nature, the views and opinions expressed are the \
> author's own and do not necessarily reflect those of Cefas. Communications on \
> Cefas' computer systems may be monitored and/or recorded to secure the effective \
> operation of the system and for other lawful purposes. 
> 
> 
> 

-- 
____________________________
           Alban Thomas
Géomatique & Environnement
http://www.alban-thomas.azerttyu.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

[prev in list] [next in list] [prev in thread] [next in thread] 