Hi Jesse, Edzer, and R-SIG-GEO!

I was able to make Kriging work, following the suggestions of Jesse and
Edzer. However, the results don't turn out to be right - the Kriging and
IDW done on a grid is not the same (in fact, it seems to be a mirror image)
of the block kriging using a separate map. The grid kriging is in fact
closer to the small area estimates (which is the actual sample).

You may take a look at the images here:
https://www.dropbox.com/sh/q5sev7n2pi6tq8f/tiqiFMjJS5

The R Code is here:
https://www.dropbox.com/sh/q5sev7n2pi6tq8f/oM87f6axv5/kriging_pov.r

The R Data is here:
https://www.dropbox.com/sh/q5sev7n2pi6tq8f/8jC9qmlBZB/POVKRIGE.RData

Would you know why block Kriging and block IDW have such a weird results?
Thank you for your time and patience.

James


On Fri, Sep 13, 2013 at 12:26 AM, James Matthew Miraflor <
james.miraflor@gmail.com> wrote:

> Thanks Jesse and Edzer! It finally worked.
>
>
> On Thu, Sep 12, 2013 at 11:44 PM, Jesse Berman <berman.jesse@gmail.com>wrote:
>
>> robin<-"+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84
>> +units=m +no_defs"
>> muni.sp2<-spTransform(muni.sp, CRS(robin))
>> brgy_poly2<-spTransform(brgy_poly, CRS(robin))
>>
>
>
>
>
> --
> *
>
> *
> *James Matthew B. Miraflor*
> MS Computer Science
> College of Engineering
> University of the Philippines
>
> http://politicsforbreakfast.blogspot.com
> ________
>
> "No problem can withstand the assault of sustained thinking." - Voltaire
> ________
>



-- 
*

*
*James Matthew B. Miraflor*
MS Computer Science
College of Engineering
University of the Philippines

http://politicsforbreakfast.blogspot.com
________

"No problem can withstand the assault of sustained thinking." - Voltaire
________

	[[alternative HTML version deleted]]

