Определить положение клика на растре в буклете, в R

Я готовлю большой латинский NetCDFraster над Rleaflet карта с использованиемshinydashboard, Когда я нажимаю на карту, появляется всплывающее окно, в котором отображаются строка, столбец, положение по ширине и значение выбранной растровой точки. (См. Воспроизводимый код ниже)

Проблема в том, что я испытываю сдвиг в растре, если растр достаточно большой. Например, здесь я щелкнул точку, которая должна иметь значение, но в результате мы определили точку, указанную выше.

Я считаю, что это связано с тем, что растр, используемыйleaflet проецируется, в то время как исходные данные, которые я использую для определения точек, являются широтно-долгосрочными, поскольку точка щелчка возвращается как широковещательнаяleaflet, Я не могу использовать проецируемый файл (depth) так как его единицы в метрах, а не в градусах! Даже если бы я попытался перепроецировать эти метры в градусы, я получил сдвиг.

Вот основной исполняемый пример кода:

#Libraries
library(leaflet)
library(raster)
library(shinydashboard)
library(shiny)

#Input data
download.file("https://www.dropbox.com/s/y9ekjod2pt09rvv/test.nc?dl=0", destfile="test.nc")
inputFile    = "test.nc"
inputVarName = "Depth"
lldepth <- raster(inputFile, varname=inputVarName)
lldepth[Which(lldepth<=0, cells=T)] <- NA #Set all cells <=0 to NA
ext <- extent(lldepth)
resol <- res(lldepth)
projection(lldepth) <- "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"

#Project for leaflet
depth <- projectRasterForLeaflet(lldepth)

#Prepare UI
sbwidth=200
sidebar <- dashboardSidebar(width=sbwidth)
body <- dashboardBody(
          box( #https://stackoverflow.com/questions/31278938/how-can-i-make-my-shiny-leafletoutput-have-height-100-while-inside-a-navbarpa
            div(class="outer",width = NULL, solidHeader = TRUE, tags$style(type = "text/css", paste0(".outer {position: fixed; top: 50px; left: ", sbwidth, "px; right: 0; bottom: 0px; overflow: hidden; padding: 0}")),
            leafletOutput("map", width = "100%", height = "100%")
            )
          )
        )
ui <- dashboardPage(
  dashboardHeader(title = "A title"),
  sidebar,
  body
)
#
#Server instance
server <- function(input, output, session) {
  output$map <- renderLeaflet({#Set extent
    leaflet()  %>%
      fitBounds(ext[1], ext[3], ext[2], ext[4])
  })

  observe({#Observer to show Popups on click
    click <- input$map_click
    if (!is.null(click)) {
      showpos(x=click$lng, y=click$lat)
    }
  })

  showpos <- function(x=NULL, y=NULL) {#Show popup on clicks
    #Translate Lat-Lon to cell number using the unprojected raster
    #This is because the projected raster is not in degrees, we cannot use it!
    cell <- cellFromXY(lldepth, c(x, y))
    if (!is.na(cell)) {#If the click is inside the raster...
      xy <- xyFromCell(lldepth, cell) #Get the center of the cell
      x <- xy[1]
      y <- xy[2]
      #Get row and column, to print later
      rc <- rowColFromCell(lldepth, cell)
      #Get value of the given cell
      val = depth[cell]
      content <- paste0("X=",rc[2],
                        "; Y=",rc[1],
                        "; Lon=", round(x, 5),
                        "; Lat=", round(y, 5),
                        "; Depth=", round(val, 1), " m")
      proxy <- leafletProxy("map")
      #add Popup
      proxy %>% clearPopups() %>% addPopups(x, y, popup = content)
      #add rectangles for testing
      proxy %>% clearShapes() %>% addRectangles(x-resol[1]/2, y-resol[2]/2, x+resol[1]/2, y+resol[2]/2)
    }
  }

  #Plot the raster
  leafletProxy("map") %>%
    addRasterImage(depth, opacity=0.8, project=FALSE, group="Example", layerId="Example", colors=colorNumeric(terrain.colors(10), values(depth), na.color = "black"))
}


print(shinyApp(ui, server))

Как правильно определить точки, если растр большой?

РЕДАКТИРОВАТЬ: Я также хотел предоставить некоторые дополнительные ссылки на (возможно) связанную документацию или вопросы:

Растровое изображение, кажется, смещено с помощью листовки для RR для перенаправления листовки при нажатии на растровое изображениеhttps://gis.stackexchange.com/questions/183918/is-it-possible-to-use-a-rasterclick-event-within-an-interactive-leaflet-mapмаркер мыши, событие клика в буклете R для блеска
 AF730 мая 2016 г., 14:09
Ах! Я думаю, что сделал это, перепроектировав данные о кликах снова. Я буду экспериментировать, а затем отправлю ответ.

Ответы на вопрос(1)

Решение Вопроса

что могу перепроецировать обратно позицию X-Y (lon-lat), заданнуюinput$map_click, В этом случае я предполагал, что входная проекция будет Lon-Lat, но я думаю, что это не обязательно должно быть. Это просто нужно иметь единицы Lat-Lon.

#Set projections
inputProj <- "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"
leafletProj <- "+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137 +towgs84=0,0,0,0,0,0,0 +units=m [email protected] +wktext +no_defs"
#Note that for some reason "[email protected] +wktext" is very important
  #as hinted to by other questions and answers linked in my question.
xy <- SpatialPoints(data.frame(x,y))
proj4string(xy) <- inputProj
xy <- as.data.frame(spTransform(xy, leafletProj))
#Get the cell number from the newly transformed metric X and Y.
cell <- cellFromXY(depth, c(xy$x, xy$y))

#At this point, you can also retrace back the center of the cell in
  #leaflet coordinates, starting from the cell number!
xy <- SpatialPoints(xyFromCell(depth, cell))
proj4string(xy) <- leafletProj
xy <- as.data.frame(spTransform(xy, inputProj))
#Here XY will again be in lat-lon, if you projection says so,
  #indicating the center of the clicked cell

Ваш ответ на вопрос