mapprofile used on a projected map

I am trying to get a data profile along an image band in projection UTM and datum WGS84. I am surprised that mapprofile does not work on a projected image. I'll be grateful to know about alternatives to use. Reprojecting to Geographical coordinates is rather time-consuming. Thanks! BM

Answer 1

[Chad Greene]

If you have a 2D projected map and its corresponding x,y coordinates, you can use interp2 to interpolate along any path you wish. You can create a straight-line path using linspace to make x and y (or lat and lon) arrays from the starting point to the ending point.

If your map is a geotiff, you can use geotiffinterp.

[BM]

Dear Chad, thank you very much for your answer. Actually, I don't think interp2 or geotiffinterp are the commands I'm looking for. Below is an illustration of what I'd like to do. This is straightforward using mapprofile, only if the image is in Lat/Long, but not in UTM. Strange... I guess I'm missing something. BM </matlabcentral/answers/uploaded_files/69411/Capture.PNG>

[Chad Greene]

BM-the image didn't get uploaded.

[BM]

[BM]

I think it is uploaded now. BM

[Chad Greene]

This is like pulling teeth. Exactly what information do you have? Do you have the x and y or lat and lon coordinates of the radiance map? Do you have x and y or lat and lon coordinates of the transect?

[BM]

Hi Chad, The images are in X,Y (UTM projection). The transect is a shape file also in UTM. Everything would be easy if the data was in Lat/Long, but it seems there are not equivalent commands for projected data (?). Strange, specially when huge amounts of Landsat data are distributed in UTM. I now I could "unproject" images to Lat/Long, but that's computationally quite inefficient, and a bit too much for my PC. I'm really looking forward for suggestions/assistance on this. Thanks! BM