Vincent Sarago

try with : ``` s3://naip-analytic/ca/2020/60cm/rgbir_cog/32114/m_3211419_nw_11_060_20200524.tif s3://naip-analytic/ca/2020/60cm/rgbir_cog/32114/m_3211420_ne_11_060_20200524.tif s3://naip-analytic/ca/2020/60cm/rgbir_cog/32114/m_3211420_nw_11_060_20200524.tif ```

🤦 my bad it should be ``` /vsis3/naip-analytic/ca/2020/60cm/rgbir_cog/32114/m_3211419_nw_11_060_20200524.tif /vsis3/naip-analytic/ca/2020/60cm/rgbir_cog/32114/m_3211420_ne_11_060_20200524.tif /vsis3/naip-analytic/ca/2020/60cm/rgbir_cog/32114/m_3211420_nw_11_060_20200524.tif ``` GDAL uses `/vsis3/` prefix and not s3://

it worked for me 🤷‍♂️

``` $ cat test.txt /vsis3/naip-analytic/ca/2020/60cm/rgbir_cog/32114/m_3211419_nw_11_060_20200524.tif /vsis3/naip-analytic/ca/2020/60cm/rgbir_cog/32114/m_3211420_ne_11_060_20200524.tif /vsis3/naip-analytic/ca/2020/60cm/rgbir_cog/32114/m_3211420_nw_11_060_20200524.tif $ gdalbuildvrt -input_file_list test.txt ca_32114_urls.vrt --config AWS_REGION "us-west-2" --config AWS_REQUEST_PAYER "requester" 0...10...20...30...40...50...60...70...80...90...100 - done. $ cat ca_32114_urls.vrt PROJCS["NAD83 / UTM zone 11N",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4269"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-117],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","26911"]]...

> I'm guessing they are set as environment variables. yes they are in my env by default so I omitted those

update: ```python from rasterio._features import _shapes from rasterio.rio.helpers import coords from rasterio.transform import IDENTITY from shapely import geometry cpdef rgb_to_hex(rgba): return '#{:02x}{:02x}{:02x}'.format(*rgba) cpdef bytes shapes_encoder( data, mask, str layer_name =...

merge #2 in master but I'm getting shapely warnings `Self-intersection at or near point`! I guess I should really try to use wagyu as mentioned in https://github.com/cogeotiff/rio-tiler-mvt/pull/2#issuecomment-764897711

maybe we should just use pygeos instead of shapely 🤷‍♂️

pygeos 0.10 has some nice improvement and may be worth to try https://github.com/pygeos/pygeos/releases/tag/0.10

you can try this link: https://join.slack.com/t/cogeotiff/shared_invite/zt-17by75auf-39YWJEM08MkwEyM67TqeDg