Hello, First sorry for my bad English. Maybe I can help you and other people with the same doubt. Currently I' ve been working on a project that I need this functionality. def proximity(in_raster,

Hi all. I'm using Rasterio for a raster analysis but I need to calculate proximity as is done in gdal_proximity. https://gdal.org/programs/gdal_proximity.html Is this implemented somewhere in

Another update, for posterity. I tried doing the same reprojection using the gdal.Warp and it worked perfectly with far less effort around setting up transforms and whatnot. I filed an issue with

Just an update: I set the env variable CHECK_WITH_INVERT_PROJ to YES and was able to generate the output using the vertical slice but again, when I try to do this with the global raster I get a blank

I'll try my best to describe the issue although it's really confusing me so I'm not sure how well I can explain it. I'm trying to reproject some satellite imagery in the form of a raster that covers

https://rasterio.readthedocs.io/en/latest/cli.html#shapes is a command line program for extracting GeoJSON shapes from a raster. It's a thin wrapper around rasterio.features.dataset_features() (see

Hi, Rasterio's CRS.to_dict and .to_prog4 *do* return PROJ4 mappings and strings but we reduce to +init=epsg:xxxx when allowed. In the new world where there is no longer an epsg file to init from, it

Given a binary numpy array (array of 0s and 1s where 0=background), what would be the best way to polygonise this raster and write it to a geojson?

Thanks! This is exactly what I was looking for. BTW GDAL has an appropriate method to export CRS to proj4: >>> srs.ExportToProj4() '+proj=utm +zone=5 +datum=WGS84 +units=m +no_defs ' In my opinion it

pyproj.CRS has a to_proj4 method that should be able to do that: https://pyproj4.github.io/pyproj/stable/examples.html

Hi folks! Is there a way to get a PROJ.4 representation of CRS object? I've tried to use to_proj4() method:from rasterio.crs import CRSwkt = 'PROJCS["WGS 84 / UTM zone 5N",GEOGCS["WGS

Hi Eduardo, I'd suggest taking a look at the method argument that merge accepts. There are a few pre-defined methods, and you can also define your own function for selecting which pixel to select from

I'm creating a mosaic from many sentinel images with a bunch of overlapping areas, although I'm trying to understand what is criteria. For example: This is a small piece totally clear of clouds And

Hi all I am trying to figure out how to pass the height and width integer parameters to the from_bounds method. In my case they look not to be honored and the window returned has slightly different

Hi all, Version 1.0.1 of rio-color is on PyPI now: https://pypi.org/project/rio-color/1.0.1/#files. Thank you, Kyle Barron, for showing how to get started with cibuildwheel. Please note that there are

I think this is the problem: https://pyproj4.github.io/pyproj/stable/gotchas.html#internal-proj-error-sqlite-error-on-select

I just installed rasterio==1.2.0b1 on a fresh virtual env (with no GDAL/PROJ env set) and I'm getting a proj error Not sure what's going on :-(

Hi all, I got reports that the 1.0.1.dev0 wheels install and work, so I've made a more stable set and a 1.0.1a1 release. Target date for 1.0.1 is 2020-12-21. See

<sean@...> wrote: Hi Sean! Thanks a lot for your hint. My problem was caused by the fact that I was merging the polygons returned by rasterio.features.shape using unary_union. I

Hi all, Lack of an up to date wheel building system has been blocking us from making a 1.0.1 release. See https://github.com/mapbox/rio-color/issues/65. Kyle Barron suggested using cibuildwheel and