<jupyter_start><jupyter_text>Mini challenge: create a pop density map for PHL<jupyter_code>php_array, clipped_transform = msk.mask(tif_file, [mapping(geom) for geom in php.geometry.tolist()], crop=True)

fig, ax = plt.subplots(facecolor='#FCF6F5FF')
fig.set_size_inches(21, 14)

ax.imshow(php_array[0], norm=norm, cmap=our_cmap)
ax.axis('off')
ax.set_title('Population Density Heatmap of The Philippines', {'fontsize':15})

#map the data to clors
#imaj is now in RGBA (512x512x4)
imaj = our_cmap(norm(php_array[0]))

plt.show()
#plt.imsave('./resources/php_popdensity.png', imaj, dpi=2000)
php_array[0]
php_array[0].shape<jupyter_output><empty_output>