MATLAB: How to plot a latitude / longitude point on a geotiff? (Without mapping toolbox.)

coordinate conversioncoordinate-transformationgeotiffimage processinglatitudelongitudeMapping Toolboxprojections

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!

Best Answer

[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.

Related Solutions

MATLAB: How to get the points of [lat, lon] of a full great circle

I would use a combination of azimuth, gc2sc, and scircle1:

az = azimuth(lat1, lon1, lat2, lon2)
[lat, lon, r] = gc2sc(lat1, lon1, az)
[latc, lonc] = scircle1(lat, lon, r)

MATLAB: Writing multiple Lat, lon points to KML File !

The syntax for kmlwrite is:

kmlwriteline(filename,latitude,longitude,altitude);

The fourth variable passed is considered as altitude data. In your case you can create an array with multiple lat and lon points and send that array to the function.

Try this:

lat_points = [lat1;lat2;lat3;lat4];
lon_points = [lon1;lon2;lon3;lon4];
kmlwriteline(filename,lat_points,lon_points);

Hope this helps.

Best Answer

Related Solutions

MATLAB: How to get the points of [lat, lon] of a full great circle

MATLAB: Writing multiple Lat, lon points to KML File !

Related Question