diff --git a/matlab/GetJK.m b/matlab/GetJK.m index fc91ae959189b78dcddfe8c3bc49b32d6f716dfb..e05008449d8b6e66d3fba0b40a67661c0c69783f 100644 --- a/matlab/GetJK.m +++ b/matlab/GetJK.m @@ -1,6 +1,6 @@ function [J,K] = GetJK(Lat,Lon) -global mstruct LOv LAv jg kg +global LOv LAv jg kg X1 = -2600000; Y1 = 0; @@ -8,11 +8,6 @@ X2 = 2300000; Y2 = 3200000; nx = 4901; ny = 3201; -if length(which('projfwd')) > 0 - [X,Y] = projfwd(mstruct,double(Lat),double(Lon)); - J = (ny - 1)*(Y - Y1)/(Y2 - Y1) + 1; - K = (nx - 1)*(X - X1)/(X2 - X1) + 1; -else - J = interp2(LOv,LAv,jg,Lon,Lat,'linear'); - K = interp2(LOv,LAv,kg,Lon,Lat,'linear'); -end + +J = interp2(LOv,LAv,jg,Lon,Lat,'linear'); +K = interp2(LOv,LAv,kg,Lon,Lat,'linear');