improve accordance of basegeo with envi

IREA-CNR-MI · Jun 9, 2020 · ad1e642 · ad1e642
1 parent 00e59bd
commit ad1e642
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/R/prisma_basegeo.R b/R/prisma_basegeo.R
@@ -48,13 +48,15 @@ prisma_basegeo <- function(band, lon, lat, fill_gaps = TRUE) {
     minlat  <- min(lat, na.rm = TRUE)
     vals    <- raster::values(band)
 
-    columns <- round((lon - minlon) / psize_x)
+    columns <- round((lon - minlon) / psize_x) + 1
     rows    <- nrows - round((lat - minlat) / psize_y)
 
     # columns <- round((lon - minlon) / psize_x) + 1
     # rows    <- nrows - round((lat - minlat) / psize_y) - 1
-    # columns[columns > ncols] <- 0
-    # rows[rows > nrows] <- 0
+
+    #♫ remove data if out of bounds to avoid potential crashes
+    columns[columns > ncols] <- 0
+    rows[rows > nrows] <- 0
 
     for (indpix in 1:(numlines*numcols)) {
         out_grd[rows[indpix], columns[indpix]] <- vals[indpix]