diff --git a/DESCRIPTION b/DESCRIPTION
index fc2f52a..01eb9b0 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -11,6 +11,6 @@ URL: https://github.com/nplatonov/ursa
 BugReports: https://github.com/nplatonov/ursa/issues
 Depends: R (>= 4.1.0)
 Imports: utils, graphics, grDevices, stats, rgdal, png, jpeg
-Suggests: jsonlite, proj4, sf (>= 0.6-1), raster, ncdf4, locfit, knitr, rmarkdown, tcltk, sp, methods, fasterize, IRdisplay, caTools, shiny, rgeos, tools, webp, htmlwidgets, htmltools, leaflet, leafem, leafpop, RColorBrewer, ragg, widgetframe, geojsonsf (>= 2.0.0), leaflet.providers, magick, terra, stars, vapour, sys
+Suggests: jsonlite, proj4, sf (>= 0.6-1), raster, ncdf4, locfit, knitr, rmarkdown, tcltk, sp, methods, fasterize, IRdisplay, caTools, shiny, rgeos, tools, webp, htmlwidgets, htmltools, leaflet, leafem, leafpop, RColorBrewer, ragg, widgetframe, geojsonsf (>= 2.0.0), leaflet.providers, magick, terra, stars, vapour, sys, RSQLite
 NeedsCompilation: yes
 ByteCompile: no
diff --git a/NEWS.md b/NEWS.md
index c608d6f..7baa98f 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -1,4 +1,4 @@
-2023-03-24
+2023-08-24
 
 <!--
 + Depends on R (>= 4.1.0) instead of R (>= 3.0.0) due to using native pipes `|>` in code
@@ -8,11 +8,19 @@
 
 ### version 3.9.11
 
--   ongoing…
+-   Repaired parsing of Nominatim responses.
+
+-   Package **`RSQLite`** is added as suggested for using map tiles from
+    specified directory of SAS Planet cache
+
+-   In `segmentize()` 1) improved response of `by` argument, 2) order of
+    `type` argument is reversed to `c("united","conseqwent")`
+
+-   Coastline is updated to version 2023-08-24T03:31.
 
 ### version 3.9.10
 
