hide NaNs in geoshow

35 views (last 30 days)
Sagar
Sagar on 5 Nov 2015
Commented: Walter Roberson on 22 Sep 2023
I am plotting some data using geoshow as below: geoshow (yy, xx, data, 'displaytype', 'texturemap'); The 'data' has a lot of NaNs in it and I do not want to show them (or display them as white) in the resulting image, neither in the colorbar (colorbar should not show the white color). The range of data is between 0 and 1.5. I tried using alphadata as below: t = geoshow (yy, xx, data, 'displaytype', 'texturemap'); set(t, 'AlphaData', ~isnan(data)); but I get the error 'value must be numeric'. I looked in the forums and there are not really clean solution to this problem. Could you please provide a nice, clean solution?
Thank you,
-Sagar

Answers (3)

Amy Haskins
Amy Haskins on 10 Nov 2015
Try using 'DisplayType', 'surface' instead of 'texturemap'. Also try looking at the fourth example on the geoshow doc page "Create Map of Korea and Display NaNs as Transparent".
  4 Comments
Darcy Cordell
Darcy Cordell on 24 Nov 2018
This does not work if you are plotting other types of data on the map.
For example, if you are plotting points using
p = geopoint(x,y);
geoshow(p)
The points will not show up because they will be "below" the surface or something. They are still there, but hiding. Not sure if there is a nice work-around...
Chien-Han Su
Chien-Han Su on 27 Jul 2020
I agree with @Darcy Cordell. My application needs me to show my both data and coastlines, but using 'surface' as the 'DisplayType' makes the coastline covered by the data in some portion and therefore imcomplete.
Currently, I fixed this through customizing the colormap (adding one additional color to represent the NaN value). But still, I wonder if there is a better solution ...

Sign in to comment.


Zelong
Zelong on 2 Mar 2023
Edited: Walter Roberson on 22 Sep 2023
I also encountered this problem today and found the answer on a web page :( https://gis.stackexchange.com/questions/192577/set-nan-values-to-transparent-in-matlab-geoshow-function )
The sample code is as follows.
[lat,lon]= meshgrid(0:10);% sample lat/lon
data = peaks(11); % sample data
data(:,2:3) = NaN; % some NaN values
geoimg = geoshow(lat,lon,data,'DisplayType', 'texturemap') % save object
geoimg =
Surface with properties: EdgeColor: [0 0 0] LineStyle: 'none' FaceColor: 'texturemap' FaceLighting: 'flat' FaceAlpha: 1 XData: [11×11 double] YData: [11×11 double] ZData: [11×11 double] CData: [11×11 double] Show all properties
geoimg.AlphaDataMapping = 'none'; % interpet alpha values as transparency values
geoimg.FaceAlpha = 'texturemap'; % Indicate that the transparency can be different each pixel
alpha(geoimg,double(~isnan(data))) % Change transparency to matrix where if data==NaN --> transparency = 1, else 0.
  1 Comment
Walter Roberson
Walter Roberson on 22 Sep 2023
@Yan Tong comments
This answer has solved the question perferctly!

Sign in to comment.


Chad Greene
Chad Greene on 6 Nov 2015
Hey Sagar,
Yeah, that set AlphaData trick works for imagesc but not for geoshow. You can try to work around it by converting the logical to numeric like this:
set(t, 'AlphaData',double(~isnan(data)));
But I doubt it will work the way you want it to. I tend to avoid using geoshow because I'm never quite sure how it will behave. It decides which plotting function to call based on what you enter, so why not call the plotting function directly? Can you use pcolor, pcolorm or imagesc?
  3 Comments
Matthew Cooper
Matthew Cooper on 9 Jul 2019
For this to work you need to also set 'FaceAlpha' to 'texturemap' which tells geoshow to map the AlphaData values to the FaceAlpha values.
set(t,'FaceAlpha','texturemap','AlphaData',double(~isnan(data)));
SUBHASMITA DASH
SUBHASMITA DASH on 29 Jun 2023
Great! It worked for pcolorm/surfacem. Thanks a lot!

Sign in to comment.

Tags

Products

Community Treasure Hunt

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

Start Hunting!