Skip to content

Commit

Permalink
3DEP and NASADEM examples (#11)
Browse files Browse the repository at this point in the history
* update notebook, readme, lockfiles

* 3dep and nasadem notebooks

* add contributing doc

* nasadem files

* fix badge, use admonitions

* fix background header, pixi run cmd
  • Loading branch information
scottyhq authored Feb 13, 2025
1 parent 65e68a0 commit acf0787
Show file tree
Hide file tree
Showing 16 changed files with 3,030 additions and 239 deletions.
1,073 changes: 1,073 additions & 0 deletions 3dep/10m-seamless.ipynb

Large diffs are not rendered by default.

46 changes: 46 additions & 0 deletions 3dep/USGS_Seamless_DEM_13_7912.vrt
Original file line number Diff line number Diff line change
@@ -0,0 +1,46 @@
<VRTDataset rasterXSize="3888006" rasterYSize="939612" subClass="VRTWarpedDataset">
<SRS dataAxisToSRSAxisMapping="2,1,3">GEOGCRS["ITRF2014",DYNAMIC[FRAMEEPOCH[2010]],DATUM["International Terrestrial Reference Frame 2014",ELLIPSOID["GRS 1980",6378137,298.257222101,LENGTHUNIT["metre",1]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],CS[ellipsoidal,3],AXIS["geodetic latitude (Lat)",north,ORDER[1],ANGLEUNIT["degree",0.0174532925199433]],AXIS["geodetic longitude (Lon)",east,ORDER[2],ANGLEUNIT["degree",0.0174532925199433]],AXIS["ellipsoidal height (h)",up,ORDER[3],LENGTHUNIT["metre",1]],USAGE[SCOPE["Geodesy."],AREA["World."],BBOX[-90,-180,90,180]],ID["EPSG",7912]]</SRS>
<GeoTransform> -1.8000055560000001e+02, 9.2592592659999998e-05, 0.0000000000000000e+00, 7.2000555599999998e+01, 0.0000000000000000e+00, -9.2592592659999998e-05</GeoTransform>
<VRTRasterBand dataType="Float32" band="1" subClass="VRTWarpedRasterBand">
<NoDataValue>-999999</NoDataValue>
<ColorInterp>Gray</ColorInterp>
</VRTRasterBand>
<BlockXSize>256</BlockXSize>
<BlockYSize>256</BlockYSize>
<GDALWarpOptions>
<WarpMemoryLimit>5.24288e+08</WarpMemoryLimit>
<ResampleAlg>Bilinear</ResampleAlg>
<WorkingDataType>Float32</WorkingDataType>
<Option name="INIT_DEST">NO_DATA</Option>
<Option name="APPLY_VERTICAL_SHIFT">YES</Option>
<Option name="MULT_FACTOR_VERTICAL_SHIFT">1</Option>
<Option name="MULT_FACTOR_VERTICAL_SHIFT_PIPELINE">1</Option>
<Option name="ERROR_OUT_IF_EMPTY_SOURCE_WINDOW">FALSE</Option>
<SourceDataset relativeToVRT="0">/vsis3/prd-tnm/StagedProducts/Elevation/13/TIFF/USGS_Seamless_DEM_13.vrt</SourceDataset>
<Transformer>
<GenImgProjTransformer>
<SrcGeoTransform>-180.00055555609359,9.2592592660220417e-05,0,72.000555556294898,0,-9.25925926599094e-05</SrcGeoTransform>
<SrcInvGeoTransform>1944005.9985859466,10799.999992111891,0,777605.99944264861,0,-10799.999992148167</SrcInvGeoTransform>
<DstGeoTransform>-180.00055560000001,9.2592592659999998e-05,0,72.000555599999998,0,-9.2592592659999998e-05</DstGeoTransform>
<DstInvGeoTransform>1944005.9990647638,10799.9999921376,0,777605.99991390284,0,-10799.9999921376</DstInvGeoTransform>
<ReprojectTransformer>
<ReprojectionTransformer>
<SourceSRS>COMPD_CS["WGS 84 (G1150) + EGM96 height",GEOGCS["WGS 84 (G1150)",DATUM["World_Geodetic_System_1984_G1150",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","1154"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","9055"]],VERT_CS["EGM96 height",VERT_DATUM["EGM96 geoid",2005,AUTHORITY["EPSG","5171"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Gravity-related height",UP],AUTHORITY["EPSG","5773"]]]</SourceSRS>
<TargetSRS>GEOGCRS["ITRF2014",DYNAMIC[FRAMEEPOCH[2010]],DATUM["International Terrestrial Reference Frame 2014",ELLIPSOID["GRS 1980",6378137,298.257222101,LENGTHUNIT["metre",1]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],CS[ellipsoidal,3],AXIS["geodetic latitude (Lat)",north,ORDER[1],ANGLEUNIT["degree",0.0174532925199433]],AXIS["geodetic longitude (Lon)",east,ORDER[2],ANGLEUNIT["degree",0.0174532925199433]],AXIS["ellipsoidal height (h)",up,ORDER[3],LENGTHUNIT["metre",1]],USAGE[SCOPE["Geodesy."],AREA["World."],BBOX[-90,-180,90,180]],ID["EPSG",7912]]</TargetSRS>
<Options>
<Option key="COORDINATE_OPERATION">+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=vgridshift +grids=us_noaa_g2018u0.tif +multiplier=1 +step +proj=cart +ellps=GRS80 +step +inv +proj=helmert +x=1.0053 +y=-1.90921 +z=-0.54157 +rx=0.02678138 +ry=-0.00042027 +rz=0.01093206 +s=0.00036891 +dx=0.00079 +dy=-0.0006 +dz=-0.00144 +drx=6.667e-05 +dry=-0.00075744 +drz=-5.133e-05 +ds=-7.201e-05 +t_epoch=2010 +convention=coordinate_frame +step +inv +proj=cart +ellps=GRS80 +step +proj=unitconvert +xy_in=rad +xy_out=deg +step +proj=axisswap +order=2,1</Option>
</Options>
</ReprojectionTransformer>
</ReprojectTransformer>
</GenImgProjTransformer>
</Transformer>
<BandList>
<BandMapping src="1" dst="1">
<SrcNoDataReal>-999999</SrcNoDataReal>
<SrcNoDataImag>0</SrcNoDataImag>
<DstNoDataReal>-999999</DstNoDataReal>
<DstNoDataImag>0</DstNoDataImag>
</BandMapping>
</BandList>
</GDALWarpOptions>
</VRTDataset>
43 changes: 43 additions & 0 deletions CONTRIBUTING.md
Original file line number Diff line number Diff line change
@@ -0,0 +1,43 @@
# Contributing Guide

Contributions are welcome here!

- Report bugs, request features or submit feedback as a [GitHub Issue](https://docs.github.com/en/issues/tracking-your-work-with-issues/about-issues).


## Fork this repository

We recommend first forking this repository and creating a local copy:

```
git clone https://github.com/YOURACCOUNT/3D_CRS_Transformation_Resources.git
cd xarray-tutorial
```

## Add content

Develop your new content on a branch. See [JupyterBook Docs](https://jupyterbook.org/en/stable/intro.html) for guides on adding `.md`, `.ipynb` and other content.

```
git checkout -b newcontent
```

## Preview your changes

Jupyter Book will execute notebooks and render HTML pages for the website. Be sure to fix any execution errors and preview the website in your web browser to make sure everything looks good!

```
pixi run website
```

Happy with how things look? Finally, open a pull request to merge you changes:

## Open a pull request

```
git add .
git commit -m "added pages x,y and improved z"
git push
```

Follow the link reported in a terminal to open a pull request!
Loading

0 comments on commit acf0787

Please sign in to comment.