[R-sig-Geo] interp.old and write.ascii.grid
Paul Hiemstra
p.hiemstra at geo.uu.nl
Fri Oct 10 10:40:09 CEST 2008
Hi,
What is the exact error you get? Can you provide a example that
reproduces the problem? What kind of system are you running
(sessionInfo() output).
An idea: The documentation of write.asciigrid says that the variable x
(in this case sn8.interpold) should be a SpatialGridDataFrame. Does
interp.old give that output (try class(sn8.interpold))? If this is not
the case you need to convert the output to a SpatialGridDataFrame.
cheers,
Paul
Mikael Carlsson wrote:
> Hello!
>
> I have a text file that looks like
>
> ew;ns;cures
> 1394180;6200295;48
> 1398190;6201020;52
> 1399690;6200420;59
> 1395780;6205105;84
> 1397820;6206590;50
> 1395290;6208305;46
> 1394995;6210370;114
> 1399475;6213715;92
> 1410825;6214870;48
> and so on
> I manage to interpolate with
>
> sn8.interpold <- interp.old(ew, ns, cures,
> xo= seq(min(ew), max(ew), length = 400),
> yo=seq(min(ns), max(ns), length = 300),
> ncp = 4, extrap=FALSE, duplicate = "error", dupfun = NULL)
> then I want to crate a asciigrid...
> write.asciigrid(sn8.interpold,"asciisn8", na.value = -9999)
> and here I got problems...
> Any tips?
>
> Thanks Mikael
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Drs. Paul Hiemstra
Department of Physical Geography
Faculty of Geosciences
University of Utrecht
Heidelberglaan 2
P.O. Box 80.115
3508 TC Utrecht
Phone: +31302535773
Fax: +31302531145
http://intamap.geo.uu.nl/~paul
More information about the R-sig-Geo
mailing list