[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