On jeudi 25 juillet 2019 19:12:02 CEST Joe Lee wrote: > I'm glad that this motion has passed! > > I got one question after reviewing Even's latest code and RFC. > > Does it consider a case that subsets 1d dataset from 2d dataset [1]? > The organization of MOP02J (1-D SWATH) and its shape is like below: > > dset[2][n_points] > lat[n_points] > lon[n_points] > > I wish I could be able to extract 1D from 2D using dset[n_points][0] and > make a GeoTIFF grid from the scattered dset/lat/lon using a tool like [2] > or gdalwarp with 3 VRTs that correspond to 1D lat/lon/dset. > > To help you understand the case better, > the equivalent Python code is available in [3]. > > If this workflow can be already done, > can anyone give me a series of commands to be executed?
As far as I understand your need, this could probably be done with a gdalmdimtranslate invokation (or possibly several ones and assemble manually a VRT) https://gdal.org/programs/gdalmdimtranslate.html#gdalmdimtranslate to extract 1D from 2D, you would do something like gdalmdimtranslate in out -array name=dset,view=[:,0] The view argument can be used to use numpy-style indexing & slicing: https://gdal.org/api/gdalmdarray_cpp.html#_CPPv4NK11GDALMDArray7GetViewERKNSt6stringE Even -- Spatialys - Geospatial professional services http://www.spatialys.com _______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
