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 questio
and do not check email on Fri
On 25 November 2014 at 07:44, Simen Langseth 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 wrote:
Dear Marius Jigmond:
Thanks for your attempt.
Unfortunately, it is produc
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 pr
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 stre
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
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: manoj
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 you
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 "C:\Python2
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
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
NotImplementedErr
On 07/27/2011 04:16 PM, Stephen Mather wrote:
Hi All,
I'm using gdal_translate (version 1.6.3 on Ubuntu) with a scale
factor to stretch float values to 16-bit integer values. It works
fine with tif output, but with PNG it stretches the values to 8-bit:
gdal_translate -of PNG -scale 0 15
Thanks Frank, excellent tip to be aware of. It's working great now.
-marius
On Wed, 2011-07-13 at 22:19 -0700, Frank Warmerdam wrote:
> On Wed, Jul 13, 2011 at 4:23 PM, Marius Jigmond
> wrote:
> > aqfeat = aqLayer.GetNextFeature()
> > aqgeom = aqfeat.GetGeometryRef()
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 g
feat.SetField('IBOUND', 1)
gridLayer.SetFeature(feat)
-marius
On Mon, 2011-07-11 at 22:01 -0700, Frank Warmerdam wrote:
> Marius,
>
> I think a bug is called for.
>
> On Mon, Jul 11, 2011 at 8:00 PM, Marius Jigmond
> wrote:
> > After some more investigation tha
...@hotmail.com
CC: gdal-dev@lists.osgeo.org
On Tue, Jul 12, 2011 at 11:56 AM, Marius Jigmond
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(&#
Date: Tue, 12 Jul 2011 11:36:42 -0400
Subject: Re: [gdal-dev] OGR Geometry methods
From: luke.peter...@gmail.com
To: mariusjigm...@hotmail.com
does executing
sql = 'select FID from judeteEPSG3844 where DENJUD = "CLUJ"'
return anything?-
Luke Peterson
312-533-1051(m)
ile a bug?
-marius
On Mon, 2011-07-11 at 20: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()
>
= 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,
&
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
27;t affected.
>
> Any ideas how to get rid of the black and have transparent pixels
> instead, anyone?
>
> Thanks again,
>
> Michael Corey
>
>
> On 7/8/11 6:12 PM, Marius Jigmond wrote:
>
> > 1.8 compiled from source
> >
> > -marius
&
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/
On Fri, 2011-07-08 at 08:14 -0700, Michael Corey wrote:
> OK, I did a little more work on this, and I've narrowed down what's
> going on, but I could still use some help in figuring out how to solve it.
>
> Here's my original image:
>
> http://mikejcorey.com/spatial/diablo-orig-5pct.tif
>
>
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 s
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
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 results['north']
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
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 wrote
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 Imagin
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 wa
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
,25],
PARAMETER["scale_factor",0.99975],
PARAMETER["false_easting",50],
PARAMETER["false_northing",50],
UNIT["Meter",1]]
INDICATIV: String (25.0)
OGRFeature(caroiaje_100k_EPSG3844):0
POLYGON : 65 points
Is the disconnect simpl
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 :
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
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 Stereog
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
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
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 do
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 wrot
,
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
wrote:
Hi Everyone,
I am using the -crop_to_cutline option of gdalwarp but the resulting
raster is
Hi Everyone,
I am using the -crop_to_cutline option of gdalwarp but the resulting
raster is slightly shifted. Both cutline shapefile and raster are using
the same projection. On the sample raster dataset (pixel=28.5m) the
shift is about 5m to the east and 3.5m to the north. Pixel values remain
the
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 pur
-29 at 22:25 +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 -
ng 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
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
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 tha
47 matches
Mail list logo
