使用R语言mosaic合并多份GeoTIFF文件报错，求优化方案

Hey there! Let's fix that GeoTIFF mosaic issue you're running into. Your current script has a few key missteps that are triggering that error in as.data.frame(y) message—let's break them down and get your raster merge working smoothly:

1. You're using the wrong function for raster data #

The merge() function in base R is built to combine data frames, not spatial raster files like GeoTIFFs. For merging elevation rasters, you need dedicated functions from raster-processing packages (either the classic raster package or the newer, more efficient terra package).

2. Typo in list.files() parameter #

Your script uses file.name=TRUE—that's a typo! The correct parameter to get full file paths (so R can locate your GeoTIFFs properly) is full.names=TRUE. Without this, list.files() only returns filenames, not the full path to your folder, which R needs to load the files.

3. Missing package loading #

You haven't loaded a raster-processing package, which is mandatory for handling GeoTIFFs and running mosaic functions.

Optimized Script (Using terra Package) #

The terra package is the modern replacement for raster and performs better with large datasets like elevation rasters. Here's a clean, working version:

# Load the terra package (install first if needed: install.packages("terra"))
library(terra)

# Set your working directory (replace with your actual folder path)
setwd("/path/to/your/selected-tiffs-folder")

# List all GeoTIFF files with full file paths
tif_files <- list.files(pattern = "\\.tiff$", full.names = TRUE)

# Load all rasters into a SpatRaster collection
raster_collection <- rast(tif_files)

# Mosaic the rasters, using mean for overlapping areas
merged_elevation <- mosaic(raster_collection, fun = mean)

# Optional: Save the merged raster to a new GeoTIFF
writeRaster(merged_elevation, filename = "global_selected_elevation.tiff", overwrite = TRUE)

Alternative (Using raster Package) #

If you prefer the older raster package, this version will also work:

# Load the raster package (install first if needed: install.packages("raster"))
library(raster)

# Set working directory
setwd("/path/to/your/selected-tiffs-folder")

# List files with full paths
tif_files <- list.files(pattern = "\\.tiff$", full.names = TRUE)

# Create a list of raster objects
raster_list <- lapply(tif_files, raster)

# Define a helper function for mosaic (handles overlapping pixels with mean)
mosaic_mean <- function(...) {
  mean(c(...), na.rm = TRUE)
}

# Merge the rasters
merged_elevation <- do.call(mosaic, c(raster_list, fun = mosaic_mean))

# Save the result
writeRaster(merged_elevation, filename = "global_selected_elevation_raster.tiff", overwrite = TRUE)

Quick Notes: #

• Replace /path/to/your/selected-tiffs-folder with the actual path to your folder of target GeoTIFFs.
• The pattern = "\\.tiff$" ensures we only match files ending with .tiff—adjust to \\.tif$ if your files use the shorter extension.
• Adding na.rm = TRUE to the mean function ensures missing values in overlapping areas don't break the merge (remove it if you want NA where any pixel is missing).

内容的提问来源于stack exchange，提问作者Bing-Hong Huang