search

ruithnadsteud / yt_georaster

A package for handling geotiff files and georeferenced datasets within yt.
1 stars 2 forks source link

Pass desired projection at load to be used as working projection for dataset #48

Closed deastwoo closed 3 years ago

deastwoo commented 3 years ago

Rather than purely basing the dataset's projection off the data we should add the option to pass a desired dataset CRS to reproject all the data bands to.

I think this will involve adding a keyword argument to yt.load (or maybe later in the process when we define a container?) taking a CRS string (e.g. "epsg:****") which would be readable by rasterio.

This could then be set as the dataset's crs and all bands would be reprojected as normal.

deastwoo commented 3 years ago

This can now be done using the crs keyword argument at load. ds = yt.load(filename, crs="epsg:****")

brittonsmith commented 3 years ago

@deastwoo we should perhaps add some documentation for this. Can you add a note in the sphinx docs about this?