process_region function

Process regional estimate