Michael,
can you paste the definition of such STAC Json item ? Guess:
"proj:epsg" is present, but "proj:bbox" or "proj:transform" are using
WGS84 coordinates ?
Even
Le 11/09/2024 à 15:59, Michael Smith via gdal-dev a écrit :
Using gdal 3.9.2:
I’m using pystac to create stac json item collections. I’m setting the
srs of using the projection extension and gdalinfo returns the correct
srs that I am setting.
gdalinfo STACIT:13459.json
Driver: VRT/Virtual Raster
Files:
/vsicurl/https://noaa-nos-coastal-lidar-pds.s3.us-east-1.amazonaws.com/dem/WA_King_DEM_2016_8589/./kingcounty_delivery4_block2_be.tif
/vsicurl/https://noaa-nos-coastal-lidar-pds.s3.us-east-1.amazonaws.com/dem/WA_King_DEM_2016_8589/./kingcounty_delivery4_block1_be.tif
Size is 74182, 79467
Coordinate System is:
COMPOUNDCRS["NAD83(HARN) / Washington North (ftUS) + NAVD88 height
(ftUS)",
PROJCRS["NAD83(HARN) / Washington North (ftUS)",
But when I look at the origin / pixel size / bounding box, its in
geodetic coordinates. Is this how it should be reported? I was also
doing some conversions from stac to vrt to support arcpro users and
the vrt then inherits the bbox/pixel size/origin of the geodetic
values from stac. Trying to see how I can get these values in actual
projected coordinates.
Mike
--
Michael Smith
Remote Sensing/GIS Center
US Army Corps of Engineers
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
--
http://www.spatialys.com
My software is free, but my time generally not.
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev