[gdal-dev] Re: how to say gdalwarp to use nadgrid by default

Eloi Ribeiro eloi.ribeiro at gmail.com
Mon Oct 18 05:26:30 EDT 2010


Thanks Chaitanya,

I intended to avoid all this writing to make sure the nadgrid is used:

*gdalwarp -s_srs "+proj=utm +zone=30 +ellps=intl +units=m +nadgrids=R2009V9.gsb
+wktext" -t_srs "+proj=longlat +ellps=GRS80 +towgs84=0,0,0" rast_23030.tif
rast_4258.tif*

But after reading this[1] I just realise that there is not such a
configuration file to search for in GDAL.

After all is not that big issue.

Thanks anyway.

By,

[1] http://www.osgeo.org/pipermail/gdal-dev/2007-April/012561.html

<http://www.osgeo.org/pipermail/gdal-dev/2007-April/012561.html>
Eloi Ribeiro
GIS Analyst
39,45º -4,40º
http://eloiribeiro.wordpress.com


On Sun, Oct 17, 2010 at 14:46, Chaitanya kumar CH <chaitanya.ch at gmail.com>wrote:

> Eloi,
>
> You can find the files at gdal/data/pcs.csv and gdal/data/gcs.csv
> You may also want to refer to http://www.gdal.org/ogr/osr_tutorial.html
>
> On Sun, Oct 17, 2010 at 5:18 PM, Eloi Ribeiro <eloi.ribeiro at gmail.com>wrote:
>
>> In other words,
>>
>> When I ask GDAL for:
>> *gdalwarp -s_srs EPSG:23030 -t_srs EPSG:4258 rast_23030.tif rast_4258.tif
>> *
>>
>> GDAL knows that 'EPSG:23030' means '+proj=utm +zone=30 +ellps=intl
>> +units=m +no_defs'.
>> So my question is: From witch file GDAL gets that information from?
>>
>> Cheers,
>>
>> Eloi Ribeiro
>> GIS Analyst
>> 39,45º -4,40º
>> http://eloiribeiro.wordpress.com
>>
>>
>> On Mon, Oct 11, 2010 at 14:55, Eloi Ribeiro <eloi.ribeiro at gmail.com>wrote:
>>
>>> Hi all,
>>>
>>> I would like that every time a I use gdalwarp for some EPSG codes [1]
>>> transformation it would use by default the nadgrids [2] (for Portugal and
>>> Spain) that I have under the folder '/usr/share/proj/'.
>>>
>>> For that purpose in PostGIS I needed to add
>>> '+nadgrids=nadgrid_file_name.gsb' at the end of the line of the 'proj4text'
>>> field text in table 'spatial_ref_sys' for each EPSG code I wonted to use the
>>> nadgrid.
>>>
>>> So I suppose there must be a similar way of doing this with GDAL. How can
>>> I achieve this with GDAL?
>>>
>>> Thank's!
>>>
>>> [1] 27492,27493,20790,20791,23029,23030,23031,3763,4258,4326
>>> [2] pt73_e89.gsb, ptLX_e89.gsb, ptLB_e89.gsb, ptED_e89.gsb, R2009V9.gsb,
>>> BALR2009.gsb
>>>
>>>
>>> Eloi Ribeiro
>>> GIS Analyst
>>> 39,45º -4,40º
>>> http://eloiribeiro.wordpress.com
>>>
>>
>>
>> _______________________________________________
>> gdal-dev mailing list
>> gdal-dev at lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>>
>
>
>
> --
> Best regards,
> Chaitanya kumar CH.
> /tʃaɪθənjə/ /kʊmɑr/
> +91-9494447584
> 17.2416N 80.1426E
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20101018/1cd6e1a8/attachment.html


More information about the gdal-dev mailing list