alys.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
--
Edzer Pebesma (he/him)
Universität Münster, Institute for Geoinformatics
Heisenbergstrasse 2, 48149 Münster, Germany
Phone: +49 251 8333081
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
otest-3.10.1rc1.zip
The NEWS file is here:
https://github.com/OSGeo/gdal/blob/v3.10.1RC1/NEWS.md
Best regards,
Even
--
Edzer Pebesma (he/him)
Universität Münster, Institute for Geoinformatics
Heisenbergstrasse 2, 48149 Münster, Germany
Phone:
Thanks a lot, Even, I missed those. Many regards,
On 14/09/2023 14:30, Even Rouault wrote:
Edzer,
Le 14/09/2023 à 13:52, Edzer Pebesma a écrit :
This came up here: https://github.com/r-spatial/sf/issues/2166
It seems to me that GDAL queries PROJ's ability to access the network
for
e case? If this is the case, is there a way applications can ask
GDAL what this state is?
R package sf interacts with PROJ through GDAL as well as directly. This
sounds messy, but some users want to get NA coordinates back for points
that are not transformable.
--
Edzer Pebesma (he/him)
Universi
On 07/03/2023 17:11, Even Rouault wrote:
Le 07/03/2023 à 17:08, Edzer Pebesma a écrit :
Thanks, Even!
I was after the units of the horizontal coordinates; GetLinearUnits()
returns "unknown" for EPSG:4326 and, according to the docs returns the
units of the vertical coordinate
ly
outside) queries with WKT1 node names. The issue here is that eqearth is
a new projection method that didn't exist in GDAL <= 2.4, so GDAL
fallbacks to WKT2 representation, and it's not the UNIT keyword in WKT2
but LENGTHUNIT.
Using GetLinearUnits() is more robust, as independent
std::cout << proj4s << ":" << (u ? u : "NULL") << std::endl;
return 0;
}
which gives
3031:metre
+proj=eqearth:NULL
whereas the proj4string expansion of "+proj=eqearth" is "+proj=eqearth
+lon_0=0 +x_0=0 +y_0=0 +datum=WGS84
or get the maximum
temperature over a set of given regions? I think there's a big gap now
between the GIS community and the modelling communities; your MDA API
helps closing it, and we can do better.
I'll talk more about this next week Thu at our FOSS4G presentation, and
hope to meet y
ps://cfconventions.org/Data/cf-conventions/cf-conventions-1.9/cf-conventions.html#geometries
--
Edzer Pebesma
Institute for Geoinformatics
Heisenbergstrasse 2, 48151 Muenster, Germany
Phone: +49 251 8333081
First example (written by OGR NetCDF driver):
netcdf t2 {
dimensions:
t2_no
RC2
Even
Le 02/11/2021 à 23:02, Edzer Pebesma a écrit :
This may have nothing to do with this RC but more with PROJ, but I
found it while doing checks with R spatial packages against this RC.
With a CRS shown below, calling OGRSpatialReference->GetName()
resulted with PROJ 8.0.1 in
repared a GDAL/OGR 3.4.0 release candidate.
Pick up an archive among the following ones (by ascending size):
https://download.osgeo.org/gdal/3.4.0/gdal-3.4.0rc1.tar.xz
https://download.osgeo.org/gdal/3.4.0/gdal-3.4.0rc1.tar.gz
https://download.osgeo.org/gdal/3.4.0/gdal340rc1.zip
A s
2020/04/28
>>but rgdal build and GDAL runtime not in sync:
>>... consider re-installing rgdal!!
>> """
>> I'd assume everything to match perfectly in a Debian build env.
>
> Nothing has been (re)built with GDAL 3.1.0 yet.
>
> When packages
ormation *ct =
> OGRCreateCoordinateTransformation(aSRS, bSRS);
>if (ct == NULL) {
> printf("ct NULL\n");
> exit(1);
>}
>exit(0);
> }
>
> So the error is linked to having importFromProj4() and then importFromWkt()
> This doesn'
On 1/23/20 5:29 PM, Sean Gillies wrote:
> Hi,
>
> On Thu, Jan 23, 2020 at 8:01 AM Edzer Pebesma
> mailto:edzer.pebe...@uni-muenster.de>>
> wrote:
>
> The combination of GDAL 3.0.3 and PROJ 6.3.0 arriving on debian
> platforms is causing some havoc for so
ot;Popular
Visualisation Pseudo Mercator",ID["EPSG",1024]],PARAMETER["Latitude of
natural
origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude
of natural
origin",0,ANGLEUNIT["degree",0.0174532925
(should start soon hopefully), so API
> details will likely change,
> but early feedback is welcome.
>
> Even
>
--
Edzer Pebesma
Institute for Geoinformatics
Heisenbergstrasse 2, 48151 Muenster, Germany
Phone: +49 251 8333081
pEpkey.asc
Description: application/pgp-keys
_
ackages sf, lwgeom and rgdal have been updated; github
versions work, we'll be submitting CRAN version in the next few days.
I updated the wiki status pages.
--
Edzer Pebesma
Institute for Geoinformatics
Heisenbergstrasse 2, 48151 Muenster, Germany
Phone: +49 251 8333081
pEpkey.asc
Descriptio
he same name that this dimension.
>
> Can you create a ticket about that ?
Sure: https://trac.osgeo.org/gdal/ticket/7165
>
>
>
> Even
>
>
>
> --
>
> Spatialys - Geospatial professional services
>
> http://www.sp
as a sigma(kz)
variable with pressure values which vary between 0 and 1.
Is this known behaviour?
--
Edzer Pebesma
Institute for Geoinformatics
Heisenbergstrasse 2, 48151 Muenster, Germany
Phone: +49 251 8333081
___
gdal-dev mailing list
gdal-dev@lists
On 11/25/2017 05:49 PM, Edzer Pebesma wrote:
> On the system I work on (Ubuntu 17.04 with ubuntugis-unstable GDAL),
> haveGEOS() seems unreliable (giving false, where it is true).
I'm sorry, that was a wrong diagnosis of a self-caused problem.
--
Edzer Pebesma
Institute for Geo
> Even
>
>
>
> --
>
> Spatialys - Geospatial professional services
>
> http://www.spatialys.com
>
>
>
> _______
> gdal-dev mailing list
> gdal-dev@lists.osgeo.org
> https://lists.os
cept is quite old, and not restricted to
spatially gridded data. But I'm happy to restrict the discussion here
(first) to spatially gridded data.
--
Edzer Pebesma
Institute for Geoinformatics
Heisenbergstrasse 2, 48151 Muenster, Germany
Phone: +49 251 8333081
g" should never come out of a spatial filter.
I never understood the value of having all these typed EMPTYs in SFA,
and whether/how they'd differ from NULL geometries in OGR.
>
> I would say such things should disappear, i.e., the PostGIS behavior is
> logical.
>
> Ari
>
On 10/26/2017 09:39 AM, Ari Jolma wrote:
> Edzer Pebesma kirjoitti 26.10.2017 klo 10:26:
>> More in general, we have multidimensional arrays where horizontal space
>> is not an x y image (GDAL raster band) but for instance a sequence of
>> feature geometries (e.g., in
pszAxisUnit[nAxisCount] or
>> GDALGridAxis::pszUnit in my scetch
>>
>>
>>
>> > and arbitrary value (double, date, time, ...)?
>>
>>
>>
>> In my scetch, I had indeed some hesitation about the appropriate type
>> for adfAxisCoordinates
opriate type
> for adfAxisCoordinates array. I think that for the cases mentionned
> above (4D spatio-temporal) we can always represent the coordinates along
> each axis as a double value. The case where that's the less obvious is
> for date/time, but using some reference like the Unix Epoch
just
>> concatenate any amount of GRIB messages into a single file and this
>> will be seen as a multi-band raster by GDAL.
>>
>>
>>
>> Even
>>
>>
>>
>> --
>>
>> Spatialys - Geospatial professional services
>>
>
s) or require the user to set
> the values of the n > 2 dimensions. This was already the case in the current
> driver with the time dimension.
The netcdf files I looked at through the GDAL driver had time and/or
depth offered as different gdal layers, and offered different variables
(i.e. dif
s://www.unidata.ucar.edu/software/udunits/
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal of Statistical Software: http://www.jstatsoft.org/
Computers & Geosciences: http://elsevier.com/locat
29 matches
Mail list logo
