Re: [gdal-dev] ogr2ogr - double free or corruption (out)

2020-03-16 Thread Andre Joost
working with GDAL 3.0.4. HTH, Andre Joost ___ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] Gdalwarp with geolocation Array and nodata value

2019-04-21 Thread Andre Joost
remaining GCPs). Is there a config option for GDALwarp to omit those nodata values in the geolocation array? Or is it worth to open a ticket for it? Greetings, Andre Joost ___ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org

Re: [gdal-dev] convert india-latest.osm to india-latest.tab

2019-03-07 Thread Andre Joost
27; does not already exist in the output dataset, and cannot be created by the output driver. Looks like a limitation of the OGR MapInfo driver. Are you bound to use that format? I guess Mapinfo can handle others as well. HTH, Andre Joost ___ gdal-d

Re: [gdal-dev] gdalwarp cutline shift between raster and vector in GDAL 2.3.2

2018-12-17 Thread Andre Joost
Hi Pedro, if you compare gdalsrsinfo epsg:20790 and gdalsrsinfo EPSG:102164 from ESRI, you see that ESRI omits the datum shift between Lisbon Hayford and WGS84. So you should use -s_srs EPSG:20790 to get the warping correctly. HTH, Andre Joost Am 17.12.18 um 11:09 schrieb Pedro Venâncio

Re: [gdal-dev] HDF data mis-applied sinusoidal

2018-11-30 Thread Andre Joost
much about it. You can sanitize it with gdal_translate -a_ullr -180 90 180 -90 -a_srs EPSG:4326 src dst HTH, Andre Joost Am 29.11.18 um 21:31 schrieb Michael Sumner: Hello, these data in HDF4 format are interpreted as being in sinusoidal projection, I guess some conflation with the MODIS

Re: [gdal-dev] How to find correct 2D bounding box to subset in different CRS with gdal (different outputs for gdalwarp)?

2018-11-21 Thread Andre Joost
:32632 rectangle. The difference will be filled with NODATA pixels. Try to vsualize your problem in a GIS software like QGIS, and you see what I mean. HTH, Andre Joost Am 21.11.18 um 16:33 schrieb Pham Huu Bang: Hello, I'm using GDAL version 1.11.4 and I have a problem how to use gd

[gdal-dev] Netcdf coordinates swapped

2018-10-12 Thread Andre Joost
Dear list, can anybody give some insight on this netcdf file: https://github.com/Andre-J/testing/raw/master/3A-DAY.nc ncdump reports: netcdf \3a-day { dimensions: AD = 2 ; chn = 2 ; nlon = 1440 ; nlat = 536 ; and Panoply displays it fine, but gdal_translate

Re: [gdal-dev] Warping EUMetsat MSG-4 MPE dataset to webmercator projection

2018-04-25 Thread Andre Joost
Am 25.04.18 um 15:50 schrieb Kristian Evers: Andre, The PROJ 5.0.1 OSGeo4W installation is faulty so that may be the problem. Specifically the null grid is not being installed in the 5.0.1 package as reported on the mailing list [0]. This has not been resolved yet. Try copying in the null grid

Re: [gdal-dev] Warping EUMetsat MSG-4 MPE dataset to webmercator projection

2018-04-25 Thread Andre Joost
Am 25.04.18 um 14:45 schrieb Andre Joost: Am 25.04.18 um 14:29 schrieb Andre Joost: Just checked here: It worked up to GDAL 2.2.3 (from Gisinternals) on the full target extent, and fails with the current 2.2.4. ... and gdal 2.4.0dev from Gisinternals works the same as 2.2.3. So only the

Re: [gdal-dev] Warping EUMetsat MSG-4 MPE dataset to webmercator projection

2018-04-25 Thread Andre Joost
Am 25.04.18 um 14:29 schrieb Andre Joost: Just checked here: It worked up to GDAL 2.2.3 (from Gisinternals) on the full target extent, and fails with the current 2.2.4. ... and gdal 2.4.0dev from Gisinternals works the same as 2.2.3. So only the OSGEO4W version 2.2.4 fails. Maybe because of

Re: [gdal-dev] Warping EUMetsat MSG-4 MPE dataset to webmercator projection

2018-04-25 Thread Andre Joost
up to GDAL 2.2.3 (from Gisinternals) on the full target extent, and fails with the current 2.2.4. And the result was not all black. HTH, Andre Joost ___ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Warping EUMetsat MSG-4 MPE dataset to webmercator projection

2018-04-21 Thread Andre Joost
he visible sector. Greetings, Andre Joost ___ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] ogr2ogr using grid shift files

2018-03-07 Thread Andre Joost
inition files, and does normally not fall back to the proj4 folder definitions. C:\OSGeo4W64\share\gdal\datum_shift.csv is the file to look for, containg a definition for EPSG:4314 DHDN, which will be used for all DHDN based projected CRS. So try to adjust that file. HTH, Andre Joost

Re: [gdal-dev] ogrInfo error - Entry Point Not Found

2018-01-29 Thread Andre Joost
compression_methods could not be located in the dynamic link library SSLEAY32.dll Most probably you have an outdated version of SSLEAY32.dll, or different versions in different places. HTH, Andre Joost ___ gdal-dev mailing list gdal-dev@lists.osgeo.or

Re: [gdal-dev] GDALWARP Error

2018-01-16 Thread Andre Joost
er to https://gis.stackexchange.com/questions/154339/unable-to-warp-hdf5-files: lon.vrt 1 lat.vrt 1 0 0 1 1 I removed the GCP section of the mir.vrt file, and there is no error message anymore. HTH, Andre Joost

Re: [gdal-dev] Vertical datum conversion for SRTM data

2017-11-22 Thread Andre Joost
Am 21.11.2017 um 12:27 schrieb Michael Smith: Egm96 grid shift file available from http://download.osgeo.org/proj/vdatum/egm96_15/egm96_15.gtx Note that OSGEO4W64 still ships an incorrect version of the file from 2010. According to http://download.osgeo.org/proj/vdatum/egm96_15/ the file s

Re: [gdal-dev] transformation failure from WGS 84 to GRS 1980

2017-11-11 Thread Andre Joost
Am 11.11.2017 um 08:45 schrieb Saulteau Don: I'm running Arch Linux :) I don't have that running, but looking at https://www.archlinux.org/packages/?sort=&q=gdal&maintainer=&flagged= I see a package 2.2.2-2 from 2017-11-09 and a package 2.2.2-3 from today. Maybe the packaging went wrong on

Re: [gdal-dev] transformation failure from WGS 84 to GRS 1980

2017-11-10 Thread Andre Joost
Am 09.11.2017 um 20:52 schrieb Donovan Cameron: Afternoon! I recently upgraded to gdal 2.2.2 from 2.1 on linux and when I try to use ogr2ogr to append/update an existing table in PostgreSQL/PostGIS the command fails because it can't apply a transformation. The command works fine in gdal 2.1.

Re: [gdal-dev] GDAL 2.2.3 planning

2017-11-02 Thread Andre Joost
Am 01.11.2017 um 11:32 schrieb Even Rouault: Hi, I've add an exchange with Sebastiaan Couwenberg about a schedule that would give GDAL 2.2.3 a chance to enter the next Ubuntu LTS 18.04. While the date for the import freeze from Debian to Ubuntu is not know yet, based on the last release cycle, i

Re: [gdal-dev] Which Proj.4 transforms are available in GDAL?

2017-09-05 Thread Andre Joost
Am 05.09.2017 um 20:50 schrieb Even Rouault: Apparently they also offer the datasets in WGS 84, so you can perhaps check which hypothesis is right from that. Both datasets align with the Custom CRS including 3-parms datum shift I included earlier in this thread. Greetings, André Joost __

Re: [gdal-dev] Which Proj.4 transforms are available in GDAL?

2017-09-05 Thread Andre Joost
Am 05.09.2017 um 19:43 schrieb Even Rouault: The shapefile clearly states DATUM["D_North_American_1927", and the download page (from the official agency) http://www.cnr.gob.sv/geoportal-cnr/ tells the same. Hum, actually looking at the full .prj, it looks similar to EPSG:5460 Yes, projection

Re: [gdal-dev] Which Proj.4 transforms are available in GDAL?

2017-09-05 Thread Andre Joost
Am 05.09.2017 um 16:56 schrieb Even Rouault: The problem is that there is no EPSG code for the projection, Ah, the -s_srs syntax is not limited to EPSG:. You can use a proj.4 string, inlined WKT, a filename that has WKT or proj.4 string in it, etc. See http://gdal.org/gdal_utilities.html

Re: [gdal-dev] Which Proj.4 transforms are available in GDAL?

2017-09-05 Thread Andre Joost
Am 01.09.2017 um 21:05 schrieb Even Rouault: On vendredi 1 septembre 2017 20:42:34 CEST Andre Joost wrote: Am 01.09.2017 um 20:14 schrieb Even Rouault: ogr2ogr -s_srs EPSG:4326 -t_srs "+proj=lcc +lat_1=13.316667 +lat_2=14.25 +lat_0=13.78 +lon_0=-89 +x_0=50 +y_0=295809.184

[gdal-dev] Adding missing EPSG codes to PROJ.4 and GDAL

2017-09-02 Thread Andre Joost
Dear all, I have recently compared the downloadable SQL files of the EPSG database (v9.1) with the output of gdalsrsinfo (v2.2.1). There are currently 368 CRS_codes missing, and I would like to have them reduced. You can find a csv list at https://github.com/Andre-J/EPSG/blob/master

Re: [gdal-dev] Which Proj.4 transforms are available in GDAL?

2017-09-01 Thread Andre Joost
Am 01.09.2017 um 20:14 schrieb Even Rouault: ogr2ogr -s_srs EPSG:4326 -t_srs "+proj=lcc +lat_1=13.316667 +lat_2=14.25 +lat_0=13.78 +lon_0=-89 +x_0=50 +y_0=295809.184 +k_0=0.6704 +ellps=clrk66 +units=m +no_defs +towgs84=0,125,194,0,0,0,0 +wktext" LimitesNAD27.shp dptoA_WGS_

Re: [gdal-dev] Which Proj.4 transforms are available in GDAL?

2017-09-01 Thread Andre Joost
Am 25.08.2017 um 13:52 schrieb Even Rouault: For unsupported projections, add +wktext. A dummy WKT representation is then built that captures the full proj.4 string in a EXTENSION node. Can also serve for supported projections if you want to overrite the WKT -> proj.4 string conversion done by

Re: [gdal-dev] Which Proj.4 transforms are available in GDAL?

2017-08-25 Thread Andre Joost
Am 25.08.2017 um 13:52 schrieb Even Rouault: And also found some that do no agree between GDAL and proj. For example echo 4.897 52.371 | proj +proj=vandg +ellps=WGS84 +units=m 489005.93 6431581.02 gdaltransform -s_srs EPSG:4326 -t_srs "+proj=vandg +ellps=WGS84 +units=m" 4.897 52.371 4884

Re: [gdal-dev] stere

2017-07-14 Thread Andre Joost
Am 14.07.2017 um 19:06 schrieb Patrick L. Francis: I work with NWS grids a lot for nexrad data. They are in a 1/40 lfm grid. Your data looks like it might be an a 1/1 LFM grid. Here is the WKT I would use for a 1/1 LFM grid. Each cell is 190.5 km. thank you for replying brian the output

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

2017-06-06 Thread Andre Joost
Am 07.06.2017 um 02:35 schrieb 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

Re: [gdal-dev] Question in orthographic projection usage with PROJ4

2017-02-01 Thread Andre Joost
Am 31.01.2017 um 17:46 schrieb niels.vonfestenberg-packisch@rohde-schwarz.com: Hi all, I would like to transform local cartesian coordinates into longitude and latitude. According to all the documentation I could find, I tried doing this (in C++, specifying a center for the orthographic proje

Re: [gdal-dev] OSM - delineation of rivers

2017-01-02 Thread Andre Joost
Am 01.01.2017 um 22:52 schrieb Simon Lyngby Kokkendorff: And indeed - worked like a charm when I ran ogr2ogr on a pbf file containing the whole of England. The River Thames is represented nicely, as was the case with osm2pgsql. Was surprisingly fast also compared to osm2pgsql! So definitely a "t

Re: [gdal-dev] identifying the fields from the .osm file

2016-11-20 Thread Andre Joost
Am 19.11.2016 um 20:45 schrieb Djordje Spasic: Thank you Etienne. I am very interested in your reply, but I haven't quite understand you. The ogr2ogr.exe file is using a particular python script to read the .osm file and identify its unique fields? Can you be a bit more specific please? Solving

Re: [gdal-dev] splitting other_tags from .osm file

2016-10-30 Thread Andre Joost
Am 30.10.2016 um 13:55 schrieb Djordje Spasic: Jukka mentioned all_tags, of which I haven't been aware of. So my issue now is: what is the difference between other_tags and all_tags? other_tags includes all tags that are **not** covered by the attribute fields list in osmconf.ini. all_ta

Re: [gdal-dev] Error in gdalwarp

2016-10-22 Thread Andre Joost
Am 22.10.2016 um 18:30 schrieb Paolo Cavallini: gdalwarp -ot Float32 -q -of GTiff -cutline watersheds.shp -co COMPRESS=DEFLATE -co PREDICTOR=1 -co ZLEVEL=6 -wo OPTIMIZE_SIZE=TRUE dem25.tif OUTPUT.tif Works for me with GDAL 1.11.3 and 2.0-dev, but fails on GDAL 2.1.0 and 2.1.2. All Gisintenal b

[gdal-dev] unable to create jpeg200 files with self compiled gdal and openjpeg

2016-07-12 Thread Andre Hollstein
Dear gdal experts, I'm compiling gdal and a whole list of other libraries myself. Currently I'm using gdal-2.1.0 and openjepeg v2.1.1 and use the libs from python: import gdal import numpy as np gdal.AllRegister() #driver = 'JP2KAK' driver = "JP2OpenJPEG" #driver = "GTiff" #driver = "JPEG2000"

Re: [gdal-dev] prebuilt GDAL 2.0.2 windows binary

2016-05-09 Thread Andre Joost
Am 06.05.2016 um 08:54 schrieb Gane R: I see a link to download the source http://download.osgeo.org/gdal/2.0.2/gdal202.zip when can I find a prebuilt version of GDAL 2.0.2 OSGEO4W offers GDAL 2.0.2, see http://download.osgeo.org/osgeo4w/x86_64/versions.html Gisinternals has GDAL 1.11.3 a

Re: [gdal-dev] Setting up GDAL Custom Transformation from ESRI GTF File

2016-04-03 Thread Andre Joost
Am 03.04.2016 um 14:05 schrieb Derek Stuart: Andre: The only Gulshan coordinate system I see with EPSG 4682 is a geographic coordinate system. I am using a projected coordinate system called BTM. This is a fairly new coordinate system, but it is the standard here in Bangladesh. PROJCS[&quo

Re: [gdal-dev] Setting up GDAL Custom Transformation from ESRI GTF File

2016-03-31 Thread Andre Joost
Am 31.03.2016 um 17:16 schrieb Derek Stuart: I am using GDAL routines as part of RAS Mapper within HEC-RAS. I am coming from ArcMap where I am using a custom transformation GTF file to translate from the D_Everest_Bangladesh datum to WGS_1984. The GTF file includes the following: GEOGTRAN["ev

Re: [gdal-dev] Transformation methods & parameters (OGRCoordinateTransformation)

2016-03-19 Thread Andre Joost
Am 18.03.2016 um 18:57 schrieb César Martínez: Thanks André, datum_shift.csv is (partially) what I was looking for. I see however that EPSG:15933 transformation operation is missing on that list (which in this Spanish example is usually the most accurate one), maybe because it uses the NTv2 meth

Re: [gdal-dev] Transformation methods & parameters (OGRCoordinateTransformation)

2016-03-19 Thread Andre Joost
Am 18.03.2016 um 17:34 schrieb César Martínez: This is relevant because sometimes there are several transformation operations available between 2 CRSs, and we'd like our users being able to choose which operation they'd like to use. For instance, between EPSG:23030 and EPSG:4326, at least the f

[gdal-dev] Feature ID propagation in SQLite data source

2016-02-29 Thread Andre Vautour
Hi all, When I execute an SQL query on a SQLite source, I am observing that the feature identifiers are changing. This was observed using 1.11.1. ogrinfo -geom=NO -fid 1227 test2.sqlite OGRFeature(roads):1227 gml_id (String) = Roads.371c00110746 arc_flag (String) = T line_com

Re: [gdal-dev] osmconf.ini

2016-02-20 Thread Andre Joost
Am 20.02.2016 um 16:18 schrieb Even Rouault: Frank, I clicked on the "Export" button on the view centered on http://www.openstreetmap.org/relation/2827792 , and ran $ ogrinfo ~/Téléchargements/map.osm multilinestrings \ -where "osm_id = '2827792'" -geom=no -al -q \ --config OSM_USE_CUSTOM_IND

Re: [gdal-dev] GDAL 2.0.2 doesn't compile with Kakadu 7.7.1

2016-02-12 Thread Andre Hollstein
a lot again, André On 02/12/2016 10:03 AM, Even Rouault wrote: Andre, The preferred way to include all the Kakadu API's is to link against >both libkdu_v77R.o and libkdu_a77R.o, where the latter is built by the "managed" makefiles. >Anyway, it may be that

[gdal-dev] GDAL 2.0.2 doesn't compile with Kakadu 7.7.1

2016-02-12 Thread Andre Hollstein
cw=no \ --with-libtiff=no --without-libtool \ --with-hdf4=no --with-hdf5=no --with-netcdf=no \ --with-geotiff=internal \ --with-grass=no --with-sqlite3=no \ --with-kakadu=/misc/fluo6/andre/projekte/install_py_et_al/python_GFZ/opt --with-pic \ --with-python=$prefix_python/bin/python \ Output: GDAL i

Re: [gdal-dev] convert shapefile from California State Plane to EPSG:4326

2016-02-09 Thread Andre Joost
Am 09.02.2016 um 03:35 schrieb maning sambale: Thanks Daniel and Even, Using following Even's advice, I get better results: https://gist.github.com/maning/093a866fd465ea880032#better-conversion Wondering if it might be appropriate to change proj.4 defaults to the above ? +1 The datum shift

Re: [gdal-dev] pdf file - raster and vector differences

2016-02-02 Thread Andre Joost
Am 02.02.2016 um 15:37 schrieb Piotr Pachół: Hello, I would like to ask If there is a way (by means of gdal tools) to recognize that pdf has vector features ? I have hundreds of geospatial pdf files (raster and vector) and I don't want to open each of file in Acrobat Reader to verify if it is r

Re: [gdal-dev] GDAL_WMS/TMS with local CRS and bbox problem

2016-01-15 Thread Andre Joost
Am 15.01.2016 um 12:46 schrieb Stefan Keller: Hi I have a GDAL .vrt from a GDAL_WMS/TMS which was tweaked to projection/CRS EPSG:21781 [see below *]. The problem here is, that the bbox is way too large. It should be Switzerland but it's larger than europe. (in 21781 it's lowerleftX 48 lover

Re: [gdal-dev] unable to open DGN

2015-11-23 Thread Andre Vautour
On 2015-11-23 12:00, Even Rouault wrote: Le lundi 23 novembre 2015 16:48:48, Andre Vautour a écrit : Right, thanks, I wasn't thinking of that when I wrote my reply. Guess it makes sense for a closed format such as DWG. Might not make as much sense to use it for DGN which is more open D

Re: [gdal-dev] unable to open DGN

2015-11-23 Thread Andre Vautour
nt on that: In order to build the DWG driver, you or your company, need to be part of some kind of alliance to be able to download a full version of the Teigha SDK. AFAIK From: gdal-dev on behalf of Andre Vautour Sent: Monday, November 23, 2015 9:58

Re: [gdal-dev] unable to open DGN

2015-11-23 Thread Andre Vautour
On 2015-11-20 17:40, Brad Hards wrote: On Fri, 20 Nov 2015 02:20:55 PM Martin Landa wrote: I am trying to open sample DGN file, but it fails. I guess that the reason is unsupported version (man pages says: Microstation DGN files from Microstation versions predating version 8.0 are supported for

[gdal-dev] OGR error handling

2015-10-20 Thread Andre Vautour
Hi all, We are using OGRSFDriver::CopyDataSource() to copy a GML data source to a memory data source on v1.11.x. We have encountered a problem where one of the feature tags of a particular GML file we are trying to read is malformed. What we are seeing is that an error occurs in OGRGMLLayer::

Re: [gdal-dev] "Nested" vrts?

2015-09-24 Thread Andre Joost
Am 24.09.2015 um 16:43 schrieb Jonathan Greenberg: Folks: Is it possible to create a vrt comprised of other vrts? For a mosaicking problem, say I want to first trim the margins of the individual tiles by a certain number of pixels -- I can, of course, gdal_translate this but this adds a bunch o

Re: [gdal-dev] Question

2015-08-28 Thread Andre Joost
Am 26.08.2015 um 19:22 schrieb Even Rouault: Something that would be possible on Windows, but not on Linux?, come on ;-) ECW and MrSid support out-of-the-box ? ;-) SCNR, André Joost ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.o

Re: [gdal-dev] Question

2015-08-26 Thread Andre Joost
Am 26.08.2015 um 17:27 schrieb Rahkonen Jukka (MML): I noticed right now while I was testing that I have several GDAL installations in my PATH, perhaps you have the same situation. Write PATH into your command line for checking it. To avoid those PATH assignments, I use the zip versions from G

Re: [gdal-dev] guidance on NSIDC Sea Ice concentration data

2015-07-26 Thread Andre Joost
Can this VRT approach be used for MODIS sea ice data too? http://gis.stackexchange.com/questions/112547/reprojecting-modis-ice-sea-temperature-products-with-gdalwarp http://gis.stackexchange.com/questions/155326/reprojecting-modis-sea-ice-extent-and-ist-data-using-gdal note that especially olde

Re: [gdal-dev] Getting warning when opening NetCDF file using gdal?

2015-07-14 Thread Andre Joost
Have you tried simply gdalinfo A201504506.L2_LAC_OC.nc ? It works for me with GDAL 1.11.2 from gisinternals. HTH, Andre Joost Am 14.07.2015 um 05:58 schrieb oz: Do you have any idea on how will I batch process gdalwarp these HDF5 files? -- View this message in context: http://osgeo

Re: [gdal-dev] Determining GDAL WMS download limit?

2015-07-07 Thread Andre Vautour
On 2015-07-07 15:34, Carl Godkin wrote: Hi Andre, Thank you very much for the reply. I hadn't thought of querying the levels as you suggest but I'll remember next time that it makes a good test. However, I'm not sure what I would do with the 7 suggestion you make. Eve

Re: [gdal-dev] Determining GDAL WMS download limit?

2015-07-07 Thread Andre Vautour
Looks like that server only has 7 levels for that satellite layer: http://mt.google.com/vt/lyrs=s&x=0&y=0&z=7 http://mt.google.com/vt/lyrs=s&x=0&y=0&z=8 The terrain and map layers have 22 levels: http://mt.google.com/vt/lyrs=t&x=0&y=0&z=22 http://mt.google.com/vt/lyrs=t&x=0&y=0&z=23 http://mt.goo

Re: [gdal-dev] Validating a projection

2015-06-07 Thread Andre Joost
Am 08.06.2015 um 01:38 schrieb Doug McCorkle: Andre, I have updated the projection information with the recommended names. When I process GeoTIFF files with this projection everything works fine while all of the data is in memory and none of the Datasets in question are written to disk. If

Re: [gdal-dev] Validating a projection

2015-06-07 Thread Andre Joost
Am 07.06.2015 um 18:50 schrieb Doug McCorkle: Hello, I am trying to add this projection to our runtime environment: PROJCS["USA_Contiguous_Albers_Equal_Area_Conic_USGS_version", GEOGCS["GCS_North_American_1983", DATUM["D_North_American_1983", SPHEROID["GRS_1980",6378137.0,298.257222101]], PRIME

Re: [gdal-dev] New OGC standard about WKT for projections

2015-05-11 Thread Andre Vautour
The way I see it, the new spec does clear up some big holes like: 1. Being able to specify an unambiguous operation/projection method (by identifier) for projections instead of loosely relying on the projection's name (9.3.3) 2. Knowing the units of projection parameters (9.3.4) On 2015-05

Re: [gdal-dev] Problems with WMS minidriver

2015-04-28 Thread Andre Vautour
On 2015-04-28 11:19, Federico wrote: Hi guys, i'm trying to get a wmts layer but i'm getting an error. I think that i'm doing something wrong but i don't what's wrong. URL: http://opencache.statkart.no/gatekeeper/gk/gk.open_wmts GetCapabilities: http://opencache.statkart.no/gatekeeper/gk/gk.open

Re: [gdal-dev] Encoding EPSG:3857 (WebMercator) in GeoTIFF, and ArcGIS interoperability

2015-04-16 Thread Andre Vautour
On 2015-04-16 12:07, Even Rouault wrote: - All formulations that try to expand the definition with ProjCoordTransGeoKey = CT_Mercator, its projection parameters and GCS parameters aren't really appropriate, since there's no way of capturing that its a spherical mercator that must be used. Well

Re: [gdal-dev] Encoding EPSG:3857 (WebMercator) in GeoTIFF, and ArcGIS interoperability

2015-04-16 Thread Andre Vautour
On 2015-04-16 07:08, Even Rouault wrote: Le mercredi 15 avril 2015 22:48:45, Even Rouault a écrit : The way I see it, there are two different ways to model Google Mercator: 1. As a WGS84 datum/ellipsoid with a custom Mercator-based projection which uses only the semi-major axis, as seems to be

Re: [gdal-dev] Encoding EPSG:3857 (WebMercator) in GeoTIFF, and ArcGIS interoperability

2015-04-15 Thread Andre Vautour
On 2015-04-15 14:59, Even Rouault wrote: Le mercredi 15 avril 2015 18:21:51, xavier lhomme a écrit : Dis you try with 102113 or 102100 instead of 3857 ? I had once to use those code in order to be discovered as 3857 in arcgis. when i created geotiff with gdal warp/transform. Xavier, thanks f

Re: [gdal-dev] Trouble with French Ecw map

2015-04-14 Thread Andre Joost
Am 15.04.2015 um 01:23 schrieb Nicolas Ragg: Hello all I have a french map in ECW format, i guess it is from IGN. Srs seems ok but when i import the thing in qgis or try to warp the file it is all shifted to north (Around Sweden). The output from gdalinfo 1.11 is PROJ.4 : '+proj=lcc +lat_1=46.8

Re: [gdal-dev] Add Mercator_variant_A method?

2015-04-13 Thread Andre Vautour
On 2015-04-13 11:01, Pepijn Van Eeckhoudt wrote: I am not sure I follow what this has to do with the projection name in WKT. The newer WKT specification (ISO 19162) relies solely on the identifier to determine which operation method to use to project coordinates, but the older WKT tends to f

Re: [gdal-dev] Add Mercator_variant_A method?

2015-04-13 Thread Andre Vautour
On 2015-04-13 10:42, Even Rouault wrote: André, I am not sure I follow what this has to do with the projection name in WKT. The newer WKT specification (ISO 19162) relies solely on the identifier to determine which operation method to use to project coordinates, Do you know if there's a publi

Re: [gdal-dev] Add Mercator_variant_A method?

2015-04-13 Thread Andre Vautour
On 2015-04-12 19:51, Even Rouault wrote: Le samedi 11 avril 2015 15:19:48, Brad Hards a écrit : On Tue, 7 Apr 2015 07:12:20 AM Jukka Rahkonen wrote: I felt that document somehow heavy to read but it appears in many places that "Mercator variant A" is the primary name since October 2010. This s

Re: [gdal-dev] UFO format / GDAL 3.0

2015-04-01 Thread Andre Vautour
Did you happen to look at the date on your calendar this morning? :p On 2015-04-01 11:17, Daniel Morissette wrote: Good morning Even, That sounds like an ambitious world (universe?) domination plan, but before going too deep into the implementation details, can you comment on the potential fo

Re: [gdal-dev] Change the min and max value of raster image in gdal

2015-02-04 Thread Andre Joost
Am 04.02.2015 um 02:23 schrieb user_name: How will I do that in gdal? http://www.gdal.org/gdal_calc.html should be the right tool. HTH, André Joost ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-d

Re: [gdal-dev] Change the min and max value of raster image in gdal

2015-02-03 Thread Andre Joost
Am 03.02.2015 um 09:52 schrieb user_name: I want to change the min and max value of my raster images since its min and max values are always negative(e,g. -1, -1; -3.2, -3.2). What if I want to change that into -1 as the min value and 1 as the max value? Any idea? Thanks! Add 2.1 to all your va

Re: [gdal-dev] gdalwarp in MapTiler Pro Command Line

2015-01-29 Thread Andre Joost
Am 29.01.2015 um 19:13 schrieb GeospatialInfoTechSolutions: It's not working in OSGEO4W shell when I try to run it. gdalwarp -overwrite -r bilinear -q -cutline S:/cropfiles/ENR_L01.shp -crop_to_cutline -of GTiff S:/GEO-GIS/1501/d-EC_IFR/SOURCE/LOW/ENR_L01.TIF s:/geo-gis/1501/d-EC_IFR/decollar

Re: [gdal-dev] Closed polyline from DXF as line or area?

2015-01-26 Thread Andre Joost
Am 26.01.2015 um 16:36 schrieb Even Rouault: Le lundi 26 janvier 2015 16:25:04, Andre Joost a écrit : BTW: creating polygons from closed lines is still not available in GDAL or other osgeo command line tools? <http://osgeo-org.1560.x6.nabble.com/gdal-dev-MULTILINESTRING-to-MULTIPOLYG

Re: [gdal-dev] Closed polyline from DXF as line or area?

2015-01-26 Thread Andre Joost
Am 26.01.2015 um 15:29 schrieb Even Rouault: It seems the documentation matches the initial implementation but has been outdated by later changes. Indeed, ogr/ogrsf_frmts/dxf/ogrdxflayer.cpp contains the following (disabled) code : The behaviour might have changed with the implementation of i

Re: [gdal-dev] Garmin .img data?

2014-12-19 Thread Andre Joost
Am 19.12.2014 um 17:39 schrieb Paolo Cavallini: -BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi Even, Il 19/12/2014 17:01, Even Rouault ha scritto: Probably perhaps nobody has written it. The link you gave doesn't link much information about the format. Is it the vector format described at ht

Re: [gdal-dev] Can't get gdal_rasterize to work... with anything.

2014-12-14 Thread Andre Joost
I guess you are getting an empty image because you have not told gdal which field includes the height data. I was successful with gdal_rasterize -a HAEDDYPIME -ts 5000 3000 -l is50v_haedarlinur_24122013 D:/Download/NEU/IS50V_HAEDARGOGN_24122013_ISN2004/IS50V_HAEDARGOGN_SHP/is50v_haedarlinur_241

Re: [gdal-dev] Can't get gdal_rasterize to work... with anything.

2014-12-14 Thread Andre Joost
Am 14.12.2014 um 01:02 schrieb Karen Pease: And there's also the whole different dataset HAEDARLIKAN_LMI which contains an img and a 1 gig ige file, but it looks like they're not supported by gdal. You are right that the .ige file is not supported, but take the .img and you get the DEM data

Re: [gdal-dev] shp - inconsistent extent

2014-12-09 Thread Andre Joost
Am 09.12.2014 um 10:33 schrieb Martin Landa: Broken shapefile can be found [1]. [1] http://geo1.fsv.cvut.cz/landa/parcely.zip The file gets loaded into QGIS without complaints. Does ARCGIS like the result of an ogr2ogr transformation of the file? Greetings, André Joost ___

[gdal-dev] Original FID for SQLite dialect

2014-11-14 Thread Andre Vautour
Hi all, I was wondering if there is a way to get at the original feature id when executing an SQL query for the SQLite dialect. Something similar to what is described in http://www.gdal.org/ogr_sql.html for FID for the OGR dialect. That is, is there something equivalent to the following for

Re: [gdal-dev] Sample script idea: gdal_srsguess

2014-11-03 Thread Andre Joost
Am 03.11.2014 um 08:46 schrieb Paolo Cavallini: Thanks Aaron for this. IMHO an even better option would be to explicitly add to EPSG/proj.4 db the bounding boxes of validity for each projection: in this way, client GIS could show only valid projections for a given map area, and this would be a

Re: [gdal-dev] Ecw format

2014-10-25 Thread Andre Joost
Am 25.10.2014 um 11:05 schrieb Marco Casella: How Can I transform some tiff file in ecw format? Any tips is well appreciated. You have to buy a license from ERDAS to do that. Or you may find an old ECW driver with SDK v3.3 which allowed a limited filesize for file freation as well. The cur

Re: [gdal-dev] Transform pixel coordinates

2014-10-25 Thread Andre Joost
Am 25.10.2014 um 09:29 schrieb Artur Bercik: sorry, yes the pixel size are in degrees,0.047 deg.(5 km approx. at equator), in GCS and WGS84. Ok, then it is easy. Using the values from your other post: upperLeftX, upperLeftY = 150.1234, 50.1234 (They are in GCS with WGS84 datum) pixelWidth,

Re: [gdal-dev] Transform pixel coordinates

2014-10-25 Thread Andre Joost
Am 25.10.2014 um 07:05 schrieb Artur Bercik: On Sat, Oct 25, 2014 at 2:00 PM, Andre Joost wrote: Am 25.10.2014 um 06:03 schrieb Artur Bercik: Dear GDAL users I have a raster image of 406 lines and 270 samples with 5km pixel size. The center coordinates (lats and lons in degrees in GCS) of

Re: [gdal-dev] Transform pixel coordinates

2014-10-24 Thread Andre Joost
Am 25.10.2014 um 06:03 schrieb Artur Bercik: Dear GDAL users I have a raster image of 406 lines and 270 samples with 5km pixel size. The center coordinates (lats and lons in degrees in GCS) of upper left pixel, and the center coordinates (lats and lons in degrees in GCS) of the lower right pixel

Re: [gdal-dev] Conversion from ASCII to GeoTIFF with gdalwarp: incorrect output

2014-10-13 Thread Andre Joost
Am 12.10.2014 22:43, schrieb Lukasz Tracewski: Hi, I have a few thousands of ESRI ASCII rasters in World_Mollweide (EPSG:54009) projection, which I need to convert to GeoTIFF. Usually it goes fine, but there are certain special cases that resist peaceful conversion. They are all placed on Fiji, w

Re: [gdal-dev] 16bit Pixel Arrangement?

2014-10-10 Thread Andre Joost
Am 10.10.2014 23:27, schrieb Even Rouault: On Friday 10 October 2014 01:12:14 M Lewis wrote: QGIS seems to be able to display the image correctly. As grey level I guess ? And likely after adjusting the min/max values QGIS has support for pseudocolurs from single band. Maybe the author

Re: [gdal-dev] Binary Predicates in SQLite SQL dialect

2014-10-01 Thread Andre Vautour
On 2014-09-30 17:08, Even Rouault wrote: Le mardi 30 septembre 2014 20:20:14, Andre Vautour a écrit : Hi all, I am trying to use a geometry binary predicate (ST_Intersects) in ExecuteSQL() with the SQLite dialect and keep getting back an "ST_Intersect function does not exist" error

[gdal-dev] Binary Predicates in SQLite SQL dialect

2014-09-30 Thread Andre Vautour
Hi all, I am trying to use a geometry binary predicate (ST_Intersects) in ExecuteSQL() with the SQLite dialect and keep getting back an "ST_Intersect function does not exist" error back from SQLite. GDAL is built with SpatiaLite (HAVE_SPATIALITE is defined), but I believe the problem is that

Re: [gdal-dev] projection of (maybe broken) grib file

2014-09-26 Thread Andre Joost
Am 26.09.2014 13:51, schrieb teeschke: Hi list We have a grib file, that gets cropped after a projection from WGS84 to WebMercator. The issue occurs in UMN Mapserver and QGIS. ArcMap can project the grib file correctly - so I think its a GDAL thing.

Re: [gdal-dev] ERROR 1: Too many points (440 out of 441) failed to transform, unable to compute output bounds.

2014-09-17 Thread Andre Joost
Am 17.09.2014 04:59, schrieb Love: Hi, Thanks! But it is workig only for some while for some hdf's, error still appears. I assigned one uniform target extent for all the files since this was what I specified before I downloaded the files from the oceancolor nasa. The extent of the hdf file(s)

Re: [gdal-dev] ERROR 1: Too many points (440 out of 441) failed to transform, unable to compute output bounds.

2014-09-15 Thread Andre Joost
Am 15.09.2014 17:02, schrieb Simon Shak: I have had this error message show up if the s_srs is not correctly auto detected or not present in the input file. I have used a combination of looking through other nonstandard meta data files present with the image sets or using a program like global ma

Re: [gdal-dev] ERROR 1: Too many points (440 out of 441) failed to transform, unable to compute output bounds.

2014-09-15 Thread Andre Joost
Am 15.09.2014 04:22, schrieb Love: Thanks andre, will try this as soon as I can and will give you a update. I just want to ask if will it be possible to batch process my hdf files and gdalwarp everything if I'm going to specify the -te of each files? I am not sure if the -te coordi

Re: [gdal-dev] Extract subdataset from netCDF in a given directory

2014-09-12 Thread Andre Joost
Am 11.09.2014 16:58, schrieb Andre Joost: According to http://www.gdal.org/frmt_hdf5.html, the syntax should be: HDF5:file_name:subdataset So it has little use to put the path before the HDF5:. I'm not sure if GDAL accepts the path before the file_name. Running the command in the same f

Re: [gdal-dev] ERROR 1: Too many points (440 out of 441) failed to transform, unable to compute output bounds.

2014-09-12 Thread Andre Joost
Am 12.09.2014 04:54, schrieb Love: Hi, The gdalinfo in my post is the gdalinfo of the subdataset 37. I have also tried gdal_translate but when I tried to use gdalwarp using the gdal_translate output, it has error. I really need to use gdalwarp so that I could project the image into epsg:4326(+pr

Re: [gdal-dev] Extract subdataset from netCDF in a given directory

2014-09-11 Thread Andre Joost
Am 10.09.2014 20:02, schrieb VictoriaH: Hello, I am using gdal_translate to extract a subdataset from a netCDF file to a geotiff. If I run the command from the same directory where the netCDF file is located, there is no problem: Filename: PM25_2001_02_average.nc Ex. gdal_translate -of GTiff HD

Re: [gdal-dev] ERROR 1: Too many points (440 out of 441) failed to transform, unable to compute output bounds.

2014-09-11 Thread Andre Joost
Am 11.09.2014 10:45, schrieb Love: I know I'm not the first to post this topic and encountered this problem. But I've search and read some posts related to this and tried their solutions but no luck. I really can't figure out what is causing this problem. How am I going to solve it? Badly need he

Re: [gdal-dev] Extract subdataset from netCDF in a given directory

2014-09-11 Thread Andre Joost
Am 10.09.2014 20:02, schrieb VictoriaH: Hello, I am using gdal_translate to extract a subdataset from a netCDF file to a geotiff. If I run the command from the same directory where the netCDF file is located, there is no problem: Filename: PM25_2001_02_average.nc Ex. gdal_translate -of GTiff HD

[gdal-dev] OGRStyleTool::GetStyleString() loosing quotes after modification

2014-09-04 Thread Andre Vautour
I ran into a problem where a quoted text string in a OGR Style String gets lost while interacting with it via an OGRStyleTool. Using the same style string afterwords causes OGRStyleTool::Parse() to fail because the quote string contained parentheses which would have typically been ignored becau

Re: [gdal-dev] Problem encountered in stitching/merging the raster files

2014-08-29 Thread Andre Joost
Am 29.08.2014 09:04, schrieb user_name: I've tried your suggestion but then the result was this: Note: I used OpenEV to open the file. Perhaps you should switch to QGIS. OpenEV is quite old, and may not work with all output of GDAL. The grey parts between the images should have NODATA value.

  1   2   >