Is there any good, complete API for the GDAL Python bindings?
I am interpolating and rasterizing points with elevation data using the
output = gdal.Grid('outcome.tif','newpoints.vrt')
In its most basic sense, it's working perfectly. But now, I'd like to try the linear approach;
output = gdal.Grid('outcome.tif','newpoints.vrt', algorithm = 'linear:radius=0')
That does generate a new .tif file (512kb), but is just one big black image. Also, not all options seem to be accepted.. for the
invdist:power=4 for example it also results in a useless image. What is going wrong, am I missing something? Below is the .vrt I'm using;
<OGRVRTDataSource> <OGRVRTLayer name="newpoints"> <SrcDataSource>newpoints.csv</SrcDataSource> <GeometryType>wkbPoint</GeometryType> <GeometryField encoding="PointFromColumns" x="long" y="lat" z="depth"/> </OGRVRTLayer>
As a next example I used the
gdal.GridOptions(), putting the parameters in that one. It is resulting in a correct file, but still ignores the options. How to actually pass the
GridOptions onto the
I found this question also; GDAL grid produces NoData Is the linear interpolation just a bug in GDAL?
See my example workfiles; https://gist.github.com/willemvanopstal/1eae4999fcfa62e89cca39a350b504d8