# -*- coding: utf-8 -*-
"""
Created on Thu Aug 27 07:58:12 2026

@author: chuffard
"""



import numpy as np
import matplotlib.pyplot as plt
from datetime import datetime
import pandas as pd
import sys
import os.path
import warnings
import glob
from pathlib import Path
import math


MARS2025atn = r"W:/Analyses/2025 MARS deployment/SES_out_atn_Red_5only.csv"
MARS2025atn = pd.read_csv(MARS2025atn)

# Extract the pattern into a new column
MARS2025atn['file_name'] = MARS2025atn['file_name'].astype(str).str.strip()

MARS2025atn[['date', 'time']] = MARS2025atn['file_name'].str.extract(r'(\d{6})-(\d{6})')

# Convert to actual datetime objects for analysis
MARS2025atn['timestamp'] = pd.to_datetime(MARS2025atn['date'] + MARS2025atn['time'], format='%y%m%d%H%M%S')

MARS2025atn.columns



MARS2025atn["DateTime"] = pd.to_datetime(MARS2025atn['timestamp'])
#convert atn to POC flux

#from Estapa et al. 2025 optical sediment trap calibration. also accounting for image and collect area since those aren't in our og atn calculation

SEStraparea = 0.5;
#imgarea = length(find(isfinite(bgr(:))))./(0.07194.^2)./(1000000.^2);  % pix / (pix/um^2) / (um^2/m^2)#this was in Meg's example code but it's not for SES

#This is adapted to pyton from matlab code Meg gave us
SESimgarea = math.pi*(25/1000/2)**2  # diameter of lower end of funnel is 25 mm, compute area in m^2
SEScolltime = 5/(60*24) #in days
MARS2025atn['F_atn']= MARS2025atn['atn']/SEScolltime*SESimgarea/SEStraparea


#These are from SES at Station M in 2025 paper
MARS2025atn['POCmodel_mg']= 10**3.4 * MARS2025atn['F_atn']**0.94
MARS2025atn['POCmodel_g']= MARS2025atn['POCmodel_mg']/1000

MARS2025atn['POCmodel_g'].mean()#(0.2753247346815628)
MARS2025atn['POCmodel_g'].plot()
plt.show()

MARS2025atn.to_csv('//Thalassa/ProjectLibrary/902305_Event_Detection_with_SES/Analyses/2025 MARS deployment/MARS2025atn_POC_red.csv')