diff --git a/src/transport_performance/population/rasterpop.py b/src/transport_performance/population/rasterpop.py
index fdfe984a..a1419bfb 100644
--- a/src/transport_performance/population/rasterpop.py
+++ b/src/transport_performance/population/rasterpop.py
@@ -282,8 +282,8 @@ def _read_and_clip(
         # cell touched and `from_disk` for improved reading speeds.
         self._xds = (
             rioxarray.open_rasterio(self.__filepath, masked=True)
-            .sel(band=band)
             .rio.clip([aoi_bounds], from_disk=True, all_touched=True)
+            .sel(band=band)
         )
 
         # checks that array has only two dimensions