DENMAN.mds

```R
prj <- "+proj=stere +lat_0=-90 +lon_0=100 +lat_ts=-70"

library(raadtools)
files <- read_amsr2_3k_ice(returnfiles = TRUE)
files <- files[nrow(files):1, ]
im <- vapour::gdal_raster_data(files$fullname[1], target_ext = e, target_crs = prj)
p <- reproj::reproj(cbind(99.3, -64.5), prj, source = "EPSG:4326")
e <- rep(p, each = 2) + c(-1, 1, -1, 1) * 180000
im <- vapour::gdal_raster_data(files$fullname[1], target_ext = e, target_crs = prj)
ximage::ximage(im)
library(SOmap)
SOmap_auto(c(96, 104), c(-72, -66))


```