Thank you very much Edzer.

I need to project "rich".

Best,

Manuel

On 08/11/2010 12:55 a.m., Edzer Pebesma wrote:
>
> On 11/07/2010 11:58 PM, Manuel Spínola wrote:
>> Dear list members,
>>
>> I got this error when running function krige:
>>
>> rich.vgm = variogram(N ~ 1, data = rich)
>>
>> rich.fit = fit.variogram(rich.vgm, model = vgm(8.8, "Ste", 0.003))
>>
>>   >  rich.kriged = krige(N~1, rich, as(r, "SpatialGridDataFrame"), model =
>> rich.fit)
>> Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim,  :
>>     var1 : data item in gstat object and newdata have different
>> coordinate reference systems
>>
> First thing I would do after seeing this error message is to compare
> their respective projections:
>
> proj4string(rich)
> proj4string(as(r, "SpatialGridDataFrame"))
>
> as krige apparently would like these two to be identical.


-- 
*Manuel Spínola, Ph.D.*
Instituto Internacional en Conservación y Manejo de Vida Silvestre
Universidad Nacional
Apartado 1350-3000
Heredia
COSTA RICA
mspinola@una.ac.cr
mspinola10@gmail.com
Teléfono: (506) 2277-3598
Fax: (506) 2237-7036
Personal website: Lobito de río 
<https://sites.google.com/site/lobitoderio/>
Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/>

	[[alternative HTML version deleted]]

