Re: [gdal-dev] No Spatial Reference System from UTM Ozi .map file
Le 13/11/2012 01:03, Nik Sands a écrit : Thanks for the replies so far. Unfortunately, adding a 3rd calibration point to the .map file has not helped. Of course I'm assuming that the added calibration point is valid (eg, I tried the one suggested below, which looks like it should be OK). Does anybody have any other suggestions as to why this .map file cannot be interpreted correctly by GDAL? Is it a bug, or an incomplete .map parsing implementation, or am I doing something wrong? Hi Nik, I have looked at the Ozi driver source code, and found that the UTM projection is not supported. Please file a bug (https://trac.osgeo.org/gdal/), provide some sample files, and I'll try to improve the driver. Jean-Claude ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Best way to mosaic and compress lots of tiff Black/white?
Hello all I have ~1000 tiff/tfw RGB but with only 0 or 255 color value in all bands (black/white). I want to mosaic them in a single image to be used in qgis and other gis software. Which is the best way to do this? And to compress the resulting image? May be Tiff?? Single band?? Another problem is that I must apply a formula to input images to transform 255 value to null, because of some white portion of each image that overlap to black portion of the adiacent one.. Thanks Pietro -- View this message in context: http://osgeo-org.1560.n6.nabble.com/Best-way-to-mosaic-and-compress-lots-of-tiff-Black-white-tp5016021.html Sent from the GDAL - Dev mailing list archive at Nabble.com. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Use of GDALTransformerFunc
Hi, I´ m using Gdal 1.9 and trying to realize a polynomial transformation using GDALTransformerFunc. I wanna use GDALTransformerFunchttp://www.gdal.org/gdal__alg_8h.html#a7df61123fec15deb3da3acabce19e647to perform the transformation. The documentation says that the first argument (void *pTransformerArg) is a callback data function. GDALCreateGCPTransformer()http://www.gdal.org/gdal__alg_8h.html#ac57c80f7bf6752fac1dff8c2a21a15a3 serves as callback data function for GDALTransformerFunc? Is there any way to obtain the coefficients of the polynomial transformation? Regards, --- Gabriel Fusca SUR Emprendimientos Tecnológicos Perú 345 Piso 5to Oficina B (C1067AAG) Ciudad de Buenos Aires, Argentina Tel. +54 (11) 4342-2976/84 gabrielfu...@suremptec.com.ar http://www.suremptec.com/ ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] SQLite Driver and non-spatial tables
I have a Spatialite enabled SQLite file with both spatial and non-spatial tables. Since my database has a table called geometry_columns, all non-spatial tables are not loaded as layers. If I don't have a geometry_columns OGR attempts to load geometry from all the tables, but still loads non-spatial tables. Is that correct? The behavior doesn't seem consistent in the two cases. I just want to make sure I understand the driver's behavior. kss ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Create overview tiff from a set of tiff files
I have a set f tiff files at different zoom levels/resolutions 0.tif (2048x1024) 1.tif (4096x2048) 2.tif (8192x4096) I would like to create one tif that is 8192x4096 but have 2 overviews that are the 0.tif and 1.tif. Is there a GDAL command that I can give multiple tiffs to create one tiff with overviews? Thanks ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Create overview tiff from a set of tiff files
Selon Billy Newman newman...@gmail.com: I have a set f tiff files at different zoom levels/resolutions 0.tif (2048x1024) 1.tif (4096x2048) 2.tif (8192x4096) I would like to create one tif that is 8192x4096 but have 2 overviews that are the 0.tif and 1.tif. Is there a GDAL command that I can give multiple tiffs to create one tiff with overviews? Rename 1.tif as 2.tif.ovr, and 0.tif as 2.tif.ovr.ovr then gdal_translate 2.tif 2_with_ovr.tif -co COPY_SRC_OVERVIEWS=YES ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] No Spatial Reference System from UTM Ozi .map file
Thanks very much for looking into this. I've filed the bug at https://trac.osgeo.org/gdal/ticket/4894 On 13/11/2012, at 10:02 PM, Jean-Claude Repetto jrepe...@free.fr wrote: Le 13/11/2012 01:03, Nik Sands a écrit : Thanks for the replies so far. Unfortunately, adding a 3rd calibration point to the .map file has not helped. Of course I'm assuming that the added calibration point is valid (eg, I tried the one suggested below, which looks like it should be OK). Does anybody have any other suggestions as to why this .map file cannot be interpreted correctly by GDAL? Is it a bug, or an incomplete .map parsing implementation, or am I doing something wrong? Hi Nik, I have looked at the Ozi driver source code, and found that the UTM projection is not supported. Please file a bug (https://trac.osgeo.org/gdal/), provide some sample files, and I'll try to improve the driver. Jean-Claude ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] reboxing, or 3D regridding
summary: I'd appreciate advice regarding tools and methods for transforming values from an unprojected global 3D spatial grid to a projected 3D grid with different horizontal and vertical resolution, or pointers to other resources to consult. details: I have output from a global atmospheric model that I'd like to use as initial and boundary conditions for a regional model. The global input netCDF has dimensions=2.5° lon x 1.875° lat x 56 vertical levels. The regional model runs over North America using a 12-km grid projected Lambert Comformal Conic (LCC), with 34 vertical levels. Since its top height is less than that of the global input, the extents of the output domain are fully contained within the input domain. Each box or voxel defined by the global input grid contains an estimate for the N2O concentration for that volume. From those I want to compute the concentrations for each output gridbox volume. I'd appreciate your recommendations for tools that can handle this usecase. Particularly, is this something that gdalwarp can do? FWIW, the best tool I've seen so far is R package=gstat, but (IIUC) - gstat expects projected input. I'm not sure if I can work around that for this usecase. - as the name implies, 'gstat' is doing geostatistics, e.g., variogram- and covariance-based modeling. I'm not sure either how to setup the distance weighting for my scenario, and, frankly, I remain unconvinced that a statistical approach is necessary for this application (though it may be a sufficient or the best-available approach). This may be due to my statistical ignorance, however. your assistance is appreciated, Tom Roche tom_ro...@pobox.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Inifitine recursion in ogr.py
Hi! I am using GDAL 1.9.2 x64, Python 3.3 x64 and Windows 7 and I have a problem with infinite recursion in ogr.py. My initial call is: feature = ogr.Feature(layer.GetLayerDefn()) and the call stack: __init__ Python\lib\site-packages\osgeo\ogr.py __getattr__ Python\lib\site-packages\osgeo\ogr.py GetField Python\lib\site-packages\osgeo\ogr.py GetFieldIndex Python\lib\site-packages\osgeo\ogr.py __getattr__ Python\lib\site-packages\osgeo\ogr.py GetField Python\lib\site-packages\osgeo\ogr.py ... repeating... ... I have tried using Eclipse 4.2 with Pydev 2.7 and PyScripter 2.5.3, and in both cases the debug functionality is messed up after the above initial call. When just running the script, no problems occur. What can I do to avoid this problem? Best regards Casper ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev