<br><br><div class="gmail_quote">On 19 July 2012 15:25, Jose Gomez-Dans <span dir="ltr"><<a href="mailto:jgomezdans@gmail.com" target="_blank">jgomezdans@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<div><br><br><div class="gmail_quote"><div class="im">On 19 July 2012 04:46, sigologo <span dir="ltr"><<a href="mailto:dlopezaspe@gmail.com" target="_blank">dlopezaspe@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

I'm trying import A20030012003008.L3m_8D_NSST_4, with gdal with this script:<br>
<br>
for file in A*L3m_8D_NSST_4 ; do<br>
  echo "map: $file"<br>
<br>
  gdal_translate -a_srs "+init=epsg:4326" -a_nodata 65535 \<br>
    -a_ullr -180 90 180 -90 -co "COMPRESS=PACKBITS" \<br>
    HDF4_SDS:UNKNOWN:"$file":0  ${file}_out.tif<br>
<br>
done<br>
<br>
all good, but the data not be in degree C............. I find for it another<br>
script but no found in gdal...<br></blockquote><div><br></div></div><div>The problem is that the data are linearly scaled (from the GRASS script you quoted). I am unfamiliar with the ocean products (that's the bit of the data I blank out ;-D), and you know the linear scaling (either because it's constant, or beause you can read it out from the files's metadata), I think you can use gdal_translates 'scale' option: -scale [src_min src_max [dst_min dst_max]]. Note that hte syntax will be slightly different to the slope/intercept one in your GRASS script.</div>
</div></div></blockquote><div><br></div><div>Here is your script modified to take the linear correction into account. Note that the destination spread (-2 to 45) was read from the metadata. Also note that you need to make the output format Float, otherwise, you get a quantised version. </div>
<div><br></div><div><b><font face="courier new, monospace">for file in A*L3m_8D_NSST_4 ; </font></b></div><div><b><font face="courier new, monospace">do</font></b></div><div><b><font face="courier new, monospace">  echo "map: $file"</font></b></div>
<div><div><b><font face="courier new, monospace">  gdal_translate -of GTiff -ot Float32 \</font></b></div><div><b><font face="courier new, monospace">                      -scale 0 65535 -2. 45. \</font></b></div><div><b><font face="courier new, monospace">                      -a_srs "EPSG:4326" \</font></b></div>
<div><b><font face="courier new, monospace">                      -a_nodata 65535  \</font></b></div><div><b><font face="courier new, monospace">                      -a_ullr -180 90 180 -90 \</font></b></div><div><b><font face="courier new, monospace">                      -co "COMPRESS=PACKBITS" \</font></b></div>
<div><b><font face="courier new, monospace">                      'HDF4_SDS:UNKNOWN:"A20030012003008.L3m_8D_NSST_4":0' \</font></b></div><div><b><font face="courier new, monospace">                      A20030012003008.L3m_8D_NSST_4.tif</font></b></div>
</div><div><b><font face="courier new, monospace">done</font></b></div><div><br></div><div><br></div><div>Jose</div></div>