Sun diff max and min occur at the same place 2022 and 2023. Is that correct?

[ScientiaFelis]

The sebms_minmax_sunhour() produce the same sites for 2022 and 2023.

This is highly improbable Is it a bug or a actual coincidence?

• [ ] It is a bug, fixed
• [ ] It was a mere coincidence, checked from SMHI original data.

[larspett]

Ping @ScientiaFelis

Jag har undersökt detta en hel del nu och det verkar som att det skapas NA (290 st) vid st_join, troligen pga att det är points och att minsta decimal innebär en miss när man gör join med intersect. Vi behöver skapa ett index (id-nummer) i rastret som hämtas från SMHI, dvs på det råa, fullständiga rastret, och sen joina med left join på id istället för st_join på position. NA-mönstren ser ut at ligga bakom krumelurerna vi såg.

Sedan är det fortsatt märkligt att max och min diff hamnar på i princip exakt samma positioner 2023 som 2022. Det känns som att positionerna tas baserat på data från 2022 snarare än 2023. Möjligen körs diffar mot artificiellt låga eller höga värden pga NA-problematiken. En bra början blir att köra joins mot id för att få bort NA.

[larspett]

också märkligt att max/min funktionen så anges 'total_sunH' utan '' i default, försöker man skicka in diff-namnet utan '' blir det error, skickar man med '' börjar den "söka igenom data" men egentligen skapar den nog en ny fil

[ScientiaFelis]

Ok, bra! Ja, st_join känns kanske lite osäkert, speciellt för punkter. Det var kanske inte det bästa designvalet jag gjort. Men det borde gå att fixa med ett skapat id.

Sen kan vi ju kolla igen om det fixar övriga problem som du sa.

[larspett]

En grej till - när man gör per month och per day, lägg till en räknare, t ex layer_count som räknar så det är 6 lager med i månadsdata och 183 i dagsdata