Fix ionosphere error by adding back in xarray. Update configs/burst…

… DB files. (#249)

* Fix ionosphere error by adding back in `xarray`

* Update configs: overrides, burst_db files, runconfigs pointing to these

* fix test path

* add missing symlink

conda-env.yml

@@ -19,7 +19,7 @@ dependencies:
   # Corrections
   - pysolid  # For solid earth tide corrections
   # - pyaps3
-  # - xarray
+  - xarray
   # Unwrappers:
   - snaphu>=0.4.1
   - isce3-cpu>=0.16.0 # For baseline

configs/algorithm_parameters_historical.yaml

@@ -1,11 +1,11 @@
 # JSON file containing frame-specific algorithm parameters to override the defaults passed
 #   in the `algorithm_parameters.yaml`.
 #   Type: string | null.
-algorithm_parameters_overrides_json: opera-disp-s1-algorithm-parameters-overrides-2024-12-23.json
+algorithm_parameters_overrides_json: opera-disp-s1-algorithm-parameters-overrides-2025-01-09.json
 ps_options:
   # Amplitude dispersion threshold to consider a pixel a PS.
   #   Type: number.
-  amp_dispersion_threshold: 0.20
+  amp_dispersion_threshold: 0.2
 phase_linking:
   # Size of the ministack for sequential estimator.
   #   Type: integer.
@@ -84,7 +84,7 @@ unwrap_options:
   unwrap_method: snaphu
   # Number of interferograms to unwrap in parallel.
   #   Type: integer.
-  n_parallel_jobs: 2
+  n_parallel_jobs: 4
   # Set wrapped phase/correlation to 0 where mask is 0 before unwrapping. .
   #   Type: boolean.
   zero_where_masked: false
@@ -116,7 +116,7 @@ unwrap_options:
       - 400
     # Number of tiles to unwrap in parallel for each interferogram.
     #   Type: integer.
-    n_parallel_tiles: 5
+    n_parallel_tiles: 7
     # Initialization method for SNAPHU.
     #   Type: string.
     #   Options: ['mcf', 'mst'].
@@ -156,6 +156,10 @@ unwrap_options:
     #   `similarity_threshold` *or* above the temporal coherence threshold is chosen.
     #   Type: number.
     similarity_threshold: 0.5
+    # After running spurt, interpolate the values that were masked during unwrapping (which are
+    #   otherwise left as nan).
+    #   Type: boolean.
+    run_ambiguity_interpolation: true
     general_settings:
       # Tile up data spatially.
       #   Type: boolean.
@@ -210,7 +214,6 @@ unwrap_options:
       min_overlap_points: 25
       # Currently, only 'dirichlet' is supported.
       #   Type: dirichlet.
-      #   Options: ['dirichlet'].
       method: dirichlet
       # Method used to estimate bulk offset between tiles.
       #   Type: string.

configs/opera-disp-s1-algorithm-parameters-overrides-2025-01-09.json

@@ -0,0 +1,34 @@
+{
+  "23210": {
+    "phase_linking": {
+      "compressed_slc_plan": "last_per_ministack"
+    },
+    "unwrap_options": {
+      "unwrap_method": "spurt"
+    }
+  },
+  "23211": {
+    "phase_linking": {
+      "compressed_slc_plan": "last_per_ministack"
+    },
+    "unwrap_options": {
+      "unwrap_method": "spurt"
+    }
+  },
+  "33038": {
+    "phase_linking": {
+      "compressed_slc_plan": "last_per_ministack"
+    },
+    "unwrap_options": {
+      "unwrap_method": "spurt"
+    }
+  },
+  "33039": {
+    "phase_linking": {
+      "compressed_slc_plan": "last_per_ministack"
+    },
+    "unwrap_options": {
+      "unwrap_method": "spurt"
+    }
+  }
+}

configs/runconfig_historical.yaml

@@ -32,10 +32,10 @@ dynamic_ancillary_file_group:
 static_ancillary_file_group:
   # JSON file containing the mapping from frame_id to frame/burst information.
   #   Type: string | null.
-  frame_to_burst_json: static_ancillary_files/opera-s1-disp-0.5.0-frame-to-burst.json.zip
+  frame_to_burst_json: static_ancillary_files/opera-s1-disp-0.7.0-frame-to-burst.json.zip
   # JSON file containing list of reference date changes for each frame.
   #   Type: string | null.
-  reference_date_database_json: static_ancillary_files/opera-disp-s1-reference-dates-minimal-2024-09-19.json
+  reference_date_database_json: static_ancillary_files/opera-disp-s1-reference-dates-2025-01-15-minimal.json
 primary_executable:
   # Product type of the PGE.
   #   Type: string.
@@ -64,11 +64,11 @@ worker_settings:
   # Number of threads to use per worker. This sets the OMP_NUM_THREADS environment variable in
   #   each python process.
   #   Type: integer.
-  threads_per_worker: 16
+  threads_per_worker: 2
   # If processing separate spatial bursts, number of bursts to run in parallel for wrapped-
   #   phase-estimation.
   #   Type: integer.
-  n_parallel_bursts: 4
+  n_parallel_bursts: 9
   # Size (rows, columns) of blocks of data to load at a time.
   #   Type: array.
   block_shape: