[gdal-dev] Projection aeqd and geotransform

[Paco Regodón]

Hello all,

i'm trying to create an VRT file for a raster in azimutal equidistant 
projection. I am almost sure the projection i should use is

+proj=aeqd +lat_0=YYYYY +lon_0=XXXXX +x_0=0 +y_0=0 +units=m

XXXXX and YYYYY being the raster's center coordinates.

I know the pixel width is 1000 meters, but i don't know how to write a 
geotransformation since i don't know the upper left pixel coordinates.

Anyone can help me? Can i define a geotransform with a central pixel 
instead of the upper left one?


<VRTDataset rasterXSize="480" rasterYSize="480">
<SRS>+proj=aeqd +lat_0=39.177500 +lon_0=0.249722 +x_0=0 +y_0=0 
+a=6371000 +b=6371000 +units=m</SRS>
<GeoTransform>??????, 1000, 0, ??????, 0, -1000</GeoTransform>
<VRTRasterBand dataType="Byte" band="1" subClass="VRTRawRasterBand">
<SourceFilename relativeToVRT="1">cullera.raw</SourceFilename>
</VRTRasterBand>
</VRTDataset>

Thanks,

Paco Regodón
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20101202/4377399d/attachment.html