Michael Sumner mdsumner

https://gdal.org/en/stable/user/coordinate_epoch.html

ogrinfo /vsicurl/https://github.com/OSGeo/gdal/raw/refs/heads/master/autotest/ogr/data/flatgeobuf/test_ogr_flatgeobuf_coordinate_epoch.fgb foo -so

INFO: Open of `/vsicurl/https://github.com/OSGeo/gdal/raw/refs/heads/master/autotest/ogr/data/flatgeobuf/test_ogr_flatgeobuf_coordinate_epoch.fgb'
      using driver `FlatGeobuf' successful.

Sandbox: https://app.sandbox.dea.ga.gov.au/hub/login?next=%2Fhub%2F

remotes::install_github("mdsumner/nuyina.underway")

library(nuyina.underway)
d <- nuyina_underway() ## every time this runs it is up to date
range(d$datetime)
#[1] "2021-12-23 05:00:00 UTC" "2025-11-11 23:42:00 UTC"

gdal mdim mosaic.is here, a general rechunker, virtualize to (mdim) VRT or kerchunk parq https://gdal.org/en/stable/programs/gdal_mdim_mosaic.html this is afaik the first tool to input a wide variety of sources that can then serve GDAL or xarray interfaces https://github.com/OSGeo/gdal/blob/v3.12.0/NEWS.md

this query

https://earth-search.aws.element84.com/v1/search?collections=sentinel-2-c1-l2a&bbox=110.477,-66.332,110.577,-66.232&datetime=2021-11-22T00:00:00Z/2021-11-22T23:59:59Z

has only assets coastal, nir09, preview - other bands are missing.

this appears to be working (on setonix I just had to uv venv .venv and install virtual-tiff and it runs)

import obstore
from virtualizarr.registry import ObjectStoreRegistry
from virtual_tiff import VirtualTIFF
import xarray as xr

# Configuration
endpoint = "https://projects.pawsey.org.au/"

You should update your understanding

don't need anything client side to read a single polygon or field from within a shape file at a URL (or in a zip at a URL). Same goes for parquet

I think you're confused about tools you use.

accuracy matters as does not spreading bullshit

https://bsky.app/profile/postholer.com/post/3lmb6nxh7fs2x

minute check on this PR in GDAL for GetRawBlockInfo

OSGeo/gdal#13260

from osgeo import gdal
gdal.UseExceptions()
ds = gdal.OpenEx("/vsicurl/https://thredds.nci.org.au/thredds/fileServer/gb6/BRAN/BRAN2023/daily/ocean_temp_2010_01.nc", gdal.OF_MULTIDIM_RASTER)
info  = temp.GetRawBlockInfo([0,0,0,0])
info.GetOffset()

	library(shiny)
	library(shinyglide)

	ui <- fluidPage(
	glide(
	screen(
	h2("Image 1"),
	img(src = "image1.png", width = "100%")
	),
	screen(

	<VRTDataset>
	<Group name="/">
	<Dimension name="Time" type="TEMPORAL" size="5479" indexingVariable="Time" />
	<Dimension name="st_ocean" type="VERTICAL" direction="DOWN" size="51" indexingVariable="st_ocean" />
	<Dimension name="xt_ocean" type="HORIZONTAL_X" size="3600" indexingVariable="xt_ocean" />
	<Dimension name="yt_ocean" type="HORIZONTAL_Y" size="1500" indexingVariable="yt_ocean" />
	<Array name="Time">
	<DataType>Float64</DataType>
	<DimensionRef ref="Time" />
	<RegularlySpacedValues start="11323.5" increment="1" />