Estoy extrayendo la información de la temperatura utilizando el getData de la función raster paquete.
bio1 <- getData('worldclim', var="bio", res = 0.5, lon=1, lat=-77)
bio2 <- getData('worldclim', var="bio", res = 0.5, lon=7, lat=-70)La variable bio var=bio significa que se trata de la extracción de 19 capas de diferentes variables bioclimáticas en las coordenadas lon/lat proporcionadas.
El problema es que si obtengo la media de estas capas para ambos bio1 y bio2 aunque sean de coordenadas lat/lon diferentes, son las mismas. Y estoy usando la resolución más fina res = 0.5 :
> cellStats(bio1,stat = "mean")
bio1_46 bio2_46 bio3_46 bio4_46 bio5_46 bio6_46 bio7_46 bio8_46 bio9_46 bio10_46 bio11_46
159.38445 146.23495 54.29901 4262.23057 292.49362 23.63659 268.85703 178.42543 132.43352 211.26331 102.22708
bio12_46 bio13_46 bio14_46 bio15_46 bio16_46 bio17_46 bio18_46 bio19_46
416.39054 61.81748 13.82139 45.91476 165.43986 48.51057 129.29831 72.57548
> cellStats(bio2,stat = "mean")
bio1_46 bio2_46 bio3_46 bio4_46 bio5_46 bio6_46 bio7_46 bio8_46 bio9_46 bio10_46 bio11_46
159.38445 146.23495 54.29901 4262.23057 292.49362 23.63659 268.85703 178.42543 132.43352 211.26331 102.22708
bio12_46 bio13_46 bio14_46 bio15_46 bio16_46 bio17_46 bio18_46 bio19_46
416.39054 61.81748 13.82139 45.91476 165.43986 48.51057 129.29831 72.57548 Del manual he entendido que cellStats obtendría la media de las celdas de cada una de mis capas. Pero tal vez está haciendo algo diferente?
