[gdal-dev] changing data type of a shapefile

2011-06-15 Thread ahmet temiz
hello how can I change data type of a shapefile by ogr2ogr ? regards -- Ahmet Temiz Jeoloji Müh. Afet ve Acil Durum Yönetimi Başkanlığı Planlama ve Zarar Azaltma Dairesi Başkanlığı Bilgi ve CBS grubu Eskişehir Yolu 10. km. Lodumlu / Ankara Tel : 0 312 2872680 / 1535

Re: [gdal-dev] disformed images in its positions

2011-06-15 Thread Chaitanya kumar CH
Frederick, You can use GCPs you used to create the polynomials directly in gdal utilities to correct the deformed images. I assume that your current images don't have any projection info. 1. First use gdal_translate with the -gcp option to create a dataset with your control points. You may also

[gdal-dev] get GCP from noaa L1b

2011-06-15 Thread Ramiro Gonzalez
Hi, I’m using gdal c++ API. I’d like to implement an algorithm to georeference noaa L1b images(of Argentine in particular). The number of GCPs returned by the L1b driver are not enough to use triangulation(midcourse solution until a better algorithm can be implemented). Searching gdal-dev

[gdal-dev] Re: disformed images in its positions

2011-06-15 Thread frederick keusch
Hi Chaitanya, first of all I'd like to thank you for the information. So if I understood you correctly, it is not necesseary to use the polynomial we calculated but only a dataset of the control points we used and gdal_transalate calculates a polynomial internally. Do you know how exactly this

Re: [gdal-dev] get GCP from noaa L1b

2011-06-15 Thread Andrew Brooks
Hi From memory the l1b driver doesn't use enough points to get proper geolocation because it sub-samples the 51 which are provided in the l1b file. The best solution (that recommended by NOAA) is to use lagrange interpolation across the 51 points. A lagrange algorithm has been provided in the

Re: [gdal-dev] changing data type of a shapefile

2011-06-15 Thread Chaitanya kumar CH
Ahmet, You can use OGR SQL [1] to typecast fields and other stuff. You can use ogr2ogr [2] with the -sql option and supply the sql statement. [1]: http://trac.osgeo.org/gdal/wiki/rfc21_ogrsqlcast [2]: http://www.gdal.org/ogr2ogr.html 2011/6/15 ahmet temiz ahmettemi...@gmail.com hello how

Re: [gdal-dev] CreateCopy() with same metadata, different raster data

