How can I draw all boundaries of cities in Pennsylvania State, USA

I have draw a map for Pennsylvanya but I need to divide the map by cities. So, how can I do that? that what I did:

pennsylvania1 = shaperead('usastatelo',...
    'UseGeoCoords', true,...
    'Selector', {@(name)strcmpi('pennsylvania',name), 'Name'});
geoshow(pennsylvania1)

Answer 1

[Rajesh Balagam]

You need to download USA city/county shape file data for this. You can find it (USA_adm_shp.zip) here.

http://www.gadm.org/country

After extracting the data you can plot Pennsylvania state data as shown below.

 data = shaperead('USA_adm2.shp','UseGeoCoords',true,'Selector',{@(name)strcmpi('pennsylvania',name),'NAME_1'});
 geoshow(data)

[Abdulmannan Ajabnoor]

Thank You so much

I tried to do it but this error happened

>> Untitled Error using openShapeFiles>checkSHP (line 67) Unable to open both USA_adm2.shp and USA_adm2.SHP.

Error in openShapeFiles (line 18) [basename, ext] = checkSHP(basename,shapeExtensionProvided);

Error in shaperead (line 208) = openShapeFiles(filename,'shaperead');

Error in Untitled (line 1) data = shaperead('USA_adm2.shp','UseGeoCoords',true,'Selector',{@(name)strcmpi('Pittsburgh',name),'NAME_2'});

Do I should to put all files under Matlap-->Toolbox-->map-->mapdata?

Thank You again

[Rajesh Balagam]

This error is produced when MATLAB couldn't find the specified file.

You need to change the MATLAB present working directory path to the directory where downloaded files are placed (or) copy the USA_adm2.shp file to your working directory.

Execute the above command after this.