function geoshow(filename)

latlon = geoloc(filename);
disp(latlon)
rgb = imread(filename);
sz  = max(size(rgb));
N = 600;
if sz>N
    imshow(imresize(rgb,N/sz));
else
    imshow(rgb)
end