20250312
12542000
1.0
FALSE
002
DtDEM_2022_1ft_WebMerc
-13706923.620521
-13668166.089406
4959074.293589
4979513.308067
1
Projected
GCS_WGS_1984
Linear Unit: Meter (1.000000)
WGS_1984_Web_Mercator_Auxiliary_Sphere
<ProjectedCoordinateSystem xsi:type='typens:ProjectedCoordinateSystem' xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance' xmlns:xs='http://www.w3.org/2001/XMLSchema' xmlns:typens='http://www.esri.com/schemas/ArcGIS/2.8.0'><WKT>PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0],AUTHORITY["EPSG",3857]]</WKT><XOrigin>-20037700</XOrigin><YOrigin>-30241100</YOrigin><XYScale>148923141.92838538</XYScale><ZOrigin>-100000</ZOrigin><ZScale>10000</ZScale><MOrigin>-100000</MOrigin><MScale>10000</MScale><XYTolerance>0.001</XYTolerance><ZTolerance>0.001</ZTolerance><MTolerance>0.001</MTolerance><HighPrecision>true</HighPrecision><WKID>102100</WKID><LatestWKID>3857</LatestWKID></ProjectedCoordinateSystem>
150000000
5000
ISO19139
20250505
22092400
20250505
22092400
Raster-based data viewing and analysis.
<DIV STYLE="text-align:Left;"><DIV><DIV><P><SPAN>This is a raster with 1-ft cell resolution representing a detrended version of the ArcGIS Terrain dataset, "2022 Trinity River Focal Reach Digital Terrain" (YT TSP 2023). For further details of the source data can be found in: Pryor, C., T. Grey, and B. Hocker. 2023. 2022 Trinity River focal reach digital terrain model survey report. Report for the Trinity River Restoration Program (TRRP). Yurok Tribal Fisheries Program, Arcata, California. Available: https://www.trrp.net/library/document?id=2627.</SPAN></P><P><SPAN>Detrending used a mass point file derived from wetted cells (water depth > 0) in the 2022 SRH2D modeled flow output for 450cfs, along with the 2022 DTM's "Terrain_Extent" polygon data plus a regional buffer to create a "hydro-planar" Terrain surface that represents the water surface from the source data. The regional buffer was created by buffering the "Terrain_Extent" polygon by 1,000,000 feet. Elevation for the regional buffer was set to a uniform 1596.69 ft* as previously used for the elevation mid-point of the focal reach. Both the original DTM and this hydro-planar surface were converted to 1 ft resolution rasters. The hydro-planer raster was then subtracted from the DTM-derived raster (a DEM) to yield this Detrended Digital Elevation Model (DtDEM)</SPAN></P><P><SPAN>* = This elevation is the midpoint between the water surface elevation at the lower end of the DTM and the water surface elevation immediately below Lewiston Dam. While this uniform elevation theoretically could make water flow out of the river near the dam (river is higher than the value of 1596.69) and into the river near the North Fork, the slope created by this is minuscule. At the upper and lower ends of the river where this potential slope would be greatest, it could account for only about 0.025 ft elevation per 100 feet lateral to the river.</SPAN></P><P><SPAN>Citation: Yurok Tribe Fisheries Department Design and Technical Services Program (YT TSP), 2023. Data Package for 2022 Trinity River Focal Reach Digital Terrain. </SPAN></P></DIV></DIV></DIV>
Conversion by Trinity River Restoration Program from data provided by YT TSP.
Detrended Digital Elevation Model (DtDEM) Derived from the 2022 Trinity River Focal Reach Digital Terrain
TRRP (Trinity River Restoration Program) and YT TSP (Yurok Tribe Fisheries Department Design and Technical Services Program ) 2023. Detrended Digital Elevation Model (DtDEM) Derived from the 2022 Trinity River Focal Reach Digital Terrain. Data Package. TRRP, Weaverville, California.
2023-10-09T00:00:00
2023-10-10T00:00:00
Eric Peterson
Trinity River Restoration Program
Science Coordinator
Weaverville
CA
96093
US
530-623-1800
<DIV STYLE="text-align:Left;"><DIV><DIV><P><SPAN>Data provided as-is, with no warranty or guarantee. This data is for planning purposes only and should not be used for legal or cadastral purposes. Please contact the Trinity River Restoration Program (TRRP ) for information regarding the use of this data.</SPAN></P></DIV></DIV></DIV>
-123.131376
-122.783212
40.777054
40.637871
1
Eric Peterson
Trinity River Restoration Program
Science Coordinator
Eric Peterson
Jeanne McSloy
Trinity River Restoration Program
Natural Resource Specialist
Jeanne McSloy
Weaverville
CA
96093
US
530-623-1800
Cort Pryor
Yurok Tribal Fisheries Program
Geomatics and Survey Branch Manager
180 Klamath Blvd
Klamath
CA
95548
cpryor@yuroktribe.nsn.us
US
707-482-1350
Detrended DEM
20250505
dataset
Eric Peterson
Trinity River Restoration Program
Science Coordinator
Eric Peterson
Jeanne McSloy
Trinity River Restoration Program
Natural Resource Specialist
<DIV STYLE="text-align:Left;"><DIV><P><SPAN>Data provided as-is, with no warranty or guarantee.</SPAN></P></DIV></DIV>
File Geodatabase Raster Dataset
na
LZW
File Geodatabase Raster Dataset
3.3.1
Eric Peterson
Trinity River Restoration Program
Science Coordinator
Eric Peterson
EPSG
5.3(9.0.0)
Elevation relative to the water surface as represented by water surface elevation (WSE) modeled by SRH2D at 450cfs.
A brief description of the lineage is also provided in the metadata abstract. Intermediate data layers mentioned in the process steps are available from TRRP upon request.
import arcpy
import re
# Define the list of mass points files
mass_points_files = [
"P:F23D_F2P_s450cfs"
]
#
Define the paths for hard line and soft clip
hard_line = "P:DtDEM_beyond_OG"
soft_clip
= "P:T2022_Terrain_Extent_polygon"
#
Output feature dataset to store the terrains
output_fds = "P:AllFLows_FDS"
#
Check out the 3D Analyst extension
arcpy.CheckOutExtension("3D")
# Function to extract the number before "cfs"
def extract_number(file_path):
match = re.search(r'(\d+)cfs', file_path)
return match.group(1) if match else ""
# Iterate through each mass points file to create terrains
for mass_points in mass_points_files:
try:
# Extract the relevant part of the mass points file name
number = extract_number(mass_points)
terrain_name = f"Ter_{number}cfs"
# Create a new Terrain Dataset in the output feature dataset
terrain_path = arcpy.ddd.CreateTerrain(in_feature_dataset=output_fds, out_terrain_name=terrain_name, average_point_spacing=50, pyramid_type="ZTOLERANCE")[0]
# Add terrain pyramid levels
arcpy.ddd.AddTerrainPyramidLevel(in_terrain=terrain_path, pyramid_level_definition=["1 1000"])
# Add the hard line, soft clip, and mass points file to the Terrain
arcpy.ddd.AddFeatureClassToTerrain(in_terrain=terrain_path, in_features=[
[hard_line, "SHAPE", "Hard_Clip", 3, 0, 1, True, False, "DtDEM_beyond_OG_embed", "<None>", False],
[soft_clip, "<None>", "Soft_Clip", 4, 0, 1, True, False, "T2022_Terrain_Extent_polygon_embed", "<None>", False],
[mass_points, "Shape", "Mass_Points", 5, 0, 1, True, False, f"{number}cfs_embed", "<None>", False]
])
# Build the Terrain Dataset
arcpy.ddd.BuildTerrain(in_terrain=terrain_path)
print(f"Created terrain: {terrain_name}")
except Exception as e:
print(f"Error processing {mass_points}: {e}")
# Check in the 3D Analyst extension
arcpy.CheckInExtension("3D")
Creation of "DtDEM_beyond_OG". Created from 2016TERRAIN "DtDEM_beyond": the "Focal_Reach_Terrain_Extents_Polygon" was buffered by 1,000,000 feet. A new 3 dimensional polygon Feature Class was created and the buffer polygon (2D) was copied and pasted into it. The new feature class was then edited and all vertices were given an elevation of 1596.69 ft (the midpoint between the upper and lower elevations within the reach used previously).
Detrending requires extending the water surface laterally from the river. The large buffer allows the convenience of a single elevation for the lateral perimeter without introducing substantive error. See DtDEM_1ft from 2016 TERRAIN/2016_FocalReach_DTM.gdb for more details.
Jeanne McSloy
Trinity River Restoration Program
Natural Resource Specialist
Jeanne McSloy
Weaverville
CA
96093
US
530-623-1800
Jeanne McSloy
Creation of DtDEM_HydroSurface (Terrain format). A Terrain format layer was created using SRH2D model output for 450cfs wetted (water depth>0) cell centroids with water surface elevation (WSE) assigned to Z (mass points), DtDEM_beyond_OG (hard lines), and T2022_Terrain_Extents_Polygon" (soft clip).
Provides the laterally extended water surface elevation, clipped to match the 2022 DTM.
Creation of Detrending_HydroSurface_2022_1ft (raster). The Terrain to Raster tool was used with 1 foot resolution and processing extent set to match the DEM__2022_1ft raster.
Rasterization of the extended water surface, to use for final calculation of the DtDEM.
Creation of DtDEM_2022_1ft (raster). The Spatial Analyst tool Raster Calculator was used to subtract the Detrending_Hydrosurfase_2022_1ft from DEM_2022_1ft. Positive values are relative elevations above the water surface; negative values are below the water surface (bathymetric).
Provides the final raster DtDEM.
Pixel
67057
127157
0.304801
0.304801
32
1
Upper Left
FALSE
LZ77
1
pixel codes
FGDBR
TRUE
row and column
0.304801
0.304801