NEOCP_filters repo start

NEOCP_filters/MPC.config

@@ -0,0 +1,160 @@
+repeatable
+norms
+raw
+noid
+Int
+NEO
+MC
+Hun
+Pho
+MB1
+Pal
+Han
+MB2
+MB3
+Hil
+JTr
+JFC
+obserr033=0.93
+obserr104=0.5
+obserr106=0.39
+obserr130=0.33
+obserr160=0.38
+obserr194=0.25
+obserr203=0.25
+obserr204=0.27
+obserr215=0.43
+obserr221=0.6
+obserr246=0.46
+obserr291=0.44
+obserr381=1.27
+obserr474=0.35
+obserr511=1.66
+obserr568=0.29
+obserr587=0.61
+obserr595=0.64
+obserr598=1.23
+obserr654=0.2
+obserr658=0.68
+obserr673=0.11
+obserr691=0.34
+obserr695=0.37
+obserr703=0.54
+obserr704=0.62
+obserr711=0.45
+obserr718=0.5
+obserr734=0.39
+obserr807=0.3
+obserr850=0.68
+obserr858=0.37
+obserr900=0.43
+obserr950=0.21
+obserr958=0.67
+obserrA17=0.3
+obserrA50=0.57
+obserrA77=0.57
+obserrB37=0.82
+obserrB49=0.6
+obserrB72=0.6
+obserrC23=0.4
+obserrC40=0.24
+obserrC43=0.73
+obserrC51=0.73
+obserrC57=0.74
+obserrC77=0.41
+obserrC82=0.6
+obserrC95=0.51
+obserrC97=0.29
+obserrD29=0.39
+obserrE12=0.55
+obserrE23=0.37
+obserrE94=0.68
+obserrF51=0.13
+obserrF52=0.12
+obserrF65=0.24
+obserrG34=0.4
+obserrG37=0.23
+obserrG40=0.23
+obserrG84=0.25
+obserrG96=0.29
+obserrH01=0.15
+obserrH06=0.49
+obserrH15=0.49
+obserrH21=0.26
+obserrH36=0.48
+obserrI22=0.6
+obserrI41=0.4
+obserrI47=0.34
+obserrI52=0.38
+obserrI81=0.6
+obserrI93=0.32
+obserrJ04=0.28
+obserrJ43=0.32
+obserrJ57=0.6
+obserrJ69=0.9
+obserrJ75=0.41
+obserrJ95=0.43
+obserrK19=0.49
+obserrK38=0.24
+obserrK51=0.2
+obserrK61=0.48
+obserrK63=0.48
+obserrK74=0.44
+obserrK78=0.6
+obserrK83=0.39
+obserrK87=0.44
+obserrK88=0.25
+obserrK91=0.19
+obserrK93=0.21
+obserrK96=0.6
+obserrL01=0.44
+obserrL02=0.6
+obserrL04=0.25
+obserrL06=0.6
+obserrL18=0.54
+obserrL27=0.6
+obserrL29=0.4
+obserrL34=0.46
+obserrL39=0.65
+obserrL47=0.55
+obserrL51=0.61
+obserrL87=1.0
+obserrL81=1.22
+obserrM22=0.4
+obserrN88=0.56
+obserrO18=0.20
+obserrO75=0.75
+obserrP93=0.53
+obserrQ21=0.54
+obserrQ62=1.02
+obserrR57=0.41
+obserrT05=0.35
+obserrT08=0.4
+obserrT09=0.3
+obserrT12=0.07
+obserrT14=0.07
+obserrU55=0.56
+obserrU68=0.68
+obserrU69=0.32
+obserrU74=0.40
+obserrV00=0.17
+obserrV01=0.48
+obserrV06=0.3
+obserrW16=0.22
+obserrW33=0.35
+obserrW34=0.46
+obserrW68=0.4
+obserrW84=0.12
+obserrW88=0.22
+obserrW94=0.6
+obserrW95=0.51
+obserrX31=0.58
+obserrX74=0.49
+obserrY00=0.5
+obserrZ33=0.6
+obserrZ41=0.6
+obserrZ43=0.25
+obserrZ45=0.6
+obserrZ80=0.55
+obserrZ84=0.35
+obserrZ99=0.41

NEOCP_filters/README.md

@@ -0,0 +1,20 @@
+Set of tools that analyze digest2 output and distinguishes NEOs from non-NEOs.
+
+The analysis is described in Cloete, R. and Veres, P.: Near-Earth Object Discovery Enhancement with the Machine Learning Methods
+
+Data:
+
+neocp.obs - obs80 (MPC1992) format of NEOPC data between 2019 - 2024. The designations were altered: last 6 characters are from trkid, first character is either '0' (NEO) or '1' (non-NEO).
+digest_data_19-24.csv - digest2 values for 2019 - 2023 NEOCP data
+digest_data_24.csv - digest2 values for 2024 NEOCP data
+MPC.config - configuration file for digest2
+optimal_thresholds.json - derived filter threshold (limit = 0)
+
+Tools:
+
+find_filter.py - takes digest2 output (e.g. digest_data_19-24.csv) in a csv format, creates JSON threshold model "optimal_thresholds.json" for the following tool
+example usage: python3 find_filter.py digest_data_19-24.csv
+
+neocp_filter.py - selects assumed non-NEOs from the input file (digest2 output, e.g. digest_data_24.csv) based on the JSON model
+example usage: neocp_filter.py digest_data_24.csv optimal_thresholds.json
+