0_split_tifs_refactored.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Split tifs\n",
    "\n",
    "## Inputs\n",
    "This program takes the DEM_Slope.tif files and splits them into 1024x768 tifs and jpgs for use in training the model and running predictions.\n",
    "\n",
    "DEM_Slope.tif files represent a PA State Game Land area and have an area number. The variable area_str represents the number in the program.\n",
    "\n",
    "Also, each tif is either a PAN or PAS depending on the coordinate reference system (CRS) used for it. PAN is CRS 32128. PAS is 32129. Some Game Land areas have both PAN and PAS CRS .tifs, such as area 58. The variable panorpas represents whether the area tif file being processed is a PAN or a PAS.\n",
    "\n",
    "To classify files for model training purposes, the program uses a shapefile of known charcoal hearths.\n",
    "### Current version\n",
    "/home/student/charcoalhearths/data_known_points/Charcoal Hearths 8-12-2020 BPC-WC-combo-dup_removal_area39pan.shp\n",
    "\n",
    "\n",
    "\n",
    "Old versions:\n",
    "/home/student/charcoalhearths/data_known_points/Charcoal Hearths_8-11-2020_bpc.shp\n",
    "/storage/images/charcoal_hearth_hill/downloads/Weston_Uploads/Shapefiles_5-20-2020/Charcoal-Hearths.shp\n",
    "\n",
    "\n",
    "## Function\n",
    "The program takes a large tif and then divides it into rectangles of standard size 1024x768 pixels. \n",
    "\n",
    "## Output\n",
    "The program stores .tifs in:\n",
    "/storage/images/charcoal_hearth_hill/images/split_tifs/\n",
    "\n",
    "039pan0611\n",
    "039 - State Game Land area\n",
    "pan - pan or pas - the CRS\n",
    "06 is row, starting from the bottom\n",
    "11 is column, starting from the left\n",
    "\n",
    ".jpgs in:\n",
    "/storage/images/charcoal_hearth_hill/images/jpgs/\n",
    "\n",
    "The program creates a .shp file to show the layout of the split tifs. These .shp files are in: \n",
    "/storage/images/charcoal_hearth_hill/polys/\n",
    "Files are named: area_str + panorpas \"_tile_boundary.shp\"\n",
    "\n",
    "![](http://jeffblackadar.ca/charcoal_hearths/splt_tifs_1_crop.png)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# %%pycodestyle\n",
    "import os\n",
    "# Run this to make sure all of the directories exist\n",
    "data_directory = \"/home/student/charcoal_hearth_hill/\"\n",
    "if not os.path.exists(data_directory):\n",
    "    os.makedirs(data_directory)\n",
    "print(data_directory)\n",
    "data_sub_directory = os.path.join(data_directory, \"images\")\n",
    "if not os.path.exists(data_sub_directory):\n",
    "    os.makedirs(data_sub_directory)\n",
    "print(data_sub_directory)\n",
    "data_sub_sub_directory = os.path.join(data_sub_directory, \"jpgs\")\n",
    "if not os.path.exists(data_sub_sub_directory):\n",
    "    os.makedirs(data_sub_sub_directory)\n",
    "print(data_sub_sub_directory)\n",
    "data_sub_sub_directory = os.path.join(data_sub_directory, \"split_tifs\")\n",
    "if not os.path.exists(data_sub_sub_directory):\n",
    "    os.makedirs(data_sub_sub_directory)\n",
    "print(data_sub_sub_directory)\n",
    "data_sub_sub_directory = os.path.join(data_sub_directory, \"split_tif_polys\")\n",
    "if not os.path.exists(data_sub_sub_directory):\n",
    "    os.makedirs(data_sub_sub_directory)\n",
    "print(data_sub_sub_directory)\n",
    "data_sub_directory = os.path.join(data_directory, \"polys\")\n",
    "if not os.path.exists(data_sub_directory):\n",
    "    os.makedirs(data_sub_directory)\n",
    "print(data_sub_directory)\n",
    "data_sub_directory = os.path.join(data_directory, \"area_hearths\")\n",
    "if not os.path.exists(data_sub_directory):\n",
    "    os.makedirs(data_sub_directory)\n",
    "print(data_sub_directory)\n",
    "data_sub_directory = os.path.join(data_directory, \"annots\")\n",
    "if not os.path.exists(data_sub_directory):\n",
    "    os.makedirs(data_sub_directory)\n",
    "print(data_sub_directory)\n",
    "data_sub_directory = os.path.join(data_directory, \"images_training\")\n",
    "if not os.path.exists(data_sub_directory):\n",
    "    os.makedirs(data_sub_directory)\n",
    "print(data_sub_directory)\n",
    "data_sub_directory = os.path.join(data_directory, \"predictions\")\n",
    "if not os.path.exists(data_sub_directory):\n",
    "    os.makedirs(data_sub_directory)\n",
    "print(data_sub_directory)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import geopandas as gpd\n",
    "#https://stackoverflow.com/questions/2922532/obtain-latitude-and-longitude-from-a-geotiff-file\n",
    "from osgeo import osr, gdal\n",
    "\n",
    "def get_poly_with_lat_long_from_geotif(geotif_fp):\n",
    "    # Import necessary modules\n",
    "\n",
    "\n",
    "    # get the existing coordinate system\n",
    "    ds = gdal.Open(geotif_fp)\n",
    "\n",
    "    old_cs= osr.SpatialReference()\n",
    "\n",
    "    old_cs.ImportFromWkt(ds.GetProjectionRef())\n",
    "\n",
    "    # create the new coordinate system\n",
    "    nad83_wkt = \"\"\"\n",
    "    GEOGCS[\"NAD83\",\n",
    "        DATUM[\"North_American_Datum_1983\",\n",
    "            SPHEROID[\"GRS 1980\",6378137,298.257222101,\n",
    "                AUTHORITY[\"EPSG\",\"7019\"]],\n",
    "            TOWGS84[0,0,0,0,0,0,0],\n",
    "            AUTHORITY[\"EPSG\",\"6269\"]],\n",
    "        PRIMEM[\"Greenwich\",0,\n",
    "            AUTHORITY[\"EPSG\",\"8901\"]],\n",
    "        UNIT[\"degree\",0.0174532925199433,\n",
    "            AUTHORITY[\"EPSG\",\"9122\"]],\n",
    "        AUTHORITY[\"EPSG\",\"4269\"]]\"\"\"\n",
    "\n",
    "    new_cs = osr.SpatialReference()\n",
    "\n",
    "    new_cs .ImportFromWkt(nad83_wkt)\n",
    "\n",
    "    # create a transform object to convert between coordinate systems\n",
    "    transform = osr.CoordinateTransformation(old_cs,new_cs) \n",
    "\n",
    "    #get the point to transform, pixel (0,0) in this case\n",
    "    width = ds.RasterXSize\n",
    "    height = ds.RasterYSize\n",
    "    gt = ds.GetGeoTransform()\n",
    "    minx = gt[0]\n",
    "    miny = gt[3] + width*gt[4] + height*gt[5] \n",
    "    maxx = gt[0] + width*gt[1] + height*gt[2]\n",
    "    maxy = gt[3]\n",
    "\n",
    "    #get the coordinates in lat long\n",
    "    latlong1 = transform.TransformPoint(minx,miny)\n",
    "    #print(latlong1)\n",
    "    latlong2 = transform.TransformPoint(minx,maxy)\n",
    "    #print(latlong2)\n",
    "    latlong3 = transform.TransformPoint(maxx,maxy)\n",
    "    #print(latlong3)\n",
    "    latlong4 = transform.TransformPoint(maxx,miny)\n",
    "    #print(latlong4)\n",
    "    \n",
    "    # Create an empty geopandas GeoDataFrame\n",
    "    polydata = gpd.GeoDataFrame()\n",
    "    polydata['geometry'] = None\n",
    "\n",
    "    # Create a Polygon\n",
    "    #print(latlong1[0], latlong1[1]), (latlong2[0], latlong2[1]), (latlong3[0], latlong3[1]), (latlong4[0], latlong4[1])\n",
    "    coords = [(latlong1[0], latlong1[1]), (latlong2[0], latlong2[1]), (latlong3[0], latlong3[1]), (latlong4[0], latlong4[1])]\n",
    "    poly = Polygon(coords)\n",
    "\n",
    "    return(poly)\n",
    "\n",
    "import numpy as np\n",
    "\n",
    "def tif_to_jpg(tif_file,output_file_name, output_folder):\n",
    "    #print(tif_file,output_file_name, output_folder)\n",
    "    ds = gdal.Open(tif_file)\n",
    "    #geoTrans = srcImage.GetGeoTransform()\n",
    "\n",
    "    band = ds.GetRasterBand(1)\n",
    "    width = ds.RasterXSize\n",
    "    height = ds.RasterYSize\n",
    "\n",
    "    data = band.ReadAsArray(0, 0, width, height)\n",
    "    #convert all the bad data\n",
    "    data[data==-9999.0] = 0\n",
    "    max_value = numpy.max(data)\n",
    "    color_multiplier = 255/(max_value-1)\n",
    "    #print(color_multiplier)\n",
    "    data = data*color_multiplier\n",
    "    #print(data)\n",
    "    data_int = np.array(data, dtype='int')\n",
    "    #print(data_int)\n",
    "    #clip = ds.readasarray(ds)\n",
    "    data_int = data_int.astype(gdalnumeric.numpy.uint8)\n",
    "\n",
    "    gdalnumeric.SaveArray(data_int, output_folder+output_file_name+\".jpg\", format=\"JPEG\")\n",
    "\n",
    "#https://pcjericks.github.io/py-gdalogr-cookbook/raster_layers.html?highlight=rasterize\n",
    "\n",
    "# Run this\n",
    "def write_annot(obj_type, obj_annots_dir, obj_f_num, org_f_path,org_f_width,org_f_height,org_f_depth, refPts):\n",
    "    # This generates the xml file annotation that is consumed by the MASK R-CNN process.\n",
    "    # With credit to: # https://www.geeksforgeeks.org/reading-writing-text-files-python/\n",
    "    annot_file_path = obj_annots_dir+obj_f_num+'.xml'\n",
    "    annot_file = open(annot_file_path,\"w\") \n",
    "\n",
    "    annot_file.write(\"<annotation>\\n\") \n",
    "    annot_file.write(\"\t<folder>\"+obj_type+\"</folder>\\n\") \n",
    "    annot_file.write(\"\t<filename>\"+obj_f_num+\".jpg</filename>\\n\") \n",
    "    annot_file.write(\"\t<path>\"+org_f_path+\"</path>\\n\") \n",
    "    annot_file.write(\"\t<source>\\n\") \n",
    "    annot_file.write(\"\t\t<database>Muhlenberg_charcoal_hearths</database>\\n\") \n",
    "    annot_file.write(\"\t</source>\\n\") \n",
    "    annot_file.write(\"\t<size>\\n\") \n",
    "    annot_file.write(\"\t\t<width>\"+str(org_f_width)+\"</width>\\n\") \n",
    "    annot_file.write(\"\t\t<height>\"+str(org_f_height)+\"</height>\\n\") \n",
    "    annot_file.write(\"\t\t<depth>\"+str(org_f_depth)+\"</depth>\\n\") \n",
    "    annot_file.write(\"\t</size>\\n\") \n",
    "    annot_file.write(\"\t<segmented>0</segmented>\\n\") \n",
    "    for ocn in range(0,len(refPts)):\n",
    "        refPt = refPts[ocn]\n",
    "        refPtMin = refPt[0]\n",
    "        refPtMax = refPt[1]\n",
    "        \n",
    "        #Sometimes the mouse gets dragged from the bottom to the top, etc\n",
    "        if(refPtMin[1]>refPtMax[1]):\n",
    "            ytemp = refPtMin[1]\n",
    "            refPtMin[1] = refPtMax[1]\n",
    "            refPtMax[1]=ytemp\n",
    "        if(refPtMin[0]>refPtMax[0]):\n",
    "            xtemp = refPtMin[0]\n",
    "            refPtMin[0] = refPtMax[0]\n",
    "            refPtMax[0]=xtemp \n",
    "        \n",
    "        annot_file.write(\"\t<object>\\n\") \n",
    "        annot_file.write(\"\t\t<name>\"+obj_type+\"</name>\\n\") \n",
    "        annot_file.write(\"\t\t<number>\"+str(ocn)+\"</number>\\n\") \n",
    "        annot_file.write(\"\t\t<truncated>0</truncated>\\n\") \n",
    "        annot_file.write(\"\t\t<difficult>0</difficult>\\n\") \n",
    "        annot_file.write(\"\t\t<bndbox>\\n\") \n",
    "        annot_file.write(\"\t\t\t<xmin>\"+str(refPtMin[0])+\"</xmin>\\n\") \n",
    "        annot_file.write(\"\t\t\t<ymin>\"+str(refPtMin[1])+\"</ymin>\\n\") \n",
    "        annot_file.write(\"\t\t\t<xmax>\"+str(refPtMax[0])+\"</xmax>\\n\") \n",
    "        annot_file.write(\"\t\t\t<ymax>\"+str(refPtMax[1])+\"</ymax>\\n\") \n",
    "        annot_file.write(\"\t\t</bndbox>\\n\") \n",
    "        annot_file.write(\"\t</object>\\n\") \n",
    "    annot_file.write(\"</annotation>\\n\") \n",
    "    annot_file.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import gdal\n",
    "from osgeo import ogr\n",
    "def get_pixel_with_geot(geoT, point_x, point_y):\n",
    "    # print(\"get_pixel_with_geot(geoT, point_x, point_y):\",point_x, point_y)\n",
    "    # point_x, point_y MUST BE UTM\n",
    "    \"\"\"Get the pixel for given coordinates (in the raster's convention, not\n",
    "    checked!) for a raster file.\n",
    "\n",
    "    Parameters\n",
    "    ----------\n",
    "    raster: string\n",
    "        A GDAL-friendly raster filename\n",
    "    point_x: float\n",
    "        The Easting in the same coordinates as the raster (not checked!)\n",
    "    point_y: float\n",
    "        The Northing in the same coordinates as the raster (not checked!)\n",
    "\n",
    "    Returns\n",
    "    -------\n",
    "    The row/column (or column/row, depending on how you define it)\n",
    "    \"\"\"\n",
    "    #g = gdal.Open(raster)\n",
    "    #if g is None:\n",
    "    #    raise ValueError(f\"{raster:s} cannot be opened!\")\n",
    "    #geoT = open_raster.GetGeoTransform()\n",
    "    inv_geoT = gdal.InvGeoTransform(geoT)\n",
    "    r, c = (gdal.ApplyGeoTransform(inv_geoT, point_x, point_y))\n",
    "    # print(\"_geot r\", int(r + 0.5))\n",
    "    # print(\"_geot c\", int(c + 0.5))\n",
    "    return int(r + 0.5), int(c + 0.5)\n",
    "\n",
    "def annotate_tif_if_it_has_hearths(area_hearth_data, tif_poly, tif_fp, annot_file_name, geotif_crs):\n",
    "    # Create an empty geopandas GeoDataFrame\n",
    "    # Reworked for new CRS defintions\n",
    "    hearth_data = gpd.GeoDataFrame(columns=[\"geometry\"], crs=geotif_crs)\n",
    "    hearth_data = hearth_data.to_crs(crs=geotif_crs)\n",
    "    hearth_data.geometry = hearth_data.geometry.to_crs(crs=geotif_crs)\n",
    "    # Create an empty geopandas GeoDataFrame\n",
    "    hearth_data_poly = gpd.GeoDataFrame(columns=[\"geometry\"], crs=geotif_crs)\n",
    "    hearth_data_poly = hearth_data_poly.to_crs(crs=geotif_crs)\n",
    "    hearth_data_poly.geometry = hearth_data_poly.geometry.to_crs(crs=geotif_crs)\n",
    "    \n",
    "\n",
    "    selection = area_hearth_data[0:]\n",
    "    point_count=0\n",
    "    \n",
    "    print(tif_fp)\n",
    "    src_ds=gdal.Open(tif_fp)\n",
    "    gt=src_ds.GetGeoTransform()\n",
    "    # print(\"gt\",gt)\n",
    "    \n",
    "    #split_tif = gpd.read_file(tif_fp)\n",
    "    #split_tif_4326 = split_tif.to_crs(crs=4326)\n",
    "    \n",
    "\n",
    "\n",
    "    \n",
    "    \n",
    "    # create a transform object to convert between coordinate systems\n",
    "    #transform = osr.CoordinateTransformation(old_cs,new_cs)\n",
    "    \n",
    "    target = osr.SpatialReference(wkt = src_ds.GetProjection())\n",
    "    #print(\"src_ds.GetProjection()\", src_ds.GetProjection())\n",
    "    #print(\"geotif_crs\",geotif_crs) \n",
    "    source = osr.SpatialReference()\n",
    "    #source.ImportFromEPSG(4326)\n",
    "    #source.ImportFromEPSG(32128)\n",
    "    #target.ImportFromEPSG(32128)\n",
    "    source.ImportFromEPSG(geotif_crs)\n",
    "    target.ImportFromEPSG(geotif_crs)\n",
    "    \n",
    "    transform = osr.CoordinateTransformation(source, target)\n",
    "\n",
    "    for index, row in selection.iterrows():\n",
    "        #print(\"for index, row in selection.iterrows():\")\n",
    "        pt = row['geometry']\n",
    "        #print(row)\n",
    "        #select only points inside the bounds of the geotiff AND are confirmed to be hearths\n",
    "        if(pt.within(tif_poly)==True):\n",
    "            #print(\"if(pt.within(tif_poly)==True):\")\n",
    "            point_count=point_count+1\n",
    "            #print(pt, pt.within(tif_poly))\n",
    "            hearth_data = hearth_data.append(row, ignore_index=True)\n",
    "            \n",
    "            #polygon_half_size = 0.0000895078\n",
    "            # polygon_half_size = 0.0001745402\n",
    "            polygon_half_size = 0.000175\n",
    "            # make the size of the annotated shape smaller, tighter to the hearth\n",
    "            #polygon_half_size = 0.00014\n",
    "            # This did not seem to work.  The smaller, tighter size removed the smoother surrounding, which is part of the object\n",
    "            #get the coordinates in lat long\n",
    "            # pp_ = point_poly - the polygon around the points\n",
    "            pp_minx=pt.bounds[0]-polygon_half_size\n",
    "            pp_maxx=pt.bounds[0]+polygon_half_size\n",
    "            pp_miny=pt.bounds[1]-(polygon_half_size * .8)\n",
    "            pp_maxy=pt.bounds[1]+(polygon_half_size * .8)\n",
    "            \n",
    "            pp_latlong1 = transform.TransformPoint(pp_minx,pp_miny)\n",
    "            #print(\"latlong1\",latlong1)\n",
    "            #print(\"pp_minx\",pp_minx)\n",
    "            pp_latlong2 = transform.TransformPoint(pp_minx,pp_maxy)\n",
    "            #print(latlong2)\n",
    "            pp_latlong3 = transform.TransformPoint(pp_maxx,pp_maxy)\n",
    "            #print(latlong3)\n",
    "            pp_latlong4 = transform.TransformPoint(pp_maxx,pp_miny)\n",
    "            #print(latlong4)\n",
    "            # Create a Polygon\n",
    "            pp_coords = [(pp_minx, pp_miny), (pp_minx, pp_maxy), (pp_maxx, pp_maxy), (pp_maxx, pp_miny)]\n",
    "            ppoly = Polygon(pp_coords)\n",
    "            new_pp_row = {'id':row['id'], 'geometry':ppoly}\n",
    "            #hearth_data_poly.loc[0, 'geometry'] = ppoly\n",
    "            #hearth_data_poly.loc[0, 'id'] = row['id']\n",
    "            hearth_data_poly = hearth_data_poly.append(new_pp_row, ignore_index=True)\n",
    "    #Do we have points inside this tif?  If so, annotate and copy>0\n",
    "    if(point_count)>0:\n",
    "\n",
    "        #hearth_data_poly = hearth_data_poly.to_crs(epsg=32128)\n",
    "        #hearth_data = hearth_data.to_crs(epsg=32128)\n",
    "        print(\"geotif_crs\",geotif_crs)\n",
    "        print(hearth_data_poly)\n",
    "        print(hearth_data_poly.crs)\n",
    "        \n",
    "        #hearth_data_poly.crs = (\"EPSG:\" + str(geotif_crs))\n",
    "        #hearth_data_poly.geometry = hearth_data_poly.geometry.to_crs(crs=geotif_crs)\n",
    "        #hearth_data_poly.to_crs(crs=geotif_crs)\n",
    "        #hearth_data_poly = hearth_data_poly.to_crs(epsg=geotif_crs)\n",
    "        print(\"-------\")\n",
    "        #geotif_crs_4326 = 4326\n",
    "        #hearth_data.crs = (\"EPSG:\" + str(geotif_crs))\n",
    "        #hearth_data.geometry = hearth_data.geometry.to_crs(crs=geotif_crs)\n",
    "        #hearth_data.to_crs(crs=geotif_crs)\n",
    "        #hearth_data = hearth_data.to_crs(epsg=geotif_crs)        \n",
    "\n",
    "    \n",
    "        #Write out shp files for diagnostics\n",
    "        outfp_folder = tif_fp[:-25]\n",
    "        outfp_base = tif_fp[-14:-4]\n",
    "        outfp = os.path.join(outfp_folder,\"split_tif_polys/\", (outfp_base + \"_hearth_data_poly.shp\"))\n",
    "        print(outfp_folder, outfp_base, outfp)\n",
    "        # Write the data into that Shapefile\n",
    "        hearth_data_poly.to_file(outfp)        \n",
    "        outfp = os.path.join(outfp_folder,\"split_tif_polys/\", (outfp_base + \"_hearth_data.shp\"))\n",
    "        # Write the data into that Shapefile\n",
    "        hearth_data.to_file(outfp)\n",
    "        print(\"****\")\n",
    "\n",
    "        #Get the dimentions of column and row\n",
    "        nc   = src_ds.RasterXSize\n",
    "        nr   = src_ds.RasterYSize\n",
    "        print(\"src_ds.RasterXSize\", nc,nr)\n",
    "        gdal.UseExceptions() #so it doesn't print to screen everytime point is outside grid\n",
    "        \n",
    "        #target = osr.SpatialReference(wkt=src_ds.GetProjection())\n",
    "        #source = osr.SpatialReference()\n",
    "        #source.ImportFromEPSG(4326)\n",
    "        #transform = osr.CoordinateTransformation(source, target)\n",
    "        refPts = []      \n",
    "        refPt = []\n",
    "        #I'm hoping the tifs will be 1024X768        \n",
    "        # Make annotation file for the jpeg from polygons\n",
    "        for index, row in hearth_data_poly.iterrows():\n",
    "            poly = row['geometry']\n",
    "            #print(row, row['id'])\n",
    "            #print(poly)\n",
    "            \n",
    "            boundary = poly.bounds\n",
    "            #print(\"boundary = poly.bounds\")\n",
    "            #print(boundary[0])\n",
    "            #print(boundary[1])\n",
    "            #print(boundary[2])\n",
    "            #print(boundary[3])\n",
    "\n",
    "            lon = boundary[0]\n",
    "            lat = boundary[1]\n",
    "    \n",
    "            print(\"lon,lat\",lon, lat)\n",
    "      \n",
    "            point = ogr.Geometry(ogr.wkbPoint)\n",
    "            point.AddPoint(lon, lat)\n",
    "            point.Transform(transform)\n",
    "            # Thanks to\n",
    "            # https://gis.stackexchange.com/questions/167069/simple-gdal-coordinate-conversion\n",
    "            # convert point to UTM\n",
    "            inSR = osr.SpatialReference()\n",
    "            inSR.ImportFromEPSG(4326)       # WGS84/Geographic\n",
    "            outSR = osr.SpatialReference()\n",
    "            outSR.ImportFromEPSG(geotif_crs)\n",
    "\n",
    "            point.AssignSpatialReference(inSR)    # tell the point what coordinates it's in\n",
    "            point.TransformTo(outSR)              # project it to the out spatial reference\n",
    "            print('{0},{1}'.format(point.GetX(),point.GetY())) # output projected X and Y coordinates\n",
    "            \n",
    "\n",
    "            #p1x,p1y = world_to_pixel(ds.GetGeoTransform(), point.GetX(), point.GetY())\n",
    "            p1x,p1y = get_pixel_with_geot(gt, point.GetX(), point.GetY())\n",
    "            # 2020 08 12            \n",
    "            print(\"p1\",point.GetX(), point.GetY(),p1x,p1y)\n",
    "        \n",
    "            #cx,cy=convert_coordinates(lon, lat)\n",
    "            #p1x,p1y= (get_pixel_with_geot(gt,cx,cy ))\n",
    "            lon = boundary[2]\n",
    "            lat = boundary[3]\n",
    "    \n",
    "            print(\"lon,lat\",lon, lat)\n",
    "        \n",
    "            #cx,cy=convert_coordinates(lon, lat)\n",
    "            #p2x,p2y= (get_pixel_with_geot(gt,cx,cy ))\n",
    "            point = ogr.Geometry(ogr.wkbPoint)\n",
    "            point.AddPoint(lon, lat)\n",
    "            point.Transform(transform)\n",
    "            point.AssignSpatialReference(inSR)    # tell the point what coordinates it's in\n",
    "            point.TransformTo(outSR)              # project it to the out spatial reference\n",
    "            print('{0},{1}'.format(point.GetX(),point.GetY())) # output projected X and Y coordinates\n",
    "\n",
    "            #p2x,p2y = world_to_pixel(ds.GetGeoTransform(), point.GetX(), point.GetY())\n",
    "            p2x,p2y = get_pixel_with_geot(gt, point.GetX(), point.GetY())\n",
    "            # 2020 08 12\n",
    "            \n",
    "            # print(\"p2\",point.GetX(), point.GetY(),p2x,p2y)\n",
    "\n",
    "            #p2y is smaller than p1y so we'll exchange them to keep p1 as the minimum.\n",
    "            ptemp=p1y\n",
    "            p1y=p2y\n",
    "            p2y=ptemp\n",
    "            if(p1x<0):\n",
    "                p1x=0\n",
    "            if(p1y<0):\n",
    "                p1y=0\n",
    "            if(p2x<0):\n",
    "                p2x=0\n",
    "            if(p2y<0):\n",
    "                p2y=0\n",
    "\n",
    "            if(p1x>nc):\n",
    "                p1x=nc\n",
    "            if(p1y>nr):\n",
    "                p1y=nr\n",
    "            if(p2x>nc):\n",
    "                p2x=nc\n",
    "            if(p2y>nr):\n",
    "                p2y=nr\n",
    "     \n",
    "        \n",
    "            refPt = [(p1x, p1y)]\n",
    "            refPt.append((p2x, p2y))\n",
    "            refPts.append(refPt)\n",
    "\n",
    "        object_type = 'charcoal_hearth_hill'\n",
    "        \n",
    "        object_annots_dir = '/home/student/charcoal_hearth_hill/annots/'\n",
    "        write_annot(object_type, object_annots_dir, annot_file_name, annot_file_name+\".jpg\",nc,nr,1, refPts)\n",
    "        \n",
    "        images_training_dir = '/home/student/charcoal_hearth_hill/images_training/'\n",
    "        tif_to_jpg(tif_fp,annot_file_name,images_training_dir)\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Thanks to\n",
    "#https://gis.stackexchange.com/questions/92207/split-a-large-geotiff-into-smaller-regions-with-python-and-gdal\n",
    "import os  \n",
    "import numpy\n",
    "#import gdal\n",
    "from osgeo import gdal, osr\n",
    "import math\n",
    "from itertools import chain\n",
    "#import pandas as pd\n",
    "import geopandas as gpd\n",
    "from shapely.geometry import Point, Polygon\n",
    "import numpy as np\n",
    "import gdalnumeric\n",
    "\n",
    "def get_extent(dataset):\n",
    "\n",
    "    cols = dataset.RasterXSize\n",
    "    rows = dataset.RasterYSize\n",
    "    transform = dataset.GetGeoTransform()\n",
    "    minx = transform[0]\n",
    "    maxx = transform[0] + cols * transform[1] + rows * transform[2]\n",
    "\n",
    "    miny = transform[3] + cols * transform[4] + rows * transform[5]\n",
    "    maxy = transform[3]\n",
    "\n",
    "    return {\n",
    "            \"minX\": str(minx), \"maxX\": str(maxx),\n",
    "            \"minY\": str(miny), \"maxY\": str(maxy),\n",
    "            \"cols\": str(cols), \"rows\": str(rows)\n",
    "            }\n",
    "\n",
    "#def create_tiles(minx, miny, maxx, maxy, nw, nh):\n",
    "#    width = maxx - minx\n",
    "#    height = maxy - miny\n",
    "\n",
    "#    matrix = []\n",
    "\n",
    "#    for j in range(nw, 0, -1):\n",
    "#        #changed range(0, (nh)) \n",
    "#        for i in range(1, (nh + 1)):\n",
    "\n",
    "#            ulx = minx + (width/nw) * i # 10/5 * 1\n",
    "#            uly = miny + (height/nh) * j # 10/5 * 1\n",
    "#\n",
    "#            lrx = minx + (width/nw) * (i + 1)\n",
    "#            lry = miny + (height/nh) * (j - 1)\n",
    "#            matrix.append([[ulx, uly], [lrx, lry]])\n",
    "#\n",
    "#    return matrix\n",
    "\n",
    "\n",
    "\n",
    "def create_standard_tiles(minx, miny, maxx, maxy):\n",
    "    print(\"standard\",minx, miny, maxx, maxy)\n",
    "    width = maxx - minx\n",
    "    height = maxy - miny\n",
    "    \n",
    "    standard_width = 1024\n",
    "    standard_height = 768\n",
    "   \n",
    "    nw = int(width/standard_width)+1\n",
    "    nh = int(height/standard_height)+1\n",
    "    print(nw,nh)\n",
    "    #print(nw*standard_width,nh*standard_height)\n",
    "    \n",
    "    matrix = []\n",
    "\n",
    "    # rows\n",
    "    for tifrows in range(nh, 0, -1):\n",
    "        # columns\n",
    "        #changed range(0, (nh)) \n",
    "        for tifcols in range(0, (nw)):\n",
    "\n",
    "            ulx = minx + (standard_width * tifcols) # 10/5 * 1\n",
    "            uly = miny + (standard_height * tifrows) # 10/5 * 1\n",
    "\n",
    "            lrx = minx + (standard_width * (tifcols + 1))\n",
    "            lry = miny + (standard_height * (tifrows - 1))\n",
    "            if(lrx>maxx):\n",
    "                lrx=maxx\n",
    "            if(lrx<minx):\n",
    "                lrx=minx\n",
    "              \n",
    "            if(lry>maxy):\n",
    "                lry=maxy\n",
    "            if(lry<miny):\n",
    "                lry=miny\n",
    "\n",
    "            if(ulx>maxx):\n",
    "                ulx=maxx\n",
    "            if(ulx<minx):\n",
    "                ulx=minx\n",
    "\n",
    "            if(uly>maxy):\n",
    "                uly=maxy\n",
    "            if(uly<miny):\n",
    "                uly=miny\n",
    "\n",
    "            #print(\"matrix:\",ulx, uly,lrx, lry)\n",
    "            # tifcols starts from 0\n",
    "            matrix.append([[ulx, uly], [lrx, lry], [(tifcols),tifrows]])\n",
    "    print(\"standard\",minx, miny, maxx, maxy)\n",
    "    return matrix\n",
    "\n",
    "\n",
    "def split(file_name, geotif_crs):\n",
    "    data_directory = \"/home/student/charcoal_hearth_hill/\"\n",
    "    jpgs_output_folder = os.path.join(data_directory,\"images/\",\"jpgs/\")\n",
    "    polys_output_folder = os.path.join(data_directory,\"polys/\")\n",
    "\n",
    "    #Some areas have both a PAN and a PAS, so we need to distinguish them\n",
    "    panorpas=\"\"\n",
    "    if(geotif_crs==32129):\n",
    "        panorpas=\"pas\"\n",
    "    if(geotif_crs==32128):\n",
    "        panorpas=\"pan\"\n",
    "\n",
    "    if(os.path.exists(file_name)):\n",
    "        from osgeo import osr\n",
    "        raw_file_name = os.path.splitext(os.path.basename(file_name))[0].replace(\"_downsample\", \"\")\n",
    "        print(raw_file_name)\n",
    "        area_num = raw_file_name[:3]\n",
    "        if(area_num[-1:] ==\"b\"):\n",
    "            area_num = area_num[:-1]\n",
    "        area_num = int(area_num)\n",
    "        area_str = \"0\"+str(area_num)\n",
    "        area_str = area_str[-3:]\n",
    "        print(area_str)\n",
    "        driver = gdal.GetDriverByName('GTiff')\n",
    "        dataset = gdal.Open(file_name)\n",
    "        proj_dataset = osr.SpatialReference(wkt=dataset.GetProjection())\n",
    "        #print(proj_dataset.GetAttrValue('AUTHORITY',1),geotif_crs)\n",
    "        \n",
    "        # Check that the geotif_crs matches what is in the file.\n",
    "        if(proj_dataset.GetAttrValue('AUTHORITY',1)==str(geotif_crs)):\n",
    "            print(proj_dataset)\n",
    "   \n",
    "            band = dataset.GetRasterBand(1)\n",
    "            transform = dataset.GetGeoTransform()\n",
    "\n",
    "            extent = get_extent(dataset)\n",
    "\n",
    "            cols = int(extent[\"cols\"])\n",
    "            rows = int(extent[\"rows\"])\n",
    "            print (\"Columns: \", cols)\n",
    "            print (\"Rows: \", rows)\n",
    "\n",
    "            minx = float(extent[\"minX\"])\n",
    "            maxx = float(extent[\"maxX\"])\n",
    "            miny = float(extent[\"minY\"])\n",
    "            maxy = float(extent[\"maxY\"])\n",
    "        \n",
    "            width = int(maxx - minx)\n",
    "            height = int(maxy - miny)\n",
    "\n",
    "            nw = int(width/1024)+1\n",
    "            nh = int(height/768)+1\n",
    "    \n",
    "            #******\n",
    "            # make a polygon for the area and see what points are in the whole area.\n",
    "            # get the coordinates in lat long\n",
    "\n",
    "            # Create an empty geopandas GeoDataFrame\n",
    "            area_hearth_data = gpd.GeoDataFrame()\n",
    "            area_hearth_data['geometry'] = None\n",
    "            area_poly = get_poly_with_lat_long_from_geotif(file_name)\n",
    "            # Read file using gpd.read_file()\n",
    "            #all_hearth_data = gpd.read_file(\"/storage/images/charcoal_hearth_hill/downloads/Weston_Uploads/Shapefiles_5-20-2020/Charcoal-Hearths.shp\")\n",
    "            all_hearth_data = gpd.read_file(\"/home/student/charcoalhearths/data_known_points/Charcoal Hearths 8-12-2020 BPC-WC-combo-dup_removal_area39pan.shp\")\n",
    "            selection = all_hearth_data[0:]\n",
    "            #print(selection)\n",
    "            for index, row in selection.iterrows():\n",
    "                #print(\"for index, row in selection.iterrows():\")\n",
    "                pt = row['geometry']\n",
    "                #print(row)\n",
    "                #select only points inside the bounds of the geotiff AND are confirmed to be hearths\n",
    "                if(pt.within(area_poly)==True):\n",
    "                    #print(\"if(pt.within(tif_poly)==True):\")\n",
    "                    area_hearth_data = area_hearth_data.append(row, ignore_index=True)\n",
    "    \n",
    "            if(not area_hearth_data.empty):\n",
    "            # Determine the output path for the Shapefile\n",
    "                outfp = os.path.join(data_directory, \"area_hearths\", \"area_\" + area_str + panorpas + \"_hearth.shp\")\n",
    "                # Write the data into that Shapefile\n",
    "                area_hearth_data.to_file(outfp)\n",
    "            #*************************\n",
    "            #output_path = os.path.join(\"data\", raw_file_name)\n",
    "            #if not os.path.exists(output_path):\n",
    "            #    os.makedirs(output_path)\n",
    "\n",
    "            # Polygons to reflect the tiles\n",
    "            # Create an empty geopandas GeoDataFrame\n",
    "            tile_boundary_poly = gpd.GeoDataFrame()\n",
    "            tile_boundary_poly['geometry'] = None\n",
    "            #tile_boundary_poly.crs = {'init':'epsg:4326'}\n",
    "            #geotif_crs_4326 = 4326\n",
    "            tile_boundary_poly.crs = (\"EPSG:\" + str(geotif_crs))\n",
    "            tile_boundary_poly.geometry = tile_boundary_poly.geometry.to_crs(crs=geotif_crs)\n",
    "            tile_boundary_poly.to_crs(crs=geotif_crs)\n",
    "            tile_boundary_poly = tile_boundary_poly.to_crs(epsg=geotif_crs)\n",
    "\n",
    "            #tiles = create_tiles(minx, miny, maxx, maxy, nw,nh)\n",
    "            tiles = create_standard_tiles(minx, miny, maxx, maxy)\n",
    "            transform = dataset.GetGeoTransform()\n",
    "            xOrigin = transform[0]\n",
    "            yOrigin = transform[3]\n",
    "            pixelWidth = transform[1]\n",
    "            pixelHeight = -transform[5]\n",
    "\n",
    "            srs = osr.SpatialReference()\n",
    "            #srs.ImportFromEPSG(32128)\n",
    "            srs.ImportFromEPSG(geotif_crs)\n",
    "            #srs.ImportFromEPSG(4326)\n",
    "            #dst_ds.SetProjection( srs.ExportToWkt() )\n",
    "    \n",
    "            tile_num = 0\n",
    "            for tile in tiles:\n",
    "                #print (tile)\n",
    "\n",
    "                minx = tile[0][0]\n",
    "                maxx = tile[1][0]\n",
    "                miny = tile[1][1]\n",
    "                maxy = tile[0][1]\n",
    "                tilex = \"00\"+str(tile[2][0])\n",
    "                tilex = tilex[-2:]\n",
    "                tiley = \"00\"+str(tile[2][1])\n",
    "                tiley = tiley[-2:]\n",
    "        \n",
    "                p1 = (minx, maxy)\n",
    "                p2 = (maxx, miny)\n",
    "\n",
    "                i1 = int((p1[0] - xOrigin) / pixelWidth)\n",
    "                j1 = int((yOrigin - p1[1]) / pixelHeight)\n",
    "                i2 = int((p2[0] - xOrigin) / pixelWidth)\n",
    "                j2 = int((yOrigin - p2[1]) / pixelHeight)\n",
    "\n",
    "                new_cols = i2-i1\n",
    "                new_rows = j2-j1\n",
    "\n",
    "                data = band.ReadAsArray(i1, j1, new_cols, new_rows)\n",
    "                #print(\"data\",data.mean(),data)\n",
    "                #if len(list(set(chain(*data))))<3:\n",
    "                if(np.all(data==data[0]) or (np.amin(data)==-9999 and np.amax(data)==0)):\n",
    "                    print(\"Ignore this tile. There is no usable image data for it. row:\" + tiley + \" col:\" + tilex)\n",
    "                    #print(\"data.mean\",data.mean())\n",
    "                else:\n",
    "                    #print(\"data.mean.too\",data.mean(),np.amin(data),np.amax(data) )\n",
    "                \n",
    "                    new_x = xOrigin + i1*pixelWidth\n",
    "                    new_y = yOrigin - j1*pixelHeight\n",
    "\n",
    "                    #print (new_x, new_y)\n",
    "\n",
    "                    new_transform = (new_x, transform[1], transform[2], new_y, transform[4], transform[5])\n",
    "                    output_file_name = area_str + panorpas + tiley + tilex\n",
    "                    output_file_base = output_file_name + \".tif\"\n",
    "                    output_file = os.path.join(data_directory,\"images\",\"split_tifs\",output_file_base)\n",
    "\n",
    "                    dst_ds = driver.Create(output_file,\n",
    "                                           new_cols,\n",
    "                                           new_rows,\n",
    "                                           1,\n",
    "                                           gdal.GDT_Float32)\n",
    "\n",
    "                    #writing output raster\n",
    "                    dst_ds.GetRasterBand(1).WriteArray( data )\n",
    "                    dst_ds.SetProjection( srs.ExportToWkt() )\n",
    "                    dst_ds.SetGeoTransform(new_transform)\n",
    "                    tif_metadata = {\n",
    "                        \"minX\": str(minx), \"maxX\": str(maxx),\n",
    "                        \"minY\": str(miny), \"maxY\": str(maxy)\n",
    "                    }\n",
    "                    dst_ds.SetMetadata(tif_metadata)\n",
    "                    dst_ds = None            \n",
    "\n",
    "                    #make a polygon for this tile\n",
    "                    tpoly = get_poly_with_lat_long_from_geotif(output_file)\n",
    "                    new_tp_row = {'id':tile_num, 'geometry':tpoly}\n",
    "                    tile_boundary_poly = tile_boundary_poly.append(new_tp_row, ignore_index=True)\n",
    "            \n",
    "                    # Save .jpg\n",
    "                    # close tif before re-opening it and saving as jpg\n",
    "                    # The jpg is to be consumed by the model.\n",
    "                    \n",
    "                    tif_to_jpg(output_file,output_file_name,jpgs_output_folder)\n",
    "                    # /storage/images/charcoal_hearth_hill/images_training/\n",
    "                    # Are there points inside the geotiff?\n",
    "                    annotate_tif_if_it_has_hearths(area_hearth_data,tpoly,output_file,output_file_name,geotif_crs)\n",
    "                \n",
    "                    #We need to refine to handle points on the border\n",
    "            \n",
    "                #Close output raster dataset\n",
    "                dst_ds = None\n",
    "\n",
    "                tile_num += 1\n",
    "\n",
    "            dataset = None\n",
    "    \n",
    "            # Determine the output path for the Shapefile\n",
    "            crs_4326 = 4326        \n",
    "            print(\"Tile boundary EPSG:\", str(geotif_crs))\n",
    "            tile_boundary_poly.crs = (\"EPSG:\" + str(crs_4326))\n",
    "            tile_boundary_poly.geometry = tile_boundary_poly.geometry.to_crs(crs=crs_4326)\n",
    "            tile_boundary_poly.to_crs(crs=crs_4326)\n",
    "            tile_boundary_poly = tile_boundary_poly.to_crs(epsg=crs_4326)\n",
    "        \n",
    "            # tile_boundary_poly = tile_boundary_poly.to_crs(epsg=geotif_crs)\n",
    "    \n",
    "            outfp = os.path.join(polys_output_folder, (area_str + panorpas + \"_tile_boundary.shp\"))\n",
    "            # Write the data into that Shapefile\n",
    "            tile_boundary_poly.to_file(outfp)\n",
    "        else:\n",
    "            print(proj_dataset.GetAttrValue('AUTHORITY',1),geotif_crs,\" don't match.\")\n",
    "    else:\n",
    "        print(file_name+\" does not exist.\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "split(\"/storage/images/part_1/39blast2demPAN_Slope_DEM.tif\",32128)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Instructions\n",
    "Download the Slope_DEM.tif files\n",
    "PAN = 32128\n",
    "PAS = 32129\n",
    "Run all cells above and then the one below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#split(\"/storage/images/charcoal_hearth_hill/downloads/Weston_Uploads/Lidar_and_Raster_Images_5-20-2020/39/39blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/charcoal_hearth_hill/downloads/Weston_Uploads/Lidar_and_Raster_Images_5-20-2020/41/41blast2demPAS_Slope_DEM.tif\",32129)\n",
    "#split(\"/storage/images/charcoal_hearth_hill/downloads/Weston_Uploads/Lidar_and_Raster_Images_5-20-2020/43/43blast2demPAS_Slope_DEM.tif\",32129)\n",
    "#split(\"/storage/images/charcoal_hearth_hill/downloads/Weston_Uploads/Lidar_and_Raster_Images_5-20-2020/46/46blast2demPAS_Slope_DEM.tif\",32129)\n",
    "#split(\"/storage/images/charcoal_hearth_hill/downloads/Weston_Uploads/Lidar_and_Raster_Images_5-20-2020/52/52partialDEM_Slope.tif\",26918)\n",
    "#split(\"/storage/images/charcoal_hearth_hill/downloads/Weston_Uploads/Lidar_and_Raster_Images_5-20-2020/55/55blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/charcoal_hearth_hill/downloads/Weston_Uploads/Lidar_and_Raster_Images_5-20-2020/58/58blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/12blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/13blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/14blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/24blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/25blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/26blast2demPAN_Slope_DEM.tif\",32128) - AttributeError: 'NoneType' object has no attribute 'GetProjection'\n",
    "#split(\"/storage/images/28blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/29blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/30blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/31blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#33 has PAN and PAS\n",
    "#split(\"/storage/images/part_1/33blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/part_1/33blast2demPAS_Slope_DEM.tif\",32129)\n",
    "\n",
    "#split(\"/storage/images/part_1/34blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/part_1/35blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/part_1/36blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/part_1/37blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/part_1/38blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/part_1/39blast2demPAN_Slope_DEM.tif\",32128)\n",
    "#split(\"/storage/images/part_1/40blast2demPAN_Slope_DEM.tif\",32128)\n",
    "\n",
    "\n",
    "\n",
    "#split(\"/storage/images/part_4/320blast2demPAN_DEM_Slope.tif\",32128)\n",
    "\n",
    "#52partialDEM_Slope.tif\",26918)\n",
    "#58 has PAN and PAS\n",
    "#60 has PAN and PAS\n",
    "#83 is missing - now there.\n",
    "#84 has PAN and PAS\n",
    "#95 has PAN and PAS\n",
    "#105 has PAN and PAS\n",
    "#108 has PAN and PAS\n",
    "#109 is missing\n",
    "\n",
    "# on /storage\n",
    "# 12:29:18 Current Time = 12:30:02 == 44 secs\n",
    "\n",
    "#on home\n",
    "# 12:31:29 12:31:51 == 22 secs\n",
    "        \n",
    "from datetime import datetime\n",
    "now = datetime.now()\n",
    "current_time = now.strftime(\"%H:%M:%S\")\n",
    "print(\"Current Time =\", current_time)\n",
    "\n",
    "split(\"/storage/images/part_1/12blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/13blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/14blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/24blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/25blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/26blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/28blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/29blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/30blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/31blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/33blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/33blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/34blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/35blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/36blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/37blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/38blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/39blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/40blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/41blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/42blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/43blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/44blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/45blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/46blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/47blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/48blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/49blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/50blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/51blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/52blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/53blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/54blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/55blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/56blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/57blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/58blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/58blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/59blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/60blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/60blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/61blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/62blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/63blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/64blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/65blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/66blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/67blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/68blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/69blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/70blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/71blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/72blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/73blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/74blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/75blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/76blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/77blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/78blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/79blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/80blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/81blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/82blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/83blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/84blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/84blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/85blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/86blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/87blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/88blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/89blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/90blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/91blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/92blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/93blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/94blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/95blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/95blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_1/96blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_1/97blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/100blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/101blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/102blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/103blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/104blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/105blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/105blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/106blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/107blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/108blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/108blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/109blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/110blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/111blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/112blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/113blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/114blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/115blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/116blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/117blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/118blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/119blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/120blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/120blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/121blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/122blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/123blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/124blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/126blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/127blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/128blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/129blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/130blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/131blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/132blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/132blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/133blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/134blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/135blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/136blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/137blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/137blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/138blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/139blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/140blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/141blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/142blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/143blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/144blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/145blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/146blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/147blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/148blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/149blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/150blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/151blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/151blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/152blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/153blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/154blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/155blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/156blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/157blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/158blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/158blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/159blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/160blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/161blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/162blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/163blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/164blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/165blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/166blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/167blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/168blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/168blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/169blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/170blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/171blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/172blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/173blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_2/98blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_2/99blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/174blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/175blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/176blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/178blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/179blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/180blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/181blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/182blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/183blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/184blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/185blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/186blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/187blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/188blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/189blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/190blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/191blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/192blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/193blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/193blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/194blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/195blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/196blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/197blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/198blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/199blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/200blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/201blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/201blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/202blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/203blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/204blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/205blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/206blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/207blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/208blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/209blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/210blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/211blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/212blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/212blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/213blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/214blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/215blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/216blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/217blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/217blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/218blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/219blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/220blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/221blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/222blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/223blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/224blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/225blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/226blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/227blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/228blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/229blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/230blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/231blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/232blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/233blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/233blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/234blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/235blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/236blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/237blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/238blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/239blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/242blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/243blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/244blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/245blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/246blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/247blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/248blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/249blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_3/250blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_3/251blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/252blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/253blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/254blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/255blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/256blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/257blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/257blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/258blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/258blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/259blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/260blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/261blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/262blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/263blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/264blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/265blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/266blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/267blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/268blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/269blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/270blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/271blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/272blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/273blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/274blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/275blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/276blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/277blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/278blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/278blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/279blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/280blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/281blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/282blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/283blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/284blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/284blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/285blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/286blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/287blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/288blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/290blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/291blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/292blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/293blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/294blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/295blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/296blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/297blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/298blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/299blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/300blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/301blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/302blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/303blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/304blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/305blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/306blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/307blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/309blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/310blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/311blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/312blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/313blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/314blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/315blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/316blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/317blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/317blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/318blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/320blast2demPAN_DEM_Slope.tif\",32128)\n",
    "split(\"/storage/images/part_4/321blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/322blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/323blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/324blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/325blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/326blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/327blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/328blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/329blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/329blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/330blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/331blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/332blast2demPAS_Slope_DEM.tif\",32129)\n",
    "split(\"/storage/images/part_4/333blast2demPAN_Slope_DEM.tif\",32128)\n",
    "split(\"/storage/images/part_4/335blast2demPAN_Slope_DEM.tif\",32128)\n",
    "\n",
    "now = datetime.now()\n",
    "current_time = now.strftime(\"%H:%M:%S\")\n",
    "print(\"Current Time =\", current_time)\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}