Hello everybody, I use SIrad for estimating daily extraterrestrial solar radiation for Sardinia using dem, in this way
DAY = "2019-02-13" dem$Ra <- extrat(dayOfYear(DAY), radians(dem$lat_250))$ExtraTerrestrialSolarRadiationDaily now I need to map hourly values for specific hours, but when I calculate $ExtraTerrestrialSolarRadiationHorly I obtain the whole set of 24 values for each grid point, ma question is .. how can I extract and plot only specific hourly values, for example maximum values (12 o' clock)? Thank you Regards, Michele Michele Fiori ARPAS - Agenzia Regionale per la Protezione dell'Ambiente della Sardegna Dipartimento Meteoclimatico - Servizio Meteorologico Agrometeorologico ed Ecosistemi Viale Porto Torres 119 - 07100 Sassari Tel + 39 079 258617 Fax +39 262681 <http://www.sardegnaambiente.it/arpas> www.sardegnaambiente.it/arpas [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo