Как я могу использовать и понимать crsClasses листовки (R-)?

Изменить после ввода ниже и, например. здесь: буклет должен найти компромисс между "понятным для неспециалистов по ГИС" и "правильным обращением с ГИС". Насколько я понимаю сейчас, буклет ориентирован на простоту использования, поэтому прогнозы не могут быть реализованы.

Оригинальный вопрос: Моя цель — правильно отображать «настоящую» плотность, то есть в событиях/м^2. На данный момент я борюсь с проекцией листовки. Что мне нужно сделать, чтобы правильно отображать все crsClasses? Мои данные генерируют однородные случайные точки на шере.

Думаю, я пропустил преобразование данных?

library(leaflet)
n <- 10000
z <- 2*runif (n) - 1
phi <- 2*pi*runif (n) - pi
x <- sin(phi)*sqrt(1 - z^2)
y <- cos(phi)*sqrt(1 - z^2)
theta <- acos(z)
# lat / lng
lat <- theta*180/pi - 90
lng <- phi*180/pi

# Working fine:
# plot3D::polygon3D(x, y, z)
rgl::plot3d(x, y, z)
# rgl::plot3d(sin(phi)*sin(theta), cos(phi)*sin(theta), cos(theta)) # the same...

# EPSG:3857, also known as "Google Mercator" or "Web Mercator", the first in the following list
crsClasses <- c("L.CRS.EPSG3857", "L.CRS.EPSG4326", "L.CRS.EPSG3395", "L.CRS.Simple", "L.Proj.CRS")
epsg3857 <- leafletCRS(crsClass = crsClasses[1])
leaflet(options = leafletOptions(
  crs = epsg3857, worldCopyJump = FALSE)) %>%
  addTiles() %>%
  addProviderTiles(providers$OpenStreetMap, group = "OSM") %>%
  addProviderTiles(providers$Stamen.TonerLite, group = "Toner Lite") %>%
  addLayersControl(baseGroups = c("OSM", "Toner Lite")) %>%
  addCircleMarkers(lat = lat,
                   lng = lng,
                   radius = 1, 
                   stroke = FALSE,
                   fillOpacity = 0.7)

Использование crsClasses[1]... выглядит нормально, но я не понимаю, почему это соответствует worldCopyJump = FALSE.

Использование crsClasses[2]... выглядит странно

crsClasses[3...6] тоже выглядят странно. Очень сложно найти хорошие учебники для новичка в ГИС.

Использованная литература:

• рстудия
• геокомп.

Комментарий после ввода Ивана Санчеса: Верна ли следующая стратегия?

• С помощью crsClass L.CRS.EPSG4326 карта изменилась, но данные остались без изменений. Таким образом, я также должен преобразовать данные. Это означает, что я должен использовать new <- sp::spTransform(x, CRSobj, ...). Это выглядит как
  • Мне нужно вычислить код EPSG, здесь вы мне уже сказали (и вы могли догадаться), что это EPSG:4326. Есть ли системный способ получить такой результат?
  • Мне нужно преобразовать (широта, длина) в правильный объект. Из ??spTransform: "x... объект для преобразования" и "CRSobj... объект класса CRS". Помощь приводит к sp::CRS(projargs, doCheckCRSArgs=TRUE). Но потом я застреваю...
  • Вызовите новую широту/долготу в addCircleMarkers().
• Как вы указали, я должен искать равновеликую проекцию для плотности
  • Из здесь я бы поискал азимутальный равновеликий по Ламберту (LAEA), как они рекомендуют. Я полагаю, что из здесь я беру EPSG:3575(??) для Европы, но, возможно, "WGS 84/Северный полюс LAEA Europe" неверен?
  • Я потерялся: ... найти (или построить) базовые карты для этой проекции. Я думаю, это означает для карты и данных?

Можете ли вы дать ориентир в этих джунглях или хотя бы в одном из вышеперечисленных случаев? Или у вас есть хорошая ссылка, чтобы узнать, как это сделать правильно?

Из документов листовки:

crs ... Coordinate Reference System to use. Don't change this if you're not sure what it means.

Но, к сожалению, они не говорят вам, где искать, если вам нужно изменить.