Hi Andrew But How can I create those VRT, I never used those, from my original binaries datasets (HDF)? It seems that VRT is the key for my problem. But I never used it before
Thanks for your help Luisa 2010/2/18 Andrew Clegg <a...@pml.ac.uk> > Hi Luisa, > > I've tried something similar in the past with satellite data. I had 3 > files, containing latitude, longitude and value stored as floats. I had some > very minor success with the following solution (I got out an image but it > seemed to drop to a very low resolution). I'm posting it more in the hope > that someone else may be able to spot what I did wrong, because I would > still really like to be able to do this. > > My approach was to create a virtual dataset for each file: > > sst.vrt (the values): > > <VRTDataset rasterXSize="1296" rasterYSize="2030"> > <Metadata domain="GEOLOCATION"> > <MDI key="X_DATASET">lon.vrt</MDI> > <MDI key="X_BAND">1</MDI> > <MDI key="Y_DATASET">lat.vrt</MDI> > <MDI key="Y_BAND">1</MDI> > <MDI key="PIXEL_OFFSET">0</MDI> > <MDI key="LINE_OFFSET">0</MDI> > <MDI key="PIXEL_STEP">1</MDI> > <MDI key="LINE_STEP">1</MDI> > </Metadata> > <VRTRasterBand dataType="Float32" band="1" subClass="VRTRawRasterBand"> > <SourceFilename relativetoVRT="1">sst.img</SourceFilename> > <ImageOffset>0</ImageOffset> > <PixelOffset>4</PixelOffset> > <LineOffset>5184</LineOffset> > <ByteOrder>LSB</ByteOrder> > </VRTRasterBand> > </VRTDataset> > > lon.vrt (longitude): (latitude is basically the same) > > <VRTDataset rasterXSize="1296" rasterYSize="2030"> > <VRTRasterBand dataType="Float32" band="1" subClass="VRTRawRasterBand"> > <SourceFilename relativetoVRT="1">lon.img</SourceFilename> > <ImageOffset>0</ImageOffset> > <PixelOffset>4</PixelOffset> > <LineOffset>5184</LineOffset> > <ByteOrder>LSB</ByteOrder> > </VRTRasterBand> > </VRTDataset> > > I then tried to run gdalwarp on sst.vrt, trying either '-geoloc' or '-to > METHOD=GEOLOC_ARRAY'. > > Regards, > > Andrew Clegg > Plymouth Marine Laboratory > > Luisa Peña wrote: > >> Greetings >> I've used GDAL a few times but this is the first time that I'm lost on >> what I can do with GDAL to solve my problem. I have 3 separate files, one >> with my data, other with Lattitude and another with Longitude. They have the >> same size and, as an example, it means that Pixel (1,1) of data is located >> in latitude retrieved from Pixel (1,1) from latitude and longitude (pixel >> (1,1) from longitude). All of them are in WGS84 but without GCP's points. I >> need to create a georreferenced/rectified Geotiff with this data and, to do >> that, I need lat and Long. How can I do this? >> >> Thanks >> >> Best regards, >> Luisa Peña >> S.Paulo, Brazil >> ------------------------------------------------------------------------ >> >> >> _______________________________________________ >> gdal-dev mailing list >> gdal-dev@lists.osgeo.org >> http://lists.osgeo.org/mailman/listinfo/gdal-dev >> > > > > -------------------------------------------------------------------------------- > Plymouth Marine Laboratory > > Registered Office: > Prospect Place The Hoe > Plymouth PL1 3DH > > Website: www.pml.ac.uk > Registered Charity No. 1091222 > PML is a company limited by guarantee > registered in England & Wales > company number 4178503 > > PML is a member of the Plymouth Marine Sciences Partnership > Website: > www.pmsp.org.uk-------------------------------------------------------------------------------- > This e-mail, its content and any file attachments are confidential. > > If you have received this e-mail in error please do not copy, disclose it > to any third party or use the contents or attachments in any way. Please > notify the sender by replying to this e-mail or e-mail fori...@pml.ac.ukand > then delete the email without making any copies or using it in any other > way. > > The content of this message may contain personal views which are not the > views of Plymouth Marine Laboratory unless specifically stated. > > You are reminded that e-mail communications are not secure and may contain > viruses. Plymouth Marine Laboratory accepts no liability for any loss or > damage which may be caused by viruses. > > -------------------------------------------------------------------------------- >
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev