[R-sig-Geo] problem in krige function

[temiz]

Roger Bivand wrote:
> On Thu, 5 Apr 2007, temiz wrote:
>
>   
>> following the topic GRASS and R kriging interpolation from grass6 tutorial
>>
>> I came here without problems
>>
>> But here:
>>
>> zz <- krige(z ~ 1, locations = el1 , newdata = mask_SG , model = efitted)
>> Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim,  :
>>         var1 : data item in gstat object and newdata have different 
>> coordinate reference systems
>>     
>
> If:
>
> proj4string(el1)
> proj4string(mask_SG)
>
> differ, the function will stop you predicting from data declared to be in 
> one projection to newdata points in a different projection. You can set 
> them by:
>
> proj4string(mask_SG) <- CRS(proj4string(el1))
>
> (untried). If this does not work, post the commands you used to create the 
> objects.
>
> Roger
>
>   
>> how can I solve this problem ?
>>
>> regards
>>
>>
>>     
>
>   
Thank you

I regarded some modifications in 
http://grass.gdf-hannover.de/wiki/GRASS_6_Tutorial#GRASS_and_R_kriging_interpolation

I encountered a new problem :
 > OK_pred <- krige(z ~ 1, el3, newdata=mask_SG, model=efitted)
[using ordinary kriging]

"memory.c", line 57: can't allocate memory in function m_get()
Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim,  :
        m_get



Why does this problem arise ?

I selected a great deal smaller area in R than I did in grass and 
carried out with spline tension without any problem.

is the data too dense or grid size too small ?

regards



-- 
Ahmet Temiz


-- 
This message has been scanned for viruses and\ dangerous con...{{dropped}}