Thanks Ahmadou for the answer!

I knew it was possible under GRASS but I wondered if it was also the case
under R. I have yet not found any function to accomplish that.

Your example will be of great help, though

2011/9/14 ahmadou dicko <dicko.ahmadou@gmail.com>

> Hi Mathieu
>
> I'm not aware of such package in R, but you can use GRASS trough spgrass6
> from R,  have access to almost all GRASS function.
>
> launch R inside a GRASS session (spearfish60 location for example)
>
> require(spgrass6)
>
>
> execGRASS("g.region", flags = "p", parameters = list(rast =
> "elevation.dem"))
>
> # parseGRASS("r.los")
>
> execGRASS("r.los", parameters = list(input = "elevation.dem", output =
> "los", coordinate = c("598869", "4916642"), obs_elev = 50, max_dist =
> 10000))
>
> los <- readRAST6("los")
>
> bgyr <- colorRampPalette(colors = c("blue", "green", "yellow", "red"),
> space = "rgb")
>
> point <- SpatialPoints(coord = data.frame(x = 598869, y = 4916642))
>
> image(los, col = rev(bgyr(100)))
> plot(point, add = TRUE)
>
> hope this help
>
>
>
>
> 2011/9/14 Mathieu Rajerison <mathieu.rajerison@gmail.com>
>
>> Hi,
>>
>> I'd like to know if an R package exists that provides calculation of a
>> line-of-sight/intervisibility raster given a digital elevation model and
>> on
>> xy(z) observation point.
>>
>> Mathieu
>>
>>        [[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]]

