Converts lower-left sinusoidal coordinates to lat-lon sf bounding box

mt_bbox(xllcorner, yllcorner, cellsize, nrows, ncols, transform = TRUE)

Arguments

xllcorner

lower left x coordinate as provided by mt_subset

yllcorner

lower left y coordinate as provided by mt_subset

cellsize

cell size provided by mt_subset

nrows

cell size provided by mt_subset

ncols

cell size provided by mt_subset

transform

transform the bounding box from sin to lat long coordinates, TRUE or FALSE (default = TRUE)

See also

Examples

# \donttest{ # Download some test data subset <- mt_subset(product = "MOD11A2", lat = 40, lon = -110, band = "LST_Day_1km", start = "2004-01-01", end = "2004-03-31", progress = FALSE) # convert sinusoidal to lat / lon lat_lon <- sin_to_ll(subset$xllcorner, subset$yllcorner)
#> Warning: GDAL Message 1: +init=epsg:XXXX syntax is deprecated. It might return a CRS with a non-EPSG compliant axis order.
# bind with the original dataframe subset <- cbind(subset, lat_lon) # convert to bounding box bb <- apply(subset, 1, function(x){ mt_bbox(xllcorner = x['xllcorner'], yllcorner = x['yllcorner'], cellsize = x['cellsize'], nrows = x['nrows'], ncols = x['ncols']) }) head(bb)
#> $`1.1` #> Geometry set for 1 feature #> geometry type: POLYGON #> dimension: XY #> bbox: xmin: -110.0157 ymin: 40 xmax: -109.9914 ymax: 40.00833 #> geographic CRS: WGS 84
#> POLYGON ((-110.0023 40, -110.0157 40.00833, -11...
#> #> $`2.1` #> Geometry set for 1 feature #> geometry type: POLYGON #> dimension: XY #> bbox: xmin: -110.0157 ymin: 40 xmax: -109.9914 ymax: 40.00833 #> geographic CRS: WGS 84
#> POLYGON ((-110.0023 40, -110.0157 40.00833, -11...
#> #> $`3.1` #> Geometry set for 1 feature #> geometry type: POLYGON #> dimension: XY #> bbox: xmin: -110.0157 ymin: 40 xmax: -109.9914 ymax: 40.00833 #> geographic CRS: WGS 84
#> POLYGON ((-110.0023 40, -110.0157 40.00833, -11...
#> #> $`4.1` #> Geometry set for 1 feature #> geometry type: POLYGON #> dimension: XY #> bbox: xmin: -110.0157 ymin: 40 xmax: -109.9914 ymax: 40.00833 #> geographic CRS: WGS 84
#> POLYGON ((-110.0023 40, -110.0157 40.00833, -11...
#> #> $`5.1` #> Geometry set for 1 feature #> geometry type: POLYGON #> dimension: XY #> bbox: xmin: -110.0157 ymin: 40 xmax: -109.9914 ymax: 40.00833 #> geographic CRS: WGS 84
#> POLYGON ((-110.0023 40, -110.0157 40.00833, -11...
#> #> $`6.1` #> Geometry set for 1 feature #> geometry type: POLYGON #> dimension: XY #> bbox: xmin: -110.0157 ymin: 40 xmax: -109.9914 ymax: 40.00833 #> geographic CRS: WGS 84
#> POLYGON ((-110.0023 40, -110.0157 40.00833, -11...
#>
# }