Re: [R-sig-Geo] using IDW with anisotropy

Hi, You could also opt for kriging, which supports anisotropy. The intamap package has a function to detect anisotropy, in case you would like to automate the interpolation. cheers, Paul On 03/18/2011 01:41 AM, smur...@worldbank.org wrote:

[R-sig-Geo] LMC manually fitted

Dear friends, I am doing cokriging and i could not fit the LMC to my variograms automatically, therefore i ftted the LMC manually. i have a total of 15 variograms (direct variograms and cross variograms). eg: clay.g$model$clay[1,2] clay.g$model$clay[1,2] = 0.22 clay.g$model$clay[2,2] = 0.15

[R-sig-Geo] extract cell numbers from raster along a line

I'd like to extract celle numbers from a raster using a line, but it seems the cellnumbers=TRUE option isn't working with lines in raster::extract. simple example : library(raster) r - raster(matrix(1:20, nrow=4)) l - SpatialLines(list(Lines(list(Line(list(x=c(0, 1), y=c(0, 1, ID=1))) #

[R-sig-Geo] Map digitization and classification

I have a series of scanned global maps (from a bound Atlas) of oceanographic sampling effort that I would like to classify. On 1 x1 lat/lon grids are symbols that represent sampling density. I need to read in these scanned maps, and classify the symbols (squares are classified as 1, triangles as

Re: [R-sig-Geo] extract cell numbers from raster along a line

I'd like to extract celle numbers from a raster using a line, but it seems the cellnumbers=TRUE option isn't working with lines in raster::extract. Etienne, It is not working because cellnumbers is currently only an argument for extracting values from a Raster* for polygons. That's what the

[R-sig-Geo] na.omit-method for STFDF

Hi, I wrote a na.omit-method for STFDF (attached) with the intention to extract the rows and columns of a sparse space-time grid (STSDF) that have data on all nodes via na.omit(as(stsdf, STFDF)) It turned out that I had no point in time with data at all locations in my data set. Hence, for

Re: [R-sig-Geo] Map digitization and classification

Here is an example of one of these map sets. This is the original as it was sent to me, but I would crop each map individually. https://docs.google.com/viewer?a=vpid=explorerchrome=truesrcid=0B0d3zfSSPFQsY2MxODEyZWEtZTRkZC00OTk2LTgwY2YtYTZkYzcwZGYxZDllhl=enauthkey=CLmWvWc

Re: [R-sig-Geo] Map digitization and classification

Hi Tim, wow that looks it could be rather difficult to automate. It's probably easiest just to visualize them in a map plot and then use locator() to recreate the locations for each symbol. That would not be too difficult, but there are a few options. Is the data not also published in the atlas

Re: [R-sig-Geo] Map digitization and classification

I had another look and the georegistration should be pretty accurate since there are so many grid lines. I missed that the first time. If the images need significant warping to get them regular again you could use control points with the GDAL command line tools, which is probably easier than

Re: [R-sig-Geo] extract cell numbers from raster along a line

It is not working because cellnumbers is currently only an argument for extracting values from a Raster* for polygons. That's what the docs say. But I'll also implement if for lines. I did notice that, but there was no ==line== section in the doc I have. Is my package outdated ? In your