of report categories for each categoryi filter down the data set to the coordinates labeled by that categoryplot basemap on the appropriate subplottransform the coordinatesthen plot the points using the basemap' plot methodfigaxes plt subplots(nrows= ncols= figsize=( )fig subplots_adjust(hspace= wspace= to_plot [' '' '' '' 'lllat= urlat= lllon=- urlon=- for codeax in zip(to_plotaxes flat) basic_haiti_map(axlllat=lllaturlat=urlatlllon=lllonurlon=urloncat_data data[data['category_%scode= compute map proj coordinates xy (cat_data longitudecat_data latitudem plot(xy' 'alpha= ax set_title('% % (codeenglish_mapping[code])the resulting figure can be seen in figure - it seems from the plot that most of the data is concentrated around the most populous cityport-au-prince basemap allows you to overlap additional map data which comes from what are called shapefiles first downloaded shapefile with roads in port-auprince (see conveniently has readshapefile method so thatafter extracting the road data archivei added just the following lines to my code plotting and visualization www it-ebooks info