[gdal-dev] Duplicate pixels in VRT

[jratike80]

Hi,

Check what are the extents of two adjacert images reported by GDAL 2 and 3,
as well as the sizes in pixels, and if images are of type "pixel is area" or
"pixel is point". DEM kind of images are sometimes intentionally produced to
have overlap.

-Jukka Rahkonen-



Jon Morris wrote
> Hello all,
> 
> I've come across a weird edge case and I'm wondering if anyone has seen
> something similar.
> 
> I'm updating our application from GDAL v2.2 to v3.2 and have a failing
> test where there is a 1 pixel shift in the output. We're building a VRT
> from a DTM and a flood depth raster and when you call ds.ReadAsArray() on
> the VRT there is a duplicated column of pixels at GDAL 2, which is causing
> the apparent shift. It's not worth raising a ticket, as the GDAL 3 output
> is correct, but I can't explain why this shift has happened and I'm
> wondering if this has happened anywhere else.
> 
> If I open the two VRTs in QGIS, they look the same, so maybe it's an issue
> with ReadAsArray. There is slight difference in SRS precision, which might
> be causing a rounding error, but otherwise I've no idea what might be
> causing it.
> 
> GDAL 2
> <SRS>
> PROJCS["North_America_Albers_Equal_Area_Conic",
> GEOGCS["NAD83",DATUM["North_American_Datum_1983",
> SPHEROID["GRS
> 1980",6378137,298.2572221010042,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],
> PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4269"]],
> PROJECTION["Albers_Conic_Equal_Area"],
> PARAMETER["standard_parallel_1",20],
> PARAMETER["standard_parallel_2",60],
> PARAMETER["latitude_of_center",40],
> PARAMETER["longitude_of_center",-96],
> PARAMETER["false_easting",0],
> PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]]]
> </SRS>
> GDAL 3
> <SRS dataAxisToSRSAxisMapping="1,2">
> PROJCS["North_America_Albers_Equal_Area_Conic",
> GEOGCS["NAD83",DATUM["North_American_Datum_1983",
> SPHEROID["GRS
> 1980",6378137,298.257222101004,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],
> PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4269"]],
> PROJECTION["Albers_Conic_Equal_Area"],
> PARAMETER["latitude_of_center",40],
> PARAMETER["longitude_of_center",-96],
> PARAMETER["standard_parallel_1",20],
> PARAMETER["standard_parallel_2",60],
> PARAMETER["false_easting",0],
> PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH]]
> </SRS>
> Thanks,
> 
> Jon
> 
> 
> Jon Morris
> Software Developer
> 
> COVID-19. During the current outbreak JBA remains open for business and we
> continue to deliver our services. However, we have adopted flexible/remote
> working as required. I will be receiving and reading email as normal but I
> may not always be available on the office number.
> 
> 
> [JBA COVID-19
> statement]<https://www.jbagroup.co.uk/sites/www.jbagroup.co.uk/files/documents/15-030%20JBA%20Business%20Continuity%20Briefing%20-%20Latest.pdf>
> 
> T +44 (0) 1756 799919
> www.jbarisk.com<http://www.jbarisk.com>
> 
> [Visit our website]<http://www.jbarisk.com>
> [http://www.jbagroup.co.uk/imgstore/JBA-Email-Sig-Icons-LINKEDIN.png]
> <https://www.linkedin.com/in/jon-morris-a2897b4/>  [Follow us on
> Twitter] <https://twitter.com/jbarisk>
> 
> Find out more about us here:
> www.jbarisk.com<http://www.jbarisk.com/> and follow us on Twitter
> @JBARisk<http://twitter.com/JBARisk> and
> LinkedIn<https://www.linkedin.com/company/2370847?trk=tyah&trkInfo=clickedVertical%3Acompany%2CclickedEntityId%3A2370847%2Cidx%3A2-1-2%2CtarId%3A1447414259786%2Ctas%3AJBA%20RISK%20MANAGEMENT>
> 
> The JBA Group supports the JBA Trust.
> 
> All JBA Risk Management's email messages contain confidential information
> and are intended only for the individual(s) named. If you are not the
> named addressee you should not disseminate, distribute or copy this
> e-mail.
> Please notify the sender immediately by email if you have received this
> email by mistake and delete this email from your system.
> 
> 
> JBA Risk Management Limited is registered in England, company number
> 07732946, 1 Broughton Park, Old Lane North, Broughton, Skipton, North
> Yorkshire, BD23 3FD, Telephone: +441756799919
> 
> 
> 
> _______________________________________________
> gdal-dev mailing list

> gdal-dev at .osgeo

> https://lists.osgeo.org/mailman/listinfo/gdal-dev





--
Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html