error message rgdal: dimensions do not match
On Tue, 4 Jul 2006 karl.sommer at dpi.vic.gov.au wrote:
I am encountering a error message when attempting to do "inverse distance kriging" or variogram calculations. The error is always the same ERROR: dimensions do not match: locations 1513600 and data 611452.
Indeed strange, given that 880 by 860 is 756800 (though 756800*2 is the first number). What does summary(t1) say? Is the single column named band1 or ndvi? The class of t1 will be SpatialGridDataFrame, wouldn't SpatialPointsDataFrame be more digestible? Does doing it in two steps (gstat() then predict() of the gstat object) show anything? What did traceback() say after krige()? Roger
What puzzles me is the fact that I am using the same code as I have successfully used with 2 other images. I have tried to reduce the dimensions of the imported image in using gdal_translate from the command prompt. The same error was still present but with different "address-numbers" relating to locations and data. I have tried to find some pointers re error messages under the GDAL manual but without success. I am able to plot the imported image as well as the target grid and SPDF without error. following is the code: library(rgdal) # import image
t1 <- readGDAL("images/test1.bil", half.cell=c(0.5, 0.5), silent = FALSE)
images/test1.bil has GDAL driver GTiff and has 860 rows and 880 columns Closing GDAL dataset handle 0100ACB8... destroyed ... done. Closing GDAL dataset handle 00000000... done. # get shape file of field
image(t1)
poly <- locator()
poly <- do.call("cbind", poly)
poly <- rbind(poly, poly[1,])
lines(poly)
field <- SpatialPolygons(list(Polygons(list(Polygon(poly)), ID="field")))
grd <- GridTopology(c(599040, 6228570), c(.5,.5), c(880, 860))
SG <- SpatialGrid(grd)
x1 <- overlay(SG, field)
SGDF <- SpatialGridDataFrame(grd, data=AttributeList(list(field=x1)))
SPDF <- as(SGDF, "SpatialPixelsDataFrame")
# do kriging
library(gstat) t1.idw <- krige(ndvi~1, t1, SPDF)
Error: dimensions do not match: locations 1513600 and data 611452 Regards Karl
_________________________________ Karl J Sommer, Department of Primary Industries, Catchment & Agriculture Services, PO Box 905 Mildura, VIC, Australia 3502 Tel: +61 (0)3 5051 4390 Fax +61 (0)3 5051 4534 Email: karl.sommer at dpi.vic.gov.au _______________________________________________ R-sig-Geo mailing list R-sig-Geo at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Roger Bivand Economic Geography Section, Department of Economics, Norwegian School of Economics and Business Administration, Helleveien 30, N-5045 Bergen, Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43 e-mail: Roger.Bivand at nhh.no