A full description of the sinusoidal projection is provided on the lpdaac page: https://lpdaac.usgs.gov/dataset_discovery/modis and wikipedia: https://en.wikipedia.org/wiki/Sinusoidal_projection

sin_to_ll(x, y)

Arguments

x

sinusoidal x coordinate (vector)

y

sinusoidal y coordinate (vector)

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) # bind with the original dataframe subset <- cbind(subset, lat_lon) head(subset)
#> xllcorner yllcorner cellsize nrows ncols band units #> 1.1 -9370036.39 4447802.08 926.625433055834 1 1 LST_Day_1km Kelvin #> 2.1 -9370036.39 4447802.08 926.625433055834 1 1 LST_Day_1km Kelvin #> 3.1 -9370036.39 4447802.08 926.625433055834 1 1 LST_Day_1km Kelvin #> 4.1 -9370036.39 4447802.08 926.625433055834 1 1 LST_Day_1km Kelvin #> 5.1 -9370036.39 4447802.08 926.625433055834 1 1 LST_Day_1km Kelvin #> 6.1 -9370036.39 4447802.08 926.625433055834 1 1 LST_Day_1km Kelvin #> scale latitude longitude site product start end complete #> 1.1 0.02 40 -110 sitename MOD11A2 2004-01-01 2004-03-31 TRUE #> 2.1 0.02 40 -110 sitename MOD11A2 2004-01-01 2004-03-31 TRUE #> 3.1 0.02 40 -110 sitename MOD11A2 2004-01-01 2004-03-31 TRUE #> 4.1 0.02 40 -110 sitename MOD11A2 2004-01-01 2004-03-31 TRUE #> 5.1 0.02 40 -110 sitename MOD11A2 2004-01-01 2004-03-31 TRUE #> 6.1 0.02 40 -110 sitename MOD11A2 2004-01-01 2004-03-31 TRUE #> modis_date calendar_date tile proc_date pixel value longitude_ll #> 1.1 A2004001 2004-01-01 h09v05 2015212185706 1 13098 -110.0023 #> 2.1 A2004009 2004-01-09 h09v05 2015212201022 1 13062 -110.0023 #> 3.1 A2004017 2004-01-17 h09v05 2015212213103 1 13297 -110.0023 #> 4.1 A2004025 2004-01-25 h09v05 2015213005429 1 13323 -110.0023 #> 5.1 A2004033 2004-02-02 h09v05 2015213090158 1 13315 -110.0023 #> 6.1 A2004041 2004-02-10 h09v05 2015213165253 1 13227 -110.0023 #> latitude_ll #> 1.1 40 #> 2.1 40 #> 3.1 40 #> 4.1 40 #> 5.1 40 #> 6.1 40
# }