From b0fe53801697212f49604b672bd7491ea90eb415 Mon Sep 17 00:00:00 2001 From: Michael Rudolf <rudolf@geo.tu-darmstadt.de> Date: Fri, 16 Dec 2022 14:52:05 +0100 Subject: [PATCH] Changed all paths to relative paths --- notebooks/Snippets.ipynb | 112 +++++++++--------- .../220423_template_aoi.ipynb | 41 +++---- .../area_of_interest/hohermeisner_aoi.ipynb | 41 +++---- notebooks/area_of_interest/kassel_aoi.ipynb | 41 +++---- .../area_of_interest/offenbach_aoi.ipynb | 41 +++---- .../area_of_interest/rheinmain_aoi.ipynb | 41 +++---- .../area_of_interest/werrakali_aoi.ipynb | 41 +++---- .../DataPrep_Diff_QuarryHM2119.ipynb | 31 ++--- .../DataPrep_Diff_QuarryTemplate_HM1914.ipynb | 35 +++--- .../DataPrep_Diff_Template_HM1914.ipynb | 51 ++++---- .../data_prep/DataPrep_KasselStadt.ipynb | 49 ++++---- .../DataPrep_Processing_220411.ipynb | 35 +++--- notebooks/data_prep/DataPrep_RheinMain.ipynb | 53 +++++---- ...p_TestRegionTemplate_DataPrepipynb_1.ipynb | 53 +++++---- .../DataPrep_TestRegion_Template.ipynb | 35 +++--- notebooks/data_prep/DataPrep_WerraKali.ipynb | 53 +++++---- notebooks/get_building.ipynb | 53 +++++---- notebooks/get_quarry.ipynb | 25 ++-- notebooks/maps_quarry.ipynb | 87 +++++++------- notebooks/reclass_mask_complete.ipynb | 45 +++---- 20 files changed, 491 insertions(+), 472 deletions(-) diff --git a/notebooks/Snippets.ipynb b/notebooks/Snippets.ipynb index 8946e6d..97a2ac3 100644 --- a/notebooks/Snippets.ipynb +++ b/notebooks/Snippets.ipynb @@ -16,12 +16,12 @@ "metadata": {}, "outputs": [], "source": [ - "buffer_complete = r\"F:\\newData\\hesse_shp\\boundingBox_hesse.shp\"\n", - "buffer_complete_dis = r\"F:\\newData\\hesse_shp\\bbmeged_hesse.shp\"\n", + "buffer_complete = os.path.join(base_path, r\"newData\\hesse_shp\\boundingBox_hesse.shp\")\n", + "buffer_complete_dis = os.path.join(base_path, r\"newData\\hesse_shp\\bbmeged_hesse.shp\")\n", "\n", - "inp_fishnet = r\"F:\\newData\\hesse_shp\\hesse_fishnet_50x50.shp\"\n", + "inp_fishnet = os.path.join(base_path, r\"newData\\hesse_shp\\hesse_fishnet_50x50.shp\")\n", "\n", - "clipped = r\"F:\\newData\\hesse_shp\\hesse_raster_50x50.shp\"" + "clipped = os.path.join(base_path, r\"newData\\hesse_shp\\hesse_raster_50x50.shp\")" ] }, { @@ -58,7 +58,7 @@ "with arcpy.da.UpdateCursor(clipped, [\"*\"]) as cursor:\n", " i = 0\n", " for row in cursor:\n", - " filename = r\"F:\\newData\\hesse_shp\\grid\\hesse_raster_50x50_{}.shp\".format(i)\n", + " filename = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\hesse_raster_50x50_{}.shp\").format(i)\n", " arcpy.management.CopyFeatures(clipped, filename, '', None, None, None)\n", " i = i + 1\n", "\n" @@ -70,7 +70,7 @@ "metadata": {}, "outputs": [], "source": [ - "workspace = r\"F:\\newData\\hesse_shp\\grid\"\n", + "workspace = os.path.join(base_path, r\"newData\\hesse_shp\\grid\")\n", "\n", "datalist = []\n", "\n", @@ -114,15 +114,15 @@ "outputs": [], "source": [ "# Input Files\n", - "inp_vect_water = r\"F:\\newData\\hesse_shp\\gis_osm_water_a_free_1.shp\"\n", - "inp_vect_build = r\"F:\\newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\"\n", - "inp_vect_landuse = r\"F:\\newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\"\n", - "inp_vect_roads = r\"F:\\newData\\hesse_shp\\gis_osm_roads_free_1.shp\"\n", - "inp_vect_railways = r\"F:\\newData\\hesse_shp\\gis_osm_railways_free_1.shp\"\n", + "inp_vect_water = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_water_a_free_1.shp\")\n", + "inp_vect_build = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\")\n", + "inp_vect_landuse = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\")\n", + "inp_vect_roads = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_roads_free_1.shp\")\n", + "inp_vect_railways = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_railways_free_1.shp\")\n", "\n", "# import rasters to clip\n", - "diff_1914 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "diff_2119 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", + "diff_1914 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "diff_2119 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", "\n", "# Clip analysis\n", "#arcpy.analysis.Clip(inp_vect_water, area_of_interest, out_data_water) \n", @@ -152,14 +152,14 @@ "source": [ "i = 0\n", "for file in datalist:\n", - " out_file_w = r\"F:\\newData\\hesse_shp\\grid\\osm_\" + r\"{}\\water_{}.shp\".format(i, i)\n", - " out_file_b = r\"F:\\newData\\hesse_shp\\grid\\osm_\" + r\"{}\\build_{}.shp\".format(i, i)\n", - " out_file_l = r\"F:\\newData\\hesse_shp\\grid\\osm_\" + r\"{}\\landu_{}.shp\".format(i, i)\n", - " out_file_ro = r\"F:\\newData\\hesse_shp\\grid\\osm_\" + r\"{}\\roads_{}.shp\".format(i, i)\n", - " out_file_ra = r\"F:\\newData\\hesse_shp\\grid\\osm_\" + r\"{}\\railw_{}.shp\".format(i, i)\n", + " out_file_w = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_\" + r\"{}\\water_{}.shp\").format(i, i)\n", + " out_file_b = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_\" + r\"{}\\build_{}.shp\").format(i, i)\n", + " out_file_l = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_\" + r\"{}\\landu_{}.shp\").format(i, i)\n", + " out_file_ro = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_\" + r\"{}\\roads_{}.shp\").format(i, i)\n", + " out_file_ra = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_\" + r\"{}\\railw_{}.shp\").format(i, i)\n", " \n", - " out_diff_19_14 = r\"F:\\newData\\hesse_shp\\grid\\osm_\" + r\"{}\\diff_19_14_{}.tif\".format(i, i)\n", - " out_diff_21_19 = r\"F:\\newData\\hesse_shp\\grid\\osm_\" + r\"{}\\diff_21_19_{}.tif\".format(i, i)\n", + " out_diff_19_14 = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_\" + r\"{}\\diff_19_14_{}.tif\").format(i, i)\n", + " out_diff_21_19 = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_\" + r\"{}\\diff_21_19_{}.tif\").format(i, i)\n", " \n", " #print(out_file_w, out_file_b, out_file_l,out_file_ro, out_file_ra)\n", " \n", @@ -185,7 +185,7 @@ "metadata": {}, "outputs": [], "source": [ - "workspace = r\"F:\\newData\\hesse_shp\\grid\"\n", + "workspace = os.path.join(base_path, r\"newData\\hesse_shp\\grid\")\n", "for dirpath, dirnames, filenames in arcpy.da.Walk(workspace, datatype=\"FeatureClass\"):\n", " for filename in filenames:\n", " filepath = r\"{0}\\{1}\".format(dirpath, filename)\n", @@ -200,7 +200,7 @@ "metadata": {}, "outputs": [], "source": [ - "workspace = r\"F:\\newData\\hesse_shp\\grid\"\n", + "workspace = os.path.join(base_path, r\"newData\\hesse_shp\\grid\")\n", "\n", "for dirpath, dirnames, filenames in arcpy.da.Walk(workspace, datatype=\"FeatureClass\"):\n", " for filename in filenames:\n", @@ -217,7 +217,7 @@ "outputs": [], "source": [ "# Quarrys extrahieren! -> nicht getestet wegen Landuse zu früh \"Dissolved\"\n", - "workspace = r\"F:\\newData\\hesse_shp\\grid\"\n", + "workspace = os.path.join(base_path, r\"newData\\hesse_shp\\grid\")\n", "field_name = \"fclass\"\n", "\n", "for dirpath, dirnames, filenames in arcpy.da.Walk(workspace):\n", @@ -243,9 +243,9 @@ "source": [ "# Quarries from only_quarry.shp \n", "\n", - "workspace = r\"F:\\newData\\hesse_shp\\grid\"\n", + "workspace = os.path.join(base_path, r\"newData\\hesse_shp\\grid\")\n", "distance = \"10 Meters\"\n", - "inp_vect_quarry = r\"F:\\newData\\hesse_osm_buffer\\only_quarry.shp\"\n", + "inp_vect_quarry = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\only_quarry.shp\")\n", "\n", "i = 0\n", "for dirpath, dirnames, filenames in arcpy.da.Walk(workspace):\n", @@ -287,21 +287,21 @@ "source": [ "# mergen der benutzen params\n", "\n", - "workspace = r\"F:\\newData\\hesse_shp\\grid\"\n", + "workspace = os.path.join(base_path, r\"newData\\hesse_shp\\grid\")\n", "#i = 0\n", "for i in range(0, 18): \n", - " inp_build = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_build_{}\".format(i, i)\n", - " inp_quarry = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_quarry_{}\".format(i, i)\n", - " inp_water = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_water_{}\".format(i, i)\n", + " inp_build = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_build_{}\").format(i, i)\n", + " inp_quarry = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_quarry_{}\").format(i, i)\n", + " inp_water = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_water_{}\").format(i, i)\n", " \n", - " inp_roads = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_roads_{}\".format(i, i)\n", - " inp_railways = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_railw_{}\".format(i, i)\n", + " inp_roads = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_roads_{}\").format(i, i)\n", + " inp_railways = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_railw_{}\").format(i, i)\n", " \n", - " out_merged_3 = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_merged3_{}\".format(i, i)\n", - " out_dis_merged_3 = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_merged3_{}\".format(i, i)\n", + " out_merged_3 = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_merged3_{}\").format(i, i)\n", + " out_dis_merged_3 = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_merged3_{}\").format(i, i)\n", " \n", - " out_merged_5 = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_merged5_{}\".format(i, i)\n", - " out_dis_merged_5 = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_merged5_{}\".format(i, i)\n", + " out_merged_5 = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_merged5_{}\").format(i, i)\n", + " out_dis_merged_5 = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_merged5_{}\").format(i, i)\n", " \n", " arcpy.management.Merge([inp_build,inp_quarry,inp_water], out_merged_3)\n", " arcpy.management.Dissolve(out_merged_3, out_dis_merged_3)\n", @@ -320,9 +320,9 @@ "source": [ "# b) erst mergen dann maske abziehen, dann von gesamter diff maskieren\n", "\n", - "workspace = r\"F:\\newData\\hesse_shp\\grid\"\n", - "merged_hesse3 = r\"F:\\newData\\hesse_shp\\c3merg_hes.shp\"\n", - "dis_merged_hesse3 = r\"F:\\newData\\hesse_shp\\c3dismerg_hes.shp\"\n", + "workspace = os.path.join(base_path, r\"newData\\hesse_shp\\grid\")\n", + "merged_hesse3 = os.path.join(base_path, r\"newData\\hesse_shp\\c3merg_hes.shp\")\n", + "dis_merged_hesse3 = os.path.join(base_path, r\"newData\\hesse_shp\\c3dismerg_hes.shp\")\n", "\n", "merge_list3 = []\n", "for dirpath, dirnames, filenames in arcpy.da.Walk(workspace):\n", @@ -340,7 +340,7 @@ "metadata": {}, "outputs": [], "source": [ - "erased_frame3 = r\"F:\\newData\\hesse_shp\\erased3_hes.shp\"\n", + "erased_frame3 = os.path.join(base_path, r\"newData\\hesse_shp\\erased3_hes.shp\")\n", "arcpy.analysis.Erase(\"bbmeged_hesse\", dis_merged_hesse3, erased_frame3, None)" ] }, @@ -350,11 +350,11 @@ "metadata": {}, "outputs": [], "source": [ - "diff_19_14 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "diff_21_19 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", + "diff_19_14 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "diff_21_19 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", "\n", - "diff_19_14_mask3 = r\"F:\\newData\\diff19_14_hesse_mask3.tif\"\n", - "diff_21_19_mask3 = r\"F:\\newData\\diff21_19_hesse_mask3.tif\"\n", + "diff_19_14_mask3 = os.path.join(base_path, r\"newData\\diff19_14_hesse_mask3.tif\")\n", + "diff_21_19_mask3 = os.path.join(base_path, r\"newData\\diff21_19_hesse_mask3.tif\")\n", "\n", "outExtractByMask3_1914 = arcpy.sa.ExtractByMask(diff_19_14, erased_frame3)\n", "outExtractByMask3_1914.save(diff_19_14_mask3)\n", @@ -378,9 +378,9 @@ "source": [ "# b) erst mergen dann maske abziehen, dann von gesamter diff maskieren\n", "\n", - "workspace = r\"F:\\newData\\hesse_shp\\grid\"\n", - "merged_hesse5 = r\"F:\\newData\\hesse_shp\\c5merg_hes.shp\"\n", - "dis_merged_hesse5 = r\"F:\\newData\\hesse_shp\\c5dismerg_hes.shp\"\n", + "workspace = os.path.join(base_path, r\"newData\\hesse_shp\\grid\")\n", + "merged_hesse5 = os.path.join(base_path, r\"newData\\hesse_shp\\c5merg_hes.shp\")\n", + "dis_merged_hesse5 = os.path.join(base_path, r\"newData\\hesse_shp\\c5dismerg_hes.shp\")\n", "\n", "merge_list5 = []\n", "for dirpath, dirnames, filenames in arcpy.da.Walk(workspace):\n", @@ -398,7 +398,7 @@ "metadata": {}, "outputs": [], "source": [ - "erased_frame5 = r\"F:\\newData\\hesse_shp\\erased5_hes.shp\"\n", + "erased_frame5 = os.path.join(base_path, r\"newData\\hesse_shp\\erased5_hes.shp\")\n", "arcpy.analysis.Erase(\"bbmeged_hesse\", dis_merged_hesse5, erased_frame5, None)" ] }, @@ -408,11 +408,11 @@ "metadata": {}, "outputs": [], "source": [ - "diff_19_14 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "diff_21_19 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", + "diff_19_14 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "diff_21_19 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", "\n", - "diff_19_14_mask5 = r\"F:\\newData\\diff19_14_hesse_mask5.tif\"\n", - "diff_21_19_mask5 = r\"F:\\newData\\diff21_19_hesse_mask5.tif\"\n", + "diff_19_14_mask5 = os.path.join(base_path, r\"newData\\diff19_14_hesse_mask5.tif\")\n", + "diff_21_19_mask5 = os.path.join(base_path, r\"newData\\diff21_19_hesse_mask5.tif\")\n", "\n", "outExtractByMask5_1914 = arcpy.sa.ExtractByMask(diff_19_14, erased_frame5)\n", "outExtractByMask5_1914.save(diff_19_14_mask5)\n", @@ -436,7 +436,7 @@ "source": [ "# Erase with raster_50x50 to get mask\n", "\n", - "workspace = r\"F:\\newData\\hesse_shp\\grid\"\n", + "workspace = os.path.join(base_path, r\"newData\\hesse_shp\\grid\")\n", "\n", "i = 0\n", "for dirpath, dirnames, filenames in arcpy.da.Walk(workspace):\n", @@ -444,11 +444,11 @@ " if filename.startswith(r'dis_hesse_raster_50x50'):\n", " inp_frame = r\"{0}\\{1}\".format(dirpath, filename)\n", " #inp_frame = r\"{0}\".format(filename[:25])\n", - " buffer_complete3 = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_merged3_{}\".format(i, i) \n", - " buffer_complete5 = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_merged5_{}\".format(i, i) \n", + " buffer_complete3 = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_merged3_{}\").format(i, i) \n", + " buffer_complete5 = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_dis_merged5_{}\").format(i, i) \n", " \n", - " out_erased3 = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_erased3_{}\".format(i, i)\n", - " out_erased5 = r\"F:\\newData\\hesse_shp\\grid\\osm_{}\\buf10_erased5_{}\".format(i, i)\n", + " out_erased3 = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_erased3_{}\").format(i, i)\n", + " out_erased5 = os.path.join(base_path, r\"newData\\hesse_shp\\grid\\osm_{}\\buf10_erased5_{}\").format(i, i)\n", " \n", " dis_filepath = r\"{0}\\osm_{1}\\{2}_quarry_{3}\".format(dirpath, i, filename[0:9], i)\n", " out_buf = r\"{0}\\osm_{1}\\buf10_dis_quarry_{2}\".format(dirpath, i, i)\n", diff --git a/notebooks/area_of_interest/220423_template_aoi.ipynb b/notebooks/area_of_interest/220423_template_aoi.ipynb index d8000aa..5df3926 100644 --- a/notebooks/area_of_interest/220423_template_aoi.ipynb +++ b/notebooks/area_of_interest/220423_template_aoi.ipynb @@ -29,7 +29,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -48,7 +49,7 @@ "# Paths: # toDo -> check in which grid polygone the aoi is + completely new inp paths <-> saves time and computing power \n", "\n", "### Input\n", - "aoi_poly = r\"F:\\newData\\Offenbach\\of_shp\\aoi_of.shp\"\n", + "aoi_poly = os.path.join(base_path, r\"newData\\Offenbach\\of_shp\\aoi_of.shp\")\n", "\n", "buff_dist_build = \"10 Meters\"\n", "buff_dist_quarry = \"10 Meters\"\n", @@ -57,34 +58,34 @@ "buff_dist_railw = \"10 Meters\"\n", "\n", "# Difference-Rasters / LIDAR\n", - "diff_rast_1914 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "diff_rast_2119 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", + "diff_rast_1914 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "diff_rast_2119 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", "\n", "# BBD-Data / PSI\n", - "psi_los_asce = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\"\n", - "psi_los_desc = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\"\n", - "psi_east_west = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\"\n", - "psi_vertical = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\"\n", + "psi_los_asce = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\")\n", + "psi_los_desc = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\")\n", + "psi_east_west = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\")\n", + "psi_vertical = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\")\n", "\n", "# OSM-Data\n", - "osm_build = r\"F:\\newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\"\n", - "osm_landu = r\"F:\\newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\"\n", - "osm_water = r\"F:\\newData\\hesse_shp\\gis_osm_water_a_free_1.shp\"\n", - "osm_railw = r\"F:\\newData\\hesse_shp\\gis_osm_railways_free_1.shp\"\n", - "osm_roads = r\"F:\\newData\\hesse_shp\\gis_osm_roads_free_1.shp\"\n", + "osm_build = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\")\n", + "osm_landu = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\")\n", + "osm_water = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_water_a_free_1.shp\")\n", + "osm_railw = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_railways_free_1.shp\")\n", + "osm_roads = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_roads_free_1.shp\")\n", "\n", "# Subsidence - Data\n", - "lyr_subsidence = r\"F:\\Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\"\n", - "gdb_subsidence = r\"F:\\Umwelt 4.0\\setzungen.gdb\"\n", + "lyr_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\")\n", + "gdb_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\setzungen.gdb\")\n", "\n", "### Output\n", "\n", "# Folder Paths\n", - "folder_of = r\"F:\\newData\\Offenbach\"\n", - "folder_of_rast = r\"F:\\newData\\Offenbach\\of_rast\"\n", - "folder_of_shp = r\"F:\\newData\\Offenbach\\of_shp\"\n", - "folder_of_psi = r\"F:\\newData\\Offenbach\\of_psi\"\n", - "folder_of_subs = r\"F:\\newData\\Offenbach\\of_subs\"\n", + "folder_of = os.path.join(base_path, r\"newData\\Offenbach\")\n", + "folder_of_rast = os.path.join(base_path, r\"newData\\Offenbach\\of_rast\")\n", + "folder_of_shp = os.path.join(base_path, r\"newData\\Offenbach\\of_shp\")\n", + "folder_of_psi = os.path.join(base_path, r\"newData\\Offenbach\\of_psi\")\n", + "folder_of_subs = os.path.join(base_path, r\"newData\\Offenbach\\of_subs\")\n", "\n", "\n", "# Diff files\n", diff --git a/notebooks/area_of_interest/hohermeisner_aoi.ipynb b/notebooks/area_of_interest/hohermeisner_aoi.ipynb index 37b7138..db6896e 100644 --- a/notebooks/area_of_interest/hohermeisner_aoi.ipynb +++ b/notebooks/area_of_interest/hohermeisner_aoi.ipynb @@ -29,7 +29,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline " + "base_path = os.getcwd()" +"%matplotlib inline" ] }, { @@ -48,7 +49,7 @@ "# Paths: # toDo -> check in which grid polygone the aoi is + completely new inp paths <-> saves time and computing power \n", "\n", "### Input\n", - "aoi_poly = r\"F:\\newData\\HoherMeisner\\hm_shp\\aoi_hm.shp\"\n", + "aoi_poly = os.path.join(base_path, r\"newData\\HoherMeisner\\hm_shp\\aoi_hm.shp\")\n", "\n", "buff_dist_build = \"10 Meters\"\n", "buff_dist_quarry = \"10 Meters\"\n", @@ -57,34 +58,34 @@ "buff_dist_railw = \"10 Meters\"\n", "\n", "# Difference-Rasters / LIDAR\n", - "diff_rast_1914 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "diff_rast_2119 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", + "diff_rast_1914 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "diff_rast_2119 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", "\n", "# BBD-Data / PSI\n", - "psi_los_asce = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\"\n", - "psi_los_desc = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\"\n", - "psi_east_west = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\"\n", - "psi_vertical = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\"\n", + "psi_los_asce = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\")\n", + "psi_los_desc = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\")\n", + "psi_east_west = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\")\n", + "psi_vertical = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\")\n", "\n", "# OSM-Data\n", - "osm_build = r\"F:\\newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\"\n", - "osm_landu = r\"F:\\newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\"\n", - "osm_water = r\"F:\\newData\\hesse_shp\\gis_osm_water_a_free_1.shp\"\n", - "osm_railw = r\"F:\\newData\\hesse_shp\\gis_osm_railways_free_1.shp\"\n", - "osm_roads = r\"F:\\newData\\hesse_shp\\gis_osm_roads_free_1.shp\"\n", + "osm_build = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\")\n", + "osm_landu = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\")\n", + "osm_water = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_water_a_free_1.shp\")\n", + "osm_railw = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_railways_free_1.shp\")\n", + "osm_roads = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_roads_free_1.shp\")\n", "\n", "# Subsidence - Data\n", - "lyr_subsidence = r\"F:\\Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\"\n", - "gdb_subsidence = r\"F:\\Umwelt 4.0\\setzungen.gdb\"\n", + "lyr_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\")\n", + "gdb_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\setzungen.gdb\")\n", "\n", "### Output\n", "\n", "# Folder Paths\n", - "folder_of = r\"F:\\newData\\HoherMeisner\"\n", - "folder_of_rast = r\"F:\\newData\\HoherMeisner\\hm_rast\"\n", - "folder_of_shp = r\"F:\\newData\\HoherMeisner\\hm_shp\"\n", - "folder_of_psi = r\"F:\\newData\\HoherMeisner\\hm_psi\"\n", - "folder_of_subs = r\"F:\\newData\\HoherMeisner\\hm_subs\"\n", + "folder_of = os.path.join(base_path, r\"newData\\HoherMeisner\")\n", + "folder_of_rast = os.path.join(base_path, r\"newData\\HoherMeisner\\hm_rast\")\n", + "folder_of_shp = os.path.join(base_path, r\"newData\\HoherMeisner\\hm_shp\")\n", + "folder_of_psi = os.path.join(base_path, r\"newData\\HoherMeisner\\hm_psi\")\n", + "folder_of_subs = os.path.join(base_path, r\"newData\\HoherMeisner\\hm_subs\")\n", "\n", "\n", "# Diff files\n", diff --git a/notebooks/area_of_interest/kassel_aoi.ipynb b/notebooks/area_of_interest/kassel_aoi.ipynb index 0584912..31d2fb7 100644 --- a/notebooks/area_of_interest/kassel_aoi.ipynb +++ b/notebooks/area_of_interest/kassel_aoi.ipynb @@ -29,7 +29,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -48,7 +49,7 @@ "# Paths: # toDo -> check in which grid polygone the aoi is + completely new inp paths <-> saves time and computing power \n", "\n", "### Input\n", - "aoi_poly = r\"F:\\newData\\Kassel\\ka_shp\\aoi_ka.shp\"\n", + "aoi_poly = os.path.join(base_path, r\"newData\\Kassel\\ka_shp\\aoi_ka.shp\")\n", "\n", "buff_dist_build = \"10 Meters\"\n", "buff_dist_quarry = \"10 Meters\"\n", @@ -57,34 +58,34 @@ "buff_dist_railw = \"10 Meters\"\n", "\n", "# Difference-Rasters / LIDAR\n", - "diff_rast_1914 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "diff_rast_2119 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", + "diff_rast_1914 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "diff_rast_2119 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", "\n", "# BBD-Data / PSI\n", - "psi_los_asce = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\"\n", - "psi_los_desc = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\"\n", - "psi_east_west = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\"\n", - "psi_vertical = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\"\n", + "psi_los_asce = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\")\n", + "psi_los_desc = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\")\n", + "psi_east_west = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\")\n", + "psi_vertical = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\")\n", "\n", "# OSM-Data\n", - "osm_build = r\"F:\\newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\"\n", - "osm_landu = r\"F:\\newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\"\n", - "osm_water = r\"F:\\newData\\hesse_shp\\gis_osm_water_a_free_1.shp\"\n", - "osm_railw = r\"F:\\newData\\hesse_shp\\gis_osm_railways_free_1.shp\"\n", - "osm_roads = r\"F:\\newData\\hesse_shp\\gis_osm_roads_free_1.shp\"\n", + "osm_build = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\")\n", + "osm_landu = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\")\n", + "osm_water = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_water_a_free_1.shp\")\n", + "osm_railw = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_railways_free_1.shp\")\n", + "osm_roads = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_roads_free_1.shp\")\n", "\n", "# Subsidence - Data\n", - "lyr_subsidence = r\"F:\\Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\"\n", - "gdb_subsidence = r\"F:\\Umwelt 4.0\\setzungen.gdb\"\n", + "lyr_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\")\n", + "gdb_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\setzungen.gdb\")\n", "\n", "### Output\n", "\n", "# Folder Paths\n", - "folder_of = r\"F:\\newData\\Kassel\"\n", - "folder_of_rast = r\"F:\\newData\\Kassel\\ka_rast\"\n", - "folder_of_shp = r\"F:\\newData\\Kassel\\ka_shp\"\n", - "folder_of_psi = r\"F:\\newData\\Kassel\\ka_psi\"\n", - "folder_of_subs = r\"F:\\newData\\Kassel\\ka_subs\"\n", + "folder_of = os.path.join(base_path, r\"newData\\Kassel\")\n", + "folder_of_rast = os.path.join(base_path, r\"newData\\Kassel\\ka_rast\")\n", + "folder_of_shp = os.path.join(base_path, r\"newData\\Kassel\\ka_shp\")\n", + "folder_of_psi = os.path.join(base_path, r\"newData\\Kassel\\ka_psi\")\n", + "folder_of_subs = os.path.join(base_path, r\"newData\\Kassel\\ka_subs\")\n", "\n", "\n", "# Diff files\n", diff --git a/notebooks/area_of_interest/offenbach_aoi.ipynb b/notebooks/area_of_interest/offenbach_aoi.ipynb index c88740b..3666fc1 100644 --- a/notebooks/area_of_interest/offenbach_aoi.ipynb +++ b/notebooks/area_of_interest/offenbach_aoi.ipynb @@ -29,7 +29,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -48,7 +49,7 @@ "# Paths: # toDo -> check in which grid polygone the aoi is + completely new inp paths <-> saves time and computing power \n", "\n", "### Input\n", - "aoi_poly = r\"F:\\newData\\Offenbach\\of_shp\\aoi_of.shp\"\n", + "aoi_poly = os.path.join(base_path, r\"newData\\Offenbach\\of_shp\\aoi_of.shp\")\n", "\n", "buff_dist_build = \"10 Meters\"\n", "buff_dist_quarry = \"10 Meters\"\n", @@ -57,34 +58,34 @@ "buff_dist_railw = \"10 Meters\"\n", "\n", "# Difference-Rasters / LIDAR\n", - "diff_rast_1914 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "diff_rast_2119 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", + "diff_rast_1914 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "diff_rast_2119 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", "\n", "# BBD-Data / PSI\n", - "psi_los_asce = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\"\n", - "psi_los_desc = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\"\n", - "psi_east_west = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\"\n", - "psi_vertical = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\"\n", + "psi_los_asce = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\")\n", + "psi_los_desc = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\")\n", + "psi_east_west = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\")\n", + "psi_vertical = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\")\n", "\n", "# OSM-Data\n", - "osm_build = r\"F:\\newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\"\n", - "osm_landu = r\"F:\\newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\"\n", - "osm_water = r\"F:\\newData\\hesse_shp\\gis_osm_water_a_free_1.shp\"\n", - "osm_railw = r\"F:\\newData\\hesse_shp\\gis_osm_railways_free_1.shp\"\n", - "osm_roads = r\"F:\\newData\\hesse_shp\\gis_osm_roads_free_1.shp\"\n", + "osm_build = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\")\n", + "osm_landu = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\")\n", + "osm_water = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_water_a_free_1.shp\")\n", + "osm_railw = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_railways_free_1.shp\")\n", + "osm_roads = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_roads_free_1.shp\")\n", "\n", "# Subsidence - Data\n", - "lyr_subsidence = r\"F:\\Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\"\n", - "gdb_subsidence = r\"F:\\Umwelt 4.0\\setzungen.gdb\"\n", + "lyr_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\")\n", + "gdb_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\setzungen.gdb\")\n", "\n", "### Output\n", "\n", "# Folder Paths\n", - "folder_of = r\"F:\\newData\\Offenbach\"\n", - "folder_of_rast = r\"F:\\newData\\Offenbach\\of_rast\"\n", - "folder_of_shp = r\"F:\\newData\\Offenbach\\of_shp\"\n", - "folder_of_psi = r\"F:\\newData\\Offenbach\\of_psi\"\n", - "folder_of_subs = r\"F:\\newData\\Offenbach\\of_subs\"\n", + "folder_of = os.path.join(base_path, r\"newData\\Offenbach\")\n", + "folder_of_rast = os.path.join(base_path, r\"newData\\Offenbach\\of_rast\")\n", + "folder_of_shp = os.path.join(base_path, r\"newData\\Offenbach\\of_shp\")\n", + "folder_of_psi = os.path.join(base_path, r\"newData\\Offenbach\\of_psi\")\n", + "folder_of_subs = os.path.join(base_path, r\"newData\\Offenbach\\of_subs\")\n", "\n", "\n", "# Diff files\n", diff --git a/notebooks/area_of_interest/rheinmain_aoi.ipynb b/notebooks/area_of_interest/rheinmain_aoi.ipynb index b59010f..ff1a31c 100644 --- a/notebooks/area_of_interest/rheinmain_aoi.ipynb +++ b/notebooks/area_of_interest/rheinmain_aoi.ipynb @@ -29,7 +29,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -48,7 +49,7 @@ "# Paths: # toDo -> check in which grid polygone the aoi is + completely new inp paths <-> saves time and computing power \n", "\n", "### Input\n", - "aoi_poly = r\"F:\\newData\\RheinMain\\rm_shp\\aoi_rm.shp\"\n", + "aoi_poly = os.path.join(base_path, r\"newData\\RheinMain\\rm_shp\\aoi_rm.shp\")\n", "\n", "buff_dist_build = \"10 Meters\"\n", "buff_dist_quarry = \"10 Meters\"\n", @@ -57,34 +58,34 @@ "buff_dist_railw = \"10 Meters\"\n", "\n", "# Difference-Rasters / LIDAR\n", - "diff_rast_1914 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "diff_rast_2119 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", + "diff_rast_1914 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "diff_rast_2119 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", "\n", "# BBD-Data / PSI\n", - "psi_los_asce = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\"\n", - "psi_los_desc = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\"\n", - "psi_east_west = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\"\n", - "psi_vertical = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\"\n", + "psi_los_asce = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\")\n", + "psi_los_desc = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\")\n", + "psi_east_west = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\")\n", + "psi_vertical = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\")\n", "\n", "# OSM-Data\n", - "osm_build = r\"F:\\newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\"\n", - "osm_landu = r\"F:\\newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\"\n", - "osm_water = r\"F:\\newData\\hesse_shp\\gis_osm_water_a_free_1.shp\"\n", - "osm_railw = r\"F:\\newData\\hesse_shp\\gis_osm_railways_free_1.shp\"\n", - "osm_roads = r\"F:\\newData\\hesse_shp\\gis_osm_roads_free_1.shp\"\n", + "osm_build = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\")\n", + "osm_landu = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\")\n", + "osm_water = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_water_a_free_1.shp\")\n", + "osm_railw = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_railways_free_1.shp\")\n", + "osm_roads = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_roads_free_1.shp\")\n", "\n", "# Subsidence - Data\n", - "lyr_subsidence = r\"F:\\Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\"\n", - "gdb_subsidence = r\"F:\\Umwelt 4.0\\setzungen.gdb\"\n", + "lyr_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\")\n", + "gdb_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\setzungen.gdb\")\n", "\n", "### Output\n", "\n", "# Folder Paths\n", - "folder_of = r\"F:\\newData\\RheinMain\"\n", - "folder_of_rast = r\"F:\\newData\\RheinMain\\rm_rast\"\n", - "folder_of_shp = r\"F:\\newData\\RheinMain\\rm_shp\"\n", - "folder_of_psi = r\"F:\\newData\\RheinMain\\rm_psi\"\n", - "folder_of_subs = r\"F:\\newData\\RheinMain\\rm_subs\"\n", + "folder_of = os.path.join(base_path, r\"newData\\RheinMain\")\n", + "folder_of_rast = os.path.join(base_path, r\"newData\\RheinMain\\rm_rast\")\n", + "folder_of_shp = os.path.join(base_path, r\"newData\\RheinMain\\rm_shp\")\n", + "folder_of_psi = os.path.join(base_path, r\"newData\\RheinMain\\rm_psi\")\n", + "folder_of_subs = os.path.join(base_path, r\"newData\\RheinMain\\rm_subs\")\n", "\n", "\n", "# Diff files\n", diff --git a/notebooks/area_of_interest/werrakali_aoi.ipynb b/notebooks/area_of_interest/werrakali_aoi.ipynb index 657dff0..799344e 100644 --- a/notebooks/area_of_interest/werrakali_aoi.ipynb +++ b/notebooks/area_of_interest/werrakali_aoi.ipynb @@ -29,7 +29,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -48,7 +49,7 @@ "# Paths: # toDo -> check in which grid polygone the aoi is + completely new inp paths <-> saves time and computing power \n", "\n", "### Input\n", - "aoi_poly = r\"F:\\newData\\WerraKali\\wk_shp\\aoi_wk.shp\"\n", + "aoi_poly = os.path.join(base_path, r\"newData\\WerraKali\\wk_shp\\aoi_wk.shp\")\n", "\n", "buff_dist_build = \"10 Meters\"\n", "buff_dist_quarry = \"10 Meters\"\n", @@ -57,34 +58,34 @@ "buff_dist_railw = \"10 Meters\"\n", "\n", "# Difference-Rasters / LIDAR\n", - "diff_rast_1914 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "diff_rast_2119 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", + "diff_rast_1914 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "diff_rast_2119 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", "\n", "# BBD-Data / PSI\n", - "psi_los_asce = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\"\n", - "psi_los_desc = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\"\n", - "psi_east_west = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\"\n", - "psi_vertical = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\"\n", + "psi_los_asce = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_asce_mittlere_Geschwindigkeit\")\n", + "psi_los_desc = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_LOS_desc_mittlere_Geschwindigkeit\")\n", + "psi_east_west = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Ost_West\")\n", + "psi_vertical = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\\BBD_2021_PSI_Vertikal\")\n", "\n", "# OSM-Data\n", - "osm_build = r\"F:\\newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\"\n", - "osm_landu = r\"F:\\newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\"\n", - "osm_water = r\"F:\\newData\\hesse_shp\\gis_osm_water_a_free_1.shp\"\n", - "osm_railw = r\"F:\\newData\\hesse_shp\\gis_osm_railways_free_1.shp\"\n", - "osm_roads = r\"F:\\newData\\hesse_shp\\gis_osm_roads_free_1.shp\"\n", + "osm_build = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\")\n", + "osm_landu = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\")\n", + "osm_water = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_water_a_free_1.shp\")\n", + "osm_railw = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_railways_free_1.shp\")\n", + "osm_roads = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_roads_free_1.shp\")\n", "\n", "# Subsidence - Data\n", - "lyr_subsidence = r\"F:\\Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\"\n", - "gdb_subsidence = r\"F:\\Umwelt 4.0\\setzungen.gdb\"\n", + "lyr_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\Setzungsempfindliche_Schichten.lpk\")\n", + "gdb_subsidence = os.path.join(base_path, r\"Umwelt 4.0\\setzungen.gdb\")\n", "\n", "### Output\n", "\n", "# Folder Paths\n", - "folder_of = r\"F:\\newData\\WerraKali\"\n", - "folder_of_rast = r\"F:\\newData\\WerraKali\\wk_rast\"\n", - "folder_of_shp = r\"F:\\newData\\WerraKali\\wk_shp\"\n", - "folder_of_psi = r\"F:\\newData\\WerraKali\\wk_psi\"\n", - "folder_of_subs = r\"F:\\newData\\WerraKali\\wk_subs\"\n", + "folder_of = os.path.join(base_path, r\"newData\\WerraKali\")\n", + "folder_of_rast = os.path.join(base_path, r\"newData\\WerraKali\\wk_rast\")\n", + "folder_of_shp = os.path.join(base_path, r\"newData\\WerraKali\\wk_shp\")\n", + "folder_of_psi = os.path.join(base_path, r\"newData\\WerraKali\\wk_psi\")\n", + "folder_of_subs = os.path.join(base_path, r\"newData\\WerraKali\\wk_subs\")\n", "\n", "\n", "# Diff files\n", diff --git a/notebooks/data_prep/DataPrep_Diff_QuarryHM2119.ipynb b/notebooks/data_prep/DataPrep_Diff_QuarryHM2119.ipynb index 1c6b519..a2ccbfd 100644 --- a/notebooks/data_prep/DataPrep_Diff_QuarryHM2119.ipynb +++ b/notebooks/data_prep/DataPrep_Diff_QuarryHM2119.ipynb @@ -31,7 +31,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -50,13 +51,13 @@ "metadata": {}, "outputs": [], "source": [ - "buffer_quarry = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_quarry_10m.shp\"\n", + "buffer_quarry = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_quarry_10m.shp\")\n", "\n", - "buffer_complete_dis = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\hesse_shp_quarry_hm_buffer_10m.shp\"\n", + "buffer_complete_dis = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\hesse_shp_quarry_hm_buffer_10m.shp\")\n", "\n", - "diff_21_19 = r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner.tif\"\n", + "diff_21_19 = os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner.tif\")\n", "\n", - "diff_21_19_mask = r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner_quarry_mask.tif\"\n", + "diff_21_19_mask = os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner_quarry_mask.tif\")\n", "\n", "distance = \"10 Meters\"" ] @@ -159,34 +160,34 @@ "outputs": [], "source": [ "outReclass1 = Reclassify(diff_21_19_mask, \"Value\", range_class01)\n", - "outReclass1.save(r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r01.tif\")\n", + "outReclass1.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r01.tif\"))\n", "\n", "outReclass2 = Reclassify(diff_21_19_mask, \"Value\", range_class02)\n", - "outReclass2.save(r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r02.tif\")\n", + "outReclass2.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r02.tif\"))\n", "\n", "outReclass3 = Reclassify(diff_21_19_mask, \"Value\", range_class03)\n", - "outReclass3.save(r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r03.tif\")\n", + "outReclass3.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r03.tif\"))\n", "\n", "outReclass4 = Reclassify(diff_21_19_mask, \"Value\", range_class04)\n", - "outReclass4.save(r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r04.tif\")\n", + "outReclass4.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r04.tif\"))\n", "\n", "outReclass5 = Reclassify(diff_21_19_mask, \"Value\", range_class05)\n", - "outReclass5.save(r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r05.tif\")\n", + "outReclass5.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r05.tif\"))\n", "\n", "outReclass6 = Reclassify(diff_21_19_mask, \"Value\", range_class06)\n", - "outReclass6.save(r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r06.tif\")\n", + "outReclass6.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r06.tif\"))\n", "\n", "outReclass7 = Reclassify(diff_21_19_mask, \"Value\", range_class07)\n", - "outReclass7.save(r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r07.tif\")\n", + "outReclass7.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r07.tif\"))\n", "\n", "outReclass8 = Reclassify(diff_21_19_mask, \"Value\", range_class08)\n", - "outReclass8.save(r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r08.tif\")\n", + "outReclass8.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r08.tif\"))\n", "\n", "outReclass9 = Reclassify(diff_21_19_mask, \"Value\", range_class09)\n", - "outReclass9.save(r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r09.tif\")\n", + "outReclass9.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r09.tif\"))\n", "\n", "outReclass10 = Reclassify(diff_21_19_mask, \"Value\", range_class10)\n", - "outReclass10.save(r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r10.tif\")" + "outReclass10.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_quarryHoherM_r10.tif\"))" ] }, { diff --git a/notebooks/data_prep/DataPrep_Diff_QuarryTemplate_HM1914.ipynb b/notebooks/data_prep/DataPrep_Diff_QuarryTemplate_HM1914.ipynb index 60771fd..dfeb97e 100644 --- a/notebooks/data_prep/DataPrep_Diff_QuarryTemplate_HM1914.ipynb +++ b/notebooks/data_prep/DataPrep_Diff_QuarryTemplate_HM1914.ipynb @@ -31,7 +31,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -50,15 +51,15 @@ "metadata": {}, "outputs": [], "source": [ - "buffer_quarry = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_quarry_10m.shp\"\n", + "buffer_quarry = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_quarry_10m.shp\")\n", "\n", - "buffer_complete_dis = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\hesse_shp_quarry_hm_buffer_10m.shp\"\n", + "buffer_complete_dis = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\hesse_shp_quarry_hm_buffer_10m.shp\")\n", "\n", - "diff_19_14 = r\"F:\\newData\\HoherMeisner\\diff_19_14_hohermeisner\\diff19_14_HoherMeisner.tif\"\n", - "#diff_21_19 = r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner.tif\"\n", + "diff_19_14 = os.path.join(base_path, r\"newData\\HoherMeisner\\diff_19_14_hohermeisner\\diff19_14_HoherMeisner.tif\")\n", + "#diff_21_19 = os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner.tif\")\n", "\n", - "diff_19_14_mask = r\"F:\\newData\\HoherMeisner\\diff_19_14_hohermeisner\\diff19_14_HoherMeisner_quarry_mask.tif\"\n", - "#diff_21_19_mask = r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner_quarry_mask.tif\"\n", + "diff_19_14_mask = os.path.join(base_path, r\"newData\\HoherMeisner\\diff_19_14_hohermeisner\\diff19_14_HoherMeisner_quarry_mask.tif\")\n", + "#diff_21_19_mask = os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner_quarry_mask.tif\")\n", "\n", "distance = \"10 Meters\"" ] @@ -161,34 +162,34 @@ "outputs": [], "source": [ "outReclass1 = Reclassify(diff_19_14_mask, \"Value\", range_class01)\n", - "outReclass1.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r01.tif\")\n", + "outReclass1.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r01.tif\"))\n", "\n", "outReclass2 = Reclassify(diff_19_14_mask, \"Value\", range_class02)\n", - "outReclass2.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r02.tif\")\n", + "outReclass2.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r02.tif\"))\n", "\n", "outReclass3 = Reclassify(diff_19_14_mask, \"Value\", range_class03)\n", - "outReclass3.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r03.tif\")\n", + "outReclass3.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r03.tif\"))\n", "\n", "outReclass4 = Reclassify(diff_19_14_mask, \"Value\", range_class04)\n", - "outReclass4.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r04.tif\")\n", + "outReclass4.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r04.tif\"))\n", "\n", "outReclass5 = Reclassify(diff_19_14_mask, \"Value\", range_class05)\n", - "outReclass5.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r05.tif\")\n", + "outReclass5.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r05.tif\"))\n", "\n", "outReclass6 = Reclassify(diff_19_14_mask, \"Value\", range_class06)\n", - "outReclass6.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r06.tif\")\n", + "outReclass6.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r06.tif\"))\n", "\n", "outReclass7 = Reclassify(diff_19_14_mask, \"Value\", range_class07)\n", - "outReclass7.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r07.tif\")\n", + "outReclass7.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r07.tif\"))\n", "\n", "outReclass8 = Reclassify(diff_19_14_mask, \"Value\", range_class08)\n", - "outReclass8.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r08.tif\")\n", + "outReclass8.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r08.tif\"))\n", "\n", "outReclass9 = Reclassify(diff_19_14_mask, \"Value\", range_class09)\n", - "outReclass9.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r09.tif\")\n", + "outReclass9.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r09.tif\"))\n", "\n", "outReclass10 = Reclassify(diff_19_14_mask, \"Value\", range_class10)\n", - "outReclass10.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r10.tif\")" + "outReclass10.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_quarryHoherM_r10.tif\"))" ] }, { diff --git a/notebooks/data_prep/DataPrep_Diff_Template_HM1914.ipynb b/notebooks/data_prep/DataPrep_Diff_Template_HM1914.ipynb index f4418e4..63b3146 100644 --- a/notebooks/data_prep/DataPrep_Diff_Template_HM1914.ipynb +++ b/notebooks/data_prep/DataPrep_Diff_Template_HM1914.ipynb @@ -31,7 +31,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -50,27 +51,27 @@ "metadata": {}, "outputs": [], "source": [ - "inp_fc_buildings = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\gis_osm_buildings_a_free_1.shp\"\n", - "inp_fc_landuse = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\gis_osm_landuse_a_free_1.shp\"\n", - "inp_fc_water = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\gis_osm_water_a_free_1.shp\"\n", + "inp_fc_buildings = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\gis_osm_buildings_a_free_1.shp\")\n", + "inp_fc_landuse = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\gis_osm_landuse_a_free_1.shp\")\n", + "inp_fc_water = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\gis_osm_water_a_free_1.shp\")\n", "\n", - "only_quarry = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_quarry.shp\" \n", + "only_quarry = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_quarry.shp\") \n", "\n", - "buffer_buildings = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_buildings_10m.shp\"\n", - "buffer_quarry = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_quarry_10m.shp\"\n", - "buffer_water = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_water_10m.shp\"\n", + "buffer_buildings = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_buildings_10m.shp\")\n", + "buffer_quarry = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_quarry_10m.shp\")\n", + "buffer_water = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_water_10m.shp\")\n", "\n", - "buffer_complete = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_complete_10m.shp\"\n", - "buffer_complete_dis = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\hesse_shp_hohermeisnerbuffer_co_dis_10m.shp\"\n", + "buffer_complete = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\buffer_complete_10m.shp\")\n", + "buffer_complete_dis = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\hesse_shp_hohermeisnerbuffer_co_dis_10m.shp\")\n", "\n", "inp_frame = \"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\HoherMeisner.shp\"\n", - "erased_frame = r\"F:\\newData\\HoherMeisner\\hesse_shp_hohermeisner\\erased_frame.shp\"\n", + "erased_frame = os.path.join(base_path, r\"newData\\HoherMeisner\\hesse_shp_hohermeisner\\erased_frame.shp\")\n", "\n", - "diff_19_14 = r\"F:\\newData\\HoherMeisner\\diff_19_14_hohermeisner\\diff19_14_HoherMeisner.tif\"\n", - "#diff_21_19 = r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner.tif\"\n", + "diff_19_14 = os.path.join(base_path, r\"newData\\HoherMeisner\\diff_19_14_hohermeisner\\diff19_14_HoherMeisner.tif\")\n", + "#diff_21_19 = os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner.tif\")\n", "\n", - "diff_19_14_mask = r\"F:\\newData\\HoherMeisner\\diff_19_14_hohermeisner\\diff19_14_HoherMeisner_mask.tif\"\n", - "#diff_21_19_mask = r\"F:\\newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner_mask.tif\"\n", + "diff_19_14_mask = os.path.join(base_path, r\"newData\\HoherMeisner\\diff_19_14_hohermeisner\\diff19_14_HoherMeisner_mask.tif\")\n", + "#diff_21_19_mask = os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_21_19_hohermeisner\\diff21_19_HoherMeisner_mask.tif\")\n", "\n", "distance = \"10 Meters\"" ] @@ -251,34 +252,34 @@ "outputs": [], "source": [ "outReclass1 = Reclassify(diff_19_14_mask, \"Value\", range_class01)\n", - "outReclass1.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r01.tif\")\n", + "outReclass1.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r01.tif\"))\n", "\n", "outReclass2 = Reclassify(diff_19_14_mask, \"Value\", range_class02)\n", - "outReclass2.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r02.tif\")\n", + "outReclass2.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r02.tif\"))\n", "\n", "outReclass3 = Reclassify(diff_19_14_mask, \"Value\", range_class03)\n", - "outReclass3.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r03.tif\")\n", + "outReclass3.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r03.tif\"))\n", "\n", "outReclass4 = Reclassify(diff_19_14_mask, \"Value\", range_class04)\n", - "outReclass4.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r04.tif\")\n", + "outReclass4.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r04.tif\"))\n", "\n", "outReclass5 = Reclassify(diff_19_14_mask, \"Value\", range_class05)\n", - "outReclass5.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r05.tif\")\n", + "outReclass5.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r05.tif\"))\n", "\n", "outReclass6 = Reclassify(diff_19_14_mask, \"Value\", range_class06)\n", - "outReclass6.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r06.tif\")\n", + "outReclass6.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r06.tif\"))\n", "\n", "outReclass7 = Reclassify(diff_19_14_mask, \"Value\", range_class07)\n", - "outReclass7.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r07.tif\")\n", + "outReclass7.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r07.tif\"))\n", "\n", "outReclass8 = Reclassify(diff_19_14_mask, \"Value\", range_class08)\n", - "outReclass8.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r08.tif\")\n", + "outReclass8.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r08.tif\"))\n", "\n", "outReclass9 = Reclassify(diff_19_14_mask, \"Value\", range_class09)\n", - "outReclass9.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r09.tif\")\n", + "outReclass9.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r09.tif\"))\n", "\n", "outReclass10 = Reclassify(diff_19_14_mask, \"Value\", range_class10)\n", - "outReclass10.save(r\"F:\\newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r10.tif\")" + "outReclass10.save(os.path.join(base_path, r\"newData\\HoherMeisner\\diff_ref_19_14_hohermeisner\\diff19_14_HoherM_r10.tif\"))" ] }, { diff --git a/notebooks/data_prep/DataPrep_KasselStadt.ipynb b/notebooks/data_prep/DataPrep_KasselStadt.ipynb index a559d7a..6851eb8 100644 --- a/notebooks/data_prep/DataPrep_KasselStadt.ipynb +++ b/notebooks/data_prep/DataPrep_KasselStadt.ipynb @@ -32,7 +32,8 @@ "# from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -56,9 +57,9 @@ }, "outputs": [], "source": [ - "raster_1914 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "area_of_interest = r\"F:\\newData\\Pilotregionen\\KasselStadt.shp\"\n", - "output_file = r\"F:\\newData\\diff19_14_Kassel.tif\"\n", + "raster_1914 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "area_of_interest = os.path.join(base_path, r\"newData\\Pilotregionen\\KasselStadt.shp\")\n", + "output_file = os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\")\n", "\n", "raster_1914_clip = Clip(raster_1914, area_of_interest) \n", "raster_1914_clip.save(output_file) " @@ -82,9 +83,9 @@ }, "outputs": [], "source": [ - "raster_2119 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", - "area_of_interest = r\"F:\\newData\\Pilotregionen\\KasselStadt.shp\"\n", - "output_file = r\"F:\\newData\\diff21_19_Kassel.tif\"\n", + "raster_2119 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", + "area_of_interest = os.path.join(base_path, r\"newData\\Pilotregionen\\KasselStadt.shp\")\n", + "output_file = os.path.join(base_path, r\"newData\\diff21_19_Kassel.tif\")\n", "\n", "raster_2119_clip = Clip(raster_2119, area_of_interest) \n", "raster_2119_clip.save(output_file) " @@ -132,7 +133,7 @@ }, "outputs": [], "source": [ - "#inp_raster = r\"F:\\newData\\diff19_14_Kassel.tif\"\n", + "#inp_raster = os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\")\n", "#min_value = -0.2\n", "#max_value = 0.2\n", "\n", @@ -177,7 +178,7 @@ }, "outputs": [], "source": [ - "#inp_raster = r\"F:\\newData\\diff21_19_Kassel.tif\"\n", + "#inp_raster = os.path.join(base_path, r\"newData\\diff21_19_Kassel.tif\")\n", "#min_value = -0.2\n", "#max_value = 0.2\n", "\n", @@ -213,7 +214,7 @@ "outputs": [], "source": [ "inp_raster = r\"Differenzen\\diff19_14_Kassel.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_K_re.tif\"\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_K_re.tif\")\n", "\n", "arcpy.management.CopyRaster(inp_raster, out_raster)" ] @@ -228,8 +229,8 @@ }, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff1914_K_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_K_re1.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff1914_K_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_K_re1.tif\")\n", "\n", "range_class = RemapRange([[stats_min_1914, -1,-1],[-1,1,\"NODATA\"],[1, stats_max_1914,1]])\n", "\n", @@ -247,8 +248,8 @@ }, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff1914_K_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_K_re2.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff1914_K_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_K_re2.tif\")\n", "\n", "range_class = RemapRange([[stats_min_1914, -5,-5],[-5,-2.5,-2],[-2.5,-1.5,-1],[-1.5,-0.1,0],[-0.1,0.1,\"NODATA\"],[0.1,1.5,0],[1.5,2.5,1],[2.5,5,2],[5, stats_max_1914,5]])\n", "\n", @@ -274,7 +275,7 @@ "outputs": [], "source": [ "inp_raster = r\"Differenzen\\diff21_19_Kassel.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_K_re.tif\"\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_K_re.tif\")\n", "\n", "arcpy.management.CopyRaster(inp_raster, out_raster)" ] @@ -289,8 +290,8 @@ }, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff2119_K_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_K_re1.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff2119_K_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_K_re1.tif\")\n", "\n", "range_class = RemapRange([[stats_min_2119, -1,-1],[-1,1,\"NODATA\"],[1, stats_max_2119,1]])\n", "\n", @@ -308,8 +309,8 @@ }, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff2119_K_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_K_re2.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff2119_K_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_K_re2.tif\")\n", "\n", "range_class = RemapRange([[stats_min_2119, -5,-5],[-5,-2.5,-2],[-2.5,-1.5,-1],[-1.5,-0.1,0],[-0.1,0.1,\"NODATA\"],[0.1,1.5,0],[1.5,2.5,1],[2.5,5,2],[5, stats_max_2119,5]])\n", "\n", @@ -334,9 +335,9 @@ }, "outputs": [], "source": [ - "inp_vect = r\"F:\\newData\\testworkspace.gdb\\rasterUTM1x1\"\n", - "area_of_int = r\"F:\\newData\\Pilotregionen\\KasselStadt.shp\"\n", - "out_data = r\"F:\\newData\\gridUTM1x1_Ka.shp\"\n", + "inp_vect = os.path.join(base_path, r\"newData\\testworkspace.gdb\\rasterUTM1x1\")\n", + "area_of_int = os.path.join(base_path, r\"newData\\Pilotregionen\\KasselStadt.shp\")\n", + "out_data = os.path.join(base_path, r\"newData\\gridUTM1x1_Ka.shp\")\n", "\n", "vect_clip = arcpy.analysis.Clip(inp_vect, area_of_int, out_data) " ] @@ -351,7 +352,7 @@ }, "outputs": [], "source": [ - "inp_path = r\"F:\\newData\\hesse_shp\"\n", + "inp_path = os.path.join(base_path, r\"newData\\hesse_shp\")\n", "\n", "\n", "osm_list = []\n", @@ -391,7 +392,7 @@ }, "outputs": [], "source": [ - "inp_path = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\"\n", + "inp_path = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\")\n", "\n", "bbd_list = []\n", "bbd_names = []\n", diff --git a/notebooks/data_prep/DataPrep_Processing_220411.ipynb b/notebooks/data_prep/DataPrep_Processing_220411.ipynb index 101ebcb..e3f88ac 100644 --- a/notebooks/data_prep/DataPrep_Processing_220411.ipynb +++ b/notebooks/data_prep/DataPrep_Processing_220411.ipynb @@ -29,7 +29,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -60,7 +61,7 @@ }, "outputs": [], "source": [ - "list_names = [row[0] for row in arcpy.da.SearchCursor(r\"F:\\newData\\Pilotregionen\\TEST.shp\", \"Name\")]\n", + "list_names = [row[0] for row in arcpy.da.SearchCursor(os.path.join(base_path, r\"newData\\Pilotregionen\\TEST.shp\", \"Name\"))]\n", "print(list_names)" ] }, @@ -106,7 +107,7 @@ "source": [ "for i, j in zip(list_names, expression):\n", " out_file = i + \".shp\"\n", - " arcpy.conversion.FeatureClassToFeatureClass(r\"F:\\newData\\Pilotregionen\\TEST.shp\", r\"F:\\newData\\Pilotregionen\", out_file, j)" + " arcpy.conversion.FeatureClassToFeatureClass(os.path.join(base_path, r\"newData\\Pilotregionen\\TEST.shp\", r\"F:\\newData\\Pilotregionen\"), out_file, j)" ] }, { @@ -127,9 +128,9 @@ }, "outputs": [], "source": [ - "test_r = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "r_clip = Clip(test_r, r\"F:\\newData\\Pilotregionen\\KasselStadt.shp\") \n", - "r_clip.save(r\"F:\\newData\\diff19_14_Kassel.tif\") " + "test_r = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "r_clip = Clip(test_r, os.path.join(base_path, r\"newData\\Pilotregionen\\KasselStadt.shp\")) \n", + "r_clip.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\")) " ] }, { @@ -176,7 +177,7 @@ "outputs": [], "source": [ "fig, ax1 = pyplot.subplots(1, 1, figsize=(14,7))\n", - "src = rasterio.open(r\"F:\\newData\\diff19_14_Kassel.tif\")\n", + "src = rasterio.open(os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\"))\n", "show_hist(src, bins=10000, lw=0.0, stacked=False, histtype='bar', title=\"Histogram\", ax=ax1)\n", "ax1.set_xlim(xmin=-0.2, xmax=0.2)\n", "pyplot.show()" @@ -201,7 +202,7 @@ }, "outputs": [], "source": [ - "arcpy.management.CopyRaster(r\"Differenzen\\diff19_14_Kassel.tif\", r\"F:\\newData\\diff1914_K_re.tif\")" + "arcpy.management.CopyRaster(r\"Differenzen\\diff19_14_Kassel.tif\", os.path.join(base_path, r\"newData\\diff1914_K_re.tif\"))" ] }, { @@ -212,8 +213,8 @@ "source": [ "range_class = RemapRange([[stats_min, -1,-1],[-1,1,\"NODATA\"],[1, stats_max,1]])\n", "\n", - "r_clip_reclas = Reclassify(r\"F:\\newData\\diff1914_K_re.tif\", \"Value\", range_class)\n", - "r_clip_reclas.save(r\"F:\\newData\\diff1914_K_re1.tif\")" + "r_clip_reclas = Reclassify(os.path.join(base_path, r\"newData\\diff1914_K_re.tif\", \"Value\"), range_class)\n", + "r_clip_reclas.save(os.path.join(base_path, r\"newData\\diff1914_K_re1.tif\"))" ] }, { @@ -223,8 +224,8 @@ "outputs": [], "source": [ "range_class = RemapRange([[stats_min, -5,-5],[-5,-2.5,-2],[-2.5,-1.5,-1],[-1.5,-0.1,0],[-0.1,0.1,\"NODATA\"],[0.1,1.5,0],[1.5,2.5,1],[2.5,5,2],[5, stats_max,5]])\n", - "r_clip_reclas = Reclassify(r\"F:\\newData\\diff1914_K_re.tif\", \"Value\", range_class)\n", - "r_clip_reclas.save(r\"F:\\newData\\diff1914_K_re2.tif\")" + "r_clip_reclas = Reclassify(os.path.join(base_path, r\"newData\\diff1914_K_re.tif\", \"Value\"), range_class)\n", + "r_clip_reclas.save(os.path.join(base_path, r\"newData\\diff1914_K_re2.tif\"))" ] }, { @@ -240,9 +241,9 @@ "metadata": {}, "outputs": [], "source": [ - "inp_vect = r\"F:\\newData\\testworkspace.gdb\\rasterUTM1x1\"\n", - "area_of_int = r\"F:\\newData\\Pilotregionen\\KasselStadt.shp\"\n", - "out_data = r\"F:\\newData\\gridUTM1x1_Ka.shp\"\n", + "inp_vect = os.path.join(base_path, r\"newData\\testworkspace.gdb\\rasterUTM1x1\")\n", + "area_of_int = os.path.join(base_path, r\"newData\\Pilotregionen\\KasselStadt.shp\")\n", + "out_data = os.path.join(base_path, r\"newData\\gridUTM1x1_Ka.shp\")\n", "vect_clip = arcpy.analysis.Clip(inp_vect, area_of_int, out_data) " ] }, @@ -254,7 +255,7 @@ "source": [ "osm_list = []\n", "osm_names = []\n", - "walk = arcpy.da.Walk(r\"F:\\newData\\hesse_shp\", datatype=\"FeatureClass\")\n", + "walk = arcpy.da.Walk(os.path.join(base_path, r\"newData\\hesse_shp\", datatype=\"FeatureClass\"))\n", "\n", "for dirpath, dirnames, filenames in walk:\n", " for filename in filenames:\n", @@ -282,7 +283,7 @@ "source": [ "bbd_list = []\n", "bbd_names = []\n", - "walk = arcpy.da.Walk(r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\", datatype=\"FeatureClass\")\n", + "walk = arcpy.da.Walk(os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\", datatype=\"FeatureClass\"))\n", "\n", "for dirpath, dirnames, filenames in walk:\n", " for filename in filenames:\n", diff --git a/notebooks/data_prep/DataPrep_RheinMain.ipynb b/notebooks/data_prep/DataPrep_RheinMain.ipynb index f412006..2248f5a 100644 --- a/notebooks/data_prep/DataPrep_RheinMain.ipynb +++ b/notebooks/data_prep/DataPrep_RheinMain.ipynb @@ -31,7 +31,8 @@ "\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -52,9 +53,9 @@ }, "outputs": [], "source": [ - "raster_1914 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "area_of_interest = r\"F:\\newData\\Pilotregionen\\RheinMain.shp\"\n", - "output_file = r\"F:\\newData\\diff19_14_RheinMain.tif\"\n", + "raster_1914 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "area_of_interest = os.path.join(base_path, r\"newData\\Pilotregionen\\RheinMain.shp\")\n", + "output_file = os.path.join(base_path, r\"newData\\diff19_14_RheinMain.tif\")\n", "\n", "raster_1914_clip = Clip(raster_1914, area_of_interest) \n", "raster_1914_clip.save(output_file) " @@ -75,9 +76,9 @@ }, "outputs": [], "source": [ - "raster_2119 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", - "area_of_interest = r\"F:\\newData\\Pilotregionen\\RheinMain.shp\"\n", - "output_file = r\"F:\\newData\\diff21_19_RheinMain.tif\"\n", + "raster_2119 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", + "area_of_interest = os.path.join(base_path, r\"newData\\Pilotregionen\\RheinMain.shp\")\n", + "output_file = os.path.join(base_path, r\"newData\\diff21_19_RheinMain.tif\")\n", "\n", "raster_2119_clip = Clip(raster_2119, area_of_interest) \n", "raster_2119_clip.save(output_file) " @@ -119,7 +120,7 @@ }, "outputs": [], "source": [ - "#inp_raster = r\"F:\\newData\\diff19_14_Kassel.tif\"\n", + "#inp_raster = os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\")\n", "#min_value = -0.2\n", "#max_value = 0.2\n", "\n", @@ -156,7 +157,7 @@ "metadata": {}, "outputs": [], "source": [ - "#inp_raster = r\"F:\\newData\\diff21_19_Kassel.tif\"\n", + "#inp_raster = os.path.join(base_path, r\"newData\\diff21_19_Kassel.tif\")\n", "#min_value = -0.2\n", "#max_value = 0.2\n", "\n", @@ -188,8 +189,8 @@ }, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff19_14_RheinMain.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_RM_re.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff19_14_RheinMain.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_RM_re.tif\")\n", "\n", "arcpy.management.CopyRaster(inp_raster, out_raster)" ] @@ -200,8 +201,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff1914_RM_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_RM_re1.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff1914_RM_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_RM_re1.tif\")\n", "\n", "range_class = RemapRange([[stats_min_1914, -1,-1],[-1,1,\"NODATA\"],[1, stats_max_1914,1]])\n", "\n", @@ -215,8 +216,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff1914_RM_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_RM_re2.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff1914_RM_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_RM_re2.tif\")\n", "\n", "range_class = RemapRange([[stats_min_1914, -5,-5],[-5,-2.5,-2],[-2.5,-1.5,-1],[-1.5,-0.1,0],[-0.1,0.1,\"NODATA\"],[0.1,1.5,0],[1.5,2.5,1],[2.5,5,2],[5, stats_max_1914,5]])\n", "\n", @@ -237,8 +238,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff21_19_RheinMain.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_RM_re.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff21_19_RheinMain.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_RM_re.tif\")\n", "\n", "arcpy.management.CopyRaster(inp_raster, out_raster)" ] @@ -249,8 +250,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff2119_RM_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_RM_re1.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff2119_RM_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_RM_re1.tif\")\n", "\n", "range_class = RemapRange([[stats_min_2119, -1,-1],[-1,1,\"NODATA\"],[1, stats_max_2119,1]])\n", "\n", @@ -264,8 +265,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff2119_RM_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_RM_re2.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff2119_RM_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_RM_re2.tif\")\n", "\n", "range_class = RemapRange([[stats_min_2119, -5,-5],[-5,-2.5,-2],[-2.5,-1.5,-1],[-1.5,-0.1,0],[-0.1,0.1,\"NODATA\"],[0.1,1.5,0],[1.5,2.5,1],[2.5,5,2],[5, stats_max_2119,5]])\n", "\n", @@ -286,9 +287,9 @@ "metadata": {}, "outputs": [], "source": [ - "inp_vect = r\"F:\\newData\\testworkspace.gdb\\rasterUTM1x1\"\n", - "area_of_int = r\"F:\\newData\\Pilotregionen\\RheinMain.shp\"\n", - "out_data = r\"F:\\newData\\gridUTM1x1_RM.shp\"\n", + "inp_vect = os.path.join(base_path, r\"newData\\testworkspace.gdb\\rasterUTM1x1\")\n", + "area_of_int = os.path.join(base_path, r\"newData\\Pilotregionen\\RheinMain.shp\")\n", + "out_data = os.path.join(base_path, r\"newData\\gridUTM1x1_RM.shp\")\n", "\n", "vect_clip = arcpy.analysis.Clip(inp_vect, area_of_int, out_data) " ] @@ -299,7 +300,7 @@ "metadata": {}, "outputs": [], "source": [ - "inp_path = r\"F:\\newData\\hesse_shp\"\n", + "inp_path = os.path.join(base_path, r\"newData\\hesse_shp\")\n", "\n", "\n", "osm_list = []\n", @@ -331,7 +332,7 @@ "metadata": {}, "outputs": [], "source": [ - "inp_path = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\"\n", + "inp_path = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\")\n", "\n", "bbd_list = []\n", "bbd_names = []\n", diff --git a/notebooks/data_prep/DataPrep_TestRegionTemplate_DataPrepipynb_1.ipynb b/notebooks/data_prep/DataPrep_TestRegionTemplate_DataPrepipynb_1.ipynb index 698561c..ffe12a4 100644 --- a/notebooks/data_prep/DataPrep_TestRegionTemplate_DataPrepipynb_1.ipynb +++ b/notebooks/data_prep/DataPrep_TestRegionTemplate_DataPrepipynb_1.ipynb @@ -31,7 +31,8 @@ "\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -52,9 +53,9 @@ }, "outputs": [], "source": [ - "raster_1914 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "area_of_interest = r\"F:\\newData\\Pilotregionen\\KasselStadt.shp\"\n", - "output_file = r\"F:\\newData\\diff19_14_Kassel.tif\"\n", + "raster_1914 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "area_of_interest = os.path.join(base_path, r\"newData\\Pilotregionen\\KasselStadt.shp\")\n", + "output_file = os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\")\n", "\n", "raster_1914_clip = Clip(raster_1914, area_of_interest) \n", "raster_1914_clip.save(output_file) " @@ -75,9 +76,9 @@ }, "outputs": [], "source": [ - "raster_2119 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", - "area_of_interest = r\"F:\\newData\\Pilotregionen\\KasselStadt.shp\"\n", - "output_file = r\"F:\\newData\\diff21_19_Kassel.tif\"\n", + "raster_2119 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", + "area_of_interest = os.path.join(base_path, r\"newData\\Pilotregionen\\KasselStadt.shp\")\n", + "output_file = os.path.join(base_path, r\"newData\\diff21_19_Kassel.tif\")\n", "\n", "raster_2119_clip = Clip(raster_2119, area_of_interest) \n", "raster_2119_clip.save(output_file) " @@ -119,7 +120,7 @@ }, "outputs": [], "source": [ - "#inp_raster = r\"F:\\newData\\diff19_14_Kassel.tif\"\n", + "#inp_raster = os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\")\n", "#min_value = -0.2\n", "#max_value = 0.2\n", "\n", @@ -156,7 +157,7 @@ "metadata": {}, "outputs": [], "source": [ - "#inp_raster = r\"F:\\newData\\diff21_19_Kassel.tif\"\n", + "#inp_raster = os.path.join(base_path, r\"newData\\diff21_19_Kassel.tif\")\n", "#min_value = -0.2\n", "#max_value = 0.2\n", "\n", @@ -188,8 +189,8 @@ }, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff19_14_Kassel.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_K_re.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_K_re.tif\")\n", "\n", "arcpy.management.CopyRaster(inp_raster, out_raster)" ] @@ -200,8 +201,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff1914_K_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_K_re1.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff1914_K_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_K_re1.tif\")\n", "\n", "range_class = RemapRange([[stats_min, -1,-1],[-1,1,\"NODATA\"],[1, stats_max,1]])\n", "\n", @@ -215,8 +216,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff1914_K_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_K_re2.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff1914_K_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_K_re2.tif\")\n", "\n", "range_class = RemapRange([[stats_min, -5,-5],[-5,-2.5,-2],[-2.5,-1.5,-1],[-1.5,-0.1,0],[-0.1,0.1,\"NODATA\"],[0.1,1.5,0],[1.5,2.5,1],[2.5,5,2],[5, stats_max,5]])\n", "\n", @@ -237,8 +238,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff21_19_Kassel.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_K_re.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff21_19_Kassel.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_K_re.tif\")\n", "\n", "arcpy.management.CopyRaster(inp_raster, out_raster)" ] @@ -249,8 +250,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff2119_K_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_K_re1.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff2119_K_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_K_re1.tif\")\n", "\n", "range_class = RemapRange([[stats_min, -1,-1],[-1,1,\"NODATA\"],[1, stats_max,1]])\n", "\n", @@ -264,8 +265,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff2119_K_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_K_re2.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff2119_K_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_K_re2.tif\")\n", "\n", "range_class = RemapRange([[stats_min, -5,-5],[-5,-2.5,-2],[-2.5,-1.5,-1],[-1.5,-0.1,0],[-0.1,0.1,\"NODATA\"],[0.1,1.5,0],[1.5,2.5,1],[2.5,5,2],[5, stats_max,5]])\n", "\n", @@ -286,9 +287,9 @@ "metadata": {}, "outputs": [], "source": [ - "inp_vect = r\"F:\\newData\\testworkspace.gdb\\rasterUTM1x1\"\n", - "area_of_int = r\"F:\\newData\\Pilotregionen\\KasselStadt.shp\"\n", - "out_data = r\"F:\\newData\\gridUTM1x1_Ka.shp\"\n", + "inp_vect = os.path.join(base_path, r\"newData\\testworkspace.gdb\\rasterUTM1x1\")\n", + "area_of_int = os.path.join(base_path, r\"newData\\Pilotregionen\\KasselStadt.shp\")\n", + "out_data = os.path.join(base_path, r\"newData\\gridUTM1x1_Ka.shp\")\n", "\n", "vect_clip = arcpy.analysis.Clip(inp_vect, area_of_int, out_data) " ] @@ -299,7 +300,7 @@ "metadata": {}, "outputs": [], "source": [ - "inp_path = r\"F:\\newData\\hesse_shp\"\n", + "inp_path = os.path.join(base_path, r\"newData\\hesse_shp\")\n", "\n", "\n", "osm_list = []\n", @@ -332,7 +333,7 @@ "metadata": {}, "outputs": [], "source": [ - "inp_path = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\"\n", + "inp_path = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\")\n", "\n", "bbd_list = []\n", "bbd_names = []\n", diff --git a/notebooks/data_prep/DataPrep_TestRegion_Template.ipynb b/notebooks/data_prep/DataPrep_TestRegion_Template.ipynb index 101ebcb..e3f88ac 100644 --- a/notebooks/data_prep/DataPrep_TestRegion_Template.ipynb +++ b/notebooks/data_prep/DataPrep_TestRegion_Template.ipynb @@ -29,7 +29,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -60,7 +61,7 @@ }, "outputs": [], "source": [ - "list_names = [row[0] for row in arcpy.da.SearchCursor(r\"F:\\newData\\Pilotregionen\\TEST.shp\", \"Name\")]\n", + "list_names = [row[0] for row in arcpy.da.SearchCursor(os.path.join(base_path, r\"newData\\Pilotregionen\\TEST.shp\", \"Name\"))]\n", "print(list_names)" ] }, @@ -106,7 +107,7 @@ "source": [ "for i, j in zip(list_names, expression):\n", " out_file = i + \".shp\"\n", - " arcpy.conversion.FeatureClassToFeatureClass(r\"F:\\newData\\Pilotregionen\\TEST.shp\", r\"F:\\newData\\Pilotregionen\", out_file, j)" + " arcpy.conversion.FeatureClassToFeatureClass(os.path.join(base_path, r\"newData\\Pilotregionen\\TEST.shp\", r\"F:\\newData\\Pilotregionen\"), out_file, j)" ] }, { @@ -127,9 +128,9 @@ }, "outputs": [], "source": [ - "test_r = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "r_clip = Clip(test_r, r\"F:\\newData\\Pilotregionen\\KasselStadt.shp\") \n", - "r_clip.save(r\"F:\\newData\\diff19_14_Kassel.tif\") " + "test_r = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "r_clip = Clip(test_r, os.path.join(base_path, r\"newData\\Pilotregionen\\KasselStadt.shp\")) \n", + "r_clip.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\")) " ] }, { @@ -176,7 +177,7 @@ "outputs": [], "source": [ "fig, ax1 = pyplot.subplots(1, 1, figsize=(14,7))\n", - "src = rasterio.open(r\"F:\\newData\\diff19_14_Kassel.tif\")\n", + "src = rasterio.open(os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\"))\n", "show_hist(src, bins=10000, lw=0.0, stacked=False, histtype='bar', title=\"Histogram\", ax=ax1)\n", "ax1.set_xlim(xmin=-0.2, xmax=0.2)\n", "pyplot.show()" @@ -201,7 +202,7 @@ }, "outputs": [], "source": [ - "arcpy.management.CopyRaster(r\"Differenzen\\diff19_14_Kassel.tif\", r\"F:\\newData\\diff1914_K_re.tif\")" + "arcpy.management.CopyRaster(r\"Differenzen\\diff19_14_Kassel.tif\", os.path.join(base_path, r\"newData\\diff1914_K_re.tif\"))" ] }, { @@ -212,8 +213,8 @@ "source": [ "range_class = RemapRange([[stats_min, -1,-1],[-1,1,\"NODATA\"],[1, stats_max,1]])\n", "\n", - "r_clip_reclas = Reclassify(r\"F:\\newData\\diff1914_K_re.tif\", \"Value\", range_class)\n", - "r_clip_reclas.save(r\"F:\\newData\\diff1914_K_re1.tif\")" + "r_clip_reclas = Reclassify(os.path.join(base_path, r\"newData\\diff1914_K_re.tif\", \"Value\"), range_class)\n", + "r_clip_reclas.save(os.path.join(base_path, r\"newData\\diff1914_K_re1.tif\"))" ] }, { @@ -223,8 +224,8 @@ "outputs": [], "source": [ "range_class = RemapRange([[stats_min, -5,-5],[-5,-2.5,-2],[-2.5,-1.5,-1],[-1.5,-0.1,0],[-0.1,0.1,\"NODATA\"],[0.1,1.5,0],[1.5,2.5,1],[2.5,5,2],[5, stats_max,5]])\n", - "r_clip_reclas = Reclassify(r\"F:\\newData\\diff1914_K_re.tif\", \"Value\", range_class)\n", - "r_clip_reclas.save(r\"F:\\newData\\diff1914_K_re2.tif\")" + "r_clip_reclas = Reclassify(os.path.join(base_path, r\"newData\\diff1914_K_re.tif\", \"Value\"), range_class)\n", + "r_clip_reclas.save(os.path.join(base_path, r\"newData\\diff1914_K_re2.tif\"))" ] }, { @@ -240,9 +241,9 @@ "metadata": {}, "outputs": [], "source": [ - "inp_vect = r\"F:\\newData\\testworkspace.gdb\\rasterUTM1x1\"\n", - "area_of_int = r\"F:\\newData\\Pilotregionen\\KasselStadt.shp\"\n", - "out_data = r\"F:\\newData\\gridUTM1x1_Ka.shp\"\n", + "inp_vect = os.path.join(base_path, r\"newData\\testworkspace.gdb\\rasterUTM1x1\")\n", + "area_of_int = os.path.join(base_path, r\"newData\\Pilotregionen\\KasselStadt.shp\")\n", + "out_data = os.path.join(base_path, r\"newData\\gridUTM1x1_Ka.shp\")\n", "vect_clip = arcpy.analysis.Clip(inp_vect, area_of_int, out_data) " ] }, @@ -254,7 +255,7 @@ "source": [ "osm_list = []\n", "osm_names = []\n", - "walk = arcpy.da.Walk(r\"F:\\newData\\hesse_shp\", datatype=\"FeatureClass\")\n", + "walk = arcpy.da.Walk(os.path.join(base_path, r\"newData\\hesse_shp\", datatype=\"FeatureClass\"))\n", "\n", "for dirpath, dirnames, filenames in walk:\n", " for filename in filenames:\n", @@ -282,7 +283,7 @@ "source": [ "bbd_list = []\n", "bbd_names = []\n", - "walk = arcpy.da.Walk(r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\", datatype=\"FeatureClass\")\n", + "walk = arcpy.da.Walk(os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\", datatype=\"FeatureClass\"))\n", "\n", "for dirpath, dirnames, filenames in walk:\n", " for filename in filenames:\n", diff --git a/notebooks/data_prep/DataPrep_WerraKali.ipynb b/notebooks/data_prep/DataPrep_WerraKali.ipynb index 8f20153..1473046 100644 --- a/notebooks/data_prep/DataPrep_WerraKali.ipynb +++ b/notebooks/data_prep/DataPrep_WerraKali.ipynb @@ -31,7 +31,8 @@ "\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -52,9 +53,9 @@ }, "outputs": [], "source": [ - "raster_1914 = r\"F:\\newData\\testworkspace.gdb\\copyraster1914\"\n", - "area_of_interest = r\"F:\\newData\\Pilotregionen\\WerraKali.shp\"\n", - "output_file = r\"F:\\newData\\diff19_14_WerraKali.tif\"\n", + "raster_1914 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster1914\")\n", + "area_of_interest = os.path.join(base_path, r\"newData\\Pilotregionen\\WerraKali.shp\")\n", + "output_file = os.path.join(base_path, r\"newData\\diff19_14_WerraKali.tif\")\n", "\n", "raster_1914_clip = Clip(raster_1914, area_of_interest) \n", "raster_1914_clip.save(output_file) " @@ -75,9 +76,9 @@ }, "outputs": [], "source": [ - "raster_2119 = r\"F:\\newData\\testworkspace.gdb\\copyraster2119_2\"\n", - "area_of_interest = r\"F:\\newData\\Pilotregionen\\WerraKali.shp\"\n", - "output_file = r\"F:\\newData\\diff21_19_WerraKali.tif\"\n", + "raster_2119 = os.path.join(base_path, r\"newData\\testworkspace.gdb\\copyraster2119_2\")\n", + "area_of_interest = os.path.join(base_path, r\"newData\\Pilotregionen\\WerraKali.shp\")\n", + "output_file = os.path.join(base_path, r\"newData\\diff21_19_WerraKali.tif\")\n", "\n", "raster_2119_clip = Clip(raster_2119, area_of_interest) \n", "raster_2119_clip.save(output_file) " @@ -119,7 +120,7 @@ }, "outputs": [], "source": [ - "#inp_raster = r\"F:\\newData\\diff19_14_Kassel.tif\"\n", + "#inp_raster = os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\")\n", "#min_value = -0.2\n", "#max_value = 0.2\n", "\n", @@ -156,7 +157,7 @@ "metadata": {}, "outputs": [], "source": [ - "#inp_raster = r\"F:\\newData\\diff21_19_Kassel.tif\"\n", + "#inp_raster = os.path.join(base_path, r\"newData\\diff21_19_Kassel.tif\")\n", "#min_value = -0.2\n", "#max_value = 0.2\n", "\n", @@ -188,8 +189,8 @@ }, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff19_14_WerraKali.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_WK_re.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff19_14_WerraKali.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_WK_re.tif\")\n", "\n", "arcpy.management.CopyRaster(inp_raster, out_raster)" ] @@ -200,8 +201,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff1914_WK_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_WK_re1.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff1914_WK_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_WK_re1.tif\")\n", "\n", "range_class = RemapRange([[stats_min_1914, -1,-1],[-1,1,\"NODATA\"],[1, stats_max_1914,1]])\n", "\n", @@ -215,8 +216,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff1914_WK_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff1914_WK_re2.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff1914_WK_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff1914_WK_re2.tif\")\n", "\n", "range_class = RemapRange([[stats_min_1914, -5,-5],[-5,-2.5,-2],[-2.5,-1.5,-1],[-1.5,-0.1,0],[-0.1,0.1,\"NODATA\"],[0.1,1.5,0],[1.5,2.5,1],[2.5,5,2],[5, stats_max_1914,5]])\n", "\n", @@ -237,8 +238,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff21_19_WerraKali.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_WK_re.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff21_19_WerraKali.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_WK_re.tif\")\n", "\n", "arcpy.management.CopyRaster(inp_raster, out_raster)" ] @@ -249,8 +250,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff2119_WK_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_WK_re1.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff2119_WK_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_WK_re1.tif\")\n", "\n", "range_class = RemapRange([[stats_min_2119, -1,-1],[-1,1,\"NODATA\"],[1, stats_max_2119,1]])\n", "\n", @@ -264,8 +265,8 @@ "metadata": {}, "outputs": [], "source": [ - "inp_raster = r\"F:\\newData\\diff2119_WK_re.tif\"\n", - "out_raster = r\"F:\\newData\\diff2119_WK_re2.tif\"\n", + "inp_raster = os.path.join(base_path, r\"newData\\diff2119_WK_re.tif\")\n", + "out_raster = os.path.join(base_path, r\"newData\\diff2119_WK_re2.tif\")\n", "\n", "range_class = RemapRange([[stats_min_2119, -5,-5],[-5,-2.5,-2],[-2.5,-1.5,-1],[-1.5,-0.1,0],[-0.1,0.1,\"NODATA\"],[0.1,1.5,0],[1.5,2.5,1],[2.5,5,2],[5, stats_max_2119,5]])\n", "\n", @@ -286,9 +287,9 @@ "metadata": {}, "outputs": [], "source": [ - "inp_vect = r\"F:\\newData\\testworkspace.gdb\\rasterUTM1x1\"\n", - "area_of_int = r\"F:\\newData\\Pilotregionen\\WerraKali.shp\"\n", - "out_data = r\"F:\\newData\\gridUTM1x1_WK.shp\"\n", + "inp_vect = os.path.join(base_path, r\"newData\\testworkspace.gdb\\rasterUTM1x1\")\n", + "area_of_int = os.path.join(base_path, r\"newData\\Pilotregionen\\WerraKali.shp\")\n", + "out_data = os.path.join(base_path, r\"newData\\gridUTM1x1_WK.shp\")\n", "\n", "vect_clip = arcpy.analysis.Clip(inp_vect, area_of_int, out_data) " ] @@ -299,7 +300,7 @@ "metadata": {}, "outputs": [], "source": [ - "inp_path = r\"F:\\newData\\hesse_shp\"\n", + "inp_path = os.path.join(base_path, r\"newData\\hesse_shp\")\n", "\n", "\n", "osm_list = []\n", @@ -331,7 +332,7 @@ "metadata": {}, "outputs": [], "source": [ - "inp_path = r\"F:\\Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\"\n", + "inp_path = os.path.join(base_path, r\"Umwelt 4.0\\BBD_Update_2021_Hessen.gdb\")\n", "\n", "bbd_list = []\n", "bbd_names = []\n", diff --git a/notebooks/get_building.ipynb b/notebooks/get_building.ipynb index cf499a3..e279015 100644 --- a/notebooks/get_building.ipynb +++ b/notebooks/get_building.ipynb @@ -16,7 +16,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline " + "base_path = os.getcwd()" +"%matplotlib inline" ] }, { @@ -26,35 +27,35 @@ "outputs": [], "source": [ "\n", - "inp_fc_landuse = r\"F:\\newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\"\n", - "inp_fc_water = r\"F:\\newData\\hesse_shp\\gis_osm_water_a_free_1.shp\"\n", - "inp_fc_building = r\"F:\\newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\"\n", + "inp_fc_landuse = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\")\n", + "inp_fc_water = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_water_a_free_1.shp\")\n", + "inp_fc_building = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_buildings_a_free_1.shp\")\n", "\n", - "only_quarry = r\"F:\\newData\\hesse_osm_buffer\\only_quarry.shp\"\n", + "only_quarry = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\only_quarry.shp\")\n", "\n", - "buffer_quarry_01m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_q01m.shp\"\n", - "buffer_quarry_02m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_q02m.shp\"\n", - "buffer_quarry_05m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_q05m.shp\"\n", - "buffer_quarry_10m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_q10m.shp\"\n", - "buffer_quarry_15m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_q15m.shp\"\n", - "buffer_quarry_20m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_q20m.shp\"\n", - "buffer_quarry_25m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_q25m.shp\"\n", + "buffer_quarry_01m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_q01m.shp\")\n", + "buffer_quarry_02m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_q02m.shp\")\n", + "buffer_quarry_05m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_q05m.shp\")\n", + "buffer_quarry_10m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_q10m.shp\")\n", + "buffer_quarry_15m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_q15m.shp\")\n", + "buffer_quarry_20m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_q20m.shp\")\n", + "buffer_quarry_25m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_q25m.shp\")\n", "\n", - "buffer_water_01m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_w01m.shp\"\n", - "buffer_water_02m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_w02m.shp\"\n", - "buffer_water_05m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_w05m.shp\"\n", - "buffer_water_10m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_w10m.shp\"\n", - "buffer_water_15m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_w15m.shp\"\n", - "buffer_water_20m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_w20m.shp\"\n", - "buffer_water_25m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_w25m.shp\"\n", + "buffer_water_01m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_w01m.shp\")\n", + "buffer_water_02m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_w02m.shp\")\n", + "buffer_water_05m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_w05m.shp\")\n", + "buffer_water_10m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_w10m.shp\")\n", + "buffer_water_15m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_w15m.shp\")\n", + "buffer_water_20m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_w20m.shp\")\n", + "buffer_water_25m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_w25m.shp\")\n", "\n", - "buffer_building_01m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_b01m.shp\"\n", - "buffer_building_02m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_b02m.shp\"\n", - "buffer_building_05m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_b05m.shp\"\n", - "buffer_building_10m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_b10m.shp\"\n", - "buffer_building_15m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_b15m.shp\"\n", - "buffer_building_20m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_b20m.shp\"\n", - "buffer_building_25m = r\"F:\\newData\\hesse_osm_buffer\\hesbuf_b25m.shp\"\n", + "buffer_building_01m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_b01m.shp\")\n", + "buffer_building_02m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_b02m.shp\")\n", + "buffer_building_05m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_b05m.shp\")\n", + "buffer_building_10m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_b10m.shp\")\n", + "buffer_building_15m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_b15m.shp\")\n", + "buffer_building_20m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_b20m.shp\")\n", + "buffer_building_25m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hesbuf_b25m.shp\")\n", "\n", "distance_01 = \"1 Meters\"\n", "distance_02 = \"2 Meters\"\n", diff --git a/notebooks/get_quarry.ipynb b/notebooks/get_quarry.ipynb index eeb8346..162150f 100644 --- a/notebooks/get_quarry.ipynb +++ b/notebooks/get_quarry.ipynb @@ -16,7 +16,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline " + "base_path = os.getcwd()" +"%matplotlib inline" ] }, { @@ -25,19 +26,19 @@ "metadata": {}, "outputs": [], "source": [ - "#inp_fc_buildings = r\"F:\\newData\\hesse_shp_kassel\\gis_osm_buildings_a_free_1.shp\"\n", - "inp_fc_landuse = r\"F:\\newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\"\n", + "#inp_fc_buildings = os.path.join(base_path, r\"newData\\hesse_shp_kassel\\gis_osm_buildings_a_free_1.shp\")\n", + "inp_fc_landuse = os.path.join(base_path, r\"newData\\hesse_shp\\gis_osm_landuse_a_free_1.shp\")\n", "\n", - "#buffer_buildings = r\"F:\\newData\\hesse_shp_kassel\\buffer_buildings_10m.shp\"\n", - "only_quarry = r\"F:\\newData\\hesse_osm_buffer\\only_quarry.shp\"\n", + "#buffer_buildings = os.path.join(base_path, r\"newData\\hesse_shp_kassel\\buffer_buildings_10m.shp\")\n", + "only_quarry = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\only_quarry.shp\")\n", "\n", - "buffer_quarry_01m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_01m.shp\"\n", - "buffer_quarry_02m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_02m.shp\"\n", - "buffer_quarry_05m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_05m.shp\"\n", - "buffer_quarry_10m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_10m.shp\"\n", - "buffer_quarry_15m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_15m.shp\"\n", - "buffer_quarry_20m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_20m.shp\"\n", - "buffer_quarry_25m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_25m.shp\"\n", + "buffer_quarry_01m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_01m.shp\")\n", + "buffer_quarry_02m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_02m.shp\")\n", + "buffer_quarry_05m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_05m.shp\")\n", + "buffer_quarry_10m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_10m.shp\")\n", + "buffer_quarry_15m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_15m.shp\")\n", + "buffer_quarry_20m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_20m.shp\")\n", + "buffer_quarry_25m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_25m.shp\")\n", "\n", "distance_01 = \"1 Meters\"\n", "distance_02 = \"2 Meters\"\n", diff --git a/notebooks/maps_quarry.ipynb b/notebooks/maps_quarry.ipynb index 799b9c6..79d4034 100644 --- a/notebooks/maps_quarry.ipynb +++ b/notebooks/maps_quarry.ipynb @@ -17,7 +17,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline " + "base_path = os.getcwd()" +"%matplotlib inline" ] }, { @@ -27,63 +28,63 @@ "outputs": [], "source": [ "\n", - "only_quarry = r\"F:\\newData\\hesse_osm_buffer\\only_quarry.shp\"\n", + "only_quarry = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\only_quarry.shp\")\n", "\n", - "buffer_quarry_01m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_01m.shp\"\n", - "buffer_quarry_02m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_02m.shp\"\n", - "buffer_quarry_05m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_05m.shp\"\n", - "buffer_quarry_10m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_10m.shp\"\n", - "buffer_quarry_15m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_15m.shp\"\n", - "buffer_quarry_20m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_20m.shp\"\n", - "buffer_quarry_25m = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_25m.shp\"\n", + "buffer_quarry_01m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_01m.shp\")\n", + "buffer_quarry_02m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_02m.shp\")\n", + "buffer_quarry_05m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_05m.shp\")\n", + "buffer_quarry_10m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_10m.shp\")\n", + "buffer_quarry_15m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_15m.shp\")\n", + "buffer_quarry_20m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_20m.shp\")\n", + "buffer_quarry_25m = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_25m.shp\")\n", "\n", - "buffer_c_dis_00m = r\"F:\\newData\\hesse_osm_shp\\quarry_dis_00m.shp\"\n", + "buffer_c_dis_00m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_dis_00m.shp\")\n", "\n", - "buffer_c_dis_01m = r\"F:\\newData\\hesse_osm_shp\\quarry_dis_01m.shp\"\n", - "buffer_c_dis_02m = r\"F:\\newData\\hesse_osm_shp\\quarry_dis_02m.shp\"\n", - "buffer_c_dis_05m = r\"F:\\newData\\hesse_osm_shp\\quarry_dis_05m.shp\"\n", - "buffer_c_dis_10m = r\"F:\\newData\\hesse_osm_shp\\quarry_dis_10m.shp\"\n", - "buffer_c_dis_15m = r\"F:\\newData\\hesse_osm_shp\\quarry_dis_15m.shp\"\n", - "buffer_c_dis_20m = r\"F:\\newData\\hesse_osm_shp\\quarry_dis_20m.shp\"\n", - "buffer_c_dis_25m = r\"F:\\newData\\hesse_osm_shp\\quarry_dis_25m.shp\"\n", + "buffer_c_dis_01m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_dis_01m.shp\")\n", + "buffer_c_dis_02m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_dis_02m.shp\")\n", + "buffer_c_dis_05m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_dis_05m.shp\")\n", + "buffer_c_dis_10m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_dis_10m.shp\")\n", + "buffer_c_dis_15m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_dis_15m.shp\")\n", + "buffer_c_dis_20m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_dis_20m.shp\")\n", + "buffer_c_dis_25m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_dis_25m.shp\")\n", "\n", "#inp_frame = \"F:\\newData\\Pilotregionen\\KasselStadt\"\n", - "erased_frame_00m = r\"F:\\newData\\hesse_osm_shp\\quarry_erased_00m.shp\"\n", - "erased_frame_01m = r\"F:\\newData\\hesse_osm_shp\\quarry_erased_01m.shp\"\n", - "erased_frame_02m = r\"F:\\newData\\hesse_osm_shp\\quarry_erased_02m.shp\"\n", - "erased_frame_05m = r\"F:\\newData\\hesse_osm_shp\\quarry_erased_05m.shp\"\n", - "erased_frame_10m = r\"F:\\newData\\hesse_osm_shp\\quarry_erased_10m.shp\"\n", - "erased_frame_15m = r\"F:\\newData\\hesse_osm_shp\\quarry_erased_15m.shp\"\n", - "erased_frame_20m = r\"F:\\newData\\hesse_osm_shp\\quarry_erased_20m.shp\"\n", - "erased_frame_25m = r\"F:\\newData\\hesse_osm_shp\\quarry_erased_25m.shp\"\n", + "erased_frame_00m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_erased_00m.shp\")\n", + "erased_frame_01m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_erased_01m.shp\")\n", + "erased_frame_02m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_erased_02m.shp\")\n", + "erased_frame_05m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_erased_05m.shp\")\n", + "erased_frame_10m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_erased_10m.shp\")\n", + "erased_frame_15m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_erased_15m.shp\")\n", + "erased_frame_20m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_erased_20m.shp\")\n", + "erased_frame_25m = os.path.join(base_path, r\"newData\\hesse_osm_shp\\quarry_erased_25m.shp\")\n", "\n", "\n", - "diff_19_14 = r\"F:\\Umwelt 4.0\\DGM-Differenzenplan_2014-2019\\dgm1_diff14-19.gdb\\dgm1_diff14-19\"\n", - "diff_21_19 = r\"F:\\Umwelt 4.0\\DGM-Differenzenplan_2019-2021\\dgm1_diff19-21.gdb\\dgm1_diff19-21\"\n", + "diff_19_14 = os.path.join(base_path, r\"Umwelt 4.0\\DGM-Differenzenplan_2014-2019\\dgm1_diff14-19.gdb\\dgm1_diff14-19\")\n", + "diff_21_19 = os.path.join(base_path, r\"Umwelt 4.0\\DGM-Differenzenplan_2019-2021\\dgm1_diff19-21.gdb\\dgm1_diff19-21\")\n", "\n", - "diff_19_14_mask_00m = r\"F:\\newData\\quarry_diff\\diff19_14_mask_00m.tif\"\n", - "diff_21_19_mask_00m = r\"F:\\newData\\quarry_diff\\diff21_19_mask_00m.tif\"\n", + "diff_19_14_mask_00m = os.path.join(base_path, r\"newData\\quarry_diff\\diff19_14_mask_00m.tif\")\n", + "diff_21_19_mask_00m = os.path.join(base_path, r\"newData\\quarry_diff\\diff21_19_mask_00m.tif\")\n", "\n", - "diff_19_14_mask_01m = r\"F:\\newData\\quarry_diff\\diff19_14_mask_01m.tif\"\n", - "diff_21_19_mask_01m = r\"F:\\newData\\quarry_diff\\diff21_19_mask_01m.tif\"\n", + "diff_19_14_mask_01m = os.path.join(base_path, r\"newData\\quarry_diff\\diff19_14_mask_01m.tif\")\n", + "diff_21_19_mask_01m = os.path.join(base_path, r\"newData\\quarry_diff\\diff21_19_mask_01m.tif\")\n", "\n", - "diff_19_14_mask_02m = r\"F:\\newData\\quarry_diff\\diff19_14_mask_02m.tif\"\n", - "diff_21_19_mask_02m = r\"F:\\newData\\quarry_diff\\diff21_19_mask_02m.tif\"\n", + "diff_19_14_mask_02m = os.path.join(base_path, r\"newData\\quarry_diff\\diff19_14_mask_02m.tif\")\n", + "diff_21_19_mask_02m = os.path.join(base_path, r\"newData\\quarry_diff\\diff21_19_mask_02m.tif\")\n", "\n", - "diff_19_14_mask_05m = r\"F:\\newData\\quarry_diff\\diff19_14_mask_05m.tif\"\n", - "diff_21_19_mask_05m = r\"F:\\newData\\quarry_diff\\diff21_19_mask_05m.tif\"\n", + "diff_19_14_mask_05m = os.path.join(base_path, r\"newData\\quarry_diff\\diff19_14_mask_05m.tif\")\n", + "diff_21_19_mask_05m = os.path.join(base_path, r\"newData\\quarry_diff\\diff21_19_mask_05m.tif\")\n", "\n", - "diff_19_14_mask_10m = r\"F:\\newData\\quarry_diff\\diff19_14_mask_10m.tif\"\n", - "diff_21_19_mask_10m = r\"F:\\newData\\quarry_diff\\diff21_19_mask_10m.tif\"\n", + "diff_19_14_mask_10m = os.path.join(base_path, r\"newData\\quarry_diff\\diff19_14_mask_10m.tif\")\n", + "diff_21_19_mask_10m = os.path.join(base_path, r\"newData\\quarry_diff\\diff21_19_mask_10m.tif\")\n", "\n", - "diff_19_14_mask_15m = r\"F:\\newData\\quarry_diff\\diff19_14_mask_15m.tif\"\n", - "diff_21_19_mask_15m = r\"F:\\newData\\quarry_diff\\diff21_19_mask_15m.tif\"\n", + "diff_19_14_mask_15m = os.path.join(base_path, r\"newData\\quarry_diff\\diff19_14_mask_15m.tif\")\n", + "diff_21_19_mask_15m = os.path.join(base_path, r\"newData\\quarry_diff\\diff21_19_mask_15m.tif\")\n", "\n", - "diff_19_14_mask_20m = r\"F:\\newData\\quarry_diff\\diff19_14_mask_20m.tif\"\n", - "diff_21_19_mask_20m = r\"F:\\newData\\quarry_diff\\diff21_19_mask_20m.tif\"\n", + "diff_19_14_mask_20m = os.path.join(base_path, r\"newData\\quarry_diff\\diff19_14_mask_20m.tif\")\n", + "diff_21_19_mask_20m = os.path.join(base_path, r\"newData\\quarry_diff\\diff21_19_mask_20m.tif\")\n", "\n", - "diff_19_14_mask_25m = r\"F:\\newData\\quarry_diff\\diff19_14_mask_25m.tif\"\n", - "diff_21_19_mask_25m = r\"F:\\newData\\quarry_diff\\diff21_19_mask_25m.tif\"\n", + "diff_19_14_mask_25m = os.path.join(base_path, r\"newData\\quarry_diff\\diff19_14_mask_25m.tif\")\n", + "diff_21_19_mask_25m = os.path.join(base_path, r\"newData\\quarry_diff\\diff21_19_mask_25m.tif\")\n", "\n", "\n", "#distance_01 = \"1 Meters\"\n", diff --git a/notebooks/reclass_mask_complete.ipynb b/notebooks/reclass_mask_complete.ipynb index 97c02d9..42c8e75 100644 --- a/notebooks/reclass_mask_complete.ipynb +++ b/notebooks/reclass_mask_complete.ipynb @@ -17,7 +17,8 @@ "from rasterio.plot import show_hist, show\n", "from matplotlib import pyplot\n", "\n", - "%matplotlib inline \n" + "base_path = os.getcwd()", + "%matplotlib inline\n" ] }, { @@ -26,21 +27,21 @@ "metadata": {}, "outputs": [], "source": [ - "buffer_buildings = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_building_01m.shp\"\n", - "buffer_quarry = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_quarry_01m.shp\"\n", - "buffer_water = r\"F:\\newData\\hesse_osm_buffer\\hes_buf_water_01m.shp\"\n", + "buffer_buildings = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_building_01m.shp\")\n", + "buffer_quarry = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_quarry_01m.shp\")\n", + "buffer_water = os.path.join(base_path, r\"newData\\hesse_osm_buffer\\hes_buf_water_01m.shp\")\n", "\n", - "buffer_complete = r\"F:\\newData\\hesse_shp_kassel\\buffer_complete.shp\" # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", - "buffer_complete_dis = r\"F:\\newData\\hesse_shp_kassel\\buffer_co_dis.shp\" # auf alle !!!!!!!!!!!!!!!!!!!!!!!\n", + "buffer_complete = os.path.join(base_path, r\"newData\\hesse_shp_kassel\\buffer_complete.shp\") # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", + "buffer_complete_dis = os.path.join(base_path, r\"newData\\hesse_shp_kassel\\buffer_co_dis.shp\") # auf alle !!!!!!!!!!!!!!!!!!!!!!!\n", "\n", "#inp_frame = \"F:\\newData\\hesse_shp\\boundingBox_hesse2.shp\"\n", - "erased_frame = r\"F:\\newData\\hesse_shp_kassel\\erased_frame.shp\" # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", + "erased_frame = os.path.join(base_path, r\"newData\\hesse_shp_kassel\\erased_frame.shp\") # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", "\n", - "diff_19_14 = r\"F:\\newData\\diff19_14_Kassel.tif\" # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", - "diff_21_19 = r\"F:\\newData\\diff21_19_Kassel.tif\" # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", + "diff_19_14 = os.path.join(base_path, r\"newData\\diff19_14_Kassel.tif\") # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", + "diff_21_19 = os.path.join(base_path, r\"newData\\diff21_19_Kassel.tif\") # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", "\n", - "diff_19_14_mask = r\"F:\\newData\\diff19_14_Kassel_mask.tif\" # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", - "diff_21_19_mask = r\"F:\\newData\\diff21_19_Kassel_mask.tif\" # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", + "diff_19_14_mask = os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask.tif\") # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", + "diff_21_19_mask = os.path.join(base_path, r\"newData\\diff21_19_Kassel_mask.tif\") # auf alle !!!!!!!!!!!!!!!!!!!!!!!! ganz Hessen\n", "\n", "#distance = \"10 Meters\"" ] @@ -132,7 +133,7 @@ "outputs": [], "source": [ "#outReclass1 = Reclassify(diff_19_14_mask, \"Value\", range_class)\n", - "#outReclass1.save(r\"F:\\newData\\diff19_14_Kassel_mask_r7.tif\")" + "#outReclass1.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask_r7.tif\"))" ] }, { @@ -178,34 +179,34 @@ "outputs": [], "source": [ "outReclass1 = Reclassify(diff_19_14_mask, \"Value\", range_class01)\n", - "outReclass1.save(r\"F:\\newData\\diff19_14_Kassel_mask_r01.tif\")\n", + "outReclass1.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask_r01.tif\"))\n", "\n", "outReclass2 = Reclassify(diff_19_14_mask, \"Value\", range_class02)\n", - "outReclass2.save(r\"F:\\newData\\diff19_14_Kassel_mask_r02.tif\")\n", + "outReclass2.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask_r02.tif\"))\n", "\n", "outReclass3 = Reclassify(diff_19_14_mask, \"Value\", range_class03)\n", - "outReclass3.save(r\"F:\\newData\\diff19_14_Kassel_mask_r03.tif\")\n", + "outReclass3.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask_r03.tif\"))\n", "\n", "outReclass4 = Reclassify(diff_19_14_mask, \"Value\", range_class04)\n", - "outReclass4.save(r\"F:\\newData\\diff19_14_Kassel_mask_r04.tif\")\n", + "outReclass4.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask_r04.tif\"))\n", "\n", "outReclass5 = Reclassify(diff_19_14_mask, \"Value\", range_class05)\n", - "outReclass5.save(r\"F:\\newData\\diff19_14_Kassel_mask_r05.tif\")\n", + "outReclass5.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask_r05.tif\"))\n", "\n", "outReclass6 = Reclassify(diff_19_14_mask, \"Value\", range_class06)\n", - "outReclass6.save(r\"F:\\newData\\diff19_14_Kassel_mask_r06.tif\")\n", + "outReclass6.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask_r06.tif\"))\n", "\n", "outReclass7 = Reclassify(diff_19_14_mask, \"Value\", range_class07)\n", - "outReclass7.save(r\"F:\\newData\\diff19_14_Kassel_mask_r07.tif\")\n", + "outReclass7.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask_r07.tif\"))\n", "\n", "outReclass8 = Reclassify(diff_19_14_mask, \"Value\", range_class08)\n", - "outReclass8.save(r\"F:\\newData\\diff19_14_Kassel_mask_r08.tif\")\n", + "outReclass8.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask_r08.tif\"))\n", "\n", "outReclass9 = Reclassify(diff_19_14_mask, \"Value\", range_class09)\n", - "outReclass9.save(r\"F:\\newData\\diff19_14_Kassel_mask_r09.tif\")\n", + "outReclass9.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask_r09.tif\"))\n", "\n", "outReclass10 = Reclassify(diff_19_14_mask, \"Value\", range_class10)\n", - "outReclass10.save(r\"F:\\newData\\diff19_14_Kassel_mask_r10.tif\")" + "outReclass10.save(os.path.join(base_path, r\"newData\\diff19_14_Kassel_mask_r10.tif\"))" ] }, { -- GitLab