Просьба гуру в работе с библиотекой PROJ4 пояснить в чём я ошибаюсь.
Есть две проекции:
МСК-70 зона 3
и+proj=tmerc +ellps=krass +towgs84=24,-123,-94,-0.02,0.25,0.13,0.13,1.1 +units=m +lon_0=80.73333333333 +lat_0=0 +k_0=1 +x_0=3250000 +y_0 =-5912900.566
МСК-70 зона 4
в Delphi пытаюсь перевести координаты из одной зоны в другу, вот что делаю:+proj=tmerc +ellps=krass +towgs84=24,-123,-94,-0.02,0.25,0.13,0.13,1.1 +units=m +lon_0=83.73333333333 +lat_0=0 +k_0=1 +x_0=4250000 +y_0 =-5912900.566
На входе даётся точкаprocedure TForm3.Button1Click(Sender: TObject);
var P1, P2: Pointer; i: Integer;
X, Y: PDoubleArray;
begin
P1 := pj_init_plus('+proj=tmerc +ellps=krass +towgs84=24,-123,-94,-0.02,0.25,0.13,0.13,1.1 +units=m +lon_0=83.73333333333 +lat_0=0 +k_0=1 +x_0=4250000 +y_0=-5912900.566');
P2 := pj_init_plus('+proj=tmerc +ellps=krass +towgs84=24,-123,-94,-0.02,0.25,0.13,0.13,1.1 +units=m +lon_0=80.73333333333 +lat_0=0 +k_0=1 +x_0=3250000 +y_0=-5912900.566');
New(X); New(Y);
X^[0] := StrToFloat(Edit1.Text);
Y^[0] := StrToFloat(Edit2.Text);
i := pj_transform(P1, P2, 1, 1, X, Y, nil);
Label1.Caption := String(pj_strerrno(i));
Edit3.Text := FloatToStr(X^[0]);
Edit4.Text := FloatToStr(Y^[0]);
end;
на выходеX = 355256.85
Y = 4333607.02
Вопрос, что я делаю неправильно?X = 3250000
Y = 4089236.93154376