Commit eceb8b60 authored by Daniel Bourgault's avatar Daniel Bourgault

Very minor cosmetic changes made.

parent 5813a767
......@@ -38,7 +38,7 @@ for i = 1:n_gcp
% the 'if' condition here.
if h_gcp(i) > 0
[range,a12,a21] = m_idist(LON0,LAT0,lon_gcp0(i),lat_gcp0(i));
proj_factor = H/(H-h_gcp(i));
proj_factor = H/(H - h_gcp(i));
new_range = range*proj_factor;
[lon_gcp(i),lat_gcp(i),a21] = m_fdist(LON0,LAT0,a12,new_range);
if lon_gcp(i) > 180
......@@ -57,7 +57,7 @@ for i = 1:n_gcp
dx = lon_gcp0(i) - LON0;
dy = lat_gcp0(i) - LAT0;
range = sqrt(dx^2 + dy^2);
proj_factor = H/(H-h_gcp(i));
proj_factor = H/(H - h_gcp(i));
new_range = range*proj_factor;
beta = atan2(dy,dx);
......
......@@ -108,7 +108,7 @@ display(' ');
display(' Welcome to g_rect: a package for georectifying oblique images on a flat ocean');
display(' Authors: Daniel Bourgault and Rich Pawlowicz');
display(' ');
inputFname = input(' Enter the filename for the input parameters: ','s');
inputFname = input(' Enter the name of the parameter file: ','s');
fid = fopen(inputFname);
nHeaderLine = 0;
......@@ -157,9 +157,15 @@ lat_gcp0 = gcp.data(:,4);
if n_column == 5
h_gcp = gcp.data(:,5);
else % otherwise elevation are considered to be zero
h_gcp(1:ngcp) = 0.0;
h_gcp(1:ngcp,1) = 0.0;
end
i_gcp = i_gcp';
j_gcp = j_gcp';
lon_gcp0 = lon_gcp0';
lat_gcp0 = lat_gcp0';
h_gcp = h_gcp';
%%
% Check if the elevation of the GCPs are not too high and above
% a certain fraction (gamma) of the camera height. If so, stop.
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment