Use None to encode NaNs in metadata and enable JSON dumps of results

Author: aMarcireau

README.md

@@ -599,7 +599,7 @@ git submodule update --recursive
 Format the code:
 
 ```sh
-clang-format -i astrometry_extension/astrometry_extension.c
+clang-format -i astrometry_extension/astrometry_extension.c astrometry_extension/astrometry_extension_utilities.h
 ```
 
 Build a local version:

astrometry/__init__.py

@@ -2,6 +2,7 @@
 import astrometry_extension
 import dataclasses
 import enum
+import json
 import math
 import operator
 import pathlib
@@ -154,6 +155,12 @@ def has_match(self) -> bool:
     def best_match(self) -> Match:
         return self.matches[0]
 
+    def to_json(self):
+        solution_as_dict = dataclasses.asdict(self)
+        for match in solution_as_dict["matches"]:
+            match["index_path"] = str(match["index_path"])
+        return json.dumps(solution_as_dict)
+
 
 class Solver(astrometry_extension.Solver):
     def __init__(self, index_files: list[pathlib.Path]):

astrometry_extension/astrometry_extension_utilities.h

@@ -5,6 +5,7 @@
 #define PY_SSIZE_T_CLEAN
 #include <Python.h>
 #include <libgen.h>
+#include <math.h>
 #include <string.h>
 #include <structmember.h>
 
@@ -255,6 +256,13 @@ static void error_callback(
     context->save = PyEval_SaveThread();
 }
 
+static PyObject* double_to_python_object(double value) {
+    if (isnan(value) || isinf(value)) {
+        Py_RETURN_NONE;
+    }
+    return PyFloat_FromDouble(value);
+}
+
 static PyObject*
 tagalong_to_python_object(startree_t* tree, int column_index, const char* column_name, int star_id, PyObject* logging) {
     int size = startree_get_tagalong_column_array_size(tree, column_index);
@@ -276,20 +284,20 @@ tagalong_to_python_object(startree_t* tree, int column_index, const char* column
             if (row_size > 1) {
                 result = PyTuple_New(row_size);
                 for (int index = 0; index < row_size; ++index) {
-                    PyTuple_SET_ITEM(result, index, PyFloat_FromDouble(((double*)row)[index]));
+                    PyTuple_SET_ITEM(result, index, double_to_python_object(((double*)row)[index]));
                 }
             } else {
-                result = PyFloat_FromDouble(*(double*)row);
+                result = double_to_python_object(*(double*)row);
             }
             break;
         case TFITS_BIN_TYPE_E: // float
             if (row_size > 1) {
                 result = PyTuple_New(row_size);
                 for (int index = 0; index < row_size; ++index) {
-                    PyTuple_SET_ITEM(result, index, PyFloat_FromDouble(((float*)row)[index]));
+                    PyTuple_SET_ITEM(result, index, double_to_python_object(((float*)row)[index]));
                 }
             } else {
-                result = PyFloat_FromDouble(*(float*)row);
+                result = double_to_python_object(*(float*)row);
             }
             break;
         case TFITS_BIN_TYPE_A: // char

setup.py

@@ -249,7 +249,7 @@ def run(self):
 
 setuptools.setup(
     name="astrometry",
-    version="4.0.0",
+    version="4.1.0",
     url="https://github.com/neuromorphicsystems/astrometry",
     author="ICNS, Alexandre Marcireau",
     author_email="[email protected]",