Re-ordering source/destination nodata sections

geoutils/raster/raster.py

@@ -2083,7 +2083,7 @@ def reproject(
             # Warning: this will not work for multiple bands with different dtypes
             dtype = self.dtypes[0]
 
-        # Set source nodata if provided
+        # --- Set source nodata if provided -- #
         if force_source_nodata is None:
             src_nodata = self.nodata
         else:
@@ -2097,24 +2097,14 @@ def reproject(
                     )
                 )
 
-        # Create a BoundingBox if required
-        if bounds is not None:
-            if not isinstance(bounds, rio.coords.BoundingBox):
-                bounds = rio.coords.BoundingBox(
-                    bounds["left"],
-                    bounds["bottom"],
-                    bounds["right"],
-                    bounds["top"],
-                )
-
         # --- Set destination nodata if provided -- #
         # This is needed in areas not covered by the input data.
         # If None, will use GeoUtils' default, as rasterio's default is unknown, hence cannot be handled properly.
         if nodata is None:
             nodata = self.nodata
             if nodata is None:
                 nodata = _default_nodata(dtype)
-                # if nodata is already being used, raise a warning.
+                # If nodata is already being used, raise a warning.
                 # TODO: for uint8, if all values are used, apply rio.warp to mask to identify invalid values
                 if not self.is_loaded:
                     warnings.warn(
@@ -2129,6 +2119,16 @@ def reproject(
                         f"self.set_nodata()."
                     )
 
+        # Create a BoundingBox if required
+        if bounds is not None:
+            if not isinstance(bounds, rio.coords.BoundingBox):
+                bounds = rio.coords.BoundingBox(
+                    bounds["left"],
+                    bounds["bottom"],
+                    bounds["right"],
+                    bounds["top"],
+                )
+
         from geoutils.misc import resampling_method_from_str
 
         # --- Basic reprojection options, needed in all cases. --- #