73 transect

In [ ]:

import ee
import geemap
from bqplot import pyplot as plt

import ee import geemap from bqplot import pyplot as plt

In [ ]:

Map = geemap.Map()
Map.add_basemap("TERRAIN")
Map

Map = geemap.Map() Map.add_basemap("TERRAIN") Map

In [ ]:

image = ee.Image("USGS/SRTMGL1_003")
vis_params = {
    "min": 0,
    "max": 4000,
    "palette": ["006633", "E5FFCC", "662A00", "D8D8D8", "F5F5F5"],
}
Map.addLayer(image, vis_params, "SRTM DEM", True, 0.5)

image = ee.Image("USGS/SRTMGL1_003") vis_params = { "min": 0, "max": 4000, "palette": ["006633", "E5FFCC", "662A00", "D8D8D8", "F5F5F5"], } Map.addLayer(image, vis_params, "SRTM DEM", True, 0.5)

In [ ]:

# Use the drawing tool to draw any line on the map.
line = Map.user_roi
if line is None:
    line = ee.Geometry.LineString(
        [[-120.223279, 36.314849], [-118.926969, 36.712192], [-117.202217, 36.756215]]
    )
    Map.addLayer(line, {}, "ROI")
Map.centerObject(line)

# Use the drawing tool to draw any line on the map. line = Map.user_roi if line is None: line = ee.Geometry.LineString( [[-120.223279, 36.314849], [-118.926969, 36.712192], [-117.202217, 36.756215]] ) Map.addLayer(line, {}, "ROI") Map.centerObject(line)

In [ ]:

line.getInfo()

line.getInfo()

In [ ]:

reducer = "mean"  # Any ee.Reducer, e.g., mean, median, min, max, stdDev
transect = geemap.extract_transect(
    image, line, n_segments=100, reducer=reducer, to_pandas=True
)

reducer = "mean" # Any ee.Reducer, e.g., mean, median, min, max, stdDev transect = geemap.extract_transect( image, line, n_segments=100, reducer=reducer, to_pandas=True )

In [ ]:

transect

transect

In [ ]:

fig = plt.figure()
plt.plot(transect["distance"], transect[reducer])
plt.xlabel("Distance")
plt.ylabel("Elevation")
plt.show()

fig = plt.figure() plt.plot(transect["distance"], transect[reducer]) plt.xlabel("Distance") plt.ylabel("Elevation") plt.show()