I am working with a satellite image in simple binary, pixels are stored as 2-byte integers. I've used the code below to build a matrix. I am able to display the satellite's "native" view. What I need to do is map the X,Y to a mercator projection. I have the upperleft and lowerright lat/lon coordinates, so I am thinking this could probably be done, but I am not sure how or what functions I should be looking at. Any suggestions?
Thanks.
img = multibandread(file, [3712 3712 1], 'int16', 0, 'bil', 'ieee-le'); m = img(ROWa:ROWb,COLa:COLb); [X,Y] = ind2sub(size(m), 1:numel(m)); Z = m(:);
Best Answer