Problema de coordenação com interpolação bilinear

O problema é que há um arquivo de texto com três colunas, cada coluna é as coordenadas do ponto no plano (X, Y, Z), respectivamente. As coordenadas dadas com um certo passo (D) em latitude (X) e longitude (Y), ou seja, formam uma grade regular. Em cada ponto (nó da grade) definido por essas coordenadas está a altura, ou seja, a terceira coordenada (Z).
Você precisa buscar a altura do ponto (Z′), que eu insiro pelo teclado atrás das minhas coordenadas (X′,Y′). Se meu ponto não cair no nó da grade regular, então é necessário realizar a interpolação bilinear entre os nós da célula na qual meu ponto cai, ou seja, no quadrado com vértices (X1,Y1; X2,Y1; X2;Y2; X1;Y2).
Por favor, ordene em Delphi ou Python.
E dê acesso ao código em si.