pythonrasterioshapefile
This question refers to the doc of rasterio:
https://rasterio.readthedocs.io/en/latest/topics/masking-by-shapefile.html
I used the given code on the website.
out_meta.update({"driver": "GTiff",
"height": out_image.shape[1],
"width": out_image.shape[2],
"transform": out_transform})
with rasterio.open("RGB.byte.masked.tif", "w", **out_meta) as dest:
dest.write(out_image)
How do they create the RGB.byte.masked.tif file ?
Your code tries to convert from EPSG:4326 to EPSG:4326. But your data is actually in EPSG:3857 (wild guess on my end). You need to use EPSG:3857 as "source" CRS.
Try this (untested):
project = lambda x, y: pyproj.transform(
pyproj.Proj(init='epsg:3857'),
pyproj.Proj(init='epsg:4326'),
x,
y
)
yard = shapely.ops.transform(project, features)
GeoJSON had no specification until https://www.rfc-editor.org/rfc/rfc7946 so older GeoJSON might not be in WGS84 which rasterio might assume to be the case.
Your code works fine. If you look at SVDNB_npp_d20120301.vcld.tif (src.meta) you'll see nodata = 255 and you specified filled=True which means fill out_image with nodata values. So when you print your array, you are seeing the nodata values around the masked area of actual values.
If you write out the image and load it in a GIS, you'll see:
# To mask the data
with fiona.open("tehran.shp", "r") as shapefile:
shapes = [feature["geometry"] for feature in shapefile]
with rasterio.open("SVDNB_npp_d20120301.vcld.tif") as src:
out_image, transformed = rasterio.mask.mask(src, shapes, crop=True, filled=True)
out_profile = src.profile.copy()
out_profile.update({'width': out_image.shape[2],'height': out_image.shape[1], 'transform': transformed})
with rasterio.open("tehran.tif", 'w', **out_profile) as dst:
dst.write(out_image)
If you just want a 1D array of 0-5, just filter out the nodata
out_image[out_image<=5]]
Best Answer
The line below opens
"RGB.byte.masked.tif"for writing (note the"w"argument). It creates the file if it does not exist and completely overwrites it if it does exist.