baseline = 150;
model=a(((x/baseline) + 1)((y/baseline) +1)-1)baseline + b x + c y + d;
pracmodel[i_, j_]:=3.139((i/baseline + 1)(j/baseline+1)-1)baseline +  1.453i + 2.52j +.05 + RandomReal[];

rando1 = RandomReal[{-10, 10}, 1000];
rando2 = RandomReal[{-8, 8}, 1000];
data2 = MapThread[List, { rando1, rando2, MapThread[pracmodel,{rando1, rando2}]}];
FindFit[data2, model, {{a, 1}, {b,1}, {c,1}, {d, 0}}, {x,y}]
{a->3.172854854042586`,b->1.4193610694618217`,c->2.4821994272114685`,d->0.5468515898000187`}