compute blockwise region features