Elias,
parece simples, mas não posso utilizar o centróide do pixel, pois haverá descontinuidade nas minhas amostras, conceitualmente está errado.
As quadrículas (pixels) adjacentes são contínuas, mas entre um centróide e outro há descontinuidade, conceitualmente está inconsistente.
Como mencionei, existe uma solução num software comercial FME (única que encontrei até então, nem mesmo no ArcGIS encontrei tal função) e o resultado da krigagem foi coerente, realmente estimei quadrículas (áreas) não amostradas.
Att,
Alex Santos