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
Thank you.
Cheers, Mike
On Sun, Apr 30, 2023 at 1:07 AM Even Rouault
wrote:
> Hi Michael,
>
> Fixed per https://github.com/OSGeo/gdal/pull/7659
>
> Even
> Le 24/04/2023 à 05:09, Michael Sumner a écrit :
>
>
> When invoking info on this vrt connection
&
southeast-2.amazonaws.com/images/NSW041500.jp2
gdalinfo /vsicurl/
https://herbariumnsw-pds.s3-ap-southeast-2.amazonaws.com/images/NSW041500.jp2
## GDAL version and provenance
Ubuntu 20.04
GDAL 3.7.0dev-717dcc0eed, released 2023/04/22 (debug build)
--
Michael Sumner
Software and Database Engi
build of code depending on the GDAL C++ API, as the
>> GDAL ABI generally changes between feature versions.
>>
>> Even
>> Le 01/04/2023 à 07:24, Michael Sumner a écrit :
>>
>> Hi,
>>
>> I'm here to ask for help finding where this issue is. If I use
ions.
>
> Even
> Le 01/04/2023 à 07:24, Michael Sumner a écrit :
>
> Hi,
>
> I'm here to ask for help finding where this issue is. If I use an earlier
> build it still works (these are R packages sf and terra that are well
> tested against previous rel
Actually... this seems to be a problem on just one machine, I would have
thought they are the same but fwiw this seems not to be a GDAL or R package
problem.
I'll check again carefully so probably ignore this, thanks!
Cheers, Mike
On Sat, 1 Apr 2023, 16:24 Michael Sumner, wrote:
> Hi,
/src/read_ogr.cpp#L694
In sf it is
https://github.com/r-spatial/sf/blob/main/src/gdal_read.cpp#L582
Thank you.
Cheers, Mike
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
I saw this : https://github.com/mapbox/COGDumper
On Sat, 1 Apr 2023, 02:54 Travis Kirstine, wrote:
> Sorry to post this unrelated question to the GDAL group but I was hoping
> someone could provide some recommendations. Does anybody know of a way to
> efficiently publish a WTMS / TMS from a COG
fix, you're on your own, and you will likely
> have to set the PYTHONPATH environment variable to point to
> /usr/local/lib/python3/dist-packages
>
> GDAL_PYTHON_INSTALL_PREFIX will not help you here (and I'm not sure to
> remember what the actual use case for it was)
>
&g
?
Is it GDAL_PYTHON_INSTALL_PREFIX ?
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
oh dear I missed a key element in the actual command:
gdalwarp band1.vrt tiff_test.tif -ts 4000 4000 -te -116.000
660.000 4.000 780.000 -t_srs "EPSG:3857"
Best
On Thu, Feb 23, 2023 at 6:35 PM Michael Sumner wrote:
> in short, take the grid specification you wan
=Area
> Image Structure Metadata:
> COMPRESSION=LZW
> INTERLEAVE=PIXEL
> Corner Coordinates:
> Upper Left (-116.000, 780.000) ( 10d25'13.65"W, 57d11'40.60"N)
> Lower Left (-116.000, 660.000) ( 10d25'13.65"W, 50d52'46.07"N)
> Upper Right ( 4.000, 780.000)
d it with something like: docker build . --target builder -t
> whatever_tag_name_you_want
>
> Even
> Le 19/02/2023 à 19:08, Michael Sumner a écrit :
>
> I would like to build GDAL in an image on the basis of a dockerfile used
> for CI, I'm confused by the layer process that copies
em running with docker - I can unpick the builder/COPY process,
but maybe I'm missing something easy.
Thank you
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Division
Hobart, Australia
e-mail: mdsum...@gmail.com
___
gdal-dev m
the identifier name following
> the type is for human information only
>
>>
>>
> Den ons. 8. feb. 2023 kl. 00.49 skrev Michael Sumner :
>
>>
>>
>> On Wed, Feb 8, 2023 at 10:44 AM Even Rouault
>> wrote:
>>
>>> Michael,
>>>
>>&g
. but I'm getting used to surprises learning C++ :)
PR incoming, thanks!
Cheers, Mike
> There are likely hundreds of similar occurrences in the code base. I
> believe there's a compiler warning (or at least cppcheck) to detect this,
> but it is not activated.
>
> Ev
264
https://github.com/OSGeo/gdal/blob/579693c2975cda414af75239e09b8b7b952029ca/alg/gdaltransformer.cpp#L412
The discrepancy is visible in the doc here:
https://gdal.org/api/gdal_alg.html#_CPPv424GDALSuggestedWarpOutput212GDALDatasetH19GDALTransformerFuncPvPdPiPiPdi
Cheers, Mike
--
Michael Sumner
So
e tool (or a python
> example) I could use to achieve this ?
>
> Thanks
> Karsten
>
> Karsten Vennemann
> Principal
>
> Terra GIS LTD
> 2119 Boyer Ave E
> Seattle, WA 98112
> USA
> www.terragis.net
> _______
> gdal-dev m
beautiful!
This has been low level on my mind to ask about. Thanks!
Cheers, Mike
On Fri, 6 Jan 2023, 01:42 Even Rouault, wrote:
> Hi,
>
> I've submitted RFC90: Direct access to compressed raster data for review
> and comments:
>
> https://github.com/OSGeo/gdal/pull/7020
>
> Summary: "
ugh, sorry just realized my comment (sentence begins 'Further ...') about
"scaling" is completely irrelevant - you mean spatial and not data scaling.
Just ignore that one sentence.
Cheers, Mike
On Fri, Jan 6, 2023 at 10:36 AM Michael Sumner wrote:
> in short, warp wi
in short, warp will give you exactly the target extent of a (possibly) new
grid specification in that -te subdivided by the -ts.
translate will give you approximately the extent provided to projwin, but
snapped to the source grid alignment at a particular zoom implied by the
projwin subdivided by
tput grid. This would allow the grid to be assigned an arbitrary
> orientation within the CRS frame.
>
> Thanks,
> Brendan
> --
> *From:* Michael Sumner
> *Sent:* Friday, November 4, 2022 6:57 PM
> *To:* Brendan Heberlein
> *Cc:* gdal-dev
&
on for
> the output raster? None of -te, -ts or -t_srs address this.
>
>
> ------
> *From:* Michael Sumner
> *Sent:* Friday, November 4, 2022 6:41 PM
> *To:* Brendan Heberlein
> *Cc:* gdal-dev
> *Subject:* Re: [gdal-dev] Specify Grid Rotation for
that's exactly what the warper does with geolocation arrays, set the target
extent, dimension, and crs with -te, -ts, -t_srs with gdalwarp.
python will have those analogous controls for the warper.
Cheers, Mike
On Sat, 5 Nov 2022, 01:23 Brendan Heberlein via gdal-dev, <
gdal-dev@lists.osgeo.or
1 - 100 of 180 matches
Mail list logo
