I managed to plot zip codes boundaries for an entire US state using ArcGIS data and python tools such as arcpy and Bokeh. The result:
The dataset has high resolution area boundaries so, I have like 2000 data points for each polygon and there are 1700 polygons. Thus, rendering 3.4 million data points approx. is very slow and the interface lags.
The questions is, do you know any methodology to reduce the area boundary resolution so the number of data points to be plotted is reduced? The resulting polygons should not overlap at all. Any other alternatives?