adding boundary to shapefile

Sir, I plot an Image from a nc file by using imagesc(lon,lat,a);I want to add my area shapefile.I import s=shaperead(shapefile);mapshow(s);But it is showing colour insiade the shapefile.I want only boundary.How can I do this.

[Tanmoyee Bhattacharya]

when me doing plot([s.X],[s.Y], 'k:');then map is showing.And when a=imagesc(lon,lat,a);then image is coming.Not falling each other.

[Tanmoyee Bhattacharya]

[Tanmoyee Bhattacharya]

Answer 1

[KSSV]

S = shaperead('file') is a structure. Use plot(S(i).X,S(i).Y) ,,i will be from 1:length(S). Varies according to your state. I suggest you to use pcolor or surf rather imagesc for plotting.

Answer 2

[KSSV]

clc; clear all ;
file = 'cpc_rfe_v2.0_sa_dly.bin.20020101.nc' ;
lon = ncread(file,'lon') ;
lat = ncread(file,'lat') ;
RF = ncread(file,'RF') ;
%
path = [pwd filesep 'India State wise with telangana' filesep 'india27-11.shp'] ;
S = shaperead(path) ;
n = length(S) ;
x = zeros(n,1) ;  y =zeros(n,1) ;
pcolor(lon,lat,RF') ; shading interp ;
hold on
for i = 1:n
   x = S(i).X ; y = S(i).Y ;
 plot(x,y,'k')
   % rfi = interp2(lon,lat,RF',x,y) ;

end