How can I plot a latitude / longitude point on a geotiff? (Without mapping toolbox.)

Hi. I have some geotiffs in geographic coordinates (WGS84) from various areas around the globe. I also have a list of points, defined by coordinates, that mark points of interest.

How can I plot the geotiff, then plot all the points within the geotiff domain in their correct position? Ultimately I need the x,y position of each point within the geotiff image.

To make matters more complicated, I don't have the Mapping Toolbox. I am using geoimread from the FEX to read the geotiff image and metadata.

Any ideas are appreciated. Thanks!

Answer 1

[KSSV]

[A,x,y,I]=geoimread('your tif file') ;
pcolor(x,y,A)

x,y, are latitude and longitude vectors, A is your raster data. I is a structure which gives you complete information of the tiff file.

[Jon]

Derp. I didn't notice that geoimread could return the latitude/longitude vectors as well. Thank you.