<div dir="ltr">Thanks Even!  Amazing, I don't think this is a rare case either, there's a lot of data in files "like this" that I keep encountering. <div><br></div><div>Another few lines in the NetCDF behemoth, I'll keep reading (and update my AI subscriptions). </div><div><br></div><div>Cheers, Mike</div><div><br></div></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">On Fri, Jul 4, 2025 at 8:21 PM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><u></u>

  
    
  
  <div>
    <p>Hi Michael,</p>
    <p>I guess this is finally the perfect use case we were waiting for
      to find a good reason to embed a Small Language Model within GDAL
      to figure out what are the columns and lines axis.  In the
      meantime, the boring <a href="https://github.com/OSGeo/gdal/pull/12700" target="_blank">https://github.com/OSGeo/gdal/pull/12700</a> will
      fix it</p>
    <p>Even<br>
    </p>
    <div>Le 04/07/2025 à 00:48, Michael Sumner
      via gdal-dev a écrit :<br>
    </div>
    <blockquote type="cite">
      
      <div dir="ltr">
        <div>I'm confused about the orientation of this in 2D mode
          1485x1272 is correct, with 122 bands. <br>
          <br>
          (requires earthdata auth)<br>
          <br>
          gdal raster info "/vsicurl/<a href="https://obdaac-tea.earthdatacloud.nasa.gov/ob-cumulus-prod-public/PACE_OCI.20250501T234515.L2.SFREFL.V3_0.nc" target="_blank">https://obdaac-tea.earthdatacloud.nasa.gov/ob-cumulus-prod-public/PACE_OCI.20250501T234515.L2.SFREFL.V3_0.nc</a>"
            --output-format text | grep /rhos<br>
          #  SUBDATASET_1_NAME=NETCDF:"/vsicurl/<a href="https://obdaac-tea.earthdatacloud.nasa.gov/ob-cumulus-prod-public/PACE_OCI.20250501T234515.L2.SFREFL.V3_0.nc" target="_blank">https://obdaac-tea.earthdatacloud.nasa.gov/ob-cumulus-prod-public/PACE_OCI.20250501T234515.L2.SFREFL.V3_0.nc</a>":/geophysical_data/rhos<br>
          #  SUBDATASET_1_DESC=[1485x1272x122] /geophysical_data/rhos
          (32-bit floating-point)<br>
          <br>
          If we drill into the rhos sds X becomes unrolled as the 1485
          bands: <br>
          <br>
          gdal raster info "vrt:///vsicurl/<a href="https://obdaac-tea.earthdatacloud.nasa.gov/ob-cumulus-prod-public/PACE_OCI.20250501T234515.L2.SFREFL.V3_0.nc?sd_name=/geophysical_data/rhos" target="_blank">https://obdaac-tea.earthdatacloud.nasa.gov/ob-cumulus-prod-public/PACE_OCI.20250501T234515.L2.SFREFL.V3_0.nc?sd_name=/geophysical_data/rhos</a>"
          --output-format text | grep "Size"<br>
          #Size is 122, 1272<br>
          <br>
          We can (redundantly??) cast via mdim to get 1485x1272, but
          then we've lost the geolocation arrays: <br>
          <br>
          gdal raster info "vrt:///vsicurl/<a href="https://obdaac-tea.earthdatacloud.nasa.gov/ob-cumulus-prod-public/PACE_OCI.20250501T234515.L2.SFREFL.V3_0.nc?transpose=/geophysical_data/rhos:0,1" target="_blank">https://obdaac-tea.earthdatacloud.nasa.gov/ob-cumulus-prod-public/PACE_OCI.20250501T234515.L2.SFREFL.V3_0.nc?transpose=/geophysical_data/rhos:0,1</a>"
          --output-format text | grep "Size"<br>
          #Size is 1485, 1272<br>
          <br>
          At any rate in the default case the geolocation arrays don't
          match the data orientation. </div>
        <div><br>
        </div>
        <div>Cheers, Mike</div>
        <div> </div>
        <div><br>
        </div>
        <span class="gmail_signature_prefix">-- </span><br>
        <div dir="ltr" class="gmail_signature">
          <div dir="ltr">
            <div dir="ltr" style="color:rgb(34,34,34)">
              <div dir="ltr">Michael Sumner</div>
              <div dir="ltr">Ordinary Member,  Streets People Love
                Hobart Association<br>
                Research Software Engineer<br>
                Australian Antarctic Division<br>
                Hobart, Australia</div>
              <div dir="ltr">0438489030<br>
                e-mail: <a href="mailto:mdsumner@gmail.com" style="color:rgb(17,85,204)" target="_blank">mdsumner@gmail.com</a></div>
            </div>
          </div>
        </div>
      </div>
      <br>
      <fieldset></fieldset>
      <pre>_______________________________________________
gdal-dev mailing list
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a>
</pre>
    </blockquote>
    <pre cols="72">-- 
<a href="http://www.spatialys.com" target="_blank">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
  </div>

</blockquote></div><div><br clear="all"></div><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div dir="ltr" style="color:rgb(34,34,34)"><div dir="ltr">Michael Sumner</div><div dir="ltr">Ordinary Member,  Streets People Love Hobart Association<br>Research Software Engineer<br>Australian Antarctic Division<br>Hobart, Australia</div><div dir="ltr">0438489030<br>e-mail: <a href="mailto:mdsumner@gmail.com" style="color:rgb(17,85,204)" target="_blank">mdsumner@gmail.com</a></div></div></div></div>