Plot coordinates of a tif image

29 Ansichten (letzte 30 Tage)
Madmad am 21 Feb. 2024
Kommentiert: Cris LaPierre am 9 Mär. 2024
Hello, I have say 'image.tif' which contains a greyscale information and a .json file containing metadata associated.
I can plot the image using:
[A,R] = readgeoraster('image.tif');
However, it does come with X or Y axes in a unit I do not understand, value with x10^6, say 600,000 to 605,000 for Y axis. Same fo X. See figure attached from mapshow help page:
I would like to transform it into proper longitude/latitude, so °N/°E when I plot the image. In the metadata, there are different projection informationn such as projectedCRS, polygon coordinates and projection epsg/shape/transform values.
I tried using the function worldGrid and geoshow functions, but 1st) it takes ages to interpolate lat/lon and 2nd) my computer does not have enough memory to plot the results using geoshow(lat,lon,A). My images are pretty small but high resolution (20,000 x 20,000 pixels).
So my question, how can I convert my Y and X axes into latitude/longitude degrees on the plot? Any function that would do it? In the end I just need to plot a grayscale image with coordinate.
  4 Kommentare
Dyuman Joshi
Dyuman Joshi am 22 Feb. 2024
Verschoben: Dyuman Joshi am 25 Feb. 2024
Could you share the image here? Use the paperclip button to attach.
Madmad am 26 Feb. 2024
I could share but they are quite big (600 MB). I got it from public SAR dataset so you would need to make a free account (

Melden Sie sich an, um zu kommentieren.

Antworten (1)

Cris LaPierre
Cris LaPierre am 26 Feb. 2024
Note that there is a projection associated with your SAR data. You need to know this projection in order to convert your scales to lat/lon. This should be incorporated into your geotiff file.
Follow this example to either display your projected data, or to unproject the data.
  7 Kommentare
Cris LaPierre
Cris LaPierre am 5 Mär. 2024
  1. You whould be able to determine the size of your image using any file browser. In order to load the image into memory, you will need at least that much RAM.
  2. You could use the same commands as in the example I shared above, but then plot using mapshow instead of geoshow. The main difference I discovered is that I had to plot lon as X and lat as y.
p = R.ProjectedCRS;
[x,y] = worldGrid(R);
[lat,lon] = projinv(p,x,y);
xlabel('Longitude (degrees)')
ylabel('Latitude (degrees)')
Cris LaPierre
Cris LaPierre am 9 Mär. 2024
I have more details for you.
mapshow is faster because it assumes your image is in cartesian coordinates. The one we have been using - boston.tif - is, so this works. However, it is more likely than not that you geotiff is not. This is why you should use geoshow instead.
geoshow projects the image using the projection info in the geotiff. This can take some time, especially for very large files.

Melden Sie sich an, um zu kommentieren.


Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by