3
This is a well-known python GDAL "gotcha". The data doesn't get written until the dataset is closed and dereferenced.
lin = None
#Or
del lin
Apart from that, your code works for me:
2
A test with your CSV file and gdal_command line utility print a lot of errors. It seems that your sample data is not suitable.
gdal_grid grid.vrt -a linear -zfield Elevation out.tif
Warning 1: qhull precision error: initial facet 1 is coplanar with the interior point
Warning 1: ERRONEOUS FACET:
Warning 1: - f1
Warning 1: - flags:
Warning 1: bottom
...
2
Running gdalinfo with debugging options reveals the reason for the failure:
gdalinfo /vsicurl/https://ftpprd.ncep.noaa.gov/data/nccf/com/hrrr/prod/hrrr.20210505/conus/hrrr.t16z.wrfsfcf01.grib2 --debug on --config CPL_CURL_VERBOSE YES
The end of the log:
> GET /data/nccf/com/hrrr/prod/hrrr.20210505/conus/hrrr.t16z.wrfsfcf01.grib2 HTTP/1.1
Host: ftpprd....
2
Regarding 1: if you store the data in Memory Drivers, then yes, you will hold all of the data in RAM (it might be a little more than on disk, if e.g. the files on disk have some compression scheme).
The GDAL documentation mentions:
Native object gets destroyed when
Python object goes out of scope, or when they are assigned to None. So replace
foo.Destroy() ...
1
You can use ds.update_tags(band, tag=value). See the docs for additional info.
Here's an example with an existing dataset (opened in r+ mode).
import rasterio as rio
with rio.open('test.tif', 'r+') as ds:
for i in range(1, ds.count+1):
# Here I add a non-standard tag called "DATE_TIME"
# this can be anything, like "...
1
Both the four GCP images and the six GCP images display with an affine transformation in ArcGIS Pro 2.7, QGIS 3.4, and QGIS 3.10 on Windows 10.
QGIS 3.18 Ubuntu Linux v20.04 is not displaying the six GCP images correctly. The 6 GCP images displayed with the upper left image corner at the equator/greenwich with no rotation, scaling, or resizing. The problem ...
Only top voted, non community-wiki answers of a minimum length are eligible
