The following works for the proper orientation but there's a shift. I would
consult the documentation as they usually specify conversion params.
gdal_translate -of GTiff -a_srs EPSG:4326 -a_ullr -179.9 -72.7 180 72.7 -co
COMPRESS=LZW NETCDF:GLOBCOMPLIR_nc.2023080819:data ofn.tif
We've been using Geotrellis (https://geotrellis.io/) for Spark jobs.
-marius
From: gdal-dev on behalf of Simon
Sent: Thursday, January 23, 2020 6:30 PM
To: gdal-dev@lists.osgeo.org
Subject: [gdal-dev] Speeding up gdalwarp process
Hi gdal-devs,
I have a
work Mon-Thurs and do not check email on Fri
On 25 November 2014 at 07:44, Simen Langseth simlan...@gmail.com wrote:
Dear Respected Authors:
I would be grateful if you could share any hints on it.
On Tue, Nov 25, 2014 at 4:33 PM, Simen Langseth simlan...@gmail.com wrote:
Dear Marius Jigmond
Does the following work:
gdal_calc.py -A infile --allBands=A --outfile=outfile --calc=median(A)
Date: Mon, 24 Nov 2014 19:14:06 +0900
From: simlan...@gmail.com
To: gdal-dev@lists.osgeo.org
Subject: [gdal-dev] gdal_calc.py: produce median raster
I have one GeoTiff file with 3 bands. I want to
I noticed you're using ArcMap for visualizing the results. Did you check what
histogram/gamma stretch gets applied by default when you load the layer? I run
into this issue often as results from GDAL operations appear much worse when
loaded into ArcMap compared to originals. When I turn off
Frank,
You'll need a newer PostGIS (2.x+) if you want to compile GDAL 1.10 with PG
support.
http://trac.osgeo.org/gdal/wiki/frmts_wtkraster.html
-marius
Date: Tue, 3 Sep 2013 17:02:44 +0200
From: sc...@sarvision.nl
To: gdal-dev@lists.osgeo.org
Subject: Re: [gdal-dev] GDAL 1.10 build fails on
For Texas reservoirs you may contact Jason Kemp of the Texas Water Development
Board (see http://www.twdb.texas.gov/surfacewater/surveys/index.asp for contact
info). He may have some data or may be able to point you in the right direction.
-marius
Date: Mon, 6 May 2013 20:58:27 -0700
From:
At the time I created the ticket I had created a workaround that uses
gdal_translate to clip a rectangular window and gdal_rasterize to burn an
inverse image of the vector. Interested folks can email me for the script (does
the list take attachments?). The script makes some assumptions about
On 05/04/2012 08:16 PM, Jay L. wrote:
List,
Just upgraded to 1.9.0 on my windows 7 machine and I can no longer
access gdal via the python bindings. Command line (ogr2ogr, gdalinfo,
etc. are functioning as expected). Specifically, running 'from osgeo
import gdal', returns:
File
On 09/01/2011 08:23 PM, Ole Nielsen wrote:
Many thanks for the reply, but I need the Python bindings as we are
using ogr deeply embedded in an application and creating of new layers
using the SetField command (amongst others) is critical. Hence my
worry when we get the error
From personal experience (and others') you get some odd error messages
when you start mixing different versions of libgdal. I would expect the
error to go away if you stick to one version of libgdal. Ubuntu-GIS
should have all the packages you need without having to mix GDAL versions.
-marius
David,
Revision 22725 builds just fine for me. I normally do a make clean
before every build. Maybe something went wrong in your build system.
-marius
On Fri, 2011-07-15 at 15:15 -0400, David Burken wrote:
Hi Frank, all,
Just fyi, getting build error from svn. Looks like cpl_serv.h is
...@hotmail.com
CC: gdal-dev@lists.osgeo.org
On Tue, Jul 12, 2011 at 11:56 AM, Marius Jigmond mariusjigm...@hotmail.com
wrote:
Strangely enough, ExecuteSQL works (it returns an osgeo.ogr.layer object with
one feature) with
'select FID from judeteEPSG3844 where DENJUD = CLUJ'
but later a GetField
Hi everyone,
I am trying to test whether centroids of polygons lie/intersect within
another polygon. I have tried Intersect, Within, and Contains but they
always return false. Should these methods work for my intended purpose
or do I need to implement a point in polygon function? Thanks.
-marius
= aqfeat.GetGeometryRef()
if point.Intersect(aqgeom):
xsect = True
print 'ibound = 1'
break
if xsect:
feat.SetField('IBOUND', 1)
gridLayer.SetFeature(feat)
xsect = False
-marius
On Mon, 2011-07-11 at 20:34 -0500, Marius Jigmond wrote:
Hi everyone,
I am trying to test
:46 -0500, Marius Jigmond wrote:
I suppose a piece of code speaks louder :):
xsect = False
for i in range(gridLayer.GetFeatureCount()):
feat = gridLayer.GetFeature(i)
geom = feat.GetGeometryRef()
point = geom.Centroid()
for j in range(aqLayer.GetFeatureCount()):
aqfeat
1.8 compiled from source
-marius
On Fri, 2011-07-08 at 18:02 -0700, Michael Corey wrote:
That certainly could be. I'm running GDAL 1.8 from the kyngchaos site.
Which version are you running?
Thanks,
Michael Corey
On 7/8/11 5:56 PM, Marius Jigmond wrote:
On Fri, 2011-07
?
Thanks again,
Michael Corey
On 7/8/11 6:12 PM, Marius Jigmond wrote:
1.8 compiled from source
-marius
On Fri, 2011-07-08 at 18:02 -0700, Michael Corey wrote:
That certainly could be. I'm running GDAL 1.8 from the kyngchaos
site. Which version are you running
In QGIS you can use Raster Calculator to multiply the rasters.
-marius
On Thu, 2011-04-28 at 00:30 -0700, canduc17 wrote:
Hi everyone.
Is there a way to crop a GEOTiff image using a Mask of the same dimension?
For example, a mask with a filled circle: every pixel inside the circle is 1
They are links between tables. Usually, between the attribute table of a
feature class and a simple (non-spatial, as Mike said) table. But you
can relate any kind of attribute tables.
-marius
On Thu, 2011-04-28 at 14:30 -0700, Paul Ramsey wrote:
Are they tables or links between tables? I can
Luca,
You can also loop through the features with relatively simple Python
code to get the extents. See sample below (adapt at will):
def findPoints(geometry, results):
for i in range(geometry.GetPointCount()):
x,y,z = geometry.GetPoint(i)
if results['north'] == None or
Armin,
Not sure if this would be straightforward in GDAL but you might want to
consider a combination of GDAL + Shapely
(http://gispython.org/shapely/docs/1.2/manual.html#cascading-unions).
What you're trying to do is spatial analysis not directly implemented in
GDAL. Shapely requires that you're
On Wed, 2011-03-09 at 12:35 -0500, Frank Warmerdam wrote:
On 11-03-09 11:51 AM, Wendell Turner wrote:
Do any of the gdal utilities crop an image on non-square
boundaries?
I use gdal_translate with -projwin ulx uly lrx lry to crop
an image that is square (left and right edges are
On Mon, 2011-02-21 at 13:31 +0530, Chaitanya kumar CH wrote:
Zoltan,
The creation of 0394w.tif failed. Does the process have write
permission on directory?
I note that you did not specify the -of option for nearblack. In it's
absence, nearblack creates the output file in Erdas Imagine
With gdal_translate you can do per band re-scaling. I've used it on
Landsat, it's not magic but it works pretty well. You can group your
images, find the right formula for each group and you can quickly verify
the results with a VRT.
-marius
On Mon, 2011-02-21 at 10:39 +0200, Zoltan Szecsei
Hi Zoltan,
Are you sure that the white you see is a pure 255/255/255? Albeit, the
nearblack with its near default of 15 would have likely taken care of
that.
If you load two othophotos in QGIS and specify the transparent pixels,
does the collar between the orthophotos disappear?
Regarding the
://spatialreference.org/ref/epsg/3844/html/ .
-marius
On Sat, 2011-02-05 at 00:06 +0100, Even Rouault wrote:
Le vendredi 04 février 2011 20:25:43, Marius Jigmond a écrit :
Hi Everyone,
Some of the options (I didn't test all) introduced with version 1.8.0
seem to behave strangely. Here's what I'm doing:
1
Should it then match http://spatialreference.org/ref/epsg/3844/esriwkt/
because it does not seem to? I apologize if my questions seem silly.
-marius
On Sun, 2011-02-20 at 00:38 +0100, Even Rouault wrote:
Le dimanche 20 février 2011 00:24:16, Marius Jigmond a écrit :
Is the disconnect
If that doesn't work, going back to what Even said earlier:
1. create a gdal_plugins.conf file in /etc/ld.so.conf.d/ with one
line: /usr/lib/gdal17plugins
2. sudo ldconfig
-marius
On Sun, 2011-02-06 at 19:08 +0100, Even Rouault wrote:
Le dimanche 06 février 2011 19:00:03, Matt Wilkie a écrit :
Even,
Hermann might be right. From reading several sources:
http://gge.unb.ca/Pubs/TR46.pdf
http://www.manifold.net/doc/double_stereographic.htm
it seems that double stereographic is oblique stereographic as long as
the origin is non-equatorial or non-polar.
The actual denomination Double
Did you run sudo ldconfig?
http://trac.osgeo.org/ubuntugis/wiki/UserTutorials
-marius
On Sat, 2011-02-05 at 20:24 +0100, Even Rouault wrote:
Le samedi 05 février 2011 20:18:16, Matt Wilkie a écrit :
So you should try sudo gdal-ecw-build /usr/local because you
./configure the ECW SDK
Aníbal,
Have you tried opening the .vrt file with a text editor? For each
source .tif there should be three entries in the .vrt file (one per
band). What does gdalinfo source.tif say?
-marius
On Fri, 2011-02-04 at 11:41 -0200, Aníbal Pacheco wrote:
On 02/03/2011 04:59 PM, Frank Warmerdam
Hi Everyone,
Some of the options (I didn't test all) introduced with version 1.8.0
seem to behave strangely. Here's what I'm doing:
1. gdal_rasterize -burn 0 -i -l mask mask.shp sample.tif
works fine but I get the following warning:
Warning : the output raster dataset and the input vector layer
Eli,
Both input datasets (mask and raster) exist, which is why 1. works.
Option 1 doesn't use any of version 1.8.0 flags, which is why I thought
something else was introduced behind the scenes with the new flags.
-marius
On Fri, 2011-02-04 at 12:08 -0800, Eli Adam wrote:
Marius,
I don't
On Sat, 2011-02-05 at 00:06 +0100, Even Rouault wrote:
Le vendredi 04 février 2011 20:25:43, Marius Jigmond a écrit :
Hi Everyone,
Some of the options (I didn't test all) introduced with version 1.8.0
seem to behave strangely. Here's what I'm doing:
1. gdal_rasterize -burn 0 -i -l
,
Can you provide some sample data to investigate this? You can raise a ticket at
http://trac.osgeo.org/gdal/newticket and attach the data.
On Tue, Feb 1, 2011 at 8:12 AM, Marius Jigmond mariusjigm...@hotmail.com
wrote:
Hi Everyone,
I am using the -crop_to_cutline option of gdalwarp
images when compared to JPEG. But, because it's a new format, I'll have
to wait until it's implemented in the jai-imageio-ext package.
-marius
Le samedi 29 janvier 2011 22:40:20, Marius Jigmond a écrit :
The link should work for another day or so (automatic server purge). I
can re-upload
Hi Everyone,
I am having some trouble figuring out why adding JPEG compression to a
RGB GeoTIFF results in Nodata pixel triplets becoming data pixel
triplets. Is this expected behavior due to the noisy nature of JPEG
compression or to the fact that these are border Nodata triplets? I know
certain
to the http://www.gdal.org/frmt_gtiff.html
COPY_SRC_OVERVIEWS is applied as long as no general options are passed.
-marius
--
Sent from Ubuntu
On Sat, 2011-01-29 at 21:04 +0100, Even Rouault wrote:
Le samedi 29 janvier 2011 20:28:02, Marius Jigmond a écrit :
Thanks Even. I ran some more
+0100, Even Rouault wrote:
Le samedi 29 janvier 2011 22:17:45, Marius Jigmond a écrit :
Just realized that earlier I hit the wrong reply button (apologies).
The COPY_SRC_OVERVIEWS yielded a wild result. Here's what I did:
1. gdaladdo -ro s70rgb321.tif 2 4 8 16 32 64 128 (note: input
Hi Everyone,
Has anyone experienced the problem in the subject line? I ran:
$gdalwarp -cutline shapefile -crop_to_cutline -s_srs EPSG:31700 -t_srs
EPSG:3844 source.tif dest.tif
and I ended up with a zero filled raster with the correct extent
however. Using 1.8dev. I have a rudimentary script
41 matches
Mail list logo