How could I get the spatial coordinates for each pixel of a Landsat scene?

Hi, my problem is that I need the map coordinates of each pixel in a LandSat image.

What I have:

[A, R, bbox] = geotiffread(filename);

[X,Y] = pixcenters(R, size(A), 'makegrid');

With this I get a coordinate gird which looks fine to me in UTM coordinate, but with this my object of interest is moving from scene to scene. Which means its coordinates are changing which is definitely not the case (its a big mountain). Do you have any idea where I made the mistake?

Remark: I stacked the images in ENVI and there is no coordinate shift between the images

Thank you for your help!

Answer 1

[Chad Greene]

The easiest way is to use Aslak Grinsted's geoimread. Syntax is simply

[A,x,y] = geoimread('filename');

where A is the image and x,y is pixel coordinates in meters.