Re: [gdal-dev] gdalwarp for ED50 UTM zone 31N to WGS84

2009-05-13 Thread Jorge Arévalo 
De nada :-). Y suerte


2009/5/14 Ragi Y. Burhum 

> Gracias Jorge
>
> 2009/5/13 Jorge Arévalo 
>
>  Hello,
>>
>> I'm from Spain, and I can confirm your parameters. These worked for me as
>> well.
>>
>> Regards
>> Jorge
>>
>>
>> On Tue, May 12, 2009 at 1:03 AM, Ragi Y. Burhum  wrote:
>>
>>> I was just sending this e-mail so that it can be searched through the
>>> archives. Thanks to Even Rouault for walking me through something that
>>> should have been obvious :)
>>>
>>> Originally, I was having trouble finding the right parameters for
>>> gdalwarp to be able to warp some imagery that I got from Spain. These were
>>> in ED50 / UTM zone 31N:
>>>
>>> rbur...@cuzconix:~/p_bug$ gdalinfo exported_from_original_sid_file.tif
>>> Driver: GTiff/GeoTIFF
>>> Files: exported_from_original_sid_file.tif
>>> Size is 6946, 4800
>>> Coordinate System is:
>>> PROJCS["ED50 / UTM zone 31N",
>>> GEOGCS["ED50",
>>> DATUM["European_Datum_1950",
>>> SPHEROID["International 1924",6378388,297.00014,
>>> AUTHORITY["EPSG","7022"]],
>>> AUTHORITY["EPSG","6230"]],
>>> PRIMEM["Greenwich",0],
>>> UNIT["degree",0.0174532925199433],
>>> AUTHORITY["EPSG","4230"]],
>>> PROJECTION["Transverse_Mercator"],
>>> PARAMETER["latitude_of_origin",0],
>>> PARAMETER["central_meridian",3],
>>> PARAMETER["scale_factor",0.9996],
>>> PARAMETER["false_easting",50],
>>> PARAMETER["false_northing",0],
>>> UNIT["metre",1,
>>> AUTHORITY["EPSG","9001"]],
>>> AUTHORITY["EPSG","23031"]]
>>> Origin = (327720.000,4730783.500)
>>> Pixel Size = (0.500,-0.500)
>>> Metadata:
>>>   AREA_OR_POINT=Area
>>> Image Structure Metadata:
>>>   INTERLEAVE=PIXEL
>>> Corner Coordinates:
>>> Upper Left  (  327720.000, 4730783.500) (  0d53'46.94"E, 42d42'34.10"N)
>>> Lower Left  (  327720.000, 4728383.500) (  0d53'49.56"E, 42d41'16.35"N)
>>> Upper Right (  331193.000, 4730783.500) (  0d56'19.50"E, 42d42'36.88"N)
>>> Lower Right (  331193.000, 4728383.500) (  0d56'22.08"E, 42d41'19.12"N)
>>> Center  (  329456.500, 4729583.500) (  0d55'4.52"E, 42d41'56.62"N)
>>> Band 1 Block=6946x1 Type=Byte, ColorInterp=Red
>>> Band 2 Block=6946x1 Type=Byte, ColorInterp=Green
>>> Band 3 Block=6946x1 Type=Byte, ColorInterp=Blue
>>>
>>>
>>> In order to get it to work correctly, I had to use the right +towgs84
>>> parameters from
>>> http://earth-info.nga.mil/GandG/coordsys/onlinedatum/CountryEuropeTable.html.
>>>  Although the table showed  -84 (+/-5), -107(+/-6), -120 (+/-3), for my own
>>> purposes -87, -98, -121 seemed to work best based on my source data.
>>>
>>> so the command line would be:
>>>
>>>  gdalwarp -s_srs "+proj=utm +zone=31 +ellps=intl +units=m
>>> +towgs84=-87,-101,-121 +no_defs" -t_srs EPSG:4326
>>> exported_from_original_sid_file.tif testout.tif
>>>
>>> Hope this helps someone else,
>>>
>>> - Ragi
>>>
>>>
>>>
>>>
>>>
>>>
>>> ___
>>> 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

Re: [gdal-dev] gdalwarp for ED50 UTM zone 31N to WGS84

2009-05-13 Thread Ragi Y. Burhum 
Gracias Jorge

2009/5/13 Jorge Arévalo 

> Hello,
>
> I'm from Spain, and I can confirm your parameters. These worked for me as
> well.
>
> Regards
> Jorge
>
>
> On Tue, May 12, 2009 at 1:03 AM, Ragi Y. Burhum  wrote:
>
>> I was just sending this e-mail so that it can be searched through the
>> archives. Thanks to Even Rouault for walking me through something that
>> should have been obvious :)
>>
>> Originally, I was having trouble finding the right parameters for gdalwarp
>> to be able to warp some imagery that I got from Spain. These were in ED50 /
>> UTM zone 31N:
>>
>> rbur...@cuzconix:~/p_bug$ gdalinfo exported_from_original_sid_file.tif
>> Driver: GTiff/GeoTIFF
>> Files: exported_from_original_sid_file.tif
>> Size is 6946, 4800
>> Coordinate System is:
>> PROJCS["ED50 / UTM zone 31N",
>> GEOGCS["ED50",
>> DATUM["European_Datum_1950",
>> SPHEROID["International 1924",6378388,297.00014,
>> AUTHORITY["EPSG","7022"]],
>> AUTHORITY["EPSG","6230"]],
>> PRIMEM["Greenwich",0],
>> UNIT["degree",0.0174532925199433],
>> AUTHORITY["EPSG","4230"]],
>> PROJECTION["Transverse_Mercator"],
>> PARAMETER["latitude_of_origin",0],
>> PARAMETER["central_meridian",3],
>> PARAMETER["scale_factor",0.9996],
>> PARAMETER["false_easting",50],
>> PARAMETER["false_northing",0],
>> UNIT["metre",1,
>> AUTHORITY["EPSG","9001"]],
>> AUTHORITY["EPSG","23031"]]
>> Origin = (327720.000,4730783.500)
>> Pixel Size = (0.500,-0.500)
>> Metadata:
>>   AREA_OR_POINT=Area
>> Image Structure Metadata:
>>   INTERLEAVE=PIXEL
>> Corner Coordinates:
>> Upper Left  (  327720.000, 4730783.500) (  0d53'46.94"E, 42d42'34.10"N)
>> Lower Left  (  327720.000, 4728383.500) (  0d53'49.56"E, 42d41'16.35"N)
>> Upper Right (  331193.000, 4730783.500) (  0d56'19.50"E, 42d42'36.88"N)
>> Lower Right (  331193.000, 4728383.500) (  0d56'22.08"E, 42d41'19.12"N)
>> Center  (  329456.500, 4729583.500) (  0d55'4.52"E, 42d41'56.62"N)
>> Band 1 Block=6946x1 Type=Byte, ColorInterp=Red
>> Band 2 Block=6946x1 Type=Byte, ColorInterp=Green
>> Band 3 Block=6946x1 Type=Byte, ColorInterp=Blue
>>
>>
>> In order to get it to work correctly, I had to use the right +towgs84
>> parameters from
>> http://earth-info.nga.mil/GandG/coordsys/onlinedatum/CountryEuropeTable.html.
>>  Although the table showed  -84 (+/-5), -107(+/-6), -120 (+/-3), for my own
>> purposes -87, -98, -121 seemed to work best based on my source data.
>>
>> so the command line would be:
>>
>>  gdalwarp -s_srs "+proj=utm +zone=31 +ellps=intl +units=m
>> +towgs84=-87,-101,-121 +no_defs" -t_srs EPSG:4326
>> exported_from_original_sid_file.tif testout.tif
>>
>> Hope this helps someone else,
>>
>> - Ragi
>>
>>
>>
>>
>>
>>
>> ___
>> 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

Re: [gdal-dev] gdalwarp for ED50 UTM zone 31N to WGS84

2009-05-13 Thread Jorge Arévalo 
Hello,

I'm from Spain, and I can confirm your parameters. These worked for me as
well.

Regards
Jorge


On Tue, May 12, 2009 at 1:03 AM, Ragi Y. Burhum  wrote:

> I was just sending this e-mail so that it can be searched through the
> archives. Thanks to Even Rouault for walking me through something that
> should have been obvious :)
>
> Originally, I was having trouble finding the right parameters for gdalwarp
> to be able to warp some imagery that I got from Spain. These were in ED50 /
> UTM zone 31N:
>
> rbur...@cuzconix:~/p_bug$ gdalinfo exported_from_original_sid_file.tif
> Driver: GTiff/GeoTIFF
> Files: exported_from_original_sid_file.tif
> Size is 6946, 4800
> Coordinate System is:
> PROJCS["ED50 / UTM zone 31N",
> GEOGCS["ED50",
> DATUM["European_Datum_1950",
> SPHEROID["International 1924",6378388,297.00014,
> AUTHORITY["EPSG","7022"]],
> AUTHORITY["EPSG","6230"]],
> PRIMEM["Greenwich",0],
> UNIT["degree",0.0174532925199433],
> AUTHORITY["EPSG","4230"]],
> PROJECTION["Transverse_Mercator"],
> PARAMETER["latitude_of_origin",0],
> PARAMETER["central_meridian",3],
> PARAMETER["scale_factor",0.9996],
> PARAMETER["false_easting",50],
> PARAMETER["false_northing",0],
> UNIT["metre",1,
> AUTHORITY["EPSG","9001"]],
> AUTHORITY["EPSG","23031"]]
> Origin = (327720.000,4730783.500)
> Pixel Size = (0.500,-0.500)
> Metadata:
>   AREA_OR_POINT=Area
> Image Structure Metadata:
>   INTERLEAVE=PIXEL
> Corner Coordinates:
> Upper Left  (  327720.000, 4730783.500) (  0d53'46.94"E, 42d42'34.10"N)
> Lower Left  (  327720.000, 4728383.500) (  0d53'49.56"E, 42d41'16.35"N)
> Upper Right (  331193.000, 4730783.500) (  0d56'19.50"E, 42d42'36.88"N)
> Lower Right (  331193.000, 4728383.500) (  0d56'22.08"E, 42d41'19.12"N)
> Center  (  329456.500, 4729583.500) (  0d55'4.52"E, 42d41'56.62"N)
> Band 1 Block=6946x1 Type=Byte, ColorInterp=Red
> Band 2 Block=6946x1 Type=Byte, ColorInterp=Green
> Band 3 Block=6946x1 Type=Byte, ColorInterp=Blue
>
>
> In order to get it to work correctly, I had to use the right +towgs84
> parameters from
> http://earth-info.nga.mil/GandG/coordsys/onlinedatum/CountryEuropeTable.html. 
> Although the table showed  -84 (+/-5), -107(+/-6), -120 (+/-3), for my own
> purposes -87, -98, -121 seemed to work best based on my source data.
>
> so the command line would be:
>
>  gdalwarp -s_srs "+proj=utm +zone=31 +ellps=intl +units=m
> +towgs84=-87,-101,-121 +no_defs" -t_srs EPSG:4326
> exported_from_original_sid_file.tif testout.tif
>
> Hope this helps someone else,
>
> - Ragi
>
>
>
>
>
>
> ___
> 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

Re: [gdal-dev] gdalwarp for ED50 UTM zone 31N to WGS84

2009-05-11 Thread Jean-Claude REPETTO 

Ragi Y. Burhum a écrit :


In order to get it to work correctly, I had to use the right +towgs84 
parameters from 
http://earth-info.nga.mil/GandG/coordsys/onlinedatum/CountryEuropeTable.html 
. Although the table showed  -84 (+/-5), -107(+/-6), -120 (+/-3), for my 
own purposes -87, -98, -121 seemed to work best based on my source data.


Hello,

Other parameters are used in this PROJ.4 ML thread :


I don't know if they are better than yours, but they seem to be official.

Jean-Claude

___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

[gdal-dev] gdalwarp for ED50 UTM zone 31N to WGS84

2009-05-11 Thread Ragi Y. Burhum 
I was just sending this e-mail so that it can be searched through the
archives. Thanks to Even Rouault for walking me through something that
should have been obvious :)

Originally, I was having trouble finding the right parameters for gdalwarp
to be able to warp some imagery that I got from Spain. These were in ED50 /
UTM zone 31N:

rbur...@cuzconix:~/p_bug$ gdalinfo exported_from_original_sid_file.tif
Driver: GTiff/GeoTIFF
Files: exported_from_original_sid_file.tif
Size is 6946, 4800
Coordinate System is:
PROJCS["ED50 / UTM zone 31N",
GEOGCS["ED50",
DATUM["European_Datum_1950",
SPHEROID["International 1924",6378388,297.00014,
AUTHORITY["EPSG","7022"]],
AUTHORITY["EPSG","6230"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4230"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",3],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",50],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","23031"]]
Origin = (327720.000,4730783.500)
Pixel Size = (0.500,-0.500)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  327720.000, 4730783.500) (  0d53'46.94"E, 42d42'34.10"N)
Lower Left  (  327720.000, 4728383.500) (  0d53'49.56"E, 42d41'16.35"N)
Upper Right (  331193.000, 4730783.500) (  0d56'19.50"E, 42d42'36.88"N)
Lower Right (  331193.000, 4728383.500) (  0d56'22.08"E, 42d41'19.12"N)
Center  (  329456.500, 4729583.500) (  0d55'4.52"E, 42d41'56.62"N)
Band 1 Block=6946x1 Type=Byte, ColorInterp=Red
Band 2 Block=6946x1 Type=Byte, ColorInterp=Green
Band 3 Block=6946x1 Type=Byte, ColorInterp=Blue


In order to get it to work correctly, I had to use the right +towgs84
parameters from
http://earth-info.nga.mil/GandG/coordsys/onlinedatum/CountryEuropeTable.html.
Although the table showed  -84 (+/-5), -107(+/-6), -120 (+/-3), for my
own
purposes -87, -98, -121 seemed to work best based on my source data.

so the command line would be:

 gdalwarp -s_srs "+proj=utm +zone=31 +ellps=intl +units=m
+towgs84=-87,-101,-121 +no_defs" -t_srs EPSG:4326
exported_from_original_sid_file.tif testout.tif

Hope this helps someone else,

- Ragi
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev