%Jay Chapman 2016
%Reads in a xxx.tif DEM raster file
%I tend to export small DEM areas as a tif from arcmap

clear all
close all

[filename,pathname] = uigetfile('*.tif','Select DEM file');
full = fullfile(pathname,filename);
AAA = imread(filename);
AAA = double(AAA);  % converts to double precision integer
AAA = flipud(AAA);  % flips to put North up
save('DEM.mat','AAA');
