[gdal-dev] Reading Layers In USGS Geospatial Topo Files:

2025-04-11 Thread Mike Flannigan via gdal-dev
f Do this 2nd: gdal_translate 'Mount Rose2.tif' 'Mount Rose.tif' -b 1 -b 2 -b 3 -mask 4 -co COMPRESS=JPEG -co PHOTOMETRIC=YCBCR --config GDAL_TIFF_INTERNAL_MASK YES Works on Linux, but creates a HUGE file (~88 MG): gdal_translate -of GTiff -co "COMPRESS=DEFLATE"

Re: [gdal-dev] fail to pip install gdal

2024-11-21 Thread Mike Taves via gdal-dev
On Fri, 22 Nov 2024 at 10:26, ni hao via gdal-dev wrote: > when I tried to pip install gdal, I got: > > Collecting gdal > Downloading gdal-3.10.tar.gz (848 kB) This fails because the source distribution of gdal needs to be compiled. There are no pre-compiled binary wheels on PyPI. There are se

Re: [gdal-dev] Help using gdal_grid when src_datasource is a NetCDF File

2024-08-22 Thread Mike Taves via gdal-dev
On Fri, 23 Aug 2024 at 04:03, Ewan Tarrant via gdal-dev wrote: > Thanks in advance for your help, To read point data from a netCDF file, refer to the vector driver: https://gdal.org/drivers/vector/netcdf.html The data would need to conform to CF-1.8 for Point data: https://cfconventions.org/Data

Re: [gdal-dev] Building GDAL documentation with ReadTheDocs

2024-07-24 Thread Mike Taves via gdal-dev
On Wed, 24 Jul 2024 at 22:34, Javier Jimenez Shaw via gdal-dev wrote: > Maybe the solution is that proj.org default be "stable" and not "9.4". The reason why PROJ's RTD is configured this way is that "stable" is set at one commit (currently 875a485f) for the most recent release, as described by h

Re: [gdal-dev] Issues specifying PostGIS connection details - possibly related to special characters in passwords?

2024-07-02 Thread Mike Taves via gdal-dev
On Wed, 3 Jul 2024 at 03:36, Robin Wilson via gdal-dev wrote: > Any suggestions welcome, If you routinely access a specific database, a reliable method is to create a .pgpass file: https://www.postgresql.org/docs/current/libpq-pgpass.html This way you never need to provide a password. __

[gdal-dev] Warp Reprojection A UTM NASD Projection into Missouri State Plane Western or Central leaves Gap Artifacts in the Tile Edges due to Distortion

2022-07-09 Thread Mike Stanton
Reprojecting shows a calculation error in the warping of a UTM Projection to a Missouri State Plane West or Central projection. The manifestation is a series of crosshatch marks where the the tile Corners are and/or a black grid around all the tiles in the data set, which is dependent on the z

Re: [gdal-dev] Convert to min containing bit depth?

2022-04-21 Thread Mike Taves
On Fri, 22 Apr 2022 at 07:05, wrote: > > Idea for a small but useful python tool: scan image for min/max values and > convert to smallest possible bit depth without losing values. Surely someone > has done something like this already. Any suggestions for where to look for > prior art? This is

Re: [gdal-dev] Faster alternative to GetFeatureCount?

2022-01-20 Thread Mike
In the case if FileGDB, which driver you use to read the FileGDB makes a big difference. On Thu, Jan 20, 2022 at 8:59 AM Jon Morris wrote: > Hi Sean, > > > > That seems to be the best option. I was testing on Shapefile initially, > but after I posted, I tried a memory layer and got quite differe

Re: [gdal-dev] How to fix path for gdal .py scripts in Anaconda-Windows?

2022-01-13 Thread Mike Taves
On Thu, 13 Jan 2022 at 07:25, Rafael Lima wrote: > I installed gdal on a Windows environment using anaconda and tried using one > of the .py scripts (e.g., ogrmerge.py). I expected the script to make use of > the gdal installed in such an environment, but it seems that the script finds > anothe

Re: [gdal-dev] Large GeoJSONs and aborting file opening

2021-07-29 Thread Mike Flannigan
n a geoj: https://www.reddit.com/r/QGIS/comments/mbljrw/render_visible_layers_only/ https://gis.stackexchange.com/questions/111784/displaying-only-selected-features-on-map-in-qgis Mike On 7/29/21 11:30 AM, Simon Eves wrote: The problematic one in this case is about 30GB, with ~5.9M features, of prop

Re: [gdal-dev] Large GeoJSONs and aborting file opening

2021-07-29 Thread Mike Flannigan
with. You might try that just to see if it makes a difference for you. So you could replace "],[" with "],\n[" (\n being a newline).  The file might get even bigger, but QGIS is much happier with the newlines than with the wide rows. Mike On 7/29/21 9:26 AM, Newcomb, Doug

