enable downloading data from ASF

Author: ehavazli

prep/NISAR_prep.ipynb

@@ -81,7 +81,8 @@
     "from rasterio.warp import reproject, Resampling, transform_bounds\n",
     "import h5py\n",
     "from pyproj import Transformer\n",
-    "import tile_mate"
+    "import tile_mate\n",
+    "import asf_search as asf"
    ]
   },
   {
@@ -102,12 +103,12 @@
    "outputs": [],
    "source": [
     "# === Basic Configuration ===\n",
-    "site = \"OregonA34\"  # Cal/Val location ID\n",
+    "site = \"LosAngelesA34\"  # Cal/Val location ID\n",
     "requirement = \"Transient\"  # Options: 'Secular', 'Coseismic', 'Transient'\n",
     "dataset = \"NISAR\"  # Dataset type: 'NISAR_sample', 'ARIA_S1', 'ARIA_S1_new'\n",
     "\n",
     "\n",
-    "rundate = \"20260227\"  # Date of this Cal/Val run\n",
+    "rundate = \"20260312\"  # Date of this Cal/Val run\n",
     "version = \"1\"         # Version of this Cal/Val run\n",
     "custom_sites = \"/home/jovyan/my_nisar_sites.txt\"  # Path to custom site metadata\n",
     "\n",
@@ -217,7 +218,7 @@
     "        print(f\"Earthdata credentials found for user: {earthdata_user}\")\n",
     "\n",
     "if not earthdata:\n",
-    "    print(\"\\nNEEDED to Download ARIA GUNWs\")\n",
+    "    print(\"\\nNEEDED to Download ARIA/NISAR GUNWs\")\n",
     "    print(\"Create account at: https://urs.earthdata.nasa.gov/\")\n",
     "    earthdata_user = input(\"Earthdata username: \").strip()\n",
     "    earthdata_password = input(\"Earthdata password: \").strip()\n",
@@ -266,9 +267,40 @@
     "<hr>\n",
     "\n",
     "<a id='prep_ifg'></a>\n",
-    "# 1. Download and Prepare Interferograms\n",
-    "\n",
-    "In this initial processing step, all the necessary Level-2 unwrapped interferogram products are gathered and organized for analysis with MintPy."
+    "# 1. Download and Prepare Interferograms"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "95f2582c-8291-4cdc-b4b6-c3e1b2e66c37",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "print(f\"CalVal site: {site}\")\n",
+    "\n",
+    "# Extract site-specific metadata\n",
+    "bbox = site_info.get('download_region')\n",
+    "startdate = site_info.get('download_start_date')\n",
+    "enddate = site_info.get('download_end_date')\n",
+    "track = site_info.get('nisar_track')\n",
+    "\n",
+    "s, n, w, e = map(float, bbox.strip(\"'\").split())\n",
+    "bbox_wkt = f\"POLYGON(({w} {s}, {e} {s}, {e} {n}, {w} {n}, {w} {s}))\"\n",
+    "startdate_fmt = f\"{enddate[:4]}-{enddate[4:6]}-{enddate[6:]}\"\n",
+    "enddate_fmt = f\"{enddate[:4]}-{enddate[4:6]}-{enddate[6:]}\"\n",
+    "data_list = asf.search(dataset=asf.DATASET.NISAR, processingLevel=asf.PRODUCT_TYPE.GUNW, relativeOrbit=int(track), intersectsWith=bbox_wkt)\n",
+    "prod_list = data_list.find_urls(extension='.h5', pattern=r'.*\\d.h5')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "460d1964-b369-48c0-b915-6ca8e7a717ee",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "asf.download_urls(prod_list, path=gunw_dir)"
    ]
   },
   {
@@ -505,15 +537,24 @@
    "outputs": [],
    "source": [
     "command = f'prep_nisar.py -i \"{gunw_dir}/*.h5\" -d {out_dem} -m {out_mask} -o {mintpy_dir}'\n",
-    "process = subprocess.run(command, shell=True)\n",
+    "process = subprocess.run(command, shell=True, check=True)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "d64d0033-49da-4b8f-9a44-8cd66ecd589d",
+   "metadata": {},
+   "outputs": [],
+   "source": [
     "print('Mintpy input files:')\n",
     "[x for x in os.listdir(mintpy_dir + '/inputs') if x.endswith('.h5')]"
    ]
   },
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "d64d0033-49da-4b8f-9a44-8cd66ecd589d",
+   "id": "8cd8b317-b54b-40b8-aa79-80bef0298c98",
    "metadata": {},
    "outputs": [],
    "source": []