MATLAB: Mapshow of Geotiff Landsat8 band

geotifflandsatMapping Toolboxmapshow

I read a Geotiff Landsat8 single band and try to plot it on a map axes. I got a number of errors, (eg. due to the 'postings' tag). I have tried dozens of possibilities and couldn't display it. I am currently a bit confused by so many "referencing objects, rasters, etc.". I would be very grateful if someone can let me know how. Please, try it with Landsat8 bands, UTM WGS84.

Thanks, B

Best Answer

If you're using the same projection as the Landsat image itself (which is probably the local UTM zone), you can use L8read and simply plot in x,y coordinates. It would look something like this:

[I,x,y] = L8read('LC80270392013243LGN00','band',3); 
imagesc(x,y,I); 
axis xy image

Related Solutions

MATLAB: Must I change the geoid using UTM projection

In order to work with Google Earth, you should specify the WGS 84 ellipsoid in your code by replacing lines 4 and 5 with:

    utmstruct.geoid = almanac('earth','wgs84','meters');

if you are working in meters, or:

    utmstruct.geoid = almanac('earth','wgs84','kilometers');

if you are working in kilometers, etc.

UTM was designed in the first half of the 20th century, well before artificial satellites made a global system like WGS 84 possible. So different local datums (and their ellipsoids) were used for UTM in different parts of the world. The purpose of the utmgeoid function is to do a geographic look-up of the standard ellipsoid for a given zone. But now UTM is used with datums like NAD 83 and WGS 84, which did not exist when it was originally defined. The utmgeoid function is not relevant in such cases.

The statement in the solution, "MATLAB’s default reference ellipsoid is International, not WGS 84" is misleading. It really should say something like: "defaultm uses the International ellipsoid for UTM when the ellipsoid and zone unspecified." (The International ellipsoid is the default only in this context.)

MATLAB: How can we find the non common strings from multiple columns

First lets find the strings that are common to all cell arrays:

listA={'India','Iran','Ireland','Israel','Jordan','Kuwait','Lebanon','Lesotho','Libya','Luxembourg','Malawi','Maldives','Morocco','Nepal','Netherlands','Norway','Oman','Paraguay','Peru','Philippines',};
listB={'Denmark','El Salvador','Gabon','Gambia','Lebanon','Lesotho','Libya','Luxembourg','Malawi','Maldives','Morocco','Nepal','Netherlands','Norway','Oman','Paraguay','Peru','Philippines','Qatar','Rwanda','Saudi Arabia','Senegal','Singapore','Sudan'};
listC={'Kuwait','Lebanon','Lesotho','Nepal','Netherlands','Oman','Sweden','Syria','Tanzania','Turkey','Uruguay'};
%

lists = {listA,listB,listC}; % storing data together makes working with them much simpler.
%
exc = lists{1};
for k = 2:numel(lists)
    exc = intersect(exc,lists{k});
end

giving the strings found in all cell arrays as:

 exc =
    'Lebanon'
    'Lesotho'
    'Nepal'
    'Netherlands'
    'Oman'

Now we can easily get the input cell arrays with the strings from exc removed:

out = cell(size(lists));
for k = 1:numel(lists)
    out{k} = setdiff(lists{k},exc);
end

giving:

>> out{:}
ans = 
    'India' 'Iran' 'Ireland' 'Israel' 'Jordan' 'Kuwait' 'Libya' 'Luxembourg' 'Malawi' 'Maldives' 'Morocco' 'Norway' 'Paraguay' 'Peru' 'Philippines'
ans =
    'Denmark' 'El Salvador' 'Gabon' 'Gambia' 'Libya' 'Luxembourg' 'Malawi' 'Maldives' 'Morocco' 'Norway' 'Paraguay' 'Peru' 'Philippines' 'Qatar' 'Rwanda' 'Saudi Arabia' 'Senegal' 'Singapore' 'Sudan'
ans =
    'Kuwait' 'Sweden' 'Syria' 'Tanzania' 'Turkey' 'Uruguay'

And if you want a list of all strings together, excluding those in exc:

uni = {};
for k = 1:numel(lists)
    uni = union(uni,lists{k});
end
uni = setdiff(uni,exc)

gives:

uni =
    'Denmark' 'El Salvador' 'Gabon' 'Gambia' 'India' 'Iran' 'Ireland' 'Israel' 'Jordan' 'Kuwait' 'Libya' 'Luxembourg' 'Malawi' 'Maldives' 'Morocco' 'Norway' 'Paraguay' 'Peru' 'Philippines' 'Qatar' 'Rwanda' 'Saudi Arabia' 'Senegal' 'Singapore' 'Sudan' 'Sweden' 'Syria' 'Tanzania' 'Turkey' 'Uruguay'

Best Answer

Related Solutions

MATLAB: Must I change the geoid using UTM projection

MATLAB: How can we find the non common strings from multiple columns

Related Question