* docs: 📝 clarify usage of rpc_height in RPCTransformer usage example
Detail using a constant height offset via `rpc_height` keyword argument in instantiation of transformer.
* Apply suggestions from code review
Co-authored-by: Batuhan Kavlak <batuhan.kavlak@up42.com>
---------
Co-authored-by: Sean Gillies <sean.gillies@gmail.com>
Co-authored-by: Batuhan Kavlak <batuhan.kavlak@up42.com>
* Failing test with a .msk sidecar
* POC of sidecar support for Python openers
Towards resolving #3027
* Register opener for the parent directory/folder of the given URL
Handling the error cases properly gets all our tests passing again.
* Add opener adapter classes and a factory
Still a lot of redundant code here, needs refactoring.
* Remove unused import
* Don't raise exceptions when files aren't found
Just return 1 or NULL
* Add two tests of using fsspec filesystem objs and fix a bug.
* Allow testing with tiledb.VFS
* More tests of fsspec http filesystem and .msk sidecar
* Module clean up
* Add mtime to stats, docstrings to abstract class
* Expand documentation
* Reassign name_list in loop
* Return -1 on stat failure, remove useless except clauses
* Don't set an error when OpenFile doesn't resolve
* Add eof callback
* Support writing to Python files
This depends on PR #2898.
Resolves#2906.
* Update change log
* Fix PR reference
* Use CPLError to convey opener failure out of the callback
* Call CPLError with 4 arguments
* Register opener for filename and mode
Also make sure that deregistration happens on cleanup
* Remove with from test
* Update opener documentation
---------
Co-authored-by: Sean Gillies <seangillies@Seans-MacBook-Air.local>
* Rough draft of a new VSI plugin
Resolves#2888 maybe
* Getting closer
* Working implementation!
The missing piece was keeping a reference to opened files.
* Add more tests
* Add more documentation about opener
* Improvements based on feedback from @rouault
* Update docstrings
* Update the VSI topic in docs
* Don't accept fsspec OpenFile instances as opener
Instead, we'll make a way for rasterio.open to accept them as the
"fp" positional argument
* Store contexts of vsipyopener file objects
Eliminates the fsspec OpenFile special case
* Strip file mode for the zipfile case
* Update rasterio/_vsiopener.pyx
Co-authored-by: Alan D. Snow <alansnow21@gmail.com>
---------
Co-authored-by: Sean Gillies <seangillies@Seans-MacBook-Air.local>
Co-authored-by: Alan D. Snow <alansnow21@gmail.com>
Raster data of dtype = int8 leads to an error message from pprint, which requires the data to be one of the following:
int16, int32, uint8, uint16, float32
I added a comment catching this error, which can be solved by casting the data: data = data.astype('int16')
* add checks in order to warn agaisnt incorrect source CRS when using RPCs to reproject
* use .from_string to be consistent with _calculate_default_transform
* raise error instead
* fix test
* initial
* add TransformDirection enum, fix cython compiler issues
* make RPCTransformWarning more generic
* include GDALGCP* in gdal.pxi
* add GCPTransformer
* typo
* move xy and rowcol to TransformerBase
* change existing transform methods to use new *Transformer classes
* normalize transform direction, update docstrings
* update docstrings
* make subsubheadings of Coordinate Transform subheading
* Add new docs transform topic
* add to index
* update docstrings
* update docstring
* Add rasterio.rpc to docs
* update georeferencing.rst
* update transform.srst
* relax isinstance check for RPCTransformer
* catch invalid transform input to get_transformer
* replace None with identity transform so we test the correct ValueError
* update AffineTransformer
* add tests
* rename method to map singleton coordinate to array-like and validate all coordiantes have same dimensions
* update test
* add more tests
* update docstring of GCPTransformer init
* fix copy+paste error
* relax check on xy transformed values due to floating point precision
* remove unused imports
* fix tests and address comments
Co-authored-by: Sean Gillies <sean@mapbox.com>
* Various fixes to documentation
* Remove in-line shell prompts
* Fix rasterize shapes type
* Convert todo's to comments
* Revert back to todo directive
* [warp_add_rpcs] initial update to reproject to support rpcs
* catch kwargs in _calculate_default_transform
* fix allocation of crs to wrong dataset (was src, now dst) in _calculate_default_transform. Specifically where we went from GDALCreateGenImgProjTransformer -> GDALCreateGenImgProjTransformer2
* pass kwargs to calculate_default_transform (to be passed to transformer
* update docstrings
* update documentation
* add debug messages for kwargs used to set Transformer options
* Don't capitalize RPC_DEM transformer option value, as it might be a path
* includes existing changes from https://github.com/mapbox/rasterio/pull/2030, better support for computation of calculate_default_transform with just gcps or rpcs
* adds testse
* Remove an empty line
Mainly to trigger a new build on travis-ci.com
* add new test file (Largely derived from RGB.byte.tif)
* Fix incorrectly encoding RPC metadata as transform options
Was checking for transformer option RPC_DEM in the RPC metadata by mistake
Co-authored-by: Sebastien Tremblay-Johnston <sebastien.tremblayjohnston@mdacorporation.com>
Co-authored-by: Sean Gillies <sean@mapbox.com>
Co-authored-by: Sean Gillies <sean.gillies@gmail.com>
* This is part of an effort to introduce the concept of rational polynomial coeffecients to be used for coordinate transformations using GDALCreateRPCTransformer.
- Added a new UserDict subclass named RPC to control.py. This structure contains the rationa polynomial coefficients, and can be instatiated from user input or read from dataset. GDAL stores RPCs as metadata which are all string encoded. The RPC object deserializes string values to float or list appropriately, and also provides a serialization method when passing to lower level GDAL methods.
- A new _transform.pyx method called _rpc_transform, similar to _base._transform except using GDALCreateRPCTransformer and GDALRPCTransform.
- A new RPCTransformWarning to errors.py, indicating when a RPC coordinate transformation has failed.
- New includes from gdal.h and gdal_alg.h for RPC related structs and functions
- New rpcs property to DatasetBase and subclassese with appropriate getter and setters (error when read only mode). Corresponding updates to _base.pxd and _io.pxd
- Tests
* trigger new travis build
* This is part of an effort to introduce the concept of rational polynomial coeffecients to be used for coordinate transformations using GDALCreateRPCTransformer.
- Added a new UserDict subclass named RPC to control.py. This structure contains the rationa polynomial coefficients, and can be instatiated from user input or read from dataset. GDAL stores RPCs as metadata which are all string encoded. The RPC object deserializes string values to float or list appropriately, and also provides a serialization method when passing to lower level GDAL methods.
- A new _transform.pyx method called _rpc_transform, similar to _base._transform except using GDALCreateRPCTransformer and GDALRPCTransform.
- A new RPCTransformWarning to errors.py, indicating when a RPC coordinate transformation has failed.
- New includes from gdal.h and gdal_alg.h for RPC related structs and functions
- New rpcs property to DatasetBase and subclassese with appropriate getter and setters (error when read only mode). Corresponding updates to _base.pxd and _io.pxd
- Tests
trigger new travis build
* Move RPC Class to it's own file
* fix typo
* make private
* re-implement RPC with attr
* make err_bias, err_rand attributes optional, default to empty string
* default to None, rather than empty string
* update the threadpool example
* add benchmark number, done a VM with 4 cores
* fixes threading typo and adds a little more background on what multhreading is trying to achieve
