Hello everybody, I am working with a geo-corrected hyperspectral image composed by 224 bands. In its header the UTM location (latitude / longitude) of the upper-left corner (long0,lat0), the pixel size and the rotation appear. By clicking in Tools / Cursor Location Value it is shown all the information about the pixel the mouse cursor is over, including its value and its latitude and longitude. How is the UTM calculated for each pixel?? It would be very easy because: latitude=lat0+num_row*pixel_size_y longitude=long0+num_column*pixel_size_x But this is not exact because of the rotation. How is the rotation considered when calculating the UTM coordinates of each pixel? My final aim is to get a matrix with the same spatial resolution than the image, but with two bands. The first band would store the longitude of each pixel, and the second one would store the latitude. Is there any easier way to do this than using the expressions above (including rotation)? Does any IDL routine that given one pixel calculates de UTM coordinates exist? I hope you can help me, if you need some more information I can provide it. Thank you very much!!