How to use UtmZoneSplitter for zones that span across multiple UTM zones?

I am using UtmZoneSplitter and UtmGridSplitter to divide a zone that spans across multiple UTM zones. I then create the bbox_list and bbox_list and having that I create the GeoDataFrame which I use for plotting. I am following this example How To: Land-Use-Land-Cover Prediction for Slovenia. But when the plot is created it seems like the bboxes are divided into two groups. I have tested it in zones that do not span across multiple zones and it works perfectly.

country = gpd.read_file(os.path.join(DATA_FOLDER, 'COLOMBIA.shp'))

depto = country[country.DPTO_CNMBR == "CASANARE"]

# Get the shape in polygon format
country_shape = country.geometry.values[-1]

depto_shape = depto.geometry.values[-1]


This shows the image of the area correcly.

Here is the part where I get the bbox_list and plot the bboxes

# Create the splitter to obtain a list of bboxes
bbox_splitter = UtmGridSplitter([depto_shape], CRS.WGS84, 5000)

bbox_list = np.array(bbox_splitter.get_bbox_list())
info_list = np.array(bbox_splitter.get_info_list())

# Prepare info of selected EOPatches
geometry = [Polygon(bbox.get_polygon()) for bbox in bbox_list]

idxs = [info['index'] for info in info_list]
idxs_x = [info['index_x'] for info in info_list]
idxs_y = [info['index_y'] for info in info_list]

gdf = gpd.GeoDataFrame({'index': idxs, 'index_x': idxs_x, 'index_y': idxs_y}, 

# figure
fig, ax = plt.subplots(figsize=(30, 30))
ax.set_title('Selected 5x5  tiles from Colombia', fontsize=30);
for bbox, info in zip(bbox_list, info_list):
    geo = bbox.geometry
    ax.text(geo.centroid.x, geo.centroid.y, info['index'], ha='center', va='center')

But this plots the image below. Any idea on how to plot the boxes correctly?

Answered in the GitHub issue:

1 Like