Topography.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "from geopy.distance import geodesic\n",
    "import rasterio\n",
    "from rasterio.warp import transform\n",
    "import requests"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# ERT N. 1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Coordinates"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Function to calculate distances between coordinates\n",
    "def calculate_distance(coord1, coord2):\n",
    "    return geodesic(coord1, coord2).meters\n",
    "\n",
    "# Function to get elevation from .tif file\n",
    "def get_elevation(lat, lon, src):\n",
    "    row, col = src.index(lon, lat)\n",
    "    return src.read(1)[row, col]\n",
    "\n",
    "# Coordinates\n",
    "coord1 = (6.7997777778, -72.8512777778)  # (Latitude, Longitude)\n",
    "coord2 = (6.8012777778, -72.8505)\n",
    "\n",
    "# Create dataframe\n",
    "num_points = 64\n",
    "dist_interval = 3  # in surface meters\n",
    "\n",
    "distances = [i * dist_interval for i in range(num_points)]\n",
    "points = []\n",
    "\n",
    "for distance in distances:\n",
    "    fraction = distance / calculate_distance(coord1, coord2)\n",
    "    lat = coord1[0] + fraction * (coord2[0] - coord1[0])\n",
    "    lon = coord1[1] + fraction * (coord2[1] - coord1[1])\n",
    "    elevation = 0  # Replace this with actual elevation data extraction from .tif file\n",
    "\n",
    "    points.append({\n",
    "        'Latitude': lat,\n",
    "        'Longitude': lon,\n",
    "        'Distance': distance,\n",
    "    })\n",
    "\n",
    "# Create dataframe\n",
    "df1 = pd.DataFrame(points)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Latitude</th>\n",
       "      <th>Longitude</th>\n",
       "      <th>Distance</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>6.799778</td>\n",
       "      <td>-72.851278</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>6.799802</td>\n",
       "      <td>-72.851265</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>6.799826</td>\n",
       "      <td>-72.851253</td>\n",
       "      <td>6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>6.799850</td>\n",
       "      <td>-72.851240</td>\n",
       "      <td>9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>6.799874</td>\n",
       "      <td>-72.851228</td>\n",
       "      <td>12</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59</th>\n",
       "      <td>6.801199</td>\n",
       "      <td>-72.850541</td>\n",
       "      <td>177</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>60</th>\n",
       "      <td>6.801223</td>\n",
       "      <td>-72.850528</td>\n",
       "      <td>180</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61</th>\n",
       "      <td>6.801247</td>\n",
       "      <td>-72.850516</td>\n",
       "      <td>183</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>62</th>\n",
       "      <td>6.801271</td>\n",
       "      <td>-72.850504</td>\n",
       "      <td>186</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63</th>\n",
       "      <td>6.801295</td>\n",
       "      <td>-72.850491</td>\n",
       "      <td>189</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>64 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "    Latitude  Longitude  Distance\n",
       "0   6.799778 -72.851278         0\n",
       "1   6.799802 -72.851265         3\n",
       "2   6.799826 -72.851253         6\n",
       "3   6.799850 -72.851240         9\n",
       "4   6.799874 -72.851228        12\n",
       "..       ...        ...       ...\n",
       "59  6.801199 -72.850541       177\n",
       "60  6.801223 -72.850528       180\n",
       "61  6.801247 -72.850516       183\n",
       "62  6.801271 -72.850504       186\n",
       "63  6.801295 -72.850491       189\n",
       "\n",
       "[64 rows x 3 columns]"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Elevation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Define a function to get elevation from Google Maps Elevation API\n",
    "def get_elevation_google_maps(api_key, lat, lon):\n",
    "    base_url = \"https://maps.googleapis.com/maps/api/elevation/json\"\n",
    "    params = {\n",
    "        'locations': f'{lat},{lon}',\n",
    "        'key': api_key\n",
    "    }\n",
    "    response = requests.get(base_url, params=params)\n",
    "    data = response.json()\n",
    "\n",
    "    if 'results' in data and data['results']:\n",
    "        elevation = data['results'][0]['elevation']\n",
    "        return elevation\n",
    "    else:\n",
    "        print(f\"Failed to get elevation for coordinates ({lat}, {lon}): {data.get('error_message')}\")\n",
    "        return None\n",
    "\n",
    "# Apply the function to each row in the dataframe\n",
    "api_key = 'INSERT_API_KEY'\n",
    "df1['Elevation'] = df1.apply(lambda row: get_elevation_google_maps(api_key, row['Latitude'], row['Longitude']), axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Latitude</th>\n",
       "      <th>Longitude</th>\n",
       "      <th>Distance</th>\n",
       "      <th>Elevation</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>6.799778</td>\n",
       "      <td>-72.851278</td>\n",
       "      <td>0</td>\n",
       "      <td>1509.781738</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>6.799802</td>\n",
       "      <td>-72.851265</td>\n",
       "      <td>3</td>\n",
       "      <td>1510.302124</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>6.799826</td>\n",
       "      <td>-72.851253</td>\n",
       "      <td>6</td>\n",
       "      <td>1510.832886</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>6.799850</td>\n",
       "      <td>-72.851240</td>\n",
       "      <td>9</td>\n",
       "      <td>1511.373779</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>6.799874</td>\n",
       "      <td>-72.851228</td>\n",
       "      <td>12</td>\n",
       "      <td>1511.723755</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59</th>\n",
       "      <td>6.801199</td>\n",
       "      <td>-72.850541</td>\n",
       "      <td>177</td>\n",
       "      <td>1535.642090</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>60</th>\n",
       "      <td>6.801223</td>\n",
       "      <td>-72.850528</td>\n",
       "      <td>180</td>\n",
       "      <td>1536.085083</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61</th>\n",
       "      <td>6.801247</td>\n",
       "      <td>-72.850516</td>\n",
       "      <td>183</td>\n",
       "      <td>1536.561768</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>62</th>\n",
       "      <td>6.801271</td>\n",
       "      <td>-72.850504</td>\n",
       "      <td>186</td>\n",
       "      <td>1537.039917</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63</th>\n",
       "      <td>6.801295</td>\n",
       "      <td>-72.850491</td>\n",
       "      <td>189</td>\n",
       "      <td>1537.517944</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>64 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "    Latitude  Longitude  Distance    Elevation\n",
       "0   6.799778 -72.851278         0  1509.781738\n",
       "1   6.799802 -72.851265         3  1510.302124\n",
       "2   6.799826 -72.851253         6  1510.832886\n",
       "3   6.799850 -72.851240         9  1511.373779\n",
       "4   6.799874 -72.851228        12  1511.723755\n",
       "..       ...        ...       ...          ...\n",
       "59  6.801199 -72.850541       177  1535.642090\n",
       "60  6.801223 -72.850528       180  1536.085083\n",
       "61  6.801247 -72.850516       183  1536.561768\n",
       "62  6.801271 -72.850504       186  1537.039917\n",
       "63  6.801295 -72.850491       189  1537.517944\n",
       "\n",
       "[64 rows x 4 columns]"
      ]
     },
     "execution_count": 44,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [],
   "source": [
    "df1.to_csv('Surveys/elev_data_Nov24.csv', index=False)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# ERT N. 2"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Coordinates"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Function to calculate distances between coordinates\n",
    "def calculate_distance(coord1, coord2):\n",
    "    return geodesic(coord1, coord2).meters\n",
    "\n",
    "# Function to get elevation from .tif file\n",
    "def get_elevation(lat, lon, src):\n",
    "    row, col = src.index(lon, lat)\n",
    "    return src.read(1)[row, col]\n",
    "\n",
    "# Coordinates\n",
    "coord1 = (6.7918333333, -72.8475277778)  # (Latitude, Longitude)\n",
    "coord2 = (6.7932222222, -72.8474166667)\n",
    "\n",
    "# Create dataframe\n",
    "num_points = 64\n",
    "dist_interval = 2.5  # in surface meters\n",
    "\n",
    "distances = [i * dist_interval for i in range(num_points)]\n",
    "points = []\n",
    "\n",
    "for distance in distances:\n",
    "    fraction = distance / calculate_distance(coord1, coord2)\n",
    "    lat = coord1[0] + fraction * (coord2[0] - coord1[0])\n",
    "    lon = coord1[1] + fraction * (coord2[1] - coord1[1])\n",
    "    elevation = 0  # Replace this with actual elevation data extraction from .tif file\n",
    "\n",
    "    points.append({\n",
    "        'Latitude': lat,\n",
    "        'Longitude': lon,\n",
    "        'Distance': distance,\n",
    "    })\n",
    "\n",
    "# Create dataframe\n",
    "df2 = pd.DataFrame(points)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Latitude</th>\n",
       "      <th>Longitude</th>\n",
       "      <th>Distance</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>6.791833</td>\n",
       "      <td>-72.847528</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>6.791856</td>\n",
       "      <td>-72.847526</td>\n",
       "      <td>2.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>6.791878</td>\n",
       "      <td>-72.847524</td>\n",
       "      <td>5.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>6.791901</td>\n",
       "      <td>-72.847522</td>\n",
       "      <td>7.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>6.791923</td>\n",
       "      <td>-72.847521</td>\n",
       "      <td>10.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59</th>\n",
       "      <td>6.793163</td>\n",
       "      <td>-72.847421</td>\n",
       "      <td>147.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>60</th>\n",
       "      <td>6.793185</td>\n",
       "      <td>-72.847420</td>\n",
       "      <td>150.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61</th>\n",
       "      <td>6.793208</td>\n",
       "      <td>-72.847418</td>\n",
       "      <td>152.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>62</th>\n",
       "      <td>6.793230</td>\n",
       "      <td>-72.847416</td>\n",
       "      <td>155.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63</th>\n",
       "      <td>6.793253</td>\n",
       "      <td>-72.847414</td>\n",
       "      <td>157.5</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>64 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "    Latitude  Longitude  Distance\n",
       "0   6.791833 -72.847528       0.0\n",
       "1   6.791856 -72.847526       2.5\n",
       "2   6.791878 -72.847524       5.0\n",
       "3   6.791901 -72.847522       7.5\n",
       "4   6.791923 -72.847521      10.0\n",
       "..       ...        ...       ...\n",
       "59  6.793163 -72.847421     147.5\n",
       "60  6.793185 -72.847420     150.0\n",
       "61  6.793208 -72.847418     152.5\n",
       "62  6.793230 -72.847416     155.0\n",
       "63  6.793253 -72.847414     157.5\n",
       "\n",
       "[64 rows x 3 columns]"
      ]
     },
     "execution_count": 47,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df2"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Elevation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Define a function to get elevation from Google Maps Elevation API\n",
    "def get_elevation_google_maps(api_key, lat, lon):\n",
    "    base_url = \"https://maps.googleapis.com/maps/api/elevation/json\"\n",
    "    params = {\n",
    "        'locations': f'{lat},{lon}',\n",
    "        'key': api_key\n",
    "    }\n",
    "    response = requests.get(base_url, params=params)\n",
    "    data = response.json()\n",
    "\n",
    "    if 'results' in data and data['results']:\n",
    "        elevation = data['results'][0]['elevation']\n",
    "        return elevation\n",
    "    else:\n",
    "        print(f\"Failed to get elevation for coordinates ({lat}, {lon}): {data.get('error_message')}\")\n",
    "        return None\n",
    "\n",
    "# Apply the function to each row in the dataframe\n",
    "api_key = 'INSERT_API_KEY'\n",
    "df2['Elevation'] = df2.apply(lambda row: get_elevation_google_maps(api_key, row['Latitude'], row['Longitude']), axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Latitude</th>\n",
       "      <th>Longitude</th>\n",
       "      <th>Distance</th>\n",
       "      <th>Elevation</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>6.791833</td>\n",
       "      <td>-72.847528</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1461.690063</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>6.791856</td>\n",
       "      <td>-72.847526</td>\n",
       "      <td>2.5</td>\n",
       "      <td>1462.279419</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>6.791878</td>\n",
       "      <td>-72.847524</td>\n",
       "      <td>5.0</td>\n",
       "      <td>1462.868896</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>6.791901</td>\n",
       "      <td>-72.847522</td>\n",
       "      <td>7.5</td>\n",
       "      <td>1463.458252</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>6.791923</td>\n",
       "      <td>-72.847521</td>\n",
       "      <td>10.0</td>\n",
       "      <td>1464.047729</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>59</th>\n",
       "      <td>6.793163</td>\n",
       "      <td>-72.847421</td>\n",
       "      <td>147.5</td>\n",
       "      <td>1488.817627</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>60</th>\n",
       "      <td>6.793185</td>\n",
       "      <td>-72.847420</td>\n",
       "      <td>150.0</td>\n",
       "      <td>1489.113892</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>61</th>\n",
       "      <td>6.793208</td>\n",
       "      <td>-72.847418</td>\n",
       "      <td>152.5</td>\n",
       "      <td>1489.410767</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>62</th>\n",
       "      <td>6.793230</td>\n",
       "      <td>-72.847416</td>\n",
       "      <td>155.0</td>\n",
       "      <td>1489.708130</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63</th>\n",
       "      <td>6.793253</td>\n",
       "      <td>-72.847414</td>\n",
       "      <td>157.5</td>\n",
       "      <td>1490.006104</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>64 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "    Latitude  Longitude  Distance    Elevation\n",
       "0   6.791833 -72.847528       0.0  1461.690063\n",
       "1   6.791856 -72.847526       2.5  1462.279419\n",
       "2   6.791878 -72.847524       5.0  1462.868896\n",
       "3   6.791901 -72.847522       7.5  1463.458252\n",
       "4   6.791923 -72.847521      10.0  1464.047729\n",
       "..       ...        ...       ...          ...\n",
       "59  6.793163 -72.847421     147.5  1488.817627\n",
       "60  6.793185 -72.847420     150.0  1489.113892\n",
       "61  6.793208 -72.847418     152.5  1489.410767\n",
       "62  6.793230 -72.847416     155.0  1489.708130\n",
       "63  6.793253 -72.847414     157.5  1490.006104\n",
       "\n",
       "[64 rows x 4 columns]"
      ]
     },
     "execution_count": 49,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {},
   "outputs": [],
   "source": [
    "df2.to_csv('Surveys/elev_data_Nov25_1.csv', index=False)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# ERT N. 3"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Coordinates"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Function to calculate distances between coordinates\n",
    "def calculate_distance(coord1, coord2):\n",
    "    return geodesic(coord1, coord2).meters\n",
    "\n",
    "# Function to get elevation from .tif file\n",
    "def get_elevation(lat, lon, src):\n",
    "    row, col = src.index(lon, lat)\n",
    "    return src.read(1)[row, col]\n",
    "\n",
    "# Coordinates\n",
    "coord1 = (6.8061111111, -72.8525000000)  # (Latitude, Longitude)\n",
    "coord2 = (6.8075, -72.8522222222)\n",
    "\n",
    "# Create dataframe\n",
    "num_points = 32\n",
    "dist_interval = 5  # in surface meters\n",
    "\n",
    "distances = [i * dist_interval for i in range(num_points)]\n",
    "points = []\n",
    "\n",
    "for distance in distances:\n",
    "    fraction = distance / calculate_distance(coord1, coord2)\n",
    "    lat = coord1[0] + fraction * (coord2[0] - coord1[0])\n",
    "    lon = coord1[1] + fraction * (coord2[1] - coord1[1])\n",
    "    elevation = 0  # Replace this with actual elevation data extraction from .tif file\n",
    "\n",
    "    points.append({\n",
    "        'Latitude': lat,\n",
    "        'Longitude': lon,\n",
    "        'Distance': distance,\n",
    "    })\n",
    "\n",
    "# Create dataframe\n",
    "df3 = pd.DataFrame(points)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Latitude</th>\n",
       "      <th>Longitude</th>\n",
       "      <th>Distance</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>6.806111</td>\n",
       "      <td>-72.852500</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>6.806155</td>\n",
       "      <td>-72.852491</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>6.806200</td>\n",
       "      <td>-72.852482</td>\n",
       "      <td>10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>6.806244</td>\n",
       "      <td>-72.852473</td>\n",
       "      <td>15</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>6.806288</td>\n",
       "      <td>-72.852465</td>\n",
       "      <td>20</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>6.806333</td>\n",
       "      <td>-72.852456</td>\n",
       "      <td>25</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>6.806377</td>\n",
       "      <td>-72.852447</td>\n",
       "      <td>30</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>6.806421</td>\n",
       "      <td>-72.852438</td>\n",
       "      <td>35</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>6.806466</td>\n",
       "      <td>-72.852429</td>\n",
       "      <td>40</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>6.806510</td>\n",
       "      <td>-72.852420</td>\n",
       "      <td>45</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>6.806554</td>\n",
       "      <td>-72.852411</td>\n",
       "      <td>50</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>6.806599</td>\n",
       "      <td>-72.852402</td>\n",
       "      <td>55</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>6.806643</td>\n",
       "      <td>-72.852394</td>\n",
       "      <td>60</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>6.806687</td>\n",
       "      <td>-72.852385</td>\n",
       "      <td>65</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>6.806732</td>\n",
       "      <td>-72.852376</td>\n",
       "      <td>70</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>6.806776</td>\n",
       "      <td>-72.852367</td>\n",
       "      <td>75</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>6.806820</td>\n",
       "      <td>-72.852358</td>\n",
       "      <td>80</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>6.806865</td>\n",
       "      <td>-72.852349</td>\n",
       "      <td>85</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>6.806909</td>\n",
       "      <td>-72.852340</td>\n",
       "      <td>90</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>6.806953</td>\n",
       "      <td>-72.852332</td>\n",
       "      <td>95</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>6.806998</td>\n",
       "      <td>-72.852323</td>\n",
       "      <td>100</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>6.807042</td>\n",
       "      <td>-72.852314</td>\n",
       "      <td>105</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>6.807086</td>\n",
       "      <td>-72.852305</td>\n",
       "      <td>110</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>6.807131</td>\n",
       "      <td>-72.852296</td>\n",
       "      <td>115</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>6.807175</td>\n",
       "      <td>-72.852287</td>\n",
       "      <td>120</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>6.807219</td>\n",
       "      <td>-72.852278</td>\n",
       "      <td>125</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>6.807264</td>\n",
       "      <td>-72.852269</td>\n",
       "      <td>130</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>6.807308</td>\n",
       "      <td>-72.852261</td>\n",
       "      <td>135</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>6.807352</td>\n",
       "      <td>-72.852252</td>\n",
       "      <td>140</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>6.807397</td>\n",
       "      <td>-72.852243</td>\n",
       "      <td>145</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>6.807441</td>\n",
       "      <td>-72.852234</td>\n",
       "      <td>150</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>6.807485</td>\n",
       "      <td>-72.852225</td>\n",
       "      <td>155</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Latitude  Longitude  Distance\n",
       "0   6.806111 -72.852500         0\n",
       "1   6.806155 -72.852491         5\n",
       "2   6.806200 -72.852482        10\n",
       "3   6.806244 -72.852473        15\n",
       "4   6.806288 -72.852465        20\n",
       "5   6.806333 -72.852456        25\n",
       "6   6.806377 -72.852447        30\n",
       "7   6.806421 -72.852438        35\n",
       "8   6.806466 -72.852429        40\n",
       "9   6.806510 -72.852420        45\n",
       "10  6.806554 -72.852411        50\n",
       "11  6.806599 -72.852402        55\n",
       "12  6.806643 -72.852394        60\n",
       "13  6.806687 -72.852385        65\n",
       "14  6.806732 -72.852376        70\n",
       "15  6.806776 -72.852367        75\n",
       "16  6.806820 -72.852358        80\n",
       "17  6.806865 -72.852349        85\n",
       "18  6.806909 -72.852340        90\n",
       "19  6.806953 -72.852332        95\n",
       "20  6.806998 -72.852323       100\n",
       "21  6.807042 -72.852314       105\n",
       "22  6.807086 -72.852305       110\n",
       "23  6.807131 -72.852296       115\n",
       "24  6.807175 -72.852287       120\n",
       "25  6.807219 -72.852278       125\n",
       "26  6.807264 -72.852269       130\n",
       "27  6.807308 -72.852261       135\n",
       "28  6.807352 -72.852252       140\n",
       "29  6.807397 -72.852243       145\n",
       "30  6.807441 -72.852234       150\n",
       "31  6.807485 -72.852225       155"
      ]
     },
     "execution_count": 52,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df3"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Elevation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Define a function to get elevation from Google Maps Elevation API\n",
    "def get_elevation_google_maps(api_key, lat, lon):\n",
    "    base_url = \"https://maps.googleapis.com/maps/api/elevation/json\"\n",
    "    params = {\n",
    "        'locations': f'{lat},{lon}',\n",
    "        'key': api_key\n",
    "    }\n",
    "    response = requests.get(base_url, params=params)\n",
    "    data = response.json()\n",
    "\n",
    "    if 'results' in data and data['results']:\n",
    "        elevation = data['results'][0]['elevation']\n",
    "        return elevation\n",
    "    else:\n",
    "        print(f\"Failed to get elevation for coordinates ({lat}, {lon}): {data.get('error_message')}\")\n",
    "        return None\n",
    "\n",
    "# Apply the function to each row in the dataframe\n",
    "api_key = 'INSERT_API_KEY'\n",
    "df3['Elevation'] = df3.apply(lambda row: get_elevation_google_maps(api_key, row['Latitude'], row['Longitude']), axis=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Latitude</th>\n",
       "      <th>Longitude</th>\n",
       "      <th>Distance</th>\n",
       "      <th>Elevation</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>6.806111</td>\n",
       "      <td>-72.852500</td>\n",
       "      <td>0</td>\n",
       "      <td>1562.014038</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>6.806155</td>\n",
       "      <td>-72.852491</td>\n",
       "      <td>5</td>\n",
       "      <td>1562.830078</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>6.806200</td>\n",
       "      <td>-72.852482</td>\n",
       "      <td>10</td>\n",
       "      <td>1563.636230</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>6.806244</td>\n",
       "      <td>-72.852473</td>\n",
       "      <td>15</td>\n",
       "      <td>1564.432373</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>6.806288</td>\n",
       "      <td>-72.852465</td>\n",
       "      <td>20</td>\n",
       "      <td>1565.218506</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>6.806333</td>\n",
       "      <td>-72.852456</td>\n",
       "      <td>25</td>\n",
       "      <td>1565.994507</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>6.806377</td>\n",
       "      <td>-72.852447</td>\n",
       "      <td>30</td>\n",
       "      <td>1566.760620</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>6.806421</td>\n",
       "      <td>-72.852438</td>\n",
       "      <td>35</td>\n",
       "      <td>1567.516724</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>6.806466</td>\n",
       "      <td>-72.852429</td>\n",
       "      <td>40</td>\n",
       "      <td>1568.262817</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>6.806510</td>\n",
       "      <td>-72.852420</td>\n",
       "      <td>45</td>\n",
       "      <td>1568.998901</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>6.806554</td>\n",
       "      <td>-72.852411</td>\n",
       "      <td>50</td>\n",
       "      <td>1569.724976</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>6.806599</td>\n",
       "      <td>-72.852402</td>\n",
       "      <td>55</td>\n",
       "      <td>1570.441040</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>6.806643</td>\n",
       "      <td>-72.852394</td>\n",
       "      <td>60</td>\n",
       "      <td>1571.066772</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>6.806687</td>\n",
       "      <td>-72.852385</td>\n",
       "      <td>65</td>\n",
       "      <td>1571.671997</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>6.806732</td>\n",
       "      <td>-72.852376</td>\n",
       "      <td>70</td>\n",
       "      <td>1572.101929</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>6.806776</td>\n",
       "      <td>-72.852367</td>\n",
       "      <td>75</td>\n",
       "      <td>1572.213623</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>6.806820</td>\n",
       "      <td>-72.852358</td>\n",
       "      <td>80</td>\n",
       "      <td>1572.335327</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>6.806865</td>\n",
       "      <td>-72.852349</td>\n",
       "      <td>85</td>\n",
       "      <td>1572.466919</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>6.806909</td>\n",
       "      <td>-72.852340</td>\n",
       "      <td>90</td>\n",
       "      <td>1572.608643</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>6.806953</td>\n",
       "      <td>-72.852332</td>\n",
       "      <td>95</td>\n",
       "      <td>1572.760376</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>6.806998</td>\n",
       "      <td>-72.852323</td>\n",
       "      <td>100</td>\n",
       "      <td>1572.921997</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>6.807042</td>\n",
       "      <td>-72.852314</td>\n",
       "      <td>105</td>\n",
       "      <td>1573.093750</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>6.807086</td>\n",
       "      <td>-72.852305</td>\n",
       "      <td>110</td>\n",
       "      <td>1573.275513</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>6.807131</td>\n",
       "      <td>-72.852296</td>\n",
       "      <td>115</td>\n",
       "      <td>1573.467163</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>6.807175</td>\n",
       "      <td>-72.852287</td>\n",
       "      <td>120</td>\n",
       "      <td>1573.668945</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>6.807219</td>\n",
       "      <td>-72.852278</td>\n",
       "      <td>125</td>\n",
       "      <td>1573.880615</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>6.807264</td>\n",
       "      <td>-72.852269</td>\n",
       "      <td>130</td>\n",
       "      <td>1574.102295</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>6.807308</td>\n",
       "      <td>-72.852261</td>\n",
       "      <td>135</td>\n",
       "      <td>1574.334106</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>6.807352</td>\n",
       "      <td>-72.852252</td>\n",
       "      <td>140</td>\n",
       "      <td>1574.575806</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>6.807397</td>\n",
       "      <td>-72.852243</td>\n",
       "      <td>145</td>\n",
       "      <td>1574.827637</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>6.807441</td>\n",
       "      <td>-72.852234</td>\n",
       "      <td>150</td>\n",
       "      <td>1574.979858</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>6.807485</td>\n",
       "      <td>-72.852225</td>\n",
       "      <td>155</td>\n",
       "      <td>1575.122437</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    Latitude  Longitude  Distance    Elevation\n",
       "0   6.806111 -72.852500         0  1562.014038\n",
       "1   6.806155 -72.852491         5  1562.830078\n",
       "2   6.806200 -72.852482        10  1563.636230\n",
       "3   6.806244 -72.852473        15  1564.432373\n",
       "4   6.806288 -72.852465        20  1565.218506\n",
       "5   6.806333 -72.852456        25  1565.994507\n",
       "6   6.806377 -72.852447        30  1566.760620\n",
       "7   6.806421 -72.852438        35  1567.516724\n",
       "8   6.806466 -72.852429        40  1568.262817\n",
       "9   6.806510 -72.852420        45  1568.998901\n",
       "10  6.806554 -72.852411        50  1569.724976\n",
       "11  6.806599 -72.852402        55  1570.441040\n",
       "12  6.806643 -72.852394        60  1571.066772\n",
       "13  6.806687 -72.852385        65  1571.671997\n",
       "14  6.806732 -72.852376        70  1572.101929\n",
       "15  6.806776 -72.852367        75  1572.213623\n",
       "16  6.806820 -72.852358        80  1572.335327\n",
       "17  6.806865 -72.852349        85  1572.466919\n",
       "18  6.806909 -72.852340        90  1572.608643\n",
       "19  6.806953 -72.852332        95  1572.760376\n",
       "20  6.806998 -72.852323       100  1572.921997\n",
       "21  6.807042 -72.852314       105  1573.093750\n",
       "22  6.807086 -72.852305       110  1573.275513\n",
       "23  6.807131 -72.852296       115  1573.467163\n",
       "24  6.807175 -72.852287       120  1573.668945\n",
       "25  6.807219 -72.852278       125  1573.880615\n",
       "26  6.807264 -72.852269       130  1574.102295\n",
       "27  6.807308 -72.852261       135  1574.334106\n",
       "28  6.807352 -72.852252       140  1574.575806\n",
       "29  6.807397 -72.852243       145  1574.827637\n",
       "30  6.807441 -72.852234       150  1574.979858\n",
       "31  6.807485 -72.852225       155  1575.122437"
      ]
     },
     "execution_count": 54,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df3"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {},
   "outputs": [],
   "source": [
    "df3.to_csv('Surveys/elev_data_Nov25_2.csv', index=False)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "eda",
   "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.10.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}