<br>Hi Andrew<div><br></div><div>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</div><div><br></div><div>
Thanks for your help</div><div><br></div><div>Luisa<br><div class="gmail_quote">2010/2/18 Andrew Clegg <span dir="ltr">&lt;<a href="mailto:ancl@pml.ac.uk">ancl@pml.ac.uk</a>&gt;</span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi Luisa,<br>
<br>
I&#39;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&#39;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.<br>

<br>
My approach was to create a virtual dataset for each file:<br>
<br>
sst.vrt (the values):<br>
<br>
&lt;VRTDataset rasterXSize=&quot;1296&quot; rasterYSize=&quot;2030&quot;&gt;<br>
  &lt;Metadata domain=&quot;GEOLOCATION&quot;&gt;<br>
     &lt;MDI key=&quot;X_DATASET&quot;&gt;lon.vrt&lt;/MDI&gt;<br>
     &lt;MDI key=&quot;X_BAND&quot;&gt;1&lt;/MDI&gt;<br>
     &lt;MDI key=&quot;Y_DATASET&quot;&gt;lat.vrt&lt;/MDI&gt;<br>
     &lt;MDI key=&quot;Y_BAND&quot;&gt;1&lt;/MDI&gt;<br>
     &lt;MDI key=&quot;PIXEL_OFFSET&quot;&gt;0&lt;/MDI&gt;<br>
     &lt;MDI key=&quot;LINE_OFFSET&quot;&gt;0&lt;/MDI&gt;<br>
     &lt;MDI key=&quot;PIXEL_STEP&quot;&gt;1&lt;/MDI&gt;<br>
     &lt;MDI key=&quot;LINE_STEP&quot;&gt;1&lt;/MDI&gt;<br>
  &lt;/Metadata&gt;<br>
     &lt;VRTRasterBand dataType=&quot;Float32&quot; band=&quot;1&quot; subClass=&quot;VRTRawRasterBand&quot;&gt;<br>
    &lt;SourceFilename relativetoVRT=&quot;1&quot;&gt;sst.img&lt;/SourceFilename&gt;<br>
    &lt;ImageOffset&gt;0&lt;/ImageOffset&gt;<br>
    &lt;PixelOffset&gt;4&lt;/PixelOffset&gt;<br>
    &lt;LineOffset&gt;5184&lt;/LineOffset&gt;<br>
    &lt;ByteOrder&gt;LSB&lt;/ByteOrder&gt;<br>
  &lt;/VRTRasterBand&gt;<br>
&lt;/VRTDataset&gt;<br>
<br>
lon.vrt (longitude): (latitude is basically the same)<br>
<br>
&lt;VRTDataset rasterXSize=&quot;1296&quot; rasterYSize=&quot;2030&quot;&gt;<br>
  &lt;VRTRasterBand dataType=&quot;Float32&quot; band=&quot;1&quot; subClass=&quot;VRTRawRasterBand&quot;&gt;<br>
     &lt;SourceFilename relativetoVRT=&quot;1&quot;&gt;lon.img&lt;/SourceFilename&gt;<br>
     &lt;ImageOffset&gt;0&lt;/ImageOffset&gt;<br>
     &lt;PixelOffset&gt;4&lt;/PixelOffset&gt;<br>
     &lt;LineOffset&gt;5184&lt;/LineOffset&gt;<br>
     &lt;ByteOrder&gt;LSB&lt;/ByteOrder&gt;<br>
  &lt;/VRTRasterBand&gt;<br>
&lt;/VRTDataset&gt;<br>
<br>
I then tried to run gdalwarp on sst.vrt, trying either &#39;-geoloc&#39; or &#39;-to METHOD=GEOLOC_ARRAY&#39;.<br>
<br>
Regards,<br>
<br>
Andrew Clegg<br>
Plymouth Marine Laboratory<br>
<br>
Luisa Peņa wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div></div><div class="h5">
Greetings <br>
I&#39;ve used GDAL a few times but this is the first time that I&#39;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&#39;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?<br>

<br>
Thanks<br>
<br>
Best regards,<br>
Luisa Peņa<br>
S.Paulo, Brazil<br></div></div>
------------------------------------------------------------------------<div class="im"><br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</div></blockquote>
<br>
<br>
--------------------------------------------------------------------------------<br>
Plymouth Marine Laboratory<br>
<br>
Registered Office:<br>
Prospect Place The Hoe<br>
Plymouth  PL1 3DH<br>
<br>
Website: <a href="http://www.pml.ac.uk" target="_blank">www.pml.ac.uk</a><br>
Registered Charity No. 1091222<br>
PML is a company limited by guarantee<br>
registered in England &amp; Wales<br>
company number 4178503<br>
<br>
PML is a member of the Plymouth Marine Sciences Partnership<br>
Website: <a href="http://www.pmsp.org.uk" target="_blank">www.pmsp.org.uk</a> --------------------------------------------------------------------------------<br>
This e-mail, its content and any file attachments are confidential.<br>
<br>
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 <a href="mailto:forinfo@pml.ac.uk" target="_blank">forinfo@pml.ac.uk</a> and then delete the email without making any copies or using it in any other way.<br>

<br>
The content of this message may contain personal views which are not the views of Plymouth Marine Laboratory unless specifically stated.<br>
<br>
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.<br>
--------------------------------------------------------------------------------<br>
</blockquote></div><br></div>