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
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to