bundle_adjust_plots, corr_and_dem_plots, WIP reorg

uw-cryo · Sep 19, 2023 · 69cb28e · 69cb28e
1 parent 6e4f93b
commit 69cb28e
Show file tree

Hide file tree

Showing 2 changed files with 143 additions and 209 deletions.
diff --git a/bundle_adjust_plots.ipynb b/bundle_adjust_plots.ipynb
@@ -69,6 +69,16 @@
     "topdir = '/Users/dshean/scr/BlackSky/Utqiagvik_20220425_stereo/BSG-STEREO-102-20220425-215106-22900060-stereo_reorder_20230119/'"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "efcef52d-4f58-4a39-b0f3-e6b86657b335",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "topdir = '/Users/dshean/scr/hma_glof_samples_for_asp_plot'"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -130,6 +140,16 @@
     "#ba_brefix=\"ba/ba_ip2k_nadirpinhole\""
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "0661c87b-6c1a-4761-b967-13b485051854",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "ba_prefix = 'ba/ba_csm'"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -141,6 +161,17 @@
     "#map_crs = 'EPSG:32604'"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "4959cb3a-16f0-4bb5-a350-caff336756d0",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "map_crs = 'EPSG:32645'\n",
+    "#map_crs = 'EPSG:32604'"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -153,6 +184,16 @@
     "#refdem = 'COP30_lzw-adj_proj_ref.tif'"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "09dad1c1-be21-4de7-a792-deae88438c45",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "refdem = 'GongbatongshaTsho_COP30_lzw-adj_proj.tif'"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -435,38 +476,12 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "geodiff_init = read_geodiff(geodiff_init_csv)\n",
-    "geodiff_final = read_geodiff(geodiff_final_csv)"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "062012ed-68db-410d-b61f-3cf221971140",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "geodiff_init.describe()"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "e2b41b9a-1b91-48ab-8bd6-5ee95e095162",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "geodiff_final.describe()"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "0bdbd967-a0e1-41dd-9343-a83a1d2aac13",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "resid_plot(geodiff_init, geodiff_final, col='diff', clim=(-15, 15), cmap='RdYlBu')"
+    "if os.path.exists(geodiff_init_csv) and os.path.exists(geodiff_final_csv):\n",
+    "    geodiff_init = read_geodiff(geodiff_init_csv)\n",
+    "    geodiff_final = read_geodiff(geodiff_final_csv)\n",
+    "    geodiff_init.describe()   \n",
+    "    geodiff_final.describe()\n",
+    "    resid_plot(geodiff_init, geodiff_final, col='diff', clim=(-15, 15), cmap='RdYlBu')"
    ]
   },
   {
@@ -504,31 +519,14 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "id": "508c1374-e1f3-4593-b18e-f56160808769",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "mapproj_match_offset = read_mapproj_match_offset(mapproj_match_offset_txt)"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "a3fc0938-ca16-46cc-8464-fdc04cb590e9",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "mapproj_match_offset.describe()"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "361b42b9-9cc0-480d-b8ce-c01fdea062e7",
+   "id": "f33d7a44-743f-4577-abad-cb36ab700cd1",
    "metadata": {},
    "outputs": [],
    "source": [
-    "col='mapproj_ip_dist_meters'"
+    "if os.path.exists(mapproj_match_offset_txt):\n",
+    "    mapproj_match_offset = read_mapproj_match_offset(mapproj_match_offset_txt)\n",
+    "    mapproj_match_offset.describe()\n",
+    "    col='mapproj_ip_dist_meters'\n"
    ]
   },
   {
@@ -1353,7 +1351,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.10.8"
+   "version": "3.10.12"
   }
  },
  "nbformat": 4,