[R-sig-Geo] Still error exporting ascii grid generated by DTM grid, please any suggestions?
Paul Hiemstra
p.hiemstra at geo.uu.nl
Mon Nov 30 09:53:31 CET 2009
Hi,
Did you try to export the data using writeGDAL (rgdal-package)? Does
this generate the same problem?
cheers,
Paul
Tobin Cara wrote:
> Hello,
>
> I am still receiving a bizarre error:
>
> Error in write.asciigrid(data_int_ked[1], name1) :
> Asciigrid does not support grids with non-square cells
>
> I do not understand this as I am putting interpolated data back into the cells of an ascii DTM grid with square cells. Information on the DTM I am kriging to is below (I have combined two datasets into elev, COS7 and the Z elevation):
>
> Thank you in advance for your help,
>
> Cara
>
>
>> summary(elev)
>>
> Object of class SpatialGridDataFrame
> Coordinates:
> min max
> X 550119.63 674619.6
> Y 79073.03 167073.0
> Is projected: NA
> proj4string : [NA]
> Number of points: 2
> Grid attributes:
> cellcentre.offset cellsize cells.dim
> X 550369.63 500 249
> Y 79323.03 500 176
> Data attributes:
> Z COS7
> Min. : 380 Min. : -1.510
> 1st Qu.: 1500 1st Qu.: 0.180
> Median : 2143 Median : 0.630
> Mean : 2077 Mean : 2.083
> 3rd Qu.: 2663 3rd Qu.: 4.190
> Max. : 4471 Max. : 8.730
> NA's :22337 NA's :22337.000
> _______________________________________________
> 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: +3130 274 3113 Mon-Tue
Phone: +3130 253 5773 Wed-Fri
http://intamap.geo.uu.nl/~paul
More information about the R-sig-Geo
mailing list