Hi All,
Is it possible union two code in one (R+SAGA) like this to export from .sgrd
to ascii files? (file: DCM_1_power10.sgrd, DCM_1_power11.sgrd,
DCM_1_power12.sgrd, DCM_1_power40.sgrd and DCM_2_power10.sgrd,
DCM_2_power11.sgrd, DCM_2_power12.sgrd, DCM_2_power40.sgrd)
for(method
On Fri, Oct 24, 2008 at 8:24 AM, Edzer Pebesma
<[EMAIL PROTECTED]> wrote:
> Tomislav Hengl wrote:
>>
>> In GIS terms, getting a raster from one grid topology to (any) other grid
>> topology (also subgrids)
>> is referred to 'resampling' (I do not think that this is yet implemented
>> in any R packa
Tomislav Hengl wrote:
In GIS terms, getting a raster from one grid topology to (any) other grid
topology (also subgrids)
is referred to 'resampling' (I do not think that this is yet implemented in any
R package; a
combination of akima and spTransform will do the trick but it is not really
comp
Hi,
Thank you all very much :-)
I finally managed to put it working. It was a problem in the code, and
not the coordinates.
>> A second doubt I have is if we can set apriori, in R the markers projected
>> in google earth... how do I do if instead of the pins, I want a ball? Do I
>> have to chan
Pieter,
the easiest in terms of the least mental exercise might be to form a
polygon from the four points and then overlay that with the grid. This
will likely be too slow and/or memory intensive.
The more efficient way would be to compute the row and column index
ranges from the bounding box, a
2008/10/24 Marta Rufino <[EMAIL PROTECTED]>:
> I am trying to export a image plot to google earth. I manage to do it with
> points data... ok. but when I do it with the images, they simply go to Fiji
> (which I am sure it is a nice place... but a bit far from the Algarve, where
> it should be).
>
Dear Tomislav and Virgilio,
It is funny,... I wanted to try it, because I was going through your
reports and manuals :-)
I am not sure, but I think the projection and reference system are ok,
because when I create the klm it works fine. Only with the png, does not
function.
summary(grid.bord
Dear Marta,
I am not sure what is the projection system you use for the Algarve case study?
If it is different
from CRS("+proj=longlat +datum=WGS84"), then you do need to first reproject
your maps.
You should really look at the examples from the book by Bivand et al.:
http://www.asdar-book.org
Marta.
I would say that it is a problem with the coordinates of your data. If
your original data are not in longlat you need to re-project them (with
spTransform) before exporting to GE.
You can find an example with the meuse data set on these slides prepared
by Roger Bivand:
http://www.bias-pr
Dear list,
I would like to start by acknowledging the wonderful work of this
support list :-)
I have a question, which should be prety dummy,... but I was looking
around in the net, in reports, etc. and could not sort it out.
I am trying to export a image plot to google earth. I manage to do
to get statistics of grids over polygons, you might also consider using the
Grid Statistics for
Polygons method in SAGA GIS:
> rsaga.get.modules("shapes_grid")
$shapes_grid
code name interactive
1 0Add Grid Values to Points FALSE
2 1 G
If the problem is memory handling in R (still a big headache to do much of GIS
analysis in R), you
should instead consider reading pieces of grids that you need e.g.:
info <- GDALinfo("dem50m.asc")
info
gridmaps01 = readGDAL("dem50m.asc", region.dim =
round(c(info[["rows"]]/2,info[["columns"]])
Yes, you can export sgrd to geoTIFFs via the "io_grid_gdal" module:
> library(RSAGA)
> rsaga.get.modules("io_grid_gdal")
$io_grid_gdal
code name interactive
10Import Raster via GDAL FALSE
21Export Raster via GDAL FALSE
3
Hi,
Which version of R and sp are you using? I'm asking this because I had a
similar problem as you had with a raster of 50.000+ points and some 7000
polygons (not sure if this is comparable to your dataset). Doing the
overlay took far too long, but after a question about this to the
r-sig-ge
14 matches
Mail list logo
