¿Cuál es la mejor proyección para la rasterización de datos aleatorios de latitud y longitud en el Atlántico norte?

Question

Tengo un conjunto de datos en el sistema de coordenadas geográficas WGS84 que me gustaría interpolar a lo largo de una cuadrícula usando Python.

Parece gdal_grid no está disponible a través de los enlaces Python, así que voy a utilizar otro 2D método de interpolación después de rasterizar los datos.

Para trazar los datos, tengo la intención de volver proyecto es el uso de pyproj, a continuación, escribir en un ráster utilizando el gdal enlaces python. La proyección UTM en el sistema de coordenadas parece que sería la mejor proyección, pero mis datos se extiende por el Atlántico Norte (es decir, múltiples zonas UTM).

Actualización:
Estoy abierto a sugerencias sobre el método de interpolación, pero yo había planeado utilizar una función inversa de la distancia de la interpolación (como se muestra aquí). Estos son presencia ausencia de avistamientos de datos de la nave de crucero de los temas, que no necesariamente sigue las líneas de latitud. Abarcan desde las latitudes a mediados de los 70 a mediados de los años 50. Gracias!

Qué proyección sería mejor utilizar?

Answer 1

Avistamientos de constituir una (no aleatoria) de la muestra de algún proceso o de la población. En consecuencia, la interpolación (especialmente) IDW no es una buena idea: se resuelve un problema diferente por completo.

Considere la posibilidad de hacer un mapa de densidad. Al hacerlo, es probable que sea mejor a favor de la igualdad-área de proyecciones sobre la conformación de las proyecciones (debido a los cambios de área de sesgo de la densidad, mientras que los no-conformality no, aunque es posible cambiar el aspecto cualitativo del mapa).

Hay muchos igual en el área de las proyecciones de el mundo. (Por desgracia, popular proyecciones, incluyendo Mercator Transversal de Mercator, y Estereográfica, no están entre ellos.) La mayoría de ellos son cilíndricos. Esto es bueno en una forma-las líneas de latitud y longitud de formar una cuadrícula ortogonal en el mapa -, pero es malo, cerca de los polos, debido a la escala de distorsión es todo concentrado en un (vertical) dirección donde puede llegar a ser extrema.

Estas consideraciones sugieren un compromiso: apenas alrededor de cualquier proyección adecuado para un largish área (como un gran país), cuando se desplazó hacia el centro del Atlántico Norte, es probable que haga un buen trabajo, simplemente debido a su relativa global de la distorsión debe ser limitado. Por ejemplo, el Albers equal-área de proyección cónica es muy popular para la colinde con Estados Unidos. Cambiando su central meridian a -30 grados, su referencia paralela a 55 grados de latitud norte, y sus paralelos a los 45 y los 65 grados, re-puntos en el Atlántico Norte:

El Tissot indicatrices mostrar cuán poco la distorsión de la realidad se produce dentro de esta región: la escala de los errores están a menos de 1,5% en cualquier dirección a lo largo; eso es mucho más de precisión que se necesita para la interpolación o estimaciones de densidad de kernel.

En contraste, a ver lo que realmente se ha logrado por esta opción en comparación con otros que podrían sugerir que ellos mismos--considerar la proyección de Mercator:

Las enormes discrepancias entre los círculos interior y exterior de la Tissot indicatrices testificar a los grandes cambios en la escala de sur a norte en este mapa. En efecto, cualquier interpolación o el cálculo de la densidad basados en estas coordenadas de mapas haría uso de los barrios cerca de la parte superior que son aproximadamente la mitad del tamaño de los barrios cerca de la parte inferior (en la realidad, no en el mapa). Que la distorsión apreciable; no hay sentido de introducir en los cálculos cuando es suficientemente simple como para elegir una mejor proyección.

Answer 2

Google Mercator (EPSG:3857) debería ser la mejor solución. Abarca todo el mundo de hasta 85° Norte/sur.

Así que si sus datos no se encuentra cerca del ártico, se debe trabajar.

Y usted puede utilizar Openlayers plugin para el mapa de fondo, si quieres vizualize los datos en QGIS.

Pero tenga en cuenta que las longitudes no son reales metros.

---

EDITAR

Si no te gusta eso, usted puede tratar de lambert conformal cónica. EPSG:3034 se define para pan-europea de asignación. Paralelos 35° y 65° de latitud Norte; dando este diseño del mapa:

Usted puede cambiar lon_0 de 10° Este a algo más occidental, si no te gusta la perspectiva.