Commit ea630c82 authored by Daniel Bourgault's avatar Daniel Bourgault

Small change in g_proj_GCP to make sure that the longitude go from -180 t0 180.

parent 47ace827
......@@ -41,7 +41,9 @@ for i = 1:n_gcp
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);
lon_gcp(i) = lon_gcp(i) - 360;
if lon_gcp(i) > 180
lon_gcp(i) = lon_gcp(i) - 360;
end
else
% If there's no elevation, then there's nothing to do.
lon_gcp(i) = lon_gcp0(i);
......
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