add MintPY config file generation

Author: ehavazli

prep/NISAR_prep.ipynb

@@ -34,7 +34,8 @@
     "[**1. Download and Prepare Interferograms**](#prep_ifg)\n",
     "- [1.1.  Get Data List](#prep_data)\n",
     "- [1.2.  Download DEM](#prep_dem)\n",
-    "- [1.3.  Load Data into MintPy](#prep_load_data)\n",
+    "- [1.3.  Create MintPy configuration file](#create_config)\n",
+    "- [1.4.  Load Data into MintPy](#prep_load_data)\n",
     "\n",
     "<hr/>"
    ]
@@ -85,17 +86,17 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "8b464e26-d7a9-4834-bc77-da3968e90f57",
+   "id": "18978e7b-b33d-4f6d-a38d-a87044275edb",
    "metadata": {},
    "outputs": [],
    "source": [
     "# === Basic Configuration ===\n",
-    "site = \"NISAR_ALOS1_Hawaii_0380\"  # Cal/Val location ID\n",
+    "site = \"Erta_Ale\"  # Cal/Val location ID\n",
     "requirement = \"Transient\"  # Options: 'Secular', 'Coseismic', 'Transient'\n",
     "dataset = \"NISAR_sample\"  # Dataset type: 'NISAR_sample', 'ARIA_S1', 'ARIA_S1_new'\n",
     "\n",
     "\n",
-    "rundate = \"20251008\"  # Date of this Cal/Val run\n",
+    "rundate = \"20260205\"  # Date of this Cal/Val run\n",
     "version = \"1\"         # Version of this Cal/Val run\n",
     "custom_sites = \"/home/jovyan/my_sites.txt\"  # Path to custom site metadata\n",
     "\n",
@@ -326,37 +327,87 @@
     "\n",
     "out_dem = os.path.join(dem_dir, f\"{site}_elevation.dem\")\n",
     "\n",
-    "cmd = [\n",
-    "    \"sardem\",\n",
-    "    \"--bbox\", *bbox,\n",
-    "    \"--output\", out_dem,\n",
-    "    \"--output-format\", \"ENVI\"\n",
-    "]\n",
+    "if os.path.isfile(out_dem):\n",
+    "    print(f\"{out_dem} exists, skip downloading\")\n",
+    "else:\n",
+    "    cmd = [\n",
+    "        \"sardem\",\n",
+    "        \"--bbox\", *bbox,\n",
+    "        \"--output\", out_dem,\n",
+    "        \"--output-format\", \"ENVI\"\n",
+    "    ]\n",
+    "    \n",
+    "    print(\"Running:\", \" \".join(cmd))\n",
+    "    try:\n",
+    "        res = subprocess.run(cmd, check=True, capture_output=True, text=True)\n",
+    "        if res.stdout:\n",
+    "            print(res.stdout)\n",
+    "        print(\"DEM files:\", os.listdir(dem_dir))\n",
+    "    except subprocess.CalledProcessError as e:\n",
+    "        print(\"sardem failed:\")\n",
+    "        print(e.stdout or \"\")\n",
+    "        print(e.stderr or \"\")\n",
+    "        raise"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "54c045a9-65ce-47ac-9fd3-d32d560667f8",
+   "metadata": {},
+   "source": [
+    "### 1.3. Create MintPy configuration file <a id='create_config'></a>"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "c3817649-4966-4ab9-a213-27c6dddbfd00",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "os.chdir(mintpy_dir)\n",
     "\n",
-    "print(\"Running:\", \" \".join(cmd))\n",
-    "try:\n",
-    "    res = subprocess.run(cmd, check=True, capture_output=True, text=True)\n",
-    "    if res.stdout:\n",
-    "        print(res.stdout)\n",
-    "    print(\"DEM files:\", os.listdir(dem_dir))\n",
-    "except subprocess.CalledProcessError as e:\n",
-    "    print(\"sardem failed:\")\n",
-    "    print(e.stdout or \"\")\n",
-    "    print(e.stderr or \"\")\n",
-    "    raise"
+    "# Build config as a dictionary first\n",
+    "config_file_content = {\n",
+    "    \"mintpy.load.processor\": \"nisar\",\n",
+    "    \"mintpy.compute.cluster\": \"local\",\n",
+    "    \"mintpy.compute.numWorker\": \"auto\",\n",
+    "    \"mintpy.topographicResidual.pixelwiseGeometry\": \"no\",\n",
+    "    \"mintpy.troposphericDelay.method\": \"no\",\n",
+    "    \"mintpy.topographicResidual\": \"no\",\n",
+    "    \"mintpy.network.tempBaseMax\": site_info.get('tempBaseMax'),\n",
+    "    \"mintpy.network.startDate\": site_info.get('download_start_date'),\n",
+    "    \"mintpy.network.endDate\": site_info.get('download_end_date'),\n",
+    "    \"mintpy.velocity.startDate\": site_info.get('download_start_date'),\n",
+    "    \"mintpy.velocity.endDate\": site_info.get('download_end_date'),\n",
+    "    \"mintpy.reference.lalo\": site_info.get('reference_lalo'),\n",
+    "    \"mintpy.network.excludeDate12\": site_info.get('ifgExcludePair'),\n",
+    "    \"mintpy.network.excludeDate\" : site_info.get('ifgExcludeDate'),\n",
+    "    \"mintpy.network.excludeIfgIndex\" : site_info.get('ifgExcludeIndex'),\n",
+    "}\n",
+    "\n",
+    "# Write config dictionary to text file\n",
+    "with open(config_file, \"w\") as f:\n",
+    "    f.writelines(f\"{k} = {v}\\n\" for k, v in config_file_content.items())\n",
+    "\n",
+    "# Confirm output\n",
+    "print(f\"MintPy config file written to:\\n    {config_file}\\n\")\n",
+    "with open(config_file, \"r\") as f:\n",
+    "    print(f.read())\n"
    ]
   },
   {
    "cell_type": "markdown",
    "id": "0de4a38e-60f3-421b-949f-aff1145e5881",
    "metadata": {},
    "source": [
-    "### 1.3. Load Data into MintPy Cubes <a id='prep_load_data'></a>\n",
+    "### 1.4. Load Data into MintPy Cubes <a id='prep_load_data'></a>\n",
     "\n",
     "The output of this step is an \"inputs\" directory in 'calval_directory/calval_location/MintPy/\" containing three HDF5 files:\n",
-    "- ifgramStack.h5: This file contains 6 dataset cubes (e.g. unwrapped phase, coherence, connected components etc.) and multiple metadata\n",
-    "- geometryGeo.h5: This file contains geometrical datasets (e.g., incidence/azimuth angle, masks, etc.)\n",
-    "- ionStack.h5   : This file contains pairwise ionosphere data present in the NISAR data "
+    "- ifgramStack.h5: Contains 6 dataset cubes (e.g. unwrapped phase, coherence, connected components etc.) and multiple metadata.\n",
+    "- geometryGeo.h5: Contains geometrical datasets (e.g., incidence/azimuth angle, masks, etc.).\n",
+    "- ionStack.h5   : Contains pairwise ionosphere data present in the NISAR L2 GUNW.\n",
+    "- tropoStack.h5 : Contains 3D interpolated total tropospheric delay calculated from NISAR radar grid tropospheric delay layers."
    ]
   },
   {
@@ -389,7 +440,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.13.7"
+   "version": "3.14.3"
   }
  },
  "nbformat": 4,