Indeed they are largely redundant. The meta property came first. Then we introduced a profile property to better support the pattern of using one dataset's profile as keyword arguments for creating new similar datasets. Note that the profile may contain items that don't appear in the meta property. 

On Fri, Nov 29, 2019, 9:57 AM Amine Aboufirass <amine.aboufirass@...> wrote:

Also, what is the difference between a dataset's metadata and profile? I have the impression one of these values is redundant. Please see example below:

>>> dataset.profile
{'driver': 'AAIGrid', 'dtype': 'int32', 'nodata': -9999.0, 'width': 8, 'height': 8, 'count': 1, 'crs': None, 'transform': Affine(1.0, 0.0, 0.0,
       0.0, -1.0, 8.0), 'tiled': False}
>>> dataset.meta
{'driver': 'AAIGrid', 'dtype': 'int32', 'nodata': -9999.0, 'width': 8, 'height': 8, 'count': 1, 'crs': None, 'transform': Affine(1.0, 0.0, 0.0,
       0.0, -1.0, 8.0)}

Also, what is the difference between a dataset's metadata and profile? I have the impression one of these values is redundant. Please see example below:

>>> dataset.profile
{'driver': 'AAIGrid', 'dtype': 'int32', 'nodata': -9999.0, 'width': 8, 'height': 8, 'count': 1, 'crs': None, 'transform': Affine(1.0, 0.0, 0.0,
       0.0, -1.0, 8.0), 'tiled': False}
>>> dataset.meta
{'driver': 'AAIGrid', 'dtype': 'int32', 'nodata': -9999.0, 'width': 8, 'height': 8, 'count': 1, 'crs': None, 'transform': Affine(1.0, 0.0, 0.0,
       0.0, -1.0, 8.0)}

Hi everyone,

After relying on GIS softwares for a long time with a bit of python help with rasterio, I decided to move on using rasterio (and geopandas) for AI applications.

I would like to have access to some common tools that I have been using in GIS softwares!!!. 

One important thing is the ability to store and manipulate the datasets using only one file such as the Esri FileGeodatabase (FileGDB). 

My question is: Are there any FileGDBs that we can use in Python to store and manipulate spatial data?

I found that Fiona can create, read and write FileGDB but only for storing shapefile kind of datasets.

thank you in advance, 

Armstrong Ngolo

Hi,

I think the reason you get an error is that you have misspelled the argument (you wrote maksed instead of masked).

FYI, the list of arguments of the read() function can be found in the documentation, here: https://rasterio.readthedocs.io/en/stable/api/rasterio.io.html#rasterio.io.DatasetReader.read

Best,

Guillaume Lostis

Please can someone explain how I align a tiff file of lights with another of population to create two numpy arrays the same size and aligned. I tried:

    population = rio.open(population)
    lights = rio.open(lights)
    res = rio.warp.reproject(rio.band(lights,1), rio.band(population,1))

This gives "CPLE_IllegalArgError: Destination band 1 appears to be read-only."

I was expecting it to read the meta data (crs, height, width and affine transform) from source and destination then return something with the data of source but the shape of destination.

Hi,

In the usage description of `rio bounds --help` 

```

...

 --sequence / --collection  Write a LF-delimited sequence of texts containing
                             individual objects (the default) or write a
                             single JSON text containing a feature collection
                             object.

...

 --collection               Output as GeoJSON feature collection(s).

```

confuses me. The first instance of `--collection` seems to imply that if several input files are passed like `rio bounds --collection *.tif`, then the output would be a single GeoJSON document with a single FeatureCollection containing multiple Features corresponding to each tif file expanded in the shell. The second instance of `--collection` is perhaps more honest in setting expectation that a FeatureCollection object will be created for each tif file. I'm wondering if the first I described was ever meant to be intended behavior? or has there been a regression? 

Thank you so much!

Hi all,

Simple question, but I want to make sure I'm doing this the canonical way. 

Assume two single-band rasters with identical projections/transforms and one of the rasters has NODATA pixels defining an irregular shape and 1 otherwise.  Is the best way to mask the other raster to use rio calc? e.g.:

  rio calc "(* (read 1) (read 2))" -o out.tif input.tif mask.tif

This works, but didn't know if a different utility is better suited.

thanks, matt

Dear All,

Is there a way to read and write hdf4 datasets using rasterio? I am working with DAAC data, for example https://e4ftl01.cr.usgs.gov//DP107/MOLA/MYD13Q1.006/2013.01.09/MYD13Q1.A2013009.h09v07.006.2015254175244.hdf

I found some examples of reading the data using gdal but extracting all relevant subdatasets and writing to a new GeoTiff using gdal would be a headache, I think.

Please let me know your thoughts. Thanks in advance for your help.

Regards,

Amine