Dear list,
I have a question about mosaicing multidimensional datasets (netCDF,
Zarr etc). These are typically 2 or 3D grids with a time dimension
from weather or ocean models (3 or 4D in total).
Using a multidimensional VRT:
https://gdal.org/drivers/raster/vrt_multidimensional.html
it is poss
a
larger region than needed, or perhaps doing a more advanced query than the
simple bounding box one, I'm doing now. Should really look a bit more into
the API!
Cheers,
Simon
On Mon, Jan 2, 2017 at 7:12 PM, Andre Joost
wrote:
> Am 01.01.2017 um 22:52 schrieb Simon Lyngby Kokkendorff:
>
ry for the noise!
Simon
On Sun, Jan 1, 2017 at 9:15 PM, Simon Lyngby Kokkendorff
wrote:
> Thanks Even,
>
> I'll try to see if I can fix the issue for that particular location by
> including a larger area. Was just puzzled because things were working fine
> at other locations
t;relevant" features
intersecting the bounding box.
Will test a bit more 😄
Cheers,
Simon
Den 01/01/2017 7.41 PM skrev "Even Rouault" :
> On dimanche 1 janvier 2017 18:17:58 CET Simon Lyngby Kokkendorff wrote:
>
> > Dear List,
>
> >
>
> > I hav
Dear List,
I have some issues getting river boundaries out of OSM data when using
GDAL to read e.g. an OSM xml file downloaded via the overpass API. Using
e.g.:
>ogr2ogr -t_srs EPSG:3857 -f SQLITE -dsco SPATIALITE=YES --config
OSM_USE_CUSTOM_INDEXING NO map.sqlite map.osm
>ogrinfo --version
GD
Done: https://trac.osgeo.org/gdal/ticket/6225#ticket
Thanks once again!!
Simon
On Thu, Nov 19, 2015 at 10:59 AM, Even Rouault
wrote:
> Le jeudi 19 novembre 2015 10:47:50, Simon Lyngby Kokkendorff a écrit :
> > I can confirm also that it works when I set a target srs.
> > But
I can confirm also that it works when I set a target srs.
But must admit that the current behaviour was confusing. I knew how to
interpret the input data - under the assumption that a geotransform is not
changed when I do no transformations or interpolations:
If I know that my data are "points" rat
o do that. Anyway it happens when the
> destination geotiff is created and its fields are set. The
> geotransformation of the geotiff file are set indirectly using geotiepoints
> using the geotransformation of the source (which is correct).
>
> Ari
>
>
> On 19.11.2015 10:24, S
rm as my input center.asc, but
shifted half a cellsize from what I would expect...
In [6]: ds = gdal.Open("center_back.asc")
In [7]: ds.GetGeoTransform()
Out[7]: (0.0, 1.0, 0.0, 1.0, 0.0, -1.0)
GDAL version:
GDAL 1.11.2, released 2015/02/10
Cheers,
Simon
On Wed, Nov 18, 2015 at
Hi List,
I have a geoid in AAIGRID-format, with the following header:
ncols 1001
nrows 731
xllcenter 40.0
yllcenter 6045000.0
cellsize 500.0
NODATA_value
GDAL will read this just fine, and report:
Driver: AAIGrid/Arc/Info ASCII Grid
Files: dkgeoid13b.utm32
Size is 1001, 731
Coordinat
I can also confirm that with -co FORMAT=NC4 everything seems to work as
expected, also with other netCDF readers...
Thanks again,
Simon
On Tue, Oct 27, 2015 at 9:03 PM, Simon Lyngby Kokkendorff
wrote:
> Thanks for the tips Kyle and Even and for the quick fix!
> Not a big netCDF user
Le mardi 27 octobre 2015 16:54:56, Kyle Shannon a écrit :
> > Simon,
> >
> > On Mon, Oct 26, 2015 at 5:04 PM, Simon Lyngby Kokkendorff
> >
> > wrote:
> > > Hello list,
> > >
> > >Just observed a puzzling behaviour of the netCDF driver. W
Hello list,
Just observed a puzzling behaviour of the netCDF driver. When creating a
netCDF dataset of Byte type, the netCDF driver seems to interpret the
nodata value as a signed byte.
For example a nodata value assigned with -a_nodata 200, will become -56.
Of course one can work around this
Hi Even,
Thanks a lot for the explanation - will modify as needed! Yes, the
behviour was a little surprising, as it is not present with e.g.
'ATTRIBUTE=' .
Cheers,
Simon
On Wed, Aug 12, 2015 at 11:25 AM, Even Rouault
wrote:
> Le mardi 11 août 2015 11:36:32, Simon Lyngby Kokken
Hello list,
I'm trying to rasterize some 3d-lines programtically from python. Using
gdal 1.11.2 from the osgeo4w64 environment. Seems to work fine, except that
I get a z-offset of 255 (one byte) in my output raster. Works fine if I use
gdal_rasterize with the -3d option.
The following python c
Yes, it was after a call to:
OGR_G_SetPoint
that the geometry type was altered. Should really test if it's a 2.5D or 2D
geometry and use OGR_G_SetPoint_2D in the latter case!
Thanks,
SImon
On Mon, May 6, 2013 at 7:29 PM, Simon Lyngby Kokkendorff
wrote:
> Hi Even,
>
> I think you
Hi Even,
I think you're right. The transformation seems to alter the geometry type
somehow - need to debug a little more to see exactly what's going on.
Thanks,
Simon
On Mon, May 6, 2013 at 3:50 PM, Even Rouault
wrote:
> Selon Simon Lyngby Kokkendorff :
>
> > Hi
Hi List,
I am programming an application which can transform ogr-datasources
(using our own transformation library). In that process I use ogr (version
1.10) via the c-api to read and write datasources - usually it works like
a charm. However, when I try to write an output spatialite datasourc
Hi Even,
Thanks for the pointer - might try again with the -multi option. Anyway,
the -co TILED=YES is speeding things up drastically (something like a
factor 10).
Cheers,
Simon
On Mon, Feb 25, 2013 at 1:16 PM, Even Rouault
wrote:
> Le lundi 25 février 2013 11:13:47, Simon Lyn
Hi Even,
> When you say "crashed", you mean it exited with the integer overflow error
> (to
> be opposed as the windows that is displayed by Windows when a process
> really
> crashes) ?
>
Yes, it didn't really crash, but exited with the integer overflow error
message.
The issue is rather with
Hej List,
I was using gdal1.10_dev (precompiled from gisinternals.com/sdk) to warp
a quite large ECW image, but after running for a few days the process
crashed with an integer overflow error:
C:\Users\simlk\Downloads\release-1600-x64-gdal-mapserver>gdalwarp -s_srs
"+proj=
etmerc +ellps=GRS80 +
Hi Even,
Looks great, thanks.
Cheers,
Simon
On Wed, Jan 23, 2013 at 10:18 PM, Even Rouault wrote:
> Le mercredi 23 janvier 2013 20:59:09, Simon Lyngby Kokkendorff a écrit :
> > Hello List,
> >
> >I'm using OGR to read a FileGDB database, perform some tr
Hello List,
I'm using OGR to read a FileGDB database, perform some transformations
and then write an output FileGDB database. In that process I would like to
preserve the structure of the original FileGDB, i.e. put "Feature Classes"
into the right "Feature Datasets" etc.
Hence, my question i
Hi Luis,
I have encountered a similar issue on my machine (win xp sp3 too). I
think what I saw was due to PyQt4 using a version of libmysql.dll, which
gdal is also linked against (dll included in the binaries from Tamas), and
it is crucial which version is found first by windows. Notice that wi
Hi,
It shouldn't be necessary to have gdal and python compiled with the same
msvc-version, as long as the various needed c-runtime dll's can be found in
the PATH. In my experience problems like you describe (DL'HELL) can be
related to the order in which things appear in the PATH. For example, if
Hi List,
I have been trying to transform some FileGDB databases with OGR and have
run into a problem reading a specific field called "GlobalID". The field
simply cannot be read by OGR, which emits this warning (and the field is
left as NULL):
ERROR 1: Error: Failed to determine string value for
Hi Even,
Thanks a lot :-) I'll try the workaround and keep an eye open for when/if
the patch goes into trunk.
Cheers,
Simon
On Mon, Oct 8, 2012 at 10:07 PM, Even Rouault
wrote:
> > However, if I try to run ogr2ogr like:
> >
> > ogr2ogr -s_srs "EPSG:25832" -t_srs "+proj=etmerc +lat_0=0 +
Hi List,
For some transformations of data in Greenland we need the extended
precision of the etmerc (Extended Transverse Mercator) projection, which is
built into proj4.8. I have succesfully built the proj4.8 library and have
gotten a definition like, e.g.: +proj=etmerc +lat_0=0 +lon_0=9 +k=0.99
ot have the ones I’m
> using. Does it matter, or must proj.dll contain all the datums I intend to
> use? Is the data pointed to by GDAL_DATA sufficient to translate from any
> known datum to another?
>
> ** **
>
> *From:* Simon Lyngby Kokkendorff silyko-at-gmail.com |GDAL| [mai
Hi Marc,
Just to understand your problem fully: You have data in your database in
WGS84, in geographical coordinates, or what? That would correspond to a
srs=osr.SpatialReference(), srs.SetWellKnownGeogCS("WGS84").
What kind of projection or geographical coordinate system do you want to
pro
Hi Jerl,
If you have python with matplotlib and the gdal bindings installed, you
can use gdal to read the DEM (which I assume you are getting gdal_contour to
draw contours from), and then matplotlibs contour-function to draw the
labelled contours. You can draw on top of an imported map/image-fi
o layer extent" on a modified layer you
might end up with a view, where your data is a tiny (or invisible) blob in
one of the corners.
Best regards,
Simon
On Fri, Apr 8, 2011 at 11:40 AM, Even Rouault
wrote:
> Selon Simon Lyngby Kokkendorff :
>
> Yes indeed, you've precisely des
Hi List,
I have observed a puzzling behaviour of the reported extent of an ogr
layer after geometries of features have been modified, with the file opened
in update mode. I use gdal 1.7 via the pyhton bindings, open a dataset with
ogr.Open(fname,update=True), then loop through layers and feature
Thanks a lot. Didn't know that GetGeometryCount did that on polygons :-)
Cheers,
Simon
On Tue, Mar 8, 2011 at 11:51 AM, Even Rouault
wrote:
> Selon Mike Toews :
>
> > On 8 March 2011 23:22, Even Rouault
> wrote:
> > > Not exactly. In fact you have to use the Geometry.GetGeometryCount()
> that
Hi List,
I am using ogr via the python bindings to construct various polygons.
Here's just a simple question, to which someone might have some input. Is
there anyway to determine the topology, i.e. the number of holes, of a
polygon without e.g. having to export to WKT and examining the output
st
Hi Yosuke,
If you use MultiPolygon geometries (with your polygons as the parts), geos
(via ogr or shapely) should be able to do the job in one go, which should
speed things up considerably...
Cheers,
Simon
On Sat, Nov 6, 2010 at 5:04 AM, yosuke kimura wrote:
> All,
>
> I think I found
Hi Rudi,
py2exe includes a lot of things into the distribution by default, so
you'll need to extend the "excludes" in your options in the setup script to
get a smaller distribution, see
http://www.py2exe.org/index.cgi/GeneralTipsAndTricks. I've been using gdal
with py2exe, and I guess that if yo
Hi all,
I know that Gitte installed the python bindings for gdal from
http://pypi.python.org/pypi/GDAL/ and followed the link to the minimalistic
gdal binaries at
http://download.osgeo.org/gdal/win32/1.6/gdalwin32exe160.zip. I did the same
thing, and the python bindings run fine on my machine - h
ecause I am still PhD candidate at Czech Technical University in
> Prague I can participate as a student again...
> I proposed the project (no.5.) related to GDAL2Tiles/MapTiler at:
> http://trac.osgeo.org/gdal/wiki/SummerOfCode
>
> Best,
>
> Klokan
>
> On Wed, Mar 25,
Hello All,
This may seem stupid :-) I am doing an application using the gdal bindings
for python on windows. I am trying to bundle it all up to an executable
using py2exe and the NSIS-installer.
It all seems to go well, the dll's from my installation of gdal,
C:\gdalwin32-1.6\bin, get copied int
40 matches
Mail list logo