Re: [gdal-dev] Large GeoJSONs and aborting file opening

2021-07-29 Thread Mike Flannigan
I have no problems.  In fact I converted from SHP to geojson 2 years ago because I used to have problems with SHP linear files. I use QGIS 3.16.8 on Linux Mint. Mike On 7/28/21 2:36 PM, gdal-dev-requ...@lists.osgeo.org wrote: Date: Wed, 28 Jul 2021 12:22:12 -0700 From: Simon Eves To:gdal-dev

Re: [gdal-dev] Does gdal support Military Grid Reference System (MGRS) grid?

2021-06-14 Thread Mike Taves
On Mon, 14 Jun 2021 at 20:49, Bang Pham Huu wrote: > I've questions regarding Military Grid Reference System (MGRS) grid data. If it's simple coordinate conversions to/from MGRS and UTM, a decent library is GeographicLib: https://geographiclib.sourceforge.io/html/GeoConvert.1.html https://geogra

[gdal-dev] GDAL Reproject Run as a Batch under QGIS

2021-06-03 Thread Mike Stanton
ng output that might be named the same as the input file. What it appears to expect is individual TIFF names. Does anyone have suggestions on how to correct this? Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Removing trailing decimals from contour labels of shape files

2021-05-15 Thread Mike Taves
On Sat, 15 May 2021 at 12:28, Jason Snyder via gdal-dev wrote: > I noticed that when displaying the graph in mapserver I get label values such > as 1010. whereas I want the labels to be whole numbers. How do I > tweak this code to get whole numbers for pressure contour labels in map > s

[gdal-dev] gdal_translate resampling introduces bias, gdalwarp does not

2020-10-28 Thread Mike A
f vs src_TRANSLATE: https://drive.google.com/file/d/1XeWafUlUzN6somcyUL2dfyvlwS3IOpRP/view?usp=sharing gdalwarp src.tif vs src_WARP: https://drive.google.com/file/d/1Zv5ypaRzCEZzVUmfGf01UtTTVoPD7ADL/view?usp=sharing Thanks in advance, Mike ___ gdal-

Re: [gdal-dev] ogr2ogr: debugging slow performance when converting .e00 to other formats

2020-10-12 Thread Mike Flannigan
ay be the best option for you. Mike On 10/12/20 2:00 PM, gdal-dev-requ...@lists.osgeo.org wrote: Hello all, I am new to this list, apologies if this type of question has already been asked. I’m using GDAL/OGR version 3.0.4 on a brand new server with 256GB ram. I have a moderately sized (~110M

Re: [gdal-dev] gdal2tiles PDF --config GDAL_PDF_DPI option

2020-09-19 Thread Mike Flannigan
think I see that it is part of gdal: https://gdal.org/programs/gdal2tiles.html I do: gdal -v Surprisingly I get: Command 'gdal' not found, did you mean: . . . I do: gdal_translate -of GTiff -co "COMPRESS=DEFLATE" --config GDAL_PDF_DPI 300 /home/mike/Documents/copy3/king.pdf

Re: [gdal-dev] Affine transform: Practical uses of x-skew and y-skew

2020-05-03 Thread Mike Taves
On Sat, 2 May 2020 at 08:01, Craig Delancy wrote: > I am trying to explain to a colleague how the affine transform maps image > coordinates onto a projection system. > Specifically, I am trying to explain how mathematically speaking, any > parallelogram section of the coordinate system can be re

Re: [gdal-dev] Windows CI build times

2020-04-16 Thread Mike Taves
Hi Even, Yes, I agree these times are way longer than they should take. It looks like Howard was investigating AZP: https://github.com/OSGeo/PROJ/tree/azp how is that going? It seems that many projects have shifted from AppVeyor to AZP for performance benefits. As you mentioned, some of the dep

[gdal-dev] Adding datasets to VRT at runtime

2020-03-20 Thread Mike
s I'm that loading geotiffs based on a paging system, and I'd like to be able to efficiently handle adding these new datasets. Any suggestions for where to look or maybe a different method is appreciated too. Thanks, Mike ___ gdal-dev mailin

Re: [gdal-dev] Transform produces different results from 2.2.4 and 3.0.2 or 3.0.3

2020-01-29 Thread mike Null
? Anyways, thanks again. Mike On Wed, Jan 29, 2020, 00:30 Carsten Lockenkötter < carsten.lockenkoet...@gon.de> wrote: > I ran into the same problem a while back. > My solution was to change the AxisMapStrategy to OAMS_TRADITION_GIS_ORDER. > > You can read more about this here:

[gdal-dev] Transform produces different results from 2.2.4 and 3.0.2 or 3.0.3

2020-01-28 Thread mike Null
it agreed with the project tool on esri and epsg.io. epsg.io/transform#s_srs=4326&t_srs=32617&x=-81.5&y=30.5 Anyways I desire to use python 3 and gdal 3. What might I need to do to get it working correctly? I know there was a lot of changes

Re: [gdal-dev] FileGDB & FID

2019-11-05 Thread Mike
ayer in python. On Tue, Nov 5, 2019, 12:29 mike Null wrote: > > I'm trying to create some features in a File GDB, and have that driver > loaded. > > I've been doing something along the lines just fine to gpkg, but I hit a > FileGDB issue and wondering how to get ar

[gdal-dev] FileGDB & FID

2019-11-05 Thread mike Null
I'm trying to create some features in a File GDB, and have that driver loaded. I've been doing something along the lines just fine to gpkg, but I hit a FileGDB issue and wondering how to get around it. This line in python lyr = fgdbds.CopyLayer(inmemorylayer,somename,['OVERWRITE=YES']) I get the

Re: [gdal-dev] gdal2tiles missing 'GRA_mode' resample method

2019-09-08 Thread Mike Taves
al2tiles.py too? Hi Tim, This has been addressed with https://github.com/OSGeo/gdal/pull/1817 and should appear in the next release. Cheers, Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Python RasterizeLayer and gdal_rasterize

2019-07-24 Thread Mike Taves
On Thu, 25 Jul 2019 at 01:19, Ivan Lucena wrote: > mem_ds.SetGeoTransform((xmin, xResolut, 0, ymax, 0, yResolut,)) Check your value for yResolut, as it should be negative. In my experience, a raster full of zeros is when the geometry is outside the raster window, so be sure to manually ch

Re: [gdal-dev] ogr2ogr not creating feature ID

2019-07-10 Thread Mike Flannigan
ng", "coordinates": [ [ -93.40203, 35.97839 ],  [ -93.4031, 35.97831 ],  [ -93.40372, 35.97807 ],  [ -93.40444, 35.97744 ],  [ -93.40544, 35.97578 ],  [ -93.40589, 35.97552 ],  [ -93.40636, 35.97549 ],  [ -93.407, 35.97576 ],  [ -93.6255, 37.18735 ] ] } }, So I'm thinking the id need

Re: [gdal-dev] Ellipsoidal length of a line

2019-06-16 Thread Mike Taves
On Sat, 15 Jun 2019 at 03:32, Nicolas Cadieux wrote: > I am trying to get the length of a line in python. (Not just the straight > length between the first and last nodes). Using geopandas, (therefore the > Shapely lib) I am getting the euclidien distance even though the dataframe > holdings t

Re: [gdal-dev] Convert GEOTIFF to JPG using GDAL

2019-04-14 Thread Mike Flannigan
Most browsers will open a JPG file. In the browser, go to File - Open and open the file. Mike On 4/13/2019 2:00 PM, gdal-dev-requ...@lists.osgeo.org wrote: Date: Fri, 12 Apr 2019 23:05:06 +0300 From: עוז אלהרר To:gdal-dev@lists.osgeo.org Subject: [gdal-dev] Convert GEOTIFF to JPG using

Re: [gdal-dev] How to create GeoTiff from image and UTM Coordinates

2019-01-27 Thread Mike Taves
On Sat, 26 Jan 2019 at 05:50, Adrien André wrote: > If GDAL cannot read the .emf file, try using ImageMagick: > convert image.emf my_raster.png This is a good approach, however you will also need to pay attention to convert's -density and -units options [1] to maintain some sort of meaningful sca

[gdal-dev] gdal/CURRENT lags

2018-12-20 Thread Mike Taves
The symlink to https://download.osgeo.org/gdal/CURRENT/ hasn't been updated for a few releases. If there is a "HOWTO-RELEASE" document, this step should be added to update the link. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.o

Re: [gdal-dev] Support for S-100 and S-102

2018-08-30 Thread Mike
Indv.shtml (US1AK90M). That page also tells you the update application date and issue date. I can try to submit an issue later today. On Thu, Aug 30, 2018, 8:54 AM Even Rouault wrote: > On jeudi 30 août 2018 07:29:12 CEST Mike wrote: > > The multiple field/subfields - I think gdal already

Re: [gdal-dev] Support for S-100 and S-102

2018-08-30 Thread Mike
The multiple field/subfields - I think gdal already does this for S-57 on DSID. It just prefixes the subfields with the field name. However -side note: it seems the DSID_ISDT (issue date) is incorrect, it just has the same value as DSID_UADT (update application date). On Thu, Aug 30, 2018, 4:27

[gdal-dev] ogr2ogr -sql on a dataset level

2018-08-19 Thread Mike
I'm looking to see if there is a way to do a sql statement on a dataset level instead of a layer level. Basically, I have a bash script, and wish to merge a bunch of esrifilegeodatabase.zip into a geopackage. I can do something like this ogr2ogr -F GPKG dest.gpkg -append -update $zip -sql "SELEC

[gdal-dev] Segfault when accessing vpf dqyarea

2017-12-20 Thread Mike
any layer to get further info on, except dqyarea@dqy(*)_area, where it then segfaults. I'm unsure if this is a bug or on how it was packaged. Of note, gdal 2.2.0 (from qgis 2.18.9) on Windows reads dqyarea fine. Any tips on how to read the layer?

Re: [gdal-dev] gdalwarp resampling_method option

2017-10-05 Thread Mike Toews
Hi, I requested this feature 2 years ago: https://trac.osgeo.org/gdal/ticket/6015 I probably had some workaround in Python, but I can't recall any details. Cheers, Mike On 6 October 2017 at 00:59, sujankoirala2 wrote: > Hi there, > I was wondering if it would be possible to

Re: [gdal-dev] Improve clipping of shapefiles

2017-09-11 Thread Mike Toews
On 11 September 2017 at 21:09, Paul Meems wrote: > I have a large shapefile with over 2.8 million shapes (fishnet) and I have a > border file with only 1 multipolygon. > > I'm trying to clip the fishnet with the border. > Using code is takes about 5 min. using command line it takes even longer. >

Re: [gdal-dev] GDAL, vsis3 and vsisubfile

2017-07-24 Thread Mike Pfaffenberger
ing remote windowed reads quickly on this data. However, I'm very open to suggestions if anyone has any ideas on how it might work. Cheers. On Mon, Jul 24, 2017 at 11:21 AM, Even Rouault wrote: > Mike, > > > > (note to other readers: this is the continuation of the thread > &g

Re: [gdal-dev] VSIS3 on digital globe multiview-stereo (NITF)

2017-07-24 Thread Mike Pfaffenberger
250 Corner Coordinates: Upper Left (0.0,0.0) Lower Left (0.0,44032.0) Upper Right (43008.0,0.0) Lower Right (43008.0,44032.0) Center (21504.0,22016.0) Band 1 Block=1024x1024 Type=UInt16, ColorInterp=Gray Overviews: 21504x22016, 10752x11008, 5376x5504, 2688x2752, 1344x1376

Re: [gdal-dev] GDAL, vsis3 and vsisubfile

2017-07-24 Thread Mike Pfaffenberger
I turned on some debug options that shed some light on to what's going on. It appears that the NITF driver must internally open a JPEG 2000 Driver on a virtual subfile. In my case, that virtual subfile starts at offset 4038 and continues to the end of the file, offset 901949970. While this is a ni

[gdal-dev] VSIS3 on digital globe multiview-stereo (NITF)

2017-07-23 Thread Mike Pfaffenberger
Hi all, I'm having an issue with VSIS3 on the Digital Globe multiview-stereo imagery collection. Here is an example command I'm running -- gdal_translate -srcwin 000 000 1000 1000 /vsis3/glitch253/test2.ntf from_s3.nitf So I'm simply cropping out an upper left 1000x1000 window. When I've done th

Re: [gdal-dev] Unable to obtain EPSG from OGRSpatialReference

2017-06-07 Thread Mike K
Nik Sands, thank you, you nailed it. That was my problem: the UTM coordinate system I was expecting is projected, so I needed to use those functions you listed with the "PROJCS" key to get the proper EPSG code. Once I parsed out the code, I get 32617 as expected. Thank you all! -- View th

[gdal-dev] S-57 open options with python bindings

2017-06-07 Thread Mike
Hi, I looked to see if this was asked before but couldn't find it. Are any of the S57 open options exposed in the swig bindings? If so, how do I use them? the page I refer to is here http://www.gdal.org/drv_s57.html Thanks -Mike ___ gdal-dev ma

[gdal-dev] Unable to obtain EPSG from OGRSpatialReference

2017-06-06 Thread Mike K
Hello, I have a situation in which I would like to load up a set of polygons, and determine what EPSG code they are using. I am testing with a simple square polygon in upstate New York using the UTM projection, saved into a Shapefile . I am able to create my dataset, layer, and read the geometry

[gdal-dev] issue updating shapefile using python bindings

2017-04-22 Thread Mike
teger64' vs regular int. Any pointers on how to move past this would be good. I won't be able to look at the error till Monday, so I might not be able to answer any follow up questions till then. Thanks! -Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] ogr2ogr -sql and -nln

