Hola, estoy tratando de obtener la media de precipitación para cada cuenca hidrográfica en España. Para ello estoy siguiendo una metodología que he visto en un paper. Tengo los datos de precipitación en un grid (0,5ºx0,5º), que he convertido a shape e intersectado con el shape de los polígonos que delimitan las cuencas.
Una de las indicaciones de la metodología dice que se ha de tener en cuenta que la superficie de cada celda en el raster varía en función de su latitud. Por lo que el peso relativo de cada celda sería igual a la superficie de dicha celda que se encuentra dentro de la cuenca correspondiente multiplicado por el Coseno de la latitud de la celda.
Para eso he calculado la coordenada Y del centroide (latitud) de cada celda en grados decimales, la he pasado a radianes (unidad de entrada para la función coseno) y luego en otra columna he calculado el coseno de dicha latitud.
Mi problema es que he obtenido resultados raros, como que ángulos muy similares den cosenos muy diferentes (como se puede ver en el ejemplo de abajo). Para mi sorpresa al hacer los cálculos en Excel los resultados son exactamente iguales.
Latitud (grados) Cos (latitud)
| 42,98440000000 |
0,54207059701 |
| 42,39580000000 |
-0,01570017839 |
Les agradezco de antemano su ayuda. Por más vueltas que le doy, no llego a entender que puede haber pasado.
Un saludo