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

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

sin_to_ll, mt_subset

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...
#> 
# }