geotiff reader

This code read geotiff file, subset in pixel (or map coordinates) capability

Description

GEOTIFF_READ: read geotiff using imread and assign map info from infinfo.
 
  I = GEOTIFF_READ('filename');
  Reads whole images
  I = GEOTIFF_READ('filename','pixel_subset', [minx maxx miny maxy]);
  I = GEOTIFF_READ('filename','map_subset' , [minx maxx miny maxy]);
  extract subset of the specified.

% output:
 I.z, image data
 I.x, x coordinate in map
 I.y, y coordinate in map
 I.info, misc. info

 imshow(I.z, 'xdata', I.x, 'ydata', I.y);
 shows image with map coordinate

% Version by Yushin Ahn
% Glacier Dynamics Laboratory,
% Byrd Polar Resear Center, Ohio State University
% Referenced enviread.m (Ian Howat)

Acknowledgements

The author wishes to acknowledge the following in the creation of this submission:
ENVI file reader, updated 2/9/2010

• MATLAB release: MATLAB 7.9 (2009b)

[Leslie]

I have problem to display the matrix I.z with the image function from matlab. (I do not have imshow). The matrix I.z is ordered differently from the one obtained with A=imread(filename) and displayed with image(A).
I tried to change the way the file is read without not much success. Any suggestion

Leslie