/OSGeo/gdal/pull/12700 will fix it
>
> Even
> Le 04/07/2025 à 00:48, Michael Sumner via gdal-dev a écrit :
>
> I'm confused about the orientation of this in 2D mode 1485x1272 is
> correct, with 122 bands.
>
> (requires earthdata auth)
>
> gdal raster info "/vs
output-format text | grep "Size"
#Size is 1485, 1272
At any rate in the default case the geolocation arrays don't match the data
orientation.
Cheers, Mike
--
Michael Sumner
Ordinary Member, Streets People Love Hobart Association
Research Software Engineer
Australian
https://element84.com/software-engineering/is-zarr-the-new-cog/
--
Michael Sumner
Research Software Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https
p
Thanks, Mike
--
Michael Sumner
Research Software Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Awesome, thanks so much.
Cheers, Mike.
On Sun, May 11, 2025 at 8:29 AM Even Rouault
wrote:
>
> Le 11/05/2025 à 00:20, Michael Sumner via gdal-dev a écrit :
>
> Hello, I have a question about safely re-casting a dataset to another, I
> want to use Multidim to re-orient an MDArr
++auto poArrayDS = poMdimGroup->OpenMDArrayFromFullname(osArrayName);
auto poClassicDS = poArrayDS->AsClassicDataset(index_x, index_y);
// FIXME how to safely swap poSrcDS with poClassicDS
poSrcDS = std::unique_ptr(poClassicDS);
On Sun, May 11, 2025 at 8:20 AM Michael Sumner wrote:
>
stion is about the last line where poClassic "replaces" poSrcDS: what
am I responsible for to make that safe?
Thanks!
--
Michael Sumner
Research Software Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
anks,
>
> --
> Andrew Bell
> andrew.bell...@gmail.com
> ___
> gdal-dev mailing list
> gdal-dev@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
Michael Sumner
Research Software Engineer
Australian An
Ok, thank you that's totally understandable and I will champion your
choices here
Appreciate the (fast) response very much
Cheers, Mike
On Thu, Mar 27, 2025, 04:28 Even Rouault wrote:
> Michael,
>
>
> >
> > gdalmdiminfo
> > ZARR:"/vsicurl/
> https://data.source.coop/dynamical/noaa-gefs-analysi
curl/
https://data.source.coop/dynamical/noaa-gefs-analysis/v0.0.1-dev.zarr
this does not
gdalmdiminfo /vsicurl/
https://data.source.coop/dynamical/noaa-gefs-analysis/v0.0.1-dev.zarr
It's always the last thing I try first :)
Does it need to be this way for some special zarr reason?
Thanks
with-flyte-unionai-gdal-and-xarray-march-12-2025/4937
These presentations are recorded and will be available at the above link
afterwards.
Cheers, Mike
Michael Sumner
Research Software Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...
and the lack of overviews. I'm also a bit skeptical
> about the Python ecosystem. My impression is that Zarr just tries (and
> manages) to be a better NetCDF.
>
> Laurentiu
>
>
> On Thu, Feb 27, 2025, at 06:06, Michael Sumner via gdal-dev wrote:
>
> Just clueing into why
e,
>
> That's an interesting one for sure.
>
> Even before the ESA talk of Zarr, we have a good amount of Zarr data
> around. Just one I know of:
> https://cloud.google.com/storage/docs/public-datasets/era5
>
> On Wed, Feb 26, 2025 at 8:07 PM Michael Sumner wrote:
>
>> >> gdal-dev mailing list
>> >> gdal-dev@lists.osgeo.org
>> >> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>> >
>> > ___
>> > gdal-dev mailing list
>&g
into in practice. You're probably aware of this already, but
> instead of:
>
> > export AWS_S3_ENDPOINT=...
> > gdalmdiminfo ...
> > unset AWS_S3_ENDPOINT
> > gdalmdiminfo ...
>
> you can do:
>
> > AWS_S3_ENDPOINT=... gdalmdiminfo ...
> > gdalmdimin
to unset the option? Or, is an empty value a desirable way to
mean "default"?
Thanks for your help.
Cheers, Mike
--
Michael Sumner
Research Software Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
Where? When it's a fascist state we're all at fault.
I knew this wouldn't last but I absolutely supported it, silently, from the
sidelines, frightened of its consequences...
This is the slope. Good luck everyone, be brave
On Fri, Feb 14, 2025, 20:07 Rahkonen Jukka via gdal-dev <
gdal-dev@lists.o
D
Output driver not found.
But in this case, nothing happens and all is silent. (I will PR to set
extension metadata but wondered if this is a problem).
gdalmdimtranslate gdrivers/data/netcdf/srid.nc out.zarr -co
ARRAY:COMPRESS=ZSTD
Cheers, Mike
--
Michael Sumner
Research Software Engi
Hmm, I see CSV special-case is close to what I need so I'll try modifying
this along those lines.
We could have -type "none", and a null value for "all bands" - and use the
CSV special-case logic to write to other tabular formats.
Cheers, Mike
On Thu, Dec 12, 2024 a
:
> https://github.com/OSGeo/gdal/blob/master/apps/gdal2ogr.c
>
> Even
> Le 12/12/2024 à 06:40, Michael Sumner via gdal-dev a écrit :
>
> Hi, I'd like to be able to create a tabular version efficiently directly
> from a raster. I can do this in scripting but is there a utility
1
b (Integer) = 7
OGRFeature(file):1
a (Integer) = 3
b (Integer) = 9
OGRFeature(file):2
a (Integer) = 5
b (Integer) = 11
OGRFeature(file):3
a (Integer) = 2
b (Integer) = 8
OGRFeature(file):4
a (Integer) = 4
b (Integer) = 10
OGRFeature(file):5
a (Integer) = 6
b (Integer) = 12
ppb vol
>
>
> NETCDF:"out.input.nc
> ":NO2
> 1
>DataType="Float32" BlockXSize="174" BlockYSize="1" />
>
>
>
>
>
>
> On Tue, 3 Dec 2024 at 12:02, Michael Sumner
Try
gdalinfo NETCDF:"out.input.nc":NO2
If you see those lon and lat subdatasets listed under "geolocation" then
you can use gdalwarp to create a georeferenced dataset, automatically with
gdalwarp sds.in out.tif or by specifying a grid with some of
-tr,-te,-t_srs,-ts.
If they don't appear they ca
I the first step you are a_ssigning moll but with a longlat ullr (extent )
a_ssigned.
This isn't correct, the numbers will need to be something like -A A/2 A
-A/2 where A is something like pi * radius i.e pi * 6378137
Your source might have this information stored somewhere, or a guess and a
pray
It's not RStudio or sf specific though, I think it's worth asking here
because there are many many R packages using GDAL and they essentially
share the same simple installation, after the GDAL version choice has been
made. What really matters is your operating system and the method you have
availab
Nice one, thanks! (I had a fomenting PR to implement on the fly geolocation
array config, so I'll use your hallucination for the design phase).
🤟
On Sat, 19 Oct 2024, 14:44 Even Rouault, wrote:
>
> Le 18/10/2024 à 23:32, Michael Sumner a écrit :
>
> I didn't know you
I didn't know you could do that with -to!! That's awesome
🤟
On Sat, 19 Oct 2024, 05:01 Even Rouault via gdal-dev, <
gdal-dev@lists.osgeo.org> wrote:
> Conrad,
>
> Try something like:
>
> gdal_translate input.tif imagery.vrt -b 3
>
> gdalwarp imagery.vrt imagery_warped.tif -geoloc -to X_DATASET=i
1
> Block=4865x1 Type=Int32, ColorInterp=Gray
>
> How can I take advantage of VRT? That sounds like a great potential
> solution.
> Thanks
> Conrad
>
> On Friday, October 18, 2024 at 12:01:27 PM GMT+1, Michael Sumner <
> mdsum...@gmail.com> wrote:
>
>
> C
Check if gdalinfo sees the longitude and latitude as geolocation arrays,
then just warp to grid you want.
If they aren't auto identified this can be arranged with VRT.
Cheers, Mike
On Fri, 18 Oct 2024, 21:55 Conrad Bielski via gdal-dev, <
gdal-dev@lists.osgeo.org> wrote:
> Hi Javier,
>
> it is
n
that gave a clue, it's below 4x100x100, and doesn't seem related to single
chunk or multiple.
Cheers, Mike
--
Michael Sumner
Research Software Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
gdal-dev mai
o/gdal/pull/11005
>
> Otherwise just do "gdal_translate Rlogo.png Rlogo.zarr -of Zarr" . By
> default the Zarr driver translates GDAL bands as the Z dimension, unless
> you specify the SINGLE_ARRAY=NO creation option, in which case it will
> create a Zarr array per band (so a di
ensions()]
#[4, 100, 100]
That looks like a bug?
And, is there a better way to go about this? PNG->ZARR is my target, I'm
fine with programmatic creation of NetCDF, but this could be quite handy as
a two-step conversion from pseudo-3D to actual 3D via the cli with
"gdal_translate
> Michael,
> > >
> > > my understanding of https://geoarrow.org/format.html#memory-layouts
> is that what writes OGR is supposed to be fine since they mentionned types
> like 'List[2]>>'. Perhaps I've missed something
> or nanoarrow has st
I realize I left out the INTERLEAVING, ie.
ogr2ogr ~/fromgdal.arrow
ogr/data/arrow/from_paleolimbot_geoarrow/polygon-default.ipc -lco
GEOMETRY_ENCODING=GEOARROW_INTERLEAVED
but still, I get these list>>}>>
On Tue, Oct 8, 2024 at 12:19 AM Michael Sumner wrote:
> When I investi
s of
remaining data but found 2266 bytes in buffer
Are we in-between moves regarding specifications, or something? I'm having
good results generally and this seems like a problem in the Arrow driver
for write.
Cheers, Mike
--
Michael Sumner
Research Software Engineer
Australian Antarcti
gt; them.
>
> And obviously one strong assumption of Kerchunk is that the files
> referenced by a Kerchunk index are immutable. If for some reason, tiles are
> moved internally because of updates, chaos will arise due to (offset, size)
> tuples being out of sync.
>
> Even
&
On Wed, Jul 24, 2024 at 8:37 AM Michael Sumner wrote:
> Hi, is there any effort or thought into something like Python's kerchunk
> in GDAL? (my summary of kerchunk is below)
>
> https://github.com/fsspec/kerchunk
>
> I'll be exploring the python outputs in detai
binding for Kerchunk?
> I think that will be the main blocker for GDAL driver development.
>
> [1] https://github.com/hyoklee/kerchunk/wiki
>
> ---
> Reality is hierarchical. Store scientific reality in HDF for Spatial
> Computing.
>
>
>
> ________
eral encoding
possibilities than tif/png/jpeg provide.) This scheme is extremely
general, literally a virtualized array-like abstraction on any storage,
and with kerchunk you can transcend many legacy issues with actual formats.
Cheers, Mike
--
Michael Sumner
Research Software Engineer
Australia
Just an update, I have found out how to do this numerically.Will
consider a PR for gdaldem if I can make it work at some point.
Thanks, Mike
On Tue, Jul 9, 2024 at 7:31 AM Michael Sumner wrote:
> Is it possible to "layer" one image over another, the first at full
> tran
quot;color-relief"
meant, at first).
Cheers, Mike
--
Michael Sumner
Research Software Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
27;t cover more
> than 360 degrees of longitude. This value, when computed, is passed as a
> hint OGRCoordinateTransformation so that it can post-correct longitudes to
> apply a +/- 360 degree offset, to be in the range of the source dataset.
> I've relaxed the sanity check to al
un 1, 2024 at 4:19 PM Michael Sumner wrote:
> This data source has an odd georeferencing, it's a 36000x17999 raster in
> the extent -179.995 -89.995 180.005 89.995.
>
> vrt://NETCDF:/vsicurl/
> https://archive.podaac.earthdata.nasa.gov/podaac-ops-cumulus-protected/MUR-JPL-L4-GLO
in VRT but I'd rather like this workflow to work as-is.
(Chasing down the weird grid referencing is another project, but it's a
really massive dataset so I'm pretty unconfident of that occurring).
Cheers, Mike
--
Michael Sumner
Research Software Engineer
Australian Antarctic Division
Ho
et after breaking
the COG layout?
## now open in rw, to zap band level md (but breaks COG layout)
#lds = gdal.OpenEx(tmp, gdal.GA_Update, open_options =
["IGNORE_COG_LAYOUT_BREAK=YES"])
#lds.GetRasterBand(1).SetMetadata({})
#lds.Close()
#Warning 1: tmp1.tif: The IFD has been rewritten
lue would be good enough to handle that outlier.
>
> (Someone, super bored, should author a book about with horror stories with
> netCDF and HDF georeferencing. There's a lot of material. Although likely
> not to be a best seller)
>
> Even
> Le 30/04/2024 à 09:45, M
Ok so my naive edits are clearly not enough, they still get written as Byte
so it's deeper in the target spec and spread across a few places I'm not
ready to get across yet.
Happy to pursue in the longer term though.
Cheers, Mike
On Tue, Apr 30, 2024 at 7:00 AM Michael Sumner wrote:
ve to keep notes and explore more deeply with native HDF5 tools.
Cheers, Mike
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
htt
t; ERROR 3: Cannot read 4029 bytes
>
> gdalinfo failed - unable to open '/vsicurl/
> https://login:pas...@n5eil01u.ecs.nsidc.org/AMSA/AU_SI12.001/2012.07.02/AMSR_U2_L3_SeaIce12km_B04_20120702.he5
> '.
>
>
>
>
>
> Thanks
>
>
>
> Joaquim
>
>
>
> *Fr
se nations, please visit;
> http://www.utemountainutetribe.com/
> http://www.cheyennenation.com/
> https://cheyenneandarapaho-nsn.gov/
> https://native-land.ca/
>
> _______
> gdal-dev mailing list
> gdal-dev@lists.osgeo.org
> https:
0, 435.000) ( 42d14'27.21"E, 39d11'27.54"S)
> Lower Right ( 395.000,-395.000) (135d 0' 0.00"E, 41d23'59.41"S)
> Center ( 0.000, 20.000) ( 0d 0' 0.01"E, 88d 8'51.76"S)
> Band 1 Bl
he providers and get the right
metadata baked in ... but there's pros and cons to either. I'm not sure
there's even enough information in these files to clearly detect the
situation, it would be a bit like the NSIDCbin driver with its very strict
requirements.
Cheers, Mike
--
Mich
.warp on gdal.org directly :
> https://gdal.org/search.html?q=gdal.warp&check_keywords=yes&area=default
>
> that leads to https://gdal.org/api/python/utilities.html#osgeo.gdal.Warp
>
> Even
> Le 29/04/2024 à 00:52, Michael Sumner via gdal-dev a écrit :
>
> I'm confu
or
"Warp(" only results in RFC 59 references, so I think this might be a real
problem).
If I'm just navigating the API docs incorrectly I'd appreciate a
walkthrough to find the signature and param descripts for "Warp()". Thanks!
Cheers, Mike
--
Michael Sumner
Softwar
sgeo.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/listinf
VSIFTellL(f)
gdal.VSIFSeekL(f, 0, 0) # begin
## the desired result is this
data = gdal.VSIFReadL(1, size, f)
## and clean up
gdal.VSIFCloseL(f)
gdal.Unlink(temp_name)
Cheers, Mike
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
> the 1D coordinate array with a very small relative epsilon, like 1e-8, is
> sufficient.
>
> Even
>
>
> Le 03/04/2024 à 22:56, Michael Sumner via gdal-dev a écrit :
>
> Here's an (ahem) extremely important discussion on the prospects for
> xarray to extend from t
as well, very happy to
discuss here or offline about any of this. I'm not really a technical
expert but I have a good overview of the landscape that GDAL sits in, and
I'm making similar pitches for involvement in the R communities.
Cheers, Mike
--
Michael Sumner
Software and Database
Oh, my bad. It's documented to split into SDS for multiple CRS.
Ouch that's not how I thought it was working, but it makes sense.
Thanks, Mike
On Thu, Apr 4, 2024 at 1:26 AM Michael Sumner wrote:
> this works for me, there are ten items in the filelist:
>
> gdalinfo
deepest I've been down in that file
... I hope to be able to craft these suggestions at some point.
Cheers, Mike
> Even
>
> --
> http://www.spatialys.com
> My software is free, but my time generally not.
>
>
--
Michael Sumner
Software and Database Engineer
Australia
CIT:\"
https://earth-search.aws.element84.com/v1/search?collections=sentinel-2-c1-l2a&bbox=179,-17,180,-16&datetime=2023-01-01T00:00:00Z/2023-01-03T23:59:59Z\
":asset=visual"
It feels like something awry with the limit for the server, or the driver,
or maybe interacting.
But,
well actually, I think what I'm asking for is the intended behaviour, but
there's an error.
Is it meant to detect sets of variables on 1D dimensions and present them
as layers? That's what would make sense to me.
Still exploring.
Cheers, Mike
On Tue, Apr 2, 2024 at 5:36 AM
ims in a general way, I think, but degenerate
dimensions might be worth consideration )
Appreciate any thoughts, thanks! I'd quite like to have the vector-approach
work as well as the mdim approach, I think they are nicely complementary
and provide different pros and cons.
Cheers, Mike
--
Micha
values in the raster - that is "value" in the way I was
thinking ( I had in mind more general cases like interval ranges, but
that's another level up)
I might be entirely messing up the concept ... will check my take by
experiment.
Cheers, Mike
--
Michael Sumner
Software and Da
t is 0-255 range values, but
is index 0-based or 1? And, can index be *value* in any contexts?
index red green blue
If it's documented I'll link it in PRs. I see qml and qlr here:
https://docs.qgis.org/3.34/en/docs/user_manual/appendices/qgis_file_formats.html
Cheers, Mike
--
Michael
ws
> ( I will explore how it's getting done).
>
> It looks like they might use (previous_size + 1) / 2 (or ceil rounding)
>
> -- http://www.spatialys.com
> My software is free, but my time generally not.
>
>
--
Michael Sumner
Software and Database Engineer
Australian Ant
views")
# Overviews: 18000x9000, 9000x4500, 4500x2250, 2250x1125, 1125x563
Cheers, Mike
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
p.emodnet.eu/erddap/files/biology_6640_benthos_NorthSea_e4af_0f0e_6a73/04_2021_6640_diva_benthos_erddap.nc'
not recognized as being in a supported file format.
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Austra
SSE instructions should
> be updated.
>
> @Even knowing now that the "old" hardware is "virtually old", should we
> remove AVX2 compatibility from ubuntu-full-latest? I do not know how much
> is the performance impact.
>
> Cheers,
> Javier.
>
> On Mon
deleting the corresponding .so in
>> /usr/lib/x86_64-linux-gnu/gdalplugins
>>
>> Even
>> Le 04/02/2024 à 22:51, Michael Sumner a écrit :
>>
>> indeed there's no avx2:
>>
>> cat /proc/cpuinfo|grep sse|head -n 1
>> flags : fpu
-with-external-depencies-built-as-plugin GDAL build mode, so it is
> easy to just remove a given plugin by deleting the corresponding .so in
> /usr/lib/x86_64-linux-gnu/gdalplugins
>
> Even
> Le 04/02/2024 à 22:51, Michael Sumner a écrit :
>
> indeed there's no avx2:
>
>
9% the reason of the issue.
>
> Even
> Le 04/02/2024 à 06:20, Michael Sumner a écrit :
>
> skipping TileDB does fix:
>
> ogr2ogr /tmp/newdir
> https://github.com/SymbolixAU/geojsonsf/raw/master/inst/examples/geo_melbourne.geojson
> -f
> "ESRI Shapefile"
> e
ore dumped)
>
> Cheers, Mike
>
>
>
>
> On Sat, Feb 3, 2024 at 12:46 PM Even Rouault
> wrote:
>
>> Michael,
>>
>> I'm wondering if there not might be something wrong with your build or
>> runtime environment. Or there's something subtle, because
--args ogrinfo
> /tmp/newdir/" (type "run") to get more useful information
>
> Even
> Le 03/02/2024 à 02:35, Michael Sumner via gdal-dev a écrit :
>
> I'm getting Illegal instruction / core dumped on ogrinfo of a directory:
>
> ogr2ogr /tmp/newdir
> ht
through some docker images and it wasn't a problem in
3.6.0, but I'm getting it since 3.7.0 - or I'm doing something wrong
entirely.
Cheers, Mike
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, A
s well
>
> -- http://www.spatialys.com
> My software is free, but my time generally not.
>
>
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
gdal-dev mailing li
;t think of a reason why your script wouldn't work with a FlatGeoBuf
> output (did you get an error message?), which would be much better, as you
> could also attach the XRES, YRES, BANDCOUNT, etc. metadata items that would
> help the GTI driver avoid probing one of the tiles.
>
and in
> > initial ticket in https://github.com/OSGeo/gdal/issues/8861)
> >
> > Even
> >
> --
> http://www.spatialys.com
> My software is free, but my time generally not.
>
> ___
> gdal-dev mailing l
om/
> http://www.cheyennenation.com/
> https://cheyenneandarapaho-nsn.gov/
> https://native-land.ca/
>
> ___
> gdal-dev mailing list
> gdal-dev@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
L_OF_RASTER, nullptr, nullptr, nullptr);
> CPLFree(pszSubdatasetSource);
> break;
>}
>else {
> CPLFree(pszSubdatasetSource);
>}
> }
>}
>
>poSrcDS->ReleaseRef();
>return 1;
> }
>
> Even
>
>
std::cout << pszSubdatasetSource << "\n";
poSrcDS->ReleaseRef();
poSrcDS = GDALDataset::Open(pszSubdatasetSource, GDAL_OF_RASTER,
nullptr, nullptr, nullptr);
CPLFree(pszSubdatasetSource);
GDALDestroySubdatasetInfo(info);
}
}
}
t set it at all ... (better to
use GDALwarp() for general cases IMO).
But, is there a good reason why setting that open option affects translate
for sources with no overviews?
I'm in master at ebac6b74a429fa6eeeb94edd074b6cb60072a35f
Cheers, Mike
--
Michael Sumner
Software and Database Engineer
ent instances of a
vertex really do stay unique). I guess the result is going to be pretty
good if you choose the precision appropriately and modifications aren't
extreme?
Thanks for consideration!
Cheers, Mike
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
what's the task? what about batching the geometry and or fields? can you
run on just the first feature, does that work?
how many features, how big is the task?
On Wed, 27 Sept 2023, 13:16 Scott via gdal-dev,
wrote:
> Any tips for using ogr2ogr to use only a specified amount of RAM? I'm
> not ha
e can incrementally build up the
SUM-pyramid with external overivews via CLI it would be awesome to have the
same stats for the two resampling regimes.
Thank you
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmai
cool thanks, I'll deal with that 😃
On Thu, 21 Sept 2023, 18:48 Even Rouault,
wrote:
>
> >
> > which don't seem to exist. Is that just an outdated comment?
>
> yes, as far as I can see from history those functions have never
> existed. This comment dates back to the initial commit and they didn
On Wed, Sep 20, 2023 at 6:40 PM Even Rouault
wrote:
> Sounds OK to me
> Le 20/09/2023 à 02:32, Michael Sumner a écrit :
>
> Thanks Even, could add a AttachDomainMetadata here that expected "-dmo
> DOMAIN:META-TAG=VALUE"
>
>
t:// con so could do
vrt://{dsn}?dmo=DOMAIN:META-TAG=VALUE
Cheers, Mike
On Wed, Sep 20, 2023 at 10:19 AM Even Rouault
wrote:
> Not from cli. From the API, it is dataset.SetMetadataItem(key, value,
> domain)
> Le 20/09/2023 à 02:14, Michael Sumner a écrit :
>
> When we add metad
important
Is that possible from cli?
Thank you
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org
yte), and was only honouring the buffer data type
> of the RasterIO() request from the AAIGrid driver (Int32, always used when
> translating from source bands having integer data types, but with the
> reasoning that if you query a Byte band then you would get always value in
> [0,255] r
he
input values and output types.
https://gist.github.com/mdsumner/ee4103d8616b9aa341e82c46b44a8c8c
Cheers, Mike
--
Michael Sumner
Research Software Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
gdal-dev mailing list
gdal-dev
ame=Band1,view=[::2,::4]'])
> nds.GetRootGroup().OpenMDArray("Band1").GetShape()
> # (20, 20)
>
> This was a bug. Fixed per https://github.com/OSGeo/gdal/pull/8297
>
> A way of working it around is to use scaleAxesSpecs=['x(4)', 'y(2)']
> instead
&g
quot;Band1").GetView("[::2,::4]").GetShape()
# (10, 5)
Am I specifying the view wrong in MultiDimTranslate?
Thanks, Mike
I've put fuller notes here in case the email is a problem:
https://gist.github.com/mdsumner/028d089ad9a960267347f891c058fc9a
--
Michael Sumner
Software and
just say to shove it through the warper because that is now so fast and
reliable.
Cheers, MIke
On Thu, Aug 10, 2023 at 7:17 AM Michael Sumner wrote:
> there is no georeferencing with this one, but as you say it has internal
> coordinates (which GDAL interprets as "geolocation arra
ut
>
> ___
> gdal-dev mailing list
> gdal-dev@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gma
!
On Wed, Aug 2, 2023 at 5:23 PM Michael Sumner wrote:
> I'm using a derived subdataset dsn to try to open a STAC query:
>
> gdalinfo "STACIT:\"
> https://explorer.sandbox.dea.ga.gov.au/stac/search?bbox=122.1,-18.28,122.48,-17.91&datetime=2020-01-01/2020-02-29\
>
ail problems with the strings:
https://gist.github.com/mdsumner/90f5b509c5cb32bc0f1a26c695baa0bf
I'm using gdal dev GDAL 3.8.0dev-bba8edfc13, released 2023/08/01 (debug
build)
Thanks, Mike
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
_
ntended to be part of the GDAL interface, but there's still quite a bit
> that can be added to the documentation. I can create an issue and post the
> script I used to scan for these if others are interested in filling some of
> the gaps.
>
> Dan
>
>
>
>
>
>
>
us setting resampling with the API, and I gather also
that there is not an analogous config option for the warper.
Thank you, Mike
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
__
(i) is not a Longitude/X dimension.
Warning 1: dimension #0 (j) is not a Latitude/Y dimension.
Thank you
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
1 - 100 of 203 matches
Mail list logo
