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
swaps the coordinates, returning a 536x1440 raster. Is there any NCO or GDAL option to sanitize the axes? Should I open a ticket for it? Thanks in advance, Andre Joost ___ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman

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

2018-04-25 Thread Andre Joost
to C:\OSGeo4W64\share\proj and see if that changes anything. You can use the file from your current 4.9.3 installation. That does not improve the situation with gdalwarp. I will do some tatest with cs2cs later. Greetings, Andre Joost ___ gdal-dev

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/missing_

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

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["BTM",GEO

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

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] 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] "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] 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] 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] 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 ___

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] 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 coordinates follo

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

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.

Re: [gdal-dev] Issues building ECW/JP2ECW driver as a plugin on Ubuntu

2014-08-21 Thread Andre Joost
Am 21.08.2014 07:55, schrieb Luke: I cross-posted to the UbuntuGIS mailing list (http://article.gmane.org/gmane.comp.gis.osgeo.ubuntu/941), but have had no response. Unfortunately, yes. There seems to be little activity on the list recently. Are there any other drivers in the GDAL tree that

Re: [gdal-dev] Reveal the meaning of (rw+vs) in driver capabilities

2014-08-18 Thread Andre Joost
Am 18.08.2014 12:20, schrieb Jukka Rahkonen: Hi, Gdalinfo --formats prints a list of formats and capabilities of each driver with abbreviations like BAG -raster- (ro) ECW -raster- (rov) JP2ECW -raster- (rov) FITS -raster- (rw+) GMT -raster- (rw) gdalinfo --format [format name] e.g. gdalinf

Re: [gdal-dev] Issues building ECW/JP2ECW driver as a plugin on Ubuntu

2014-08-13 Thread Andre Joost
Am 14.08.2014 02:56, schrieb Luke Pinner: Can anybody provide some guidance on how to compile the ECW/JP2ECW drivers as a plugin on Linux? I have written down my workflow here: http://gis.stackexchange.com/questions/94870/unable-to-install-ecw-support-on-lubuntu-14-04 The plugin works in

Re: [gdal-dev] Lambert 93 problem with GDAL 1.11.0

2014-07-18 Thread Andre Joost
Am 18.07.2014 08:36, schrieb orel: Thanks for spotted this, but it doesn't change anything. I've uploaded the file here : http://speedy.sh/vNEcQ/test-bati-pts2.zip If I display it, using GlobalMapper for example (or my osgEarth based application), th

Re: [gdal-dev] Lambert 93 problem with GDAL 1.11.0

2014-07-17 Thread Andre Joost
Am 17.07.2014 16:31, schrieb Frédéric Trastour: Hi, It seems that in your wkt standard_parallel_1 et standard_parallel_2 values are inverted. Take a look at : http://spatialreference.org/ref/epsg/2154/prettywkt/ That should not do any harm. Look at ESRI:102110: http://spatialreference.org/re

Re: [gdal-dev] gdal_translate 1.11 netcdf output is a flipped Y image

2014-07-01 Thread Andre Joost
Am 01.07.2014 14:51, schrieb 00darki: *gdal_translate.exe -a_ullr "-180" "90" "180" "-90" -a_srs "EPSG:3408" D:/GDAL/MapServer/mnt0.nc D:/GDAL/MapServer/mnt00_geotif00.tif I guess that the bounds must meet the CRS of the tif. Since that is not lat/lon, the bounds make no sense to me. HT

Re: [gdal-dev] gdal 1.11 unable to read netcdf files

2014-06-27 Thread Andre Joost
Am 27.06.2014 20:00, schrieb Etienne Tourigny: If these plugin files came from gisinternals and netcdf is not supported without moving the plugin files to another directory, then that is a bug in gisinternals. Perhaps Tamas can comment on this? Gdalinfo works for me, with GDAL 1.11.0 from Tama

Re: [gdal-dev] gdal2tiles from vrt with unexpected result

2014-06-18 Thread Andre Joost
Am 17.06.2014 17:30, schrieb klo: VRT extend is approximately this: ulx, uly = (7436881, 4707326) lrx, lry = (7687838, 4511902) Not sure if it helps you, but I just solved a similar issue here:

Re: [gdal-dev] gdal2tiles from vrt with unexpected result

2014-06-16 Thread Andre Joost
unnecessarily 4 times while preserving 256px canvas. Any ideas, what might be wrong? What extent does your source vrt have? Greetings, Andre Joost ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] Question to ogr2ogr

2014-06-03 Thread Andre Joost
Am 03.06.2014 13:44, schrieb Jukka Rahkonen: ogr2ogr -f "ESRI Shapefile" -s_srs "+proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=450 +y_0=0 +datum=potsdam +nadgrids=BETA2007.gsb +units=m +no_defs" -t_srs epsg:25832 output.shp input.shp the +datum=potsdam can be misleading, this includes fixed +

Re: [gdal-dev] issues with gdal2tiles.py

