awslabs / awslabs/open-data-docs
`.tif` files for some states in single year have heterogeneous projection
- Dominant language
- Jupyter Notebook
- Stars
- 180
- Forks
- 32
- PR merge metrics
- No merged PRs in 30d
Description
I am using the data in `s3://naip-visualization/` to generate `.vrt` files and later use those files to get the raster image based on `shapely` boundary with `rasterio`.
However this is working only in few cases. For instance, I was able to build this for RI all four years available 2012, 2014, 2016,2018. Also for AL for the year 2013. I am starting to see issues for few other states.
**AZ 2015**:
_Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2015/100cm/rgb/32114/m_3211435_se_11_1_20150714.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2015/100cm/rgb/32114/m_3211435_sw_11_1_20150708.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2015/100cm/rgb/32114/m_3211436_ne_11_1_20150725.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2015/100cm/rgb/32114/m_3211436_nw_11_1_20150708.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2015/100cm/rgb/32114/m_3211436_se_11_1_20150725.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2015/100cm/rgb/32114/m_3211436_sw_11_1_20150708.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2015/100cm/rgb/32114/m_3211437_ne_11_1_20150725.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2015/100cm/rgb/32114/m_3211437_nw_11_1_20150725.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2015/100cm/rgb/32114/m_3211437_se_11_1_20150725.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2015/100cm/rgb/32114/m_3211437_sw_11_1_20150725.tif_
**CO 2013**:
_Skipping /vsis3/naip-visualization/co/2013/100cm/rgb/38108/m_3810837_ne_12_1_20130920.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 13N, got NAD83 / UTM zone 12N. Skipping /vsis3/naip-visualization/co/2017/100cm/rgb/38108/m_3810837_se_12_1_20171007.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2017/60cm/rgb/32114/m_3211440_se_11_h_20170702.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 13N, got NAD83 / UTM zone 12N. Skipping /vsis3/naip-visualization/co/2013/100cm/rgb/38108/m_3810837_nw_12_1_20130920.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 13N, got NAD83 / UTM zone 12N. Skipping /vsis3/naip-visualization/co/2015/100cm/rgb/38108/m_3810833_ne_12_1_20150909.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 13N, got NAD83 / UTM zone 12N. Skipping /vsis3/naip-visualization/co/2017/100cm/rgb/38108/m_3810837_sw_12_1_20171007.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 13N, got NAD83 / UTM zone 12N. Skipping /vsis3/naip-visualization/co/2013/100cm/rgb/38108/m_3810837_se_12_1_20130920.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 13N, got NAD83 / UTM zone 12N. Skipping /vsis3/naip-visualization/co/2017/100cm/rgb/38108/m_3810838_ne_12_1_20170826.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 13N, got NAD83 / UTM zone 12N. Skipping /vsis3/naip-visualization/co/2013/100cm/rgb/38108/m_3810837_sw_12_1_20130920.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 12N, got NAD83 / UTM zone 11N. Skipping /vsis3/naip-visualization/az/2017/60cm/rgb/32114/m_3211440_sw_11_h_20170702.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 13N, got NAD83 / UTM zone 12N. Skipping /vsis3/naip-visualization/co/2017/100cm/rgb/38108/m_3810838_nw_12_1_20171007.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 13N, got NAD83 / UTM zone 12N. Skipping /vsis3/naip-visualization/co/2013/100cm/rgb/38108/m_3810838_ne_12_1_20130804.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 13N, got NAD83 / UTM zone 12N. Skipping /vsis3/naip-visualization/co/2015/100cm/rgb/38108/m_3810833_nw_12_1_20150909.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD83 / UTM zone 13N, got NAD83 / UTM zone 12N. Skipping /vsis3/naip-visualization/co/2017/100cm/rgb/38108/m_3810838_se_12_1_20171008.tif
Warning 1: gdalbuildvrt does not support heterogeneous projection: expected NAD8_
I did see this for few other states too. Now, There can be users directly indexing into files and in someway figure out which tif to access based on the boundary but this would make it tough to use `rasterio`.
VRT is really good and `rasterio` takes care of stitching images together for overlapping boundaries between quarter quad tiles.
Contributor guide
Assessment
This issue has not been assessed yet.