2017-03-08 Thread Mike
just and not +somesuffix. How can I accomplish what I want? This is ultimately run on many datasets, so it will be appending to existing tables (or create them if they didnt exist yet). -Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] gdal on freebsd

2016-11-07 Thread Mike Toews
On 7 November 2016 at 21:33, Jyrki Von Karkki wrote: > Hi all > > I'm having a bit of an issue with the postgresql/postgis driver on freebsd. > It simply does not appear installed. How did you install GDAL? The PostgreSQL driver is not always available since it not compiled by default. If you ar

Re: [gdal-dev] Simple circle in OGR

2016-11-02 Thread Mike Toews
tude + " +lon_0=" + longitude + " +x_0=0 +y_0=0" See http://gis.stackexchange.com/a/121539/1872 for this recipe in R. You can make something similar with OSR. Cheers, -Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] VRT raster for reading formatted files?

2016-06-25 Thread Mike Toews
nd, all of other metadata, including grid size, are provided elsewhere. Cheers, -Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Shapefile to SQL Server

2016-04-18 Thread Mike Colbert
uot;, "GEOM_TYPE=geography", "-lco", "GEOM_NAME=geog", "-nln", "CM_SHAPE", "--debug", "ON" ,"-t_srs", "EPSG:4269" ,"-s_srs", "EPSG:26912" }; ogr2ogr.main(cmd); Mike From: gdal-dev [mailto:gdal

[gdal-dev] Shapefile to SQL Server

2016-04-15 Thread Mike Colbert
yer UGRB_Ozone_NAA (use -skipfailures to skip errors) I'm assuming there is nothing unusual about the shapefile. The file is available here: http://deq.wyoming.gov/media/attachments/Air%20Quality/Winter%20Ozone/Nonattainment%20Information/2012_AQD_UGRB-Ozone-Nonattainment-Area-GIS-Shape-File.zip Any help is appreciated. Thanks, Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Hard coded Pi constants

2016-03-03 Thread Mike Toews
This is not just a MSVC issue. M_PI and other constants are not part of the ISO standard, so if you try: $ gcc --std=c99 small_program_with_M_PI.c you will get "error: ‘M_PI’ undeclared" I'm not sure if you can compile GDAL with `--std=c99`, but I don't see any compelling reason to remove #ifndef

Re: [gdal-dev] RFC 61: Call for vote on adoption

2016-02-10 Thread Mike Toews
On 11 February 2016 at 02:05, Even Rouault wrote: > Would someone looking at this discussion and having access to ESRI software be > willing to generate tiny (meaning just one single shape, with the smallest > number of vertices) shapefiles of type PointZ, ArcZ, PolygonZ and MultiPointZ, > and *wi

Re: [gdal-dev] RFC 61: Call for vote on adoption

2016-02-10 Thread Mike Toews
On 11 February 2016 at 01:31, Even Rouault wrote: > Le mercredi 10 février 2016 13:05:20, Peter Halls a écrit : >> ESRI handle this in a non-intuitive way: XYM is supported, but Z >> always has a Measure, so is XYZM! The formal definition is here: >> https://www.*esri*.com/library/whitepa

Re: [gdal-dev] Reading the points from shp file

2016-01-10 Thread Mike Flannigan
laces the data into the $outfile. Place both # the .shp and .shx file in the C:/Copy directory. # The output file is also created in the same directory. I can also send the files if you run Perl. They use the Geo::ShapeFile module. Mike On 1/10/2016 2:00 PM, gdal-dev-requ...@lists.osgeo.org

Re: [gdal-dev] EPSG:6933

2015-11-10 Thread Mike Toews
On 11 November 2015 at 02:33, Even Rouault wrote: > Long answer: > This CRS probably didn't exist in the EPSG database v8.5 that was used the > last time the proj/GDAL derived EPSG database was generated. And I see that > the ellipsoidal Lamber Cylindrical Equal Area projection method 9835 was als

Re: [gdal-dev] gdal2tile

2015-08-24 Thread Mike Manning
Hi Even, Fantastic, thank you very much! I'll give this a go and let you know how I go. Cheers, Mike -Original Message- From: Even Rouault [mailto:even.roua...@spatialys.com] Sent: Monday, 24 August 2015 6:04 PM To: gdal-dev@lists.osgeo.org Cc: Mike Manning Subject: Re: [gda

[gdal-dev] gdal2tile

2015-08-23 Thread Mike Manning
enerate_base_tiles() File "/usr/bin/gdal2tiles.py", line 1276, in generate_base_tiles dsquery.WriteRaster(wx, wy, wxsize, wysize, data, band_list=list(range(1,self.dataBandsCount+1))) File "/usr/lib/python2.7/dist-packages/osgeo/gdal.py", line 796, in WriteRast

Re: [gdal-dev] Why is saving an array to ENVI file format using SaveArray so slow?

2015-07-26 Thread mike c
mike > From: even.roua...@spatialys.com > To: gdal-dev@lists.osgeo.org > Subject: Re: [gdal-dev] Why is saving an array to ENVI file format using > SaveArray so slow? > Date: Fri, 24 Jul 2015 21:31:30 +0200 > CC: mikec7...@hotmail.com > > Mike, > > I've tri

[gdal-dev] Why is saving an array to ENVI file format using SaveArray so slow?

2015-07-24 Thread mike c
our thoughts on this matter would be appreciated. Thanks Regards mike Finished save 3790 function calls in 24.405 seconds Ordered by: standard name ncalls tottime percall cumtime percall filename:lineno(function) 10.0000.0

Re: [gdal-dev] Gdal_calc.py with more than 26 input files

2015-07-16 Thread Mike Toews
not implemented. I've started an enhancement ticket for a -3D option: https://trac.osgeo.org/gdal/ticket/6040 -Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Gdal_calc.py with more than 26 input files

2015-07-16 Thread Mike Toews
al_calc.py -A var_stacked.vrt --outfile=avg_var.tif --calc="A.mean(axis=0)" However, the result avg_var.tif is not as expected due to the use of blocks. Is this ticket worthy? -Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] gdal_rasterize with text attributes.

