diff --git a/geotessera/cli.py b/geotessera/cli.py
index 0feefda..1054e82 100644
--- a/geotessera/cli.py
+++ b/geotessera/cli.py
@@ -459,7 +459,7 @@ def info(
             for tile in coverage["tiles"][:10]:
                 tiles_table.add_row(
                     Path(tile["path"]).name,
-                    f"({tile['tile_lat']}, {tile['tile_lon']})",
+                    f"({tile['tile_lon']}, {tile['tile_lat']})",
                     str(tile["bands"]),
                 )
 
diff --git a/geotessera/core.py b/geotessera/core.py
index 5e605a7..3ad148d 100644
--- a/geotessera/core.py
+++ b/geotessera/core.py
@@ -445,7 +445,7 @@ def tile_progress_callback(
 
             except Exception as e:
                 self.logger.warning(
-                    f"Failed to download tile ({tile_lat:.2f}, {tile_lon:.2f}): {e}"
+                    f"Failed to download tile ({tile_lon:.2f}, {tile_lat:.2f}): {e}"
                 )
                 if progress_callback:
                     progress_callback(
@@ -697,11 +697,11 @@ def update_progress(status: str):
             tile = local_tile_map.get((tile_lon, tile_lat))
             if tile is None or not tile.is_available():
                 self.logger.warning(
-                    f"Tile ({tile_lat:.2f}, {tile_lon:.2f}) not available after download, skipping"
+                    f"Tile ({tile_lon:.2f}, {tile_lat:.2f}) not available after download, skipping"
                 )
                 continue
 
-            update_progress(f"Loading tile ({tile_lat:.2f}, {tile_lon:.2f})")
+            update_progress(f"Loading tile ({tile_lon:.2f}, {tile_lat:.2f})")
 
             try:
                 # Load embedding from Tile (handles both NPY and GeoTIFF formats)
@@ -721,7 +721,7 @@ def update_progress(status: str):
                 # Ensure dimensions are valid integers
                 if dst_width is None or dst_height is None:
                     raise ValueError(
-                        f"Failed to calculate dimensions for tile ({tile_lat}, {tile_lon})"
+                        f"Failed to calculate dimensions for tile ({tile_lon}, {tile_lat})"
                     )
                 dst_width = int(dst_width)
                 dst_height = int(dst_height)
@@ -758,11 +758,11 @@ def update_progress(status: str):
                     dataset.write(reprojected_embedding)
 
                 reprojected_memfiles.append(memfile)
-                update_progress(f"Reprojected tile ({tile_lat:.2f}, {tile_lon:.2f})")
+                update_progress(f"Reprojected tile ({tile_lon:.2f}, {tile_lat:.2f})")
 
             except Exception as e:
                 self.logger.warning(
-                    f"Failed to process tile ({tile_lat:.2f}, {tile_lon:.2f}): {e}"
+                    f"Failed to process tile ({tile_lon:.2f}, {tile_lat:.2f}): {e}"
                 )
                 continue
 
@@ -918,7 +918,7 @@ def sample_embeddings_at_points(
                 progress_callback(
                     tile_idx,
                     total_tiles,
-                    f"Processing tile {tile_idx + 1}/{total_tiles}: ({tile_lat:.2f}, {tile_lon:.2f})",
+                    f"Processing tile {tile_idx + 1}/{total_tiles}: ({tile_lon:.2f}, {tile_lat:.2f})",
                 )
 
             try:
@@ -985,7 +985,7 @@ def sample_embeddings_at_points(
             except Exception as e:
                 # If tile fetch/load fails, leave those points as NaN
                 self.logger.warning(
-                    f"Failed to process tile ({tile_lat:.2f}, {tile_lon:.2f}): {e}"
+                    f"Failed to process tile ({tile_lon:.2f}, {tile_lat:.2f}): {e}"
                 )
                 if include_metadata:
                     for original_idx in point_indices:
@@ -1208,7 +1208,7 @@ def download_tile(
             return True
 
         except Exception as e:
-            self.logger.error(f"Failed to download tile ({lat:.2f}, {lon:.2f}): {e}")
+            self.logger.error(f"Failed to download tile ({lon:.2f}, {lat:.2f}): {e}")
             return False
 
     def download_tiles_for_points(
@@ -2157,7 +2157,7 @@ def apply_pca_to_embeddings(
                 progress_callback(
                     i,
                     total_tiles,
-                    f"Processing tile {i + 1}/{total_tiles}: ({tile_lat:.2f}, {tile_lon:.2f})",
+                    f"Processing tile {i + 1}/{total_tiles}: ({tile_lon:.2f}, {tile_lat:.2f})",
                 )
 
             # Reshape for PCA: (height, width, channels) -> (pixels, channels)