Proyecciones para sendas de autorroscantes de eclipse solar

Question

La mayoría de las proyecciones de la ruta de acceso de la totalidad para el 2017 eclipse solar aspecto

(a partir de http://xjubier.free.fr/en/site_pages/solar_eclipses/TSE_2017_GoogleMapFull.html).

Pero también hay mapas en línea, utilizando diferentes proyecciones tal que la trayectoria es una línea recta o en su mayor parte en línea recta, como en

(a partir de http://www.eclipse-maps.com)

No puedo encontrar cualquier CRS para el segundo. El ideal de CRS para mí produciría un camino recto (sólo tiene que ser sobre todo recto por la parte de la trayectoria a través de los estados unidos) y, además, sería el camino hacia arriba y hacia abajo con "occidente" lateral en la parte superior. Estoy usando QGIS, pero estoy susceptibles de cambiar a otro paquete de código abierto, si es necesario.

@MichaelZeiler dijo (y me han confirmado) que para este particular, eclipse y esta parte de la carrera (el Estados Unidos), un estándar de proyección sucede a trabajar lo suficientemente bien como para producir un satisfactoriamente recta lo suficientemente umbral ruta de acceso: EPSG:102003 = USA_Contiguous_Albers_Equal_Area_conic

+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=37.5 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +unidades=m +no_defs

Ahora estoy tratando de averiguar cómo transformar eso en un CRS que produce mapas que están rotados 67.560 grados a la derecha en relación a la anterior.

En otras palabras, quiero cambiar el CRS, por lo que me sale

Observe cómo todas las etiquetas son horizontales. Esta es una ventaja de cambiar el CRS, ya que se opuso a la Lona y otras soluciones de etiquetas de trabajo normalmente para todas las capas nuevas. Todavía no entiendo cómo cambiar el CRS después de la lectura de un número de páginas.

---

Xavier Jubier dijo que una buena aproximación de línea recta también puede obtenerse con EPSG:102009 = North_America_Lambert_Conformal_conic

+proj=lcc +lat_1=20 +lat_2=60 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +unidades=m +no_defs

He aprendido después de muchas búsquedas que no son oscuros proyecciones que le dará más exacta y astronomía general y soluciones en torno a cualquier persona que quiera tomar este enfoque más difícil va a querer saber la palabra clave quincuncial.

Answer 1

Intente Oblicuo Mecator, por ejemplo, con

"+proj=omerc +lonc=-100 +gamma=90 +alfa=10 +lat_0=30 +ellps=WGS84"

He aquí un R ejemplo, sólo para mostrar lo que he intentado. No he estudiado este tema en profundidad.

library(rgdal)

prj <- "+proj=omerc +lonc=-100 +gamma=90 +alpha=10 +lat_0=30 +ellps=WGS84"

library(USAboundaries)
us <- us_states()

library(graticule)
g <- graticule(seq(165, 300, by = 15), seq(-30, 85, by = 15), proj = prj)
gl <- graticule_labels(seq(165, 300, by = 15), seq(-30, 85, by = 15), proj = prj, xline = 285, yline = 20)
x <- spTransform(us, prj)
plot(x)
plot(g, add = TRUE)
text(gl, lab = parse(text = gl$lab))

Hay algunos más ejemplos generales de girar con un mapa real de la proyección en R aquí, aunque es sólo norma PROJ.4 por lo que la costumbre GDAL, QGIS pila funcionará bien.

https://rpubs.com/geospacedman/rotatespatial