2015-06-27 Thread Mike Toews
Rasters can only store numeric types, not strings. But you could either rasterize the primary key, or reclassify the strings into a so-called lookup table of ID and string (e.g. 1="green", 2="blue", etc.). And then you could look-up the text from a raster using the common integ

Re: [gdal-dev] Design for sub-second accuracy in OGR ?

2015-04-06 Thread Mike Toews
On 6 April 2015 at 09:39, Even Rouault wrote: > I should have mentionned what currently exists indeed : > > struct { > GInt16 Year; > GByte Month; > GByte Day; > GByte Hour; > GByte Minute; > GByte Second; > GByte TZFlag; /*

Re: [gdal-dev] OGR: How to access Coded Value Domains in File Geodatabases?

2015-01-06 Thread Mike Toews
either hash the values or do an SQL JOIN. -Mike On 7 January 2015 at 14:20, Stefan Keller wrote: > Hi, > > ArcGIS offers "Coded Value Domains" (short: "coded domain") to specify > a valid set of values (code+description) for an attribute [1]. > > Any clues o

Re: [gdal-dev] gdal_calc.py: produce median raster

2014-11-26 Thread Mike Toews
ed/numpy.median.html -Mike On 24 November 2014 at 23:14, Simen Langseth wrote: > I have one GeoTiff file with 3 bands. I want to produce a raster file > computing pixel wise median value for the three raster bands. How can I do > that? > > I tried as follows: > > gdal_calc.py -A i

Re: [gdal-dev] Call for discussion on RFC 51: RasterIO() improvements : resampling and progress callback

2014-11-25 Thread Mike Toews
fer a few prepared use-cases, or somehow provide extra parameters. Or just provide a single use case, similar to Gaussian blur routines in various raster editing software. -Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Questions about polygonization

2014-08-19 Thread Mike Toews
r[2, 2] = 12 gt = [0, 1, 0, 6, 0, -1] for (sh, v) in rasterio.features.shapes(ar, transform=gt): print('%s: %s' % (v, shape(sh))) 11: POLYGON ((0 6, 0 0, 5 0, 5 6, 0 6), (2 4, 3 4, 3 3, 2 3, 2 4)) 12: POLYGON ((2 4, 2 3, 3 3, 3 4, 2 4)) -Mike

Re: [gdal-dev] [BULK] Re: Dissolve shapefile using GDAL/OGR

2014-07-14 Thread Mike Toews
The path strings are not escaped, so dissolve[13] is '\a' or '\x07' and not 'a', which makes the path invalid. For Windows paths, the best practice is to either use forward slashes (like unix), or use raw escape with r as a string prefix, e.g. dissolve = r'C:\OSGeo4W64\apps...' On 12 July 2014 0

Re: [gdal-dev] ogr2ogr.py

2014-06-27 Thread Mike Saunt
Cool, and understood. It caught me out, that is all and I'd like to make sure it doesn't for others... Sent from my phone which isn't an iPhone On 27 Jun 2014 07:40, "Chaitanya kumar CH" wrote: > Mike, > > Thanks for the update. But, I have to mention that the

[gdal-dev] ogr2ogr.py

2014-06-26 Thread Mike Saunt
know or pointing me to the relevant web page with details it would be much appreciated Many thanks Mike​ -- ​​ -- *Mike Saunt* Astun Technology Ltd, The Coach House, 17 West Street, Epsom, Surrey, KT18 7RL, UK t:+44 1372 744 009 m:+44 788 680 ***

Re: [gdal-dev] Add to PythonGotchas: saving and closing raster datasets

2014-05-27 Thread Mike Toews
", I think I'll use `del ds` as the recommended method, with mention that any copies of the references to the dataset/source also need to be dereferenced. -Mike [1] http://www.gdal.org/ogr__api_8h.html#a9d845a6cf6652756925530418905471a ___ g

Re: [gdal-dev] compiling without writing permission - problem gdal-1.10.1

2014-05-06 Thread Mike Toews
he Python bindings manually with the appropriate --prefix: $ cd swig/python $ python setup.py install --prefix=/to/where/you/want/it or edit swig/python/GNUmakefile (line 71): $(PYTHON) setup.py install --prefix=$(DESTDIR)$(prefix) -Mike ___ gdal-dev mailin

[gdal-dev] Add to PythonGotchas: saving and closing raster datasets

2014-05-03 Thread Mike Toews
es: obj = None or del obj I'm unsure of which is regarded best practice, but encourage consistency. Please edit as necessary or reply with suggestions. -Mike === Saving and closing datasets === To save and close raster datasets, the object needs to be dereferenced, such as setting

Re: [gdal-dev] Installing GDAL on Win7

2014-03-24 Thread Mike Flannigan
go with the 85% compression. A 10.1MB PDF turns into a 26.2MB TIFF, which would zip to 26.0MB if anybody cared to zip it. Thanks for all your help. You guys saved me a lot of time. Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Installing GDAL on Win7

2014-03-23 Thread Mike Flannigan
le. These files are about 15 times bigger than the ones I am used to, but they do have quite a bit of detail and don't have the blank white space that a 7.5 minute topo has. These Tiffs are very well georeferenced. I am pleased. Mike On 3/23/2014 7:31 PM, Mike Flannigan wrote: Thank

Re: [gdal-dev] Installing GDAL on Win7

2014-03-23 Thread Mike Flannigan
wpt.tif --config GDAL_PDF_LAYERS_OFF "Map_Collar","Map_Frame.Projection_and_Grids" --config GDAL_PDF_DPI300 It gave an error: ERROR 1: --config option given without a key and value argument. I'm still working with it. I don't want tiled - I want the border of the drawing included in the G

[gdal-dev] Installing GDAL on Win7

2014-03-23 Thread Mike Flannigan
a basic convert utility. I'm going to this page: http://www.gisinternals.com/sdk/ but I don't know which download I need. I would like a 64 bit version. Should I use the 2nd item under "GDAL and MapServer", which is the "MSVC2005 (Win64) -release" one? O

[gdal-dev] ogr.GeometryTypeToName oddness

2013-04-16 Thread Mike Toews
3D geometries by ogr.GeometryTypeToName(), but this is false. Furthermore, other OGR documents refer to these geometries as so called 2.5D geometries, and not true 3D geometries. - The 2.5D geometries (e.g. wkbLineString25D) are simply reported as their geometry type, and no mention of dimensions. It seems the "3D " name prepending logic is reversed. -Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] Capture stderr in interactive Python shell

2012-08-05 Thread Mike Toews
in a shell (cmd.exe or Bash): ERROR 4: `noexist.tif' does not exist in the file system, and is not recognised as a supported dataset name. It appears this output is sent to stderr. How can this output also appear in an interactive Python shell, such as PythonWin or IDLE?

[gdal-dev] gdal_translate and TIFF resolutions(?)

2011-12-15 Thread Mike O'Malley
an I force gdal_translate to set that resolution field? I'm not a TIFF expert, but it seems like it would be a standard field if tiffinfo throws an error when it's not populated. Thanks, Mike ___ gdal-dev m

[gdal-dev] GDAL, Java and UINT16

2011-10-27 Thread Mike O'Malley
#x27;t find it. The Worldwind code handles bytes, signed int16s, and signed int32s. Thanks, Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

RE: [gdal-dev] Quick question regarding image translation in Java (I think)

2011-07-28 Thread Mike O'Malley
That was it exactly. Next time I'll read the fine print a bit more closely. Thanks Even! Mike -Original Message- From: Even Rouault [mailto:even.roua...@mines-paris.org] Sent: Thursday, July 28, 2011 2:44 PM To: gdal-dev@lists.osgeo.org Cc: Mike O'Malley Subject: Re:

[gdal-dev] Quick question regarding image translation in Java (I think)

2011-07-28 Thread Mike O'Malley
the NITFs in my set. The NITFs are viewable using FWTools and look good; the TIFFs created by gdal_translate match them. It's just my Java code that's producing blank images. Am I missing something obvious? Thanks, Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] gdal_translate and gdal_warp

2011-06-15 Thread Mike Willis
neighbour, and bicubic as the gdalwarp resampling methods but have not succeeded in getting anything useable. Any other suggestions? The files are Gb in size. We have the RPCs and full metadata for each of the files. Cheers Mike Willis Postdoctoral Research Associate Earth and Atmospheric S

Re: [gdal-dev] Lengthy GDALRasterAttributeTable output

2011-06-02 Thread Mike Toews
Perfect, thanks! That worked for both gdalinfo and for derived data via CreateCopy(). -Mike On 2 June 2011 21:21, Even Rouault wrote: > Le jeudi 02 juin 2011 00:54:48, Mike Toews a écrit : >> I'm struggling to understand why both the aux.xml file and output from >> gdalinfo

[gdal-dev] Lengthy GDALRasterAttributeTable output

2011-06-01 Thread Mike Toews
I'm struggling to understand why both the aux.xml file and output from gdalinfo are really bulky in file size for my GeoTIFFs. First, output from gdalinfo on a 375 KB GeoTIFF file has 196664 lines: Driver: GTiff/GeoTIFF Files: HorizB.tif HorizB.aux Size is 317, 301 Coordinate System is: ...

Re: [gdal-dev] Adding field to Shapefile via Python

2011-03-20 Thread Mike Toews
s in it. I would recommend adding a "See also" link from AddFieldDefn to: http://www.gdal.org/ogr/ogr__api_8h.html#ab585ef1166c61c4819f7fd46ee4a275 -Mike On 18 March 2011 18:45, Chaitanya kumar CH wrote: > > Mike, > > A new field cannot be added to a feature definition while

[gdal-dev] Adding field to Shapefile via Python

2011-03-17 Thread Mike Toews
re() feature.GetField('EXIST') # an existing float field with data feature.SetField('EXIST', 111.1) # works feature.SetField('MYFLD', 123) # does nothing layer.SetFeature(feature) Which successfully updates the data in the existing field, but not

Re: [gdal-dev] Polygon topology

2011-03-08 Thread Mike Toews
o, except that assumption doesn't work for MultiPolygon geometry types. -Mike ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Polygon topology

2011-03-08 Thread Mike Toews
x27;, g.num_holes, 'hole'+('s' if g.num_holes != 1 else '') feature = layer.GetNextFeature() num += 1 -Mike On 8 March 2011 21:57, Simon Lyngby Kokkendorff wrote: > Hi List, > >   I am using ogr via the python bindings to construct various poly

[gdal-dev] Re: gdal-dev Digest, Vol 82, Issue 4

2011-03-03 Thread Mike Leahy
of other bits that will still need fixing. Mike > Date: Wed, 02 Mar 2011 10:46:02 -0400 > From: Jeff McKenna > Subject: Re: [gdal-dev] Re: PHP bindings > To: gdal-dev@lists.osgeo.org > Message-ID: <4d6e582a.5080...@gatewaygeomatics.com> > Content-Type: text/p

Re: [gdal-dev] Re: PHP bindings

2011-03-01 Thread Mike Leahy
be accounted for in some way through the php typemaps. Regards, Mike > Date: Tue, 1 Mar 2011 19:59:48 +0100 > From: Even Rouault > Subject: Re: [gdal-dev] Re: PHP bindings > To: gdal-dev@lists.osgeo.org, mgle...@alumni.uwaterloo.ca > Message-ID: <201103011959.48700.even.roua...

Re: [gdal-dev] Re: PHP bindings

2011-02-28 Thread Mike Leahy
k that's out of my league, however. Hopefully this is enough to peak the interest of someone who has a little more familarity with c/swig/php, and we could get a more fully-supported set of PHP bindings. Best regards, Mike On Thursday, February 24, 2011 19:45:16 Mike Leahy wrote: > Hell

Re: [gdal-dev] Re: PHP bindings

2011-02-24 Thread Mike Leahy
ring: So...with this in mind, does this suggest that PHP will essentially need an entirely rewritten list of functions that substitute those that are currently there now? I would hope that there's a more efficient/effective way to go about this, and that I'm just looking at this with

[gdal-dev] Re: gdal-dev Digest, Vol 81, Issue 27

2011-02-11 Thread Mike Leahy
ts with the way you've done it. Best regards, Mike On Friday, February 11, 2011 12:00:11 Joris van de Sande wrote: > Date: Fri, 11 Feb 2011 12:50:41 +0100 > From: Joris van de Sande > Subject: [gdal-dev] Reviving the PHP Bindings > To: gdal-dev@lists.osgeo.org > > Message-

Re: [gdal-dev] Re: PHP bindings

2011-02-09 Thread Mike Leahy
43295]]" So maybe I just don't know how to interact with the spatial reference object within the PHP environment. Any further suggestions here? Regards, Mike On Wednesday, February 09, 2011 14:14:55 Even Rouault wrote: > Le mercredi 09 février 2011 08:47:25, Mike Leahy a éc

Re: [gdal-dev] Re: PHP bindings

2011-02-08 Thread Mike Leahy
= At this point, there's not much more I can do to without some guidance from someone with a better understanding of where things could be going wrong. Is this perhaps a topic that would be better discussed with the SWIG users list? Best regards, Mike On Tuesday, February 08,

[gdal-dev] Re: PHP bindings

2011-02-07 Thread Mike Leahy
n the php module isn't able to access the installed gdal libraries? Any suggestions would be welcome. I'm testing this in a minimal 64-bit Ubuntu 10.11 VM. Regards, Mike On Sunday, February 06, 2011 14:50:34 Mike Leahy wrote: > Hello List, > > Does anyone have an ide

  1   2   >