ga_readonly ga_update gdaldatasetbbsize_f gdaldatasetsimpleread_f gdalrastersimpleread_f gdt_byte gdt_cfloat32 gdt_cfloat64 gdt_cint16 gdt_cint32 gdt_float32 gdt_float64 gdt_int16 gdt_int32 gdt_uint16 gdt_uint32 gdt_unknown gf_read gf_write

subroutine gdal::gdaldatasetsimpleread_f ( type(gdaldataseth), value  hds, real(kind=c_double), intent(in)  bbxmin, real(kind=c_double), intent(in)  bbymin, real(kind=c_double), intent(in)  bbxmax, real(kind=c_double), intent(in)  bbymax, real(kind=c_float), dimension(:,:,:), intent(out), allocatable  pbuffer, real(kind=c_double), intent(out)  xmin, real(kind=c_double), intent(out)  ymin, real(kind=c_double), intent(out)  xmax, real(kind=c_double), intent(out)  ymax  ) Even more simplified method for importing data from a dataset within a bounding box specified in georeferenced coordinates. This method imports all raster bands of a dataset hds keeping the data included in a rectangular bounding box specified in terms of georeferenced coordinates and clipped to the dataset domain extension. Datasets with a rotational geotransform are not supported. The result is stored in a 3-d buffer which must be declared as ALLOCATABLE and is allocated inside the method (f2003 feature). The actual bounding box of data obtained, in terms of georeferenced coordinates of centers of first and last grid boxes, and the grid step are returned in output. If an error occurs while accessing the dataset, the buffer is not allocated (check with the ALLOCATED() intrinsic function), while if the bounding box lies outside of the dataset domain one or more dimensions are allocated to zero size (check with the SIZE() intrinsic function). The data is kept in the dataset resolution, no interpolation is done. WARNING: this is an experimental method, so the interface may change in the future, use with care! Todo: add the possibility to swap/transpose the data to the array order required by the caller through an additional parameter Todo: convert to a generic interface with different types for pbuffer Todo: make a write-equivalent? Parameters hds dataset to read [in] bbxmin minimum x georeferenced coordinate of the bounding box [in] bbymin minimum y georeferenced coordinate of the bounding box [in] bbxmax maximum x georeferenced coordinate of the bounding box [in] bbymax maximum y georeferenced coordinate of the bounding box [out] pbuffer buffer containing output data, allocated by the present method, its previous contents, if any, is lost [out] xmin minimum x georeferenced coordinate of the data read into pbuffer (center of first grid box) [out] ymin minimum y georeferenced coordinate of the data read into pbuffer (center of first grid box) [out] xmax maximum x georeferenced coordinate of the data read into pbuffer (center of last grid box) [out] ymax maximum y georeferenced coordinate of the data read into pbuffer (center of last grid box) Definition at line 1050 of file gdal.F90.