2011-06-15 Thread Frank Warmerdam
On 11-06-14 09:18 PM, Cole, Derek wrote: Hrmm, that is sort of unfortunate. I have set up a tiling algorithm to read in tiles at a time of imagery data, which is already slow-ish, even doing it it with threads. Having to do a whole new read of the entire image at once (since it doesnt seem

Re: [gdal-dev] changing data type of a shapefile

2011-06-15 Thread Eli Adam
Ahmet, You can also change geometry type, it is also an option with ogr2ogr [2], -nlt. I've only used -nlt to promote to multi-geometry. Eli On 6/15/2011 at 7:27 AM, in message BANLkTik8UghqaRv8E-pES6OSv-+e=a3...@mail.gmail.com, Chaitanya kumar CH chaitanya...@gmail.com wrote: Ahmet,

[gdal-dev] SetUTM - how to say how far north Python

2011-06-15 Thread Chris English
Hi, Assuming the following: src_ds = gdal.Open('test_3.ter') dst_ds = driver.CreateCopy('test_6.hf2', src_ds, 0) dst_ds.SetGeoTransform([582495,1,0.5,4512717,0.5,-1])0 srs = osr.SpatialReference() srs.SetUTM(18,1)0 srs.SetWellKnownGeogCS('WGS84')0 dst_ds.SetProjection(

Re: [gdal-dev] Re: disformed images in its positions

2011-06-15 Thread Chaitanya kumar CH
Frederick, GCP is a control point. Ground Control Point. They are not visible in the image. Just a map of pixel-line to coordinates. I am not familiar with how the warping algorithm handles the GCPs. But if you want, you can also set the projection info after you use your own method to correct

[gdal-dev] gdal_translate and gdal_warp

2011-06-15 Thread Mike Willis
Hi Folks, Multipart question I am using gdal_translate to pull out files from JPEG2000 files wrapped by an NTF wrapper. It seems fine, but I am losing all metadata and no headers are built when I do this. I can get readable ENVI files or Geotiffs but they're often flipped back to front, or

[gdal-dev] gdal_translate compress is larger??

2011-06-15 Thread Matt Wilkie
Hi Folks, I've run into a strange thing, gdal_translate -co compress=lzw is creating files that are larger than the source uncompressed image. Input: 04/04/2006 08:39 AM 959,236,993 A1.tif Results of gdal_translate -co compress=[lzw,packbits,none] 15/06/2011 10:35 AM

RE: [gdal-dev] Building Overviews for Tiles of Large images

2011-06-15 Thread Cole, Derek
Bringing this back up, as I am curious. I have started playing around with using RasterIO and providing smaller output buffers for the read-in data, thus making GDAL do the downsampling (which apparently is nearest neighbor type?). Do you think this is faster than creating Overview imagery,

Re: [gdal-dev] gdal_translate compress is larger??

2011-06-15 Thread Adam Nowacki
Its just how compression works. If some data is reduced in size then other has to grow in size. Try -co compress=deflate -co PREDICTOR=2 On 2011-06-15 19:43, Matt Wilkie wrote: Hi Folks, I've run into a strange thing, gdal_translate -co compress=lzw is creating files that are larger than the

Re: [gdal-dev] gdal_translate compress is larger??

2011-06-15 Thread Frank Warmerdam
On 11-06-15 01:43 PM, Matt Wilkie wrote: Hi Folks, I've run into a strange thing, gdal_translate -co compress=lzw is creating files that are larger than the source uncompressed image. Input: 04/04/2006 08:39 AM 959,236,993 A1.tif Results of gdal_translate -co compress=[lzw,packbits,none]

[gdal-dev] hdf5image metadata problems

2011-06-15 Thread Alex Mantaut
Hi all, I'm trying to open the metadata from a subdataset of an HDF5 file... When I open the file by it's filename and then ask for the metadata it returns the metadata asociated to the root and all the subdatasets But when I open the file by specifying the subdataset name to

Re: [gdal-dev] hdf5image metadata problems

2011-06-15 Thread Frank Warmerdam
On 11-06-15 03:19 PM, Alex Mantaut wrote: Hi all, I'm trying to open the metadata from a subdataset of an HDF5 file... When I open the file by it's filename and then ask for the metadata it returns the metadata asociated to the root and all the subdatasets But when I open

Re: [gdal-dev] gdal_translate compress is larger??

2011-06-15 Thread Matt Wilkie
I'm curious as to why I've not encountered this before as I use lzw compression all the time. Anyway, as always thanks for the education. :) I've discovered that for my present project that jpeg-in-geotiff compression is good enough, which makes me happy because even with predictor=2 the

Re: [gdal-dev] hdf5image metadata problems

2011-06-15 Thread Alex Mantaut
Hi Frank: Thanks for the swift reply. Right now it opens the metadata asociated to the root, and also the metadata asociated with all the other groups, but it doesn't open the metadata asociated to the specified element (bug?) regards 2011/6/15 Frank Warmerdam

Re: [gdal-dev] hdf5image metadata problems

2011-06-15 Thread Frank Warmerdam
On 11-06-15 03:50 PM, Alex Mantaut wrote: Hi Frank: Thanks for the swift reply. Right now it opens the metadata asociated to the root, and also the metadata asociated with all the other groups, but it doesn't open the metadata asociated to the specified element (bug?)

Re: [gdal-dev] hdf5image metadata problems

2011-06-15 Thread Alex Mantaut
Frank: I opened a new defect ticket, and suggested a way around the problem (I'm not sure if it's the best solution) http://trac.osgeo.org/gdal/ticket/4121 Regards 2011/6/15 Frank Warmerdam warmer...@pobox.com On 11-06-15 03:50 PM, Alex Mantaut wrote: Hi Frank:

[gdal-dev] Re: disformed images in its positions

2011-06-15 Thread frederick keusch
Hi Chaitanya, thanks again for your help. I tried to find some information about how to set the projection info after our method and also how and with which command apply the modified GCP's... hoping to find a solution soon. Thank you again, frederick -- View this message in context:

RE: [gdal-dev] CreateCopy() with same metadata, different raster data

2011-06-15 Thread Cole, Derek
What about using the VRT driver? Is that for the situation I am describing? According to the CreateCopy documentation: It is intended that CreateCopy() will often be used with a source dataset which is a virtual dataset allowing configuration of band types, and other information without

Re: [gdal-dev] hdf5image metadata problems

2011-06-15 Thread Antonio Valentino
Hi Frank, hi Alex, Il 15/06/2011 22:34, Frank Warmerdam ha scritto: On 11-06-15 03:50 PM, Alex Mantaut wrote: Hi Frank: Thanks for the swift reply. Right now it opens the metadata asociated to the root, and also the metadata asociated with all the other groups,