Create get_by_province_gis.R