How can I remove the region outside of the vector data when I combine the vector and raster data on the same map?

Hi All,

I just plot a map combine the vector data and raster data with the following code,

load mapdemo.mat
ax = worldmap(latlim, lonlim);
surfm(Lat_lon,Lon_lat,China_data);
symbols=makesymbolspec('Polygon',{'default','FaceColor','none', 'LineStyle','--','LineWidth',0.2, 'EdgeColor',[0 0 0]});
geoshow('bou2_4p.shp', 'FaceColor','w','FaceAlpha',0,'EdgeColor','black','LineWidth', 1)
colorbar

However, I want to obtain a map like this,

I am wondering how can I remove the region outside of the vector data in the first figure to get the second figure. Could you help me? If Matlab mapping toolbox is not propoer to realize such function, is there other software is capable of this function?

Thank you so much.

Yours,

Ricky

Answer 1

[KSSV]

You must be having the boundary coordinates in hand........use inpolygon and get the points lying inside the closed region. Read about inpolygon

[Rui Li (AIRicky)]

Thank you so much. I've solved my problem with inpolygon and boundary coordinates.