hide NaNs in geoshow

20 次查看(过去 30 天)
Sagar
Sagar 2015-11-5
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

回答(3 个)

Amy Haskins
Amy Haskins 2015-11-10
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 个评论
Darcy Cordell
Darcy Cordell 2018-11-24
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 2020-7-27
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 ...

请先登录,再进行评论。


Zelong
Zelong 2023-3-2
编辑:Walter Roberson 2023-9-22
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.

Chad Greene
Chad Greene 2015-11-6
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 个评论
Matthew Cooper
Matthew Cooper 2019-7-9
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 2023-6-29
Great! It worked for pcolorm/surfacem. Thanks a lot!

请先登录,再进行评论。

标签

Community Treasure Hunt

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

Start Hunting!

Translated by