2014-05-29 Thread Andre Joost
Am 29.05.2014 17:23, schrieb Marcela Campo: Sure, this is the output of gdalinfo -nomd -norat on the original img file Origin = (4392504.764280953444541,7439406.971977045759559) Pixel Size = (30.015846093058180,-30.015846093058180) Corner Coordinates: Upper Left ( 4392504.764, 7439406.972) ( 64

Re: [gdal-dev] issues with gdal2tiles.py

2014-05-29 Thread Andre Joost
Am 27.05.2014 18:45, schrieb Marcela Campo: I am not sure about source projection, I found it out doing: gdalinfo -proj4 input.img Can you add the full output of gdalinfo input.img (before any reprojection)? Greetings, André Joost ___ gdal-dev m

Re: [gdal-dev] issues with gdal2tiles.py

2014-05-27 Thread Andre Joost
Am 27.05.2014 18:45, schrieb Marcela Campo: Hi André, thank you for your response. I am not using EPSG:900913 anywhere You mentioned an error message about EPSG:900913 in your first post. I did my gdal2tiling a few years ago, and at that time 900913 was still listed in proj. Now it has been r

Re: [gdal-dev] issues with gdal2tiles.py

2014-05-27 Thread Andre Joost
Am 27.05.2014 05:31, schrieb Marcela Campo: Hi, I am trying to use gdal2tiles with a set of geotiff and I get the following error ERROR 6: EPSG PCS/GCS code 900913 not found in EPSG support files. Is this a valid EPSG coordinate system? ERROR 6: No translation for an empty SRS to PROJ.4 form

Re: [gdal-dev] GDAL_DRIVER_PATH dll hell

2014-05-15 Thread Andre Joost
Am 15.05.2014 12:07, schrieb Dmitriy Baryshnikov: I found out that another program used the gdal 1.8 was installed on that PC. And this program create env variable GDAL_DRIVER_PATH which pointed on some non standard driver linked with gdal 1.8 I think some version checking needed or something

Re: [gdal-dev] Update of EPSG database to v8.4

2014-05-15 Thread Andre Joost
Am 14.05.2014 12:53, schrieb Even Rouault: Regarding PROJ.4 'espg' and PostGIS spatial_ref_sys.sql files, I've also added the list of SRS that are GEOCCS (GeoCentric) and COMPD_CS (Compound Horizontal + Vertical). Please let me know if you see issues. The GEOCGS will likely cause a lot of c

Re: [gdal-dev] OSM driver to create separate tables for tags

2014-05-15 Thread Andre Joost
Am 14.05.2014 23:56, schrieb Jukka Rahkonen: Read the driver manual http://www.gdal.org/ogr/drv_osm.html Multilinestrings are mostly route relations. Line members of these relations go also into line layer. Ok, I missed that one. I do not see any reason why it should be much if any slower

Re: [gdal-dev] OSM driver to create separate tables for tags

2014-05-14 Thread Andre Joost
Am 14.05.2014 16:41, schrieb Jukka Rahkonen: I wonder if it would make any sense in the GDAL context to have an option in the OSM driver to separate geometries and tags and write them into their own tables. The structure of the resulting tables (5 + 5 for points, lines, multilines, multipolygons

Re: [gdal-dev] how to know which proj4 version gdal uses internally ?

2014-05-02 Thread Andre Joost
Am 02.05.2014 13:23, schrieb Yasin Selçuk Berber: how can i know which proj4 version uses gdal v1.9.2 internally ? Open a command window, and type proj This will return the currently installed version number. HTH, André Joost ___ gdal-dev maili

Re: [gdal-dev] MODIS re-projection issues (gdalwarp -tps)

2014-04-16 Thread Andre Joost
Am 16.04.2014 16:27, schrieb Giuseppe Amatulli: Is there is a way to solve this problem? I think that, cut the image considering the -180 +180 transition would solve the problem but has to be cut at GCP level and i do not have any clue how to do it. There are no GCP points in the file, but

Re: [gdal-dev] ogr2ogr ERROR 6: No translation

2014-04-10 Thread Andre Joost
Am 10.04.2014 17:13, schrieb Jason Mathis: Hi All, thanks for the responses. I believe my extents are: (-120.00591219477626, -114.03946520116759, 35.00247168585889, 42.00234483532862) This sounds reasonable for Nevada. $ ogr2ogr version is: GDAL 1.10.1, released 2013/08/26 I not sure how

Re: [gdal-dev] ogr2ogr ERROR 6: No translation

2014-04-09 Thread Andre Joost
Am 10.04.2014 04:31, schrieb Jason Mathis: Hi All, On to the next problem:) I discovered today that I need to reproject the shape file before rasterizing it. I get the below error. Any ideas? Is it not possible? What extent does your shapefile have? Some projections are limited to the visib

Re: [gdal-dev] Problem using gdal_contour after gdal_merge.py

2014-03-22 Thread Andre Joost
Am 22.03.2014 18:59, schrieb Tom Beddard: Are there any tricks to generating contours from merged files or is it just going to be a case of scripting the generation of lots of .shp contour files from each of the elevation tiles? Thanks for any help, I’m still on the steep bit of the learning c

  1   2   >