--   Foreign functions are symbolic insead of character after [R-devel
+-   Foreign functions are symbolic instead of character after [R-devel
     BUG FIXES
     2023-03-20](https://developer.r-project.org/blosxom.cgi/R-devel/NEWS/2023/03/20#n2023-03-20).
 
diff --git a/R/compose_legend.R b/R/compose_legend.R
index 1dd994d..2d83f76 100644
--- a/R/compose_legend.R
+++ b/R/compose_legend.R
@@ -148,7 +148,7 @@
          arglist2 <- c(quote(obj[[i]]),arglist) ## 20180308 change 'arglist[-1]'?
         # arglist[[1]] <- quote(obj[[i]])
          if (.is.colortable(obj[[i]])) {
-            arglist2[["units"]] <- units[i]
+            arglist2[["unit"]] <- units[i]
            # str(arglist2)
             do.call("legend_colorbar",arglist2)
            # arglist[["units"]] <- NULL
diff --git a/R/conn.open_gdal.R b/R/conn.open_gdal.R
index 175c5c4..4a7741e 100644
--- a/R/conn.open_gdal.R
+++ b/R/conn.open_gdal.R
@@ -8,9 +8,9 @@
   #    .elapsedTime("rgdal has been loaded")
    opW <- options(warn=0-!verbose,show.error.messages=TRUE) ## to prevent 'GeoTransform values not available'
    on.exit(options(opW))
-   if (devel <- T) {
-      a <- vapour:::gdalinfo_internal(fname,json=TRUE,stats=FALSE)
-      a <- jsonlite::fromJSON(a)
+   if (devel <- F) {
+     # a <- vapour:::gdalinfo_internal(fname,json=TRUE,stats=FALSE)
+     # a <- jsonlite::fromJSON(a)
    }
    a <- try(rgdal::GDALinfo(fname,returnStats=FALSE,returnRAT=FALSE
                 ,returnColorTable=TRUE,returnCategoryNames=TRUE),silent=TRUE)
diff --git a/R/conn.read_gdal.R b/R/conn.read_gdal.R
index 9c46d4d..2bdeafb 100644
--- a/R/conn.read_gdal.R
+++ b/R/conn.read_gdal.R
@@ -81,7 +81,6 @@
      # str(ds <- sf::gdal_subdatasets(fname,name=TRUE))
       opW <- options(warn=ifelse(.isPackageInUse(),-1,1))
       res <- as_ursa(sf::gdal_read(fname))
-      options(opW)
       if (!is.null(band))
          res <- res[band]
    }
@@ -92,7 +91,8 @@
       res <- if (!is.null(band)) obj[band] else obj[]
       close(obj)
    }
-   if (T & length(grep("^\\d{8}\\.s1ab\\.1km\\.n\\.mos[13]d\\.jpg$",basename(fname)))) {
+   if (T & length(grep("^(\\d{8}\\.s1ab\\.1km\\.n\\.mos[13]d|.+sentinel1-n-[13]daymos)\\.jpg$"
+                      ,basename(fname)))) {
      ## patch to seaice.dk Sentinel-1 mosaic
       g0 <- ursa_grid(res)
       if ((g0$columns==4500L)&&(g0$rows==5500L)) {
diff --git a/R/glance.R b/R/glance.R
index 1b990c4..ebfddba 100644
--- a/R/glance.R
+++ b/R/glance.R
@@ -247,7 +247,19 @@
       s2 <- strsplit(style,split="\\s+")
       s <- expand.grid(s1,s2)
    }
-   isUrl <- .lgrep("http(s)*://",style)
+   canUrl <- length(unlist(regmatches(style,gregexpr("\\{[xyz]\\}",style))))==3
+   if (!canUrl) {
+      sascache <- gsub("\\s(color|grey|gray)$","",style)
+      if (!length(ind <- which(dir.exists(sascache))))
+         sascache <- file.path(getOption("SAS_Planet_cache"),sascache)
+      if (length(ind <- which(dir.exists(sascache)))>0) {
+         list1 <- dir(path=sascache[ind],pattern="\\.sqlitedb$",recursive=TRUE)
+         canUrl <- length(list1)>0
+         if (canUrl)
+            style <- basename(sascache)
+      }
+   }
+   isUrl <- .lgrep("http(s)*://",style) | canUrl
    if (isUrl) {
       ind <- .grep("http(s)*://",style)
       style[ind] <- unlist(strsplit(style[ind],split="\\s+"))[1]
@@ -282,7 +294,7 @@
          nextStyle <- .grep(ostyle[1],staticMap
                            ,invert=TRUE,value=TRUE)
       else
-         nextStyle <- .grep(ostyle[1],c("mapnik","mapsurfer","cartoDB","opentopomap")
+         nextStyle <- .grep(ostyle[1],c("mapnik","cartoDB","opentopomap")
                            ,invert=TRUE,value=TRUE)[seq(3)]
       nsize <- length(nextStyle)+1
       cache <- .getPrm(arglist,"cache",class=c("logical","character"),default=TRUE)
@@ -510,6 +522,8 @@
             panel_new(...) #fill=ifelse(isWeb,"transparent","chessboard"))
          if (before) {
             panel_plot(basemap,alpha=basemap.alpha)
+            if (toCoast)
+               panel_coastline(cline)
          }
         # if ((!length(ct))||(all(is.na(ct[[i]]$index)))) {
          if ((!length(ct))||(!hasField)) {
@@ -583,9 +597,9 @@
          }
          if (after) {
             panel_plot(basemap,alpha=basemap.alpha)
+            if (toCoast)
+               panel_coastline(cline)
          }
-         if (toCoast)
-            panel_coastline(cline)
          if ((geocodeStatus)||(!length(ct)))
             panel_graticule(gline,margin=c(T,T,F,F))
          else
diff --git a/R/legend_colorbar.R b/R/legend_colorbar.R
index 47cebea..6e7ecdf 100644
--- a/R/legend_colorbar.R
+++ b/R/legend_colorbar.R
@@ -36,7 +36,7 @@
    }
    units <- .getPrm(arglist,name="unit(s)*",kwd=kwd
                    ,class=list("character","expression"),default="",verbose=FALSE)
-   if (!nchar(units)) {
+   if (all(!nchar(units))) {
       aname <- names(arglist[1])
       if (!is.null(aname))
          units <- names(arglist[1])
@@ -799,7 +799,7 @@
       panel_box()
    }
   # options(ursaWidth=width,ursaHeight=height)
-   if (((is.character(units))&&(is.na(units)))||(!nchar(units)))
+   if (((is.character(units))&&(is.na(units[1])))||(!nchar(units[1])))
       return(invisible(10L))
   # b <- 2*width+0.5+height
   # print(b)
diff --git a/R/panel_annotation.R b/R/panel_annotation.R
index 272ecb9..786e1f2 100644
--- a/R/panel_annotation.R
+++ b/R/panel_annotation.R
@@ -75,7 +75,7 @@
       str(list(label=label,position=position,cex=cex,adjust=adjust,fg=fg,bg=bg
               ,fill=fill,buffer=buffer,vertical=vertical,verbose=verbose))
    opt <- par(family=font)
-   g1 <- session_grid()
+   g1 <- getOption("ursaPngPanelGrid") # session_grid()
    minx <- g1$minx
    miny <- g1$miny
    maxx <- g1$maxx
@@ -98,8 +98,8 @@
       pos <- c((xy[1,1]-minx)/(maxx-minx),(xy[1,2]-miny)/(maxy-miny))
    }
    else if ((!is.na(x))&&(!is.na(y))) {
-      if ((x<minx)||(x>maxx)||(y<miny)||(y>maxy))
-         return(NULL)
+     # if ((x<minx)||(x>maxx)||(y<miny)||(y>maxy))
+     #    return(NULL)
       pos <- c((x-minx)/(maxx-minx),(y-miny)/(maxy-miny))
    }
    else
@@ -117,7 +117,7 @@
       if (!isE)
          label <- paste(label,collapse="\n")
       if (!nchar(label))
-         return (10L)
+         return(invisible(NULL)) ## 10L
    }
    if (is.factor(pos))
       pos <- as.character(pos)
@@ -237,6 +237,8 @@
       if (verbose)
          .elapsedTime("label:finish")
    }
+   if (verbose)
+      str(list(x=x,y=y,adj=c(adjust,vadj),labels=label,cex=mycex,col=fg,srt=srt))
    text(x=x,y=y,adj=c(adjust,vadj),labels=label,cex=mycex,col=fg,srt=srt)
    par(opt)
    invisible(NULL)
diff --git a/R/panel_coastline.R b/R/panel_coastline.R
index a4d1d04..5ea2079 100644
--- a/R/panel_coastline.R
+++ b/R/panel_coastline.R
@@ -628,13 +628,11 @@
 
 'update_coastline' <- function(merge=TRUE) {
    missedSF <- !requireNamespace("sf",quietly=TRUE) 
-   missedLW <- !requireNamespace("lwgeom",quietly=TRUE)
-   if ((missedSF)&&(missedLW))
-      stop("Suggested packages 'sf' and 'lwgeom' are required for this operation")
-   else if (missedSF)
+   if (missedSF)
       stop("Suggested package 'sf' is required for this operation")
-   else if (missedLW)
-      stop("Suggested package 'lwgeom' is required for this operation")
+  # missedLW <- !requireNamespace("lwgeom",quietly=TRUE)
+  # if (missedLW)
+  #    stop("Suggested package 'lwgeom' is required for this operation")
    dpath <- getOption("ursaRequisite")
    ftemp <- tempfile(tmpdir=dpath)
    opW <- options(warn=-1)
diff --git a/R/regrid.R b/R/regrid.R
index d1f7e57..8cc5fc5 100644
--- a/R/regrid.R
+++ b/R/regrid.R
@@ -11,7 +11,7 @@
          session_grid(result)
       return(result)
    }
-   fun <- "resize" #as.character(match.call())[1]
+   fun <- "regrid" #as.character(match.call())[1]
    if (.is.ursa_stack(x)) {
       return(lapply(x,regrid,...))
    }
diff --git a/R/session.R b/R/session.R
index 6df6087..67c03bc 100644
--- a/R/session.R
+++ b/R/session.R
@@ -127,8 +127,11 @@
 'session_proj' <- 'session_proj4' <- 'session_crs' <- function() session_grid()$crs
 'session_cellsize' <- function() with(session_grid(),sqrt(as.numeric(resx)*as.numeric(resy)))
 'session_dim' <- function() with(session_grid(),c(lines=rows,samples=columns))
-'session_bbox' <- function() with(session_grid()
-                                 ,c(minx=minx,miny=miny,maxx=maxx,maxy=maxy))
+'session_bbox' <- function() {
+   ret <- with(session_grid(),c(minx=minx,miny=miny,maxx=maxx,maxy=maxy))
+   attr(ret,"crs") <- session_crs()
+   ret
+}
 'session_pngviewer' <- function(allow=NA) {
    opV <- getOption("ursaAllowPngViewer")
   # str(list(allow=allow,opV=opV,isRscript=.isRscript()))
diff --git a/R/spatial_engine.R b/R/spatial_engine.R
index ed190d7..638ed33 100644
--- a/R/spatial_engine.R
+++ b/R/spatial_engine.R
@@ -814,10 +814,13 @@
                                                ,recursive=FALSE,ignore.case=TRUE) {
    patt0 <- "\\.(gpkg|tab|kml|json|geojson|mif|sqlite|shp|osm)(\\.(zip|gz|bz2))*$"
    if (devel <- TRUE & all(!dir.exists(path))) {
-      dpath <- list.dirs(dirname(path),full.names=FALSE)
-      ind <- grep(basename(path),dpath,ignore.case=ignore.case)
-      if (length(ind)==1)
-         path <- file.path(dirname(path),dpath[ind])
+      dname <- dirname(path)
+      if (dname!=".") {
+         dpath <- list.dirs(dirname(path),full.names=FALSE)
+         ind <- grep(basename(path),dpath,ignore.case=ignore.case)
+         if (length(ind)==1)
+            path <- file.path(dirname(path),dpath[ind])
+      }
    }
    res <- dir(path=path,pattern=patt0,full.names=full.names
              ,recursive=recursive,ignore.case=ignore.case)
@@ -1299,12 +1302,14 @@
 }
 'spatial_bind' <- function(...) {
    arglist <- list(...)
+   if (length(ind <- which(sapply(arglist,is.null))))
+      arglist <- arglist[-ind]
+   if (!length(arglist))
+      return(NULL)
    if ((length(arglist)==1)&&
        (any(sapply(arglist[[1]],is_spatial)))&&
        (!is_spatial(arglist[[1]])))
       arglist <- arglist[[1]]
-   if (length(ind <- which(sapply(arglist,is.null))))
-      arglist <- arglist[-ind]
    res <- arglist[[1]]
    isSF <- .isSF(res)
    isSP <- .isSP(res)
diff --git a/R/ursa_as.R b/R/ursa_as.R
index 5e0e42a..f5c4be3 100644
--- a/R/ursa_as.R
+++ b/R/ursa_as.R
@@ -400,10 +400,14 @@
          ##~ str(a1[[1]]$finalize())
          ##~ q()
       }
-      bbox <- obj@ptr$extent$vector[c(1,3,2,4)]
-      res <- obj@ptr$res
-      crs <- obj@ptr$get_crs("proj4")
-      aname <- obj@ptr$names
+     # bbox <- obj@ptr$extent$vector[c(1,3,2,4)]
+     # res <- obj@ptr$res
+     # crs <- obj@ptr$get_crs("proj4")
+     # aname <- obj@ptr$names
+      bbox <- as.vector(obj)[c(1,3,2,4)]
+      res <- terra::res(obj)
+      crs <- terra::crs(obj,proj=TRUE)
+      aname <- terra::names(obj)
       g1 <- regrid(bbox=bbox,res=res,crs=crs)
       if (identical(bbox,c(0,0,1,1)))
          g1 <- regrid(bbox=c(0,0,rev(dim(g1))),res=1,crs=crs)
diff --git a/R/xxx.geocode.R b/R/xxx.geocode.R
index 184cba6..d52c9c9 100644
--- a/R/xxx.geocode.R
+++ b/R/xxx.geocode.R
@@ -52,6 +52,7 @@
                   # ,"&polygon_text=1"
                    ,"&format=xml","&bounded=0","&accept-language=en-US,ru")
      # dst <- tempfile() # "nominatim.xml" # tempfile()
+     # print(src)
       dst <- .ursaCacheDownload(src,quiet=!verbose)
       xmlstring <- scan(dst,character(),quiet=!verbose)
      # Encoding(xmlstring) <- "UTF-8"
@@ -64,11 +65,11 @@
          print(geotext)
       }
       ptype <- .grep("^type=",xmlstring,value=TRUE)
-      ptype <- .gsub(".*\'(.+)\'.*","\\1",ptype)
+      ptype <- .gsub(".*(\'|\")(.+)(\'|\").*","\\2",ptype)
       pclass <- .grep("^class=",xmlstring,value=TRUE)
-      pclass <- .gsub(".*\'(.+)\'.*","\\1",pclass)
+      pclass <- .gsub(".*(\'|\")(.+)(\'|\").*","\\2",pclass)
       prank <- .grep("^place_rank=",xmlstring,value=TRUE)
-      prank <- .gsub(".*\'(.+)\'.*","\\1",prank)
+      prank <- .gsub(".*(\'|\")(.+)(\'|\").*","\\2",prank)
       if (FALSE) {
          prank[prank=="8"] <- "state"
          prank[prank=="10"] <- "region"
@@ -77,12 +78,12 @@
         # prank[prank=="17"] <- "town|island"
         # prank[prank=="18"] <- "village|"
         # prank[prank=="19"] <- ""
-     }
+      }
       ptype <- paste(pclass,ptype,prank,sep=":")
       lon <- .grep("lon=",xmlstring,value=TRUE)
-      lon <- as.numeric(.gsub(".*\'(.+)\'.*","\\1",lon))
+      lon <- as.numeric(.gsub(".*(\'|\")(.+)(\'|\").*","\\2",lon))
       lat <- .grep("lat=",xmlstring,value=TRUE)
-      lat <- as.numeric(.gsub(".*\'(.+)\'.*","\\1",lat))
+      lat <- as.numeric(.gsub(".*(\'|\")(.+)(\'|\").*","\\2",lat))
       pt <- cbind(lon=lon,lat=lat)
       if (!nrow(pt))
          return(NULL)
@@ -99,7 +100,7 @@
       ann <- .grep("display_name=",xmlstring,value=TRUE)
       ann <- .gsub(".*\"(.+)\".*","\\1",ann)
       importance <- .grep("importance",xmlstring,value=TRUE)
-      importance <- as.numeric(.gsub(".*\'(.+)\'.*","\\1",importance))
+      importance <- as.numeric(.gsub(".*(\'|\")(.+)(\'|\").*","\\2",importance))
      # type <- NULL ## debug
       typeInd <- integer()
       if ((is.character(place))&&(nchar(place))) {
@@ -140,8 +141,8 @@
               # print(file.size(dst))
                b <- unlist(strsplit(b,split="'"))
                b <- .grep("(MULTI)*(POLYGON|LINESTRING)",b,value=TRUE,ignore.case=FALSE)
-               b <- .gsub("((MULTI)*(POLYGON|LINESTRING)|\\(|\\))"," ",b)
-               b <- .gsub("(^\\s+|\\s+$)","",b)
+               b <- .gsub("(^.+(MULTI)*(POLYGON|LINESTRING)|\\(|\\))"," ",b)
+               b <- .gsub("(^\\s+|\\s+(\".+)*$)","",b)
                b <- lapply(b,function(b2) {
                   b3 <- unlist(strsplit(b2,split=","))
                   b3 <- unlist(strsplit(b3,split="\\s+"))
diff --git a/R/xxx.geomap.R b/R/xxx.geomap.R
index 2e413db..1c55d64 100644
--- a/R/xxx.geomap.R
+++ b/R/xxx.geomap.R
@@ -33,6 +33,19 @@
      # print(art);q()
      # proj <- "merc"
    }
+   canUrl <- length(unlist(regmatches(style,gregexpr("\\{[xyz]\\}",style))))==3
+   isSAScache <- FALSE
+   if (!canUrl) {
+      sascache <- style
+      if (!length(ind <- which(dir.exists(sascache))))
+         sascache <- file.path(getOption("SAS_Planet_cache"),style)
+      if (length(ind <- which(dir.exists(sascache)))>0) {
+         list1 <- dir(path=sascache[ind],pattern="\\.sqlitedb$",recursive=TRUE)
+         canUrl <- length(list1)>0
+         isSAScache <- TRUE
+         art <- "sascache"
+      }
+   }
    if ((art=="none")&&(length(style)==1)&&
          (requireNamespace("leaflet",quietly=.isPackageInUse()))&&
          (requireNamespace("leaflet.providers",quietly=.isPackageInUse()))) {
@@ -127,9 +140,21 @@
    if (isStatic) {
       len[len>mlen] <- mlen
    }
-   isUrl <- .lgrep("^http(s)*://",style)>0
+   if (trytodeprecate20230723 <- TRUE) {
+      canUrl <- length(unlist(regmatches(style,gregexpr("\\{[xyz]\\}",style))))==3
+      if (!canUrl) {
+         sascache <- style
+         if (!length(ind <- which(dir.exists(sascache))))
+            sascache <- file.path(getOption("SAS_Planet_cache"),style)
+         if (length(ind <- which(dir.exists(sascache)))>0) {
+            list1 <- dir(path=sascache[ind],pattern="\\.sqlitedb$",recursive=TRUE)
+            canUrl <- length(list1)>0
+         }
+      }
+   }
+   isUrl <- .lgrep("^http(s)*://",style)>0 | canUrl
   # canTile <- .lgrep(art,eval(as.list(args(".tileService()"))$server))>0
-   canTile <- isUrl | .lgrep(art,.tileService())>0
+   canTile <- isUrl | .lgrep(art,.tileService())>0 | canUrl
    isTile <- .lgrep("tile",style)>0 & canTile
    if ((!isStatic)&&(!isTile)) {
       if (art %in% staticMap)
@@ -153,7 +178,8 @@
    isWeb <- .lgrep(tilePatt,art)>0 | isUrl
    if (verbose)
       print(data.frame(art=art,color=isColor,grey=isGrey,static=isStatic
-                      ,canTile=canTile,tile=isTile,web=isWeb,row.names="geomap:"))
+                      ,canTile=canTile,canUrl=canUrl,tile=isTile,web=isWeb
+                      ,row.names="geomap:"))
    geocodeStatus <- FALSE
    if (.isSP(loc)) {
       proj4 <- sp::proj4string(loc)
@@ -174,7 +200,6 @@
          loc <- c(.project(matrix(loc,ncol=2,byrow=TRUE),proj4
                           ,inv=TRUE))[c(1,3,2,4)]
    }
-   
    isWMS <- isUrl & .is.wms(style)
    notYetGrid <- TRUE
    g3 <- NULL
@@ -242,7 +267,7 @@
          lon_0[lon_0>=(-70) && lon_0<(-25)] <- -40
          lon_0[lon_0>=(-25) && lon_0<(+50)] <- 10
          lon_0[lon_0>=(50) && lon_0<(+135)] <- 90
-         proj4 <- paste("","+proj=laea +lat_0=90",paste0("+lon_0=",lon_0)
+         proj4 <- paste("+proj=laea +lat_0=90",paste0("+lon_0=",lon_0)
                        ,"+x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs")
          if (bbox[3]<bbox[1])
             bbox[3] <- bbox[3]+360
@@ -455,9 +480,12 @@
          tgr$maxy <- xy["y",]+res/2
          h <- sort(unique(tgr[,"x"]))
          v <- sort(unique(tgr[,"y"]))
-         cat("-=---------------------------------------------\n")
+         if (F & !.isPackageInUse())
+            cat("-=---------------------------------------------\n")
          g1 <- with(g0,regrid(g1,bbox=c(minx,miny,maxx,maxy),crs=proj4
-                             ,zero="keep",verbose=!FALSE))
+                             ,zero="keep",verbose=F & !.isPackageInUse()))
+         if (F & !.isPackageInUse())
+            cat("-=---------------------------------------------\n")
          if (TRUE) {
             x <- seq(g1$minx,g1$maxx,by=g1$resx)
             y <- seq(g1$miny,g1$maxy,by=g1$resy)
@@ -639,12 +667,13 @@
          tgr$lat2 <- y2
       }
       img1 <- vector("list",nrow(tgr))
-      for (i in sample(seq(nrow(tgr))))
+      for (i in sample(seq(nrow(tgr)))) {
          img1[[i]] <- .tileGet(z=zoom,x=tgr[i,"x"],y=tgr[i,"y"]
                               ,minx=tgr[i,"minx"],miny=tgr[i,"miny"]
                               ,maxx=tgr[i,"maxx"],maxy=tgr[i,"maxy"]
                               ,retina=retina,url=tile$url
                               ,fileext=tile$fileext,cache=cache,verbose=verbose)
+      }
       nb <- sapply(img1,function(x) {
          if (!is.array(x))
             return(0)
@@ -673,17 +702,20 @@
             x
          })
       }
+      ind <- sapply(img1,function(x) length(dim(x))==3)
      # str(lapply(img1,dim))
-      dimb <- apply(list2DF(lapply(img1,dim)),1,max)
+      dimb <- apply(list2DF(lapply(img1[ind],dim)),1,max)
       img <- array(0L,dim=c(dimb[1]*length(v),dimb[2]*length(h),nbmax))
       for (i in sample(seq(nrow(tgr)))) {
         # img[igr[i,"y"]*256L+seq(256),igr[i,"x"]*256+seq(256),] <- img2[,,1:3]
         # img[igr[i,"y"]*256L+seq(256),igr[i,"x"]*256+seq(256),] <- img2[,,seq(nb)]
         # img[igr[i,"y"]*256L+seq(256),igr[i,"x"]*256+seq(256),seq(nb)] <- img2[,,seq(nb)]
          img2 <- img1[[i]]
+         dima <- dim(img2)
+         if (!length(dima))
+            next
          if (inherits(img2,"try-error"))
             next
-         dima <- dim(img2)
          if (!((dima[1]==dimb[1])&&(dima[2]==dimb[1]))) {
            # .elapsedTime("everytime 0205a")
             img2 <- as.array(regrid(as.ursa(img2),res=c(dima[1]/dimb[1],dima[2]/dimb[1])))
diff --git a/R/xxx.panel_cluster.R b/R/xxx.panel_cluster.R
index b7aab10..dcc9848 100644
--- a/R/xxx.panel_cluster.R
+++ b/R/xxx.panel_cluster.R
@@ -3,15 +3,11 @@
                            ,fun=c("count","sum","mean","label"),label=fun %in% "count"
                            ,ngroup=NA,separate=FALSE,repel=20L,legend="bottomright"
                            ,title=NULL,verbose=FALSE,...) {
-   if (.skipPlot(TRUE))
-      return(NULL)
    ##~ method <- c('1'="ward.D",'2'="ward.D2",'3'="single",'4'="complete"
               ##~ ,'5'="average",'6'="mcquitty",'7'="median"
               ##~ ,'8'="centroid")[4] ## 3 4! 8 
    method <- match.arg(method)
    fun <- match.arg(fun)
-   if (isFALSE(legend))
-      legend <- NULL
    cutted <- 1.05
    da <- spatial_data(obj)
   # str(colnames(da))
@@ -81,6 +77,7 @@
       else
          bname <- ".undefined"
    }
+   bname <- bname[!is.na(bname)]
    if ((separate)&&(isCat)&&(!is.na(ngroup))&&(length(ngroup)!=length(bname))&&(ngroup[1]>2)) {
       tg <- table(xy4[[nameCat]])
      # print(tg)
@@ -102,7 +99,9 @@
      # print(ngroup)
    }
   # print(ngroup)
-   lutList <- lapply(if (separate) bname else ".+",function(sep) {
+  # print(xy4)
+   bname <- bname[bname %in% xy4[[3]]]
+   lutList <- lapply(if (separate) sample(bname) else ".+",function(sep) {
      # message(sQuote(sep),":")
      # print(ngroup[sep])
      # return(NULL)
@@ -126,7 +125,13 @@
             ng <- ngroup[sep]
          else
             ng <- ngroup
-         if (min(len[len>0])>s*0.75)
+         if (all(len==0)) {
+            cat("----------------------------------------\n")
+            str(xy5[,c("x","y")])
+            str(len)
+            ng <- NA
+         }
+         else if (min(len[len>0])>s*0.75)
             ng <- NA
         # print(data.frame(ng=ng,len=min(len[len>0])))
          if (!is.na(ng)) {
@@ -301,7 +306,7 @@
    if (length(bg)==1)
       bg <- rep(bg,length(ctInd))
    bg <- col2rgb(bg)/255
-   bg <- rgb(bg[1,],bg[2,],bg[3,],alpha=ifelse(separate,0.6,0.2))
+   bg <- rgb(bg[1,],bg[2,],bg[3,],alpha=0.2)
    s2 <- s/ifelse(label,2,1.5)/overlap
    for (i in seq(nrow(lut))) {
       x <- lut$.x[i] # <- mean(da2$x)
@@ -333,10 +338,9 @@
             lab <- lname[which(lab>0)]
          }
         # str(lab)
-         panel_annotation(x=x,y=y,label=as.character(lab),cex=cex^0.25,adj=c(0.5,0.53)
+         panel_annotation(x=x,y=y,label=as.character(lab),cex=cex,adj=c(0.5,0.53)
                         # ,fg="#FFFFFFA0",bg="#000000AF"
-                         ,buffer=2/scale ## commment it
-                        # ,buffer=2
+                        # ,buffer=2/scale
                          )
       }
       if (F) {
@@ -346,17 +350,9 @@
       }
    }
    if (!is.null(legend)) {
-      sc <- 96*getOption("ursaPngRetina")/getOption("ursaPngDpi")
-      if (T)
-         str(list('par(cex)'=par("cex"),cex=cex,scale=getOption("ursaPngScale")
-                 ,retina=getOption("ursaPngRetina"),ps=getOption("ursaPngPointsize")
-                 ,dpi=getOption("ursaPngDpi"),sc=sc))
       legend(legend,legend=bname,title=title
-            ,col=ctInd
-            ,cex=c(1,cex)[1]/par("cex")
-            ,pch=21
-            ,pt.lwd=ifelse(label,1,0)*2.4/par("cex")*sc
-            ,pt.cex=1.8/par("cex")
+            ,col=ctInd,cex=c(1,cex)[1]/par("cex")
+            ,pch=21,pt.lwd=ifelse(label,1,0)*2.4/par("cex"),pt.cex=1.8/par("cex")
             ,box.lwd=0.1,bg="#FFFFFFAF"
            # ,pt.bg=ursa_colortable(colorize(seq_along(ctInd),pal=ctInd,alpha="30"))
             ,pt.bg=if (label) bg else ctInd
@@ -439,13 +435,8 @@
    }
    if (is.null(col))
       col <- ursa_colortable(colorize(seq(nz)))
-   if (!is.null(border)) {
-      border1 <- rep_len(border,nz)
-      border2 <- rep_len("#FFFFFF",nz)
-      gscale <- colSums(col2rgb(col)*c(0.30,0.59,0.11))
-      if (all(gscale>160))
-         border2[] <- "#00000080"
-   }
+   if (!is.null(border)) 
+      border <- rep_len(border, nz)
    if (!is.null(lty)) 
       lty <- rep_len(lty, nz)
    if (!is.null(lwd)) 
@@ -490,18 +481,18 @@
    }
    for (i in seq_len(nz)) {
       n <- max(2,floor(edges*dz[i]))
-      if (!ball) { ## external circle
+      if (!ball) {
          P <- t2xy(seq.int(z[i],z[i+1],length.out=n),scale=c(1,0.65))
          ##~ polygon(c(P$x,0+x),c(P$y,0+y),border=border[i],col=col[i]
                 ##~ ,lty=lty[i],lwd=lwd[i]
                 ##~ )
-         polypath(P$x,P$y,border=border2[i],col=col1[i]
+         polypath(P$x,P$y,border=border[i],col=col1[i]
                  ,lty=lty[i],lwd=lwd[i]
                  ,rule=c("winding","evenodd")[2]
                  )
       }
       P <- t2xy(seq.int(z[i],z[i+1],length.out=n),scale=ifelse(ball,0.75,0.65))
-      polypath(P$x,P$y,border=border1[i],col=if (ball) col1[i] else col2[i]
+      polypath(P$x,P$y,border=border[i],col=if (ball) col1[i] else col2[i]
               ,lty=lty[i],lwd=lwd[i]
               ,rule=c("winding","evenodd")[2]
               )
diff --git a/R/xxx.polarmap.R b/R/xxx.polarmap.R
index 63bd1e3..97d7192 100644
--- a/R/xxx.polarmap.R
+++ b/R/xxx.polarmap.R
@@ -31,8 +31,14 @@
                       ,style=c("Arctic Connect","Arctic SDI")) {
    isSDI <- isTRUE(.lgrep("sdi",style[1])>0)
    isConnect <- !isSDI
-   if (missing(obj))
-      obj <- spatialize(c(0,50,360,50))
+   if (onlyBasemap <- missing(obj)) {
+      if (is.null(getOption("ursaSessionGrid")))
+         obj <- spatialize(c(0,50,360,50))
+      else {
+         obj <- spatial_transform(polygonize(session_bbox()),4326)
+      }
+      addFeature <- FALSE
+   }
    for (pkg in c("leaflet","leafem","leafpop")) 
       if (!requireNamespace(pkg,quietly=.isPackageInUse()))
          stop(paste("Package",sQuote(pkg),"is required for this operation"))
diff --git a/R/xxx.spatialize.R b/R/xxx.spatialize.R
index 5c97ed0..c31e5bc 100644
--- a/R/xxx.spatialize.R
+++ b/R/xxx.spatialize.R
@@ -1257,7 +1257,7 @@
          if (verbose)
             print(data.frame(lon0=lon_0,lat0=lat_0,lat_ts=lat_ts,row.names=proj))
          if (proj=="stere") {
-            t_srs <- paste("","+proj=stere"
+            t_srs <- paste("+proj=stere"
                           ,paste0("+lat_0=",lat_0)
                           ,paste0("+lat_ts=",lat_ts)
                           ,paste0("+lon_0=",lon_0)
@@ -1276,7 +1276,7 @@
                lon_0[lon_0>=(-25) && lon_0<(+50)] <- 10
                lon_0[lon_0>=(50) && lon_0<(+135)] <- 90
             }
-            t_srs <- paste("","+proj=laea"
+            t_srs <- paste("+proj=laea"
                           ,paste0("+lat_0=",lat_0)
                           ,paste0("+lon_0=",lon_0)
                           ,"+k=1","+x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs")
diff --git a/R/yyy.project.R b/R/yyy.project.R
index 49b543a..48e12ed 100644
--- a/R/yyy.project.R
+++ b/R/yyy.project.R
@@ -121,18 +121,32 @@
          if (!tryMatrix)
             a <- .try(res[-ind,] <- unclass(sf::st_transform(sf::st_sfc(
                               sf::st_multipoint(xy[-ind,]),crs=crs_s),crs_t)[[1]]))
-         else
+         else {
             a <- .try(res[-ind,] <- sf::sf_project(from=crs_s,to=crs_t,pts=xy[-ind,]))
+         }
       }
       else {
          if (!tryMatrix)
             a <- .try(res <- unclass(sf::st_transform(sf::st_sfc(
                                      sf::st_multipoint(xy),crs=crs_s),crs_t)[[1]]))
-         else
+         else {
+            if (T & !sf::sf_proj_network())
+               sf::sf_proj_network(url="",TRUE)
+            if (F & verbose) {
+               print(xy)
+               print(crs_s)
+               print(crs_t)
+               print(sf::st_crs(crs_t)$proj4string)
+            }
             a <- .try(res <- sf::sf_project(from=crs_s,to=crs_t,pts=xy))
+         }
       }
    }
    if (!a) {
+      if (verbose) {
+         print("WWW")
+         q()
+      }
       if (verbose)
          message("'rgdal' is used")
       if (!("rgdal" %in% loadedNamespaces()))
@@ -155,7 +169,7 @@
       if (!a) {
          if (verbose)
             message("patch with backing to 'proj4' via interim projection")
-         .epsg3411 <- paste("","+proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1"
+         .epsg3411 <- paste("+proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1"
                            ,"+x_0=0 +y_0=0 +a=6378273 +b=6356889.449 +units=m +no_defs")
          if (requireNamespace("proj4",quietly=.isPackageInUse())) {
             a <- .try(res <- proj4::project(xy=xy,proj=.epsg3411))
@@ -203,6 +217,8 @@
   # a <- try(as.integer(code),silent=TRUE)
    if ((is.character(code))&&(!nchar(code)))
       return(code)
+   if (is.matrix(code))
+      return(NULL)
    if (!.lgrep("\\D",code))
       p4epsg <- paste0(c("+init=epsg:","EPSG:")[2],code)
    else if (.lgrep("^epsg:\\d+",code))
diff --git a/R/yyy.tile.R b/R/yyy.tile.R
index f389471..b426700 100644
--- a/R/yyy.tile.R
+++ b/R/yyy.tile.R
@@ -137,8 +137,6 @@
                      ,paste0(osmCr,", \uA9 OpenTopoMap"))
    s$wiki <- c("https://maps.wikimedia.org/osm-intl/{z}/{x}/{y}{r}.png"
               ,paste0("Wikimedia | ",osmCr))
-   s$polarmap <- c("https://{abc}.tiles.arcticconnect.ca/osm_{l}/{z}/{x}/{y}.png"
-                  ,paste0("Map \uA9 ArcticConnect. Data ",osmCr))
    s$ArcticConnect <- c("https://{abc}.tiles.arcticconnect.ca/osm_{l}/{z}/{x}/{y}.png"
                   ,paste0("Map \uA9 ArcticConnect. Data ",osmCr))
    s$ArcticSDI <- c(paste0("http://basemap.arctic-sdi.org/mapcache?"
@@ -155,6 +153,9 @@
                    )
                   ,paste0("ArcticSDI")
                   ,"jpg")
+  # s$polarmap <- c("https://{abc}.tiles.arcticconnect.ca/osm_{l}/{z}/{x}/{y}.png"
+  #                ,paste0("Map \uA9 ArcticConnect. Data ",osmCr))
+   s$polarmap <- if (T) s$ArcticSDI else s$ArcticConnect
    s$google.h <- c(paste0("https://mt{0123}.google.com/vt/lyrs=h" ## roads only
                                ,"&x={x}&y={y}&z={z}&hl=",language),googleCr)
    s$google <- c(paste0("https://mt{0123}.google.com/vt/lyrs=m" ## standard roadmap
@@ -250,6 +251,21 @@
    if (server[1] %in% names(s))
       tile$name <- server
    indUrl <- .grep("^http(s)*://",style)
+   if (!length(indUrl)) {
+      patt <- dirname(style)
+      indUrl <- which(patt!=".")
+      if (!length(indUrl)) {
+         patt <- file.path(getOption("SAS_Planet_cache"),style)
+         indUrl <- which(dir.exists(patt))
+         if (length(indUrl)==1)
+            style[indUrl] <- patt[indUrl]
+      }
+      else
+         if (!dir.exists(patt[indUrl]))
+            indUrl <- integer()
+      if (length(style)==1)
+         style <- c(style,"from SAS Planet cache")
+   }
    if (!length(indUrl))
       return(names(s))
    indUrl <- indUrl[1]
@@ -329,7 +345,37 @@
    if (!.isPackageInUse())
       print(tile)
   # q()
-   fname <- .ursaCacheDownload(tile,mode="wb",cache=cache,quiet=!verbose)
+   if ((dir.exists(tile))&&(requireNamespace("RSQLite",quietly=!.isPackageInUse()))) {
+     ## try SAS Planet
+      list1 <- dir(path=file.path(tile,paste0("z",z+1L)),pattern="\\.sqlitedb$"
+                  ,recursive=TRUE,full.names=TRUE)
+      if (!length(list1))
+         return("missed in cache")
+      for (i in seq_along(list1)) {
+         mydb <- RSQLite::dbConnect(RSQLite::SQLite(),list1[i])
+         a <- RSQLite::dbGetQuery(mydb
+                ,paste0("SELECT * FROM t where t.x=",x," and ","t.y=",y))
+         RSQLite::dbDisconnect(mydb)
+         if (nrow(a))
+            break
+      }
+      if (!nrow(a))
+         return("missed in cache")
+      fname <- a$b[[1]]
+      d <- rawToChar(fname[1:8],multiple=TRUE) |> paste(collapse="")# |> as.character()
+      Encoding(d) <- "latin1" # |> print()
+      if (grepl("PNG",d))
+         fileext <- "png"
+      else if (grepl("JFIF",d))
+         fileext <- "jpg"
+      else
+         fileext <- "jpg"
+   }
+   else if (file.exists(tile)) {
+      fname <- tile
+   }
+   else
+      fname <- .ursaCacheDownload(tile,mode="wb",cache=cache,quiet=!verbose)
    if (inherits(fname,"try-error")) {
       return(fname)
      # message(a)
diff --git a/R/yyy.util.R b/R/yyy.util.R
index 9af9cbe..9dab362 100644
--- a/R/yyy.util.R
+++ b/R/yyy.util.R
@@ -322,9 +322,6 @@
       stop("NA values are not applicable")
    if (is.na(spherical))
       spherical <- isLonLat
-   ##~ print(summary(src))
-   ##~ print(summary(dst))
-   ##~ print(spherical)
    b1 <- .Cursa(C_dist2dist,x1=as.numeric(dst[,"x"]),y1=as.numeric(dst[,"y"])
                            ,x2=as.numeric(src[,"x"]),y2=as.numeric(src[,"y"])
                ,nxy=nrow(dst),ndf=nrow(src),positive=as.integer(positive)
diff --git a/inst/optional/sponsorship/wwf.png b/inst/optional/sponsorship/wwf.png
deleted file mode 100644
index c775e85..0000000
Binary files a/inst/optional/sponsorship/wwf.png and /dev/null differ
diff --git a/inst/requisite/coast-l.rds b/inst/requisite/coast-l.rds
index dc335bf..c1c85f6 100644
Binary files a/inst/requisite/coast-l.rds and b/inst/requisite/coast-l.rds differ
diff --git a/man/classConnection.Rd b/man/classConnection.Rd
index 158193e..bf232c3 100644
--- a/man/classConnection.Rd
+++ b/man/classConnection.Rd
@@ -76,7 +76,7 @@ Non-public function \code{.con.skeleton()} is used to generate the blank \code{u
 %%  ~~ If necessary, more details than the description above ~~
 }
 \value{
-\code{ursaConnection} is a list. The most of names have a relation to specification of \href{https://www.l3harrisgeospatial.com/docs/ENVIHeaderFiles.html}{ENVI Header Files}. Items:
+\code{ursaConnection} is a list. The most of names have a relation to specification of \href{https://www.nv5geospatialsoftware.com/docs/ENVIHeaderFiles.html}{ENVI Header Files}. Items:
 \item{driver}{Character. Keyword of supported image formats. Allowed "ENVI" or "GDAL".}
 \item{samples}{Integer. Number of image columns (samples)}
 \item{lines}{Integer. Number of image rows (lines)}
diff --git a/man/conn.create_any.Rd b/man/conn.create_any.Rd
index 70cd6ba..3295242 100644
--- a/man/conn.create_any.Rd
+++ b/man/conn.create_any.Rd
@@ -71,7 +71,7 @@ Keywords:
 \item{12, "UInt16" = Integer: 16-bit unsigned integer}
 \item{13, "UInt32" = Long: 32-bit unsigned integer}
 }
-Specification \url{https://www.l3harrisgeospatial.com/docs/ENVIHeaderFiles.html} is used.
+Specification \url{https://www.nv5geospatialsoftware.com/docs/ENVIHeaderFiles.html} is used.
 }
 \item{\code{byteorder} - numeric (integer). Byte order.}
 \item{\code{bands} - numeric( integer). Number of bands/layers}
diff --git a/man/conn.write_envi.Rd b/man/conn.write_envi.Rd
index e2a266f..a288803 100644
--- a/man/conn.write_envi.Rd
+++ b/man/conn.write_envi.Rd
@@ -28,7 +28,7 @@ close(a)
 }
 }
 \value{
-Integer code of ENVI data type. See values of the \dQuote{data type} field in description of the \href{https://www.l3harrisgeospatial.com/docs/ENVIHeaderFiles.html}{ENVI Header Format}.
+Integer code of ENVI data type. See values of the \dQuote{data type} field in description of the \href{https://www.nv5geospatialsoftware.com/docs/ENVIHeaderFiles.html}{ENVI Header Format}.
 }
 %%~ \references{
 %%~ %% ~put references to the literature/web site here ~
diff --git a/man/conn.write_gdal.Rd b/man/conn.write_gdal.Rd
index 489973f..b747846 100644
--- a/man/conn.write_gdal.Rd
+++ b/man/conn.write_gdal.Rd
@@ -36,7 +36,7 @@ close(a)
 }
 }
 \value{
-Integer code of ENVI data type. See values of the \dQuote{data type} field in description of the \href{https://www.l3harrisgeospatial.com/docs/ENVIHeaderFiles.html}{ENVI Header Format}.
+Integer code of ENVI data type. See values of the \dQuote{data type} field in description of the \href{https://www.nv5geospatialsoftware.com/docs/ENVIHeaderFiles.html}{ENVI Header Format}.
 }
 %%~ \references{
 %%~ %% ~put references to the literature/web site here ~
diff --git a/man/spatial_engine.Rd b/man/spatial_engine.Rd
index f1cc632..9ce9db7 100644
--- a/man/spatial_engine.Rd
+++ b/man/spatial_engine.Rd
@@ -271,7 +271,7 @@ Nikita Platonov \email{platonov@sevin.ru}
 %% ~Make other sections like Warning with \section{Warning }{....} ~
 
 \section{Acknowledgements }{
-   The great improvement for development of functions for manipulation with spatial objects has been reached during work in series of projects (2015-2018) for design of marine protected areas in the Arctic, which were supported by WWF Russia.
+   The great improvement for development of functions for manipulation with spatial objects has been reached during work in series of projects (2015-2018) for design of marine protected areas in the Arctic.
 }
 
 \examples{
diff --git a/src/ursa.c b/src/ursa.c
index 5103085..647d256 100644
--- a/src/ursa.c
+++ b/src/ursa.c
@@ -3688,7 +3688,7 @@ void dist2dist(double *x1,double *y1,double *x2,double *y2
          lon1[i]=x1[i]*PI/180.0;
          slat1[i]=sin(y1[i]*PI/180.0);
          clat1[i]=cos(y1[i]*PI/180.0);
-        // Rprintf("x1=%.2f cos=%.2f sin=%.2f\n",x1[j],clat1[j],slat1[j]);
+        // Rprintf("x1=%.2f cos=%.2f sin=%.2f\n",x1[i],clat1[i],slat1[i]);
       }
    }
    for (j=0;j<n2;j++)
@@ -3699,7 +3699,7 @@ void dist2dist(double *x1,double *y1,double *x2,double *y2
       {
          if (spherical) {
            // d=0;
-           // Rprintf("k=%d lon1=%.2f\n",k,lon1[i]);
+           // Rprintf("k=%d lon1=%.2f (%.2f)\n",k,lon1[i],lon1[i]*180.0/PI);
             d=6371000.0*acos(slat1[i]*slat2[j]+clat1[i]*clat2[j]*cos(fabs(lon1[i]-lon2[j])));
            // d=acos(cos(lat1[i]));
          }
@@ -3718,10 +3718,16 @@ void dist2dist(double *x1,double *y1,double *x2,double *y2
             i2=i;
          }
       }
-      dist[j]=sqrt(minD);
+      if (spherical)
+         dist[j]=minD;
+      else
+         dist[j]=sqrt(minD);
       ind[j]=i2;
       if (verbose)
          progressBar(j,n2,"");
+      if ((0)&&(spherical)) {
+         Rprintf("ind=%d dist=%.3f\n",j,minD);
+      }
    }
    free(clat1);
    free(clat2);