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.
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)
Those "Tinfo." should be changed to "Tinfo(1)." to handle images with multiple bands.
Hi, - I get:
Index exceeds matrix dimensions.
Error in geotiffread (line 120)
tmp(i-sub(3)+1,:,j) = t(R);
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
Inspired by: ENVI file reader, updated 2/9/2010