153 reputation
8
bio website
location
age
visits member for 2 years, 3 months
seen Dec 5 at 15:21

Nov
19
comment How to easily shift all features in a vector dataset?
Brilliant -- a simple one line CLI solution.
Aug
12
answered gdal_grid unable to get nodata outside of original data area
Aug
12
answered gdal_grid failing to set nodata
Aug
11
comment Writing numpy array to raster file
The "rotation is in twice" to account for the effect of a rotated bit of y on x and the rotated bit of x on y. See lists.osgeo.org/pipermail/gdal-dev/2011-July/029449.html which tries to explain the interrelationships between the "rotation" parameters.
Aug
5
revised Bilinear interpolation of point data on a raster in Python?
deleted 2 characters in body
Aug
5
answered Bilinear interpolation of point data on a raster in Python?
Aug
5
answered gdal_grid has some problem with reading vrt file
Aug
4
comment gdal_grid failing to set nodata
They might be thinking of this as a default value to store in the array, while identifying that value as flag for nodata value is a separate function. That said, one would think that this might work: ######################## gdal_grid -co PROFILE=GDALGeoTIFF -co TIFFTAG_GDAL_NODATA=-9999 .... ############# but with gdal 11.1 it fails with "Warning 6: Driver GTiff does not support TIFFTAG_GDAL_NODATA creation option"
Aug
4
comment Bilinear interpolation of point data on a raster in Python?
ax,ay might have a problem if the grid is rotated at all. It might be better to transform your interpolate-to points into pixel or data coordinates. Also, if there is an off-by-one problem with them, then you might be going beyond the size of the band.
Aug
1
comment What is the inverse of QgsGeometry.asGeometryCollection() in qgis python?
combine looks like it would merge the parts, losing their individual identities. I don't want them merged into one, I want a new multipart geometry.
Aug
1
comment Negative UTM easting
110,000m/1852/60 = 0.998 degrees wide, which is too narrow. 6 degrees would be about 6*1852*60=666,720m wide.
Jul
30
comment What is the inverse of QgsGeometry.asGeometryCollection() in qgis python?
My motivation is the same specified distance redistributing-vertices use-case as in gis.stackexchange.com/q/107999/10229 but I'd like to know how to write things to deal with polygons, and potentially 2.5D polygons and polylines as well. (Ultimately, I need to redistribute the vertices unevenly, based on a raster or some function of a raster, but that's for later.)
Jul
30
revised What is the inverse of QgsGeometry.asGeometryCollection() in qgis python?
expand the question
Jul
30
comment What is the inverse of QgsGeometry.asGeometryCollection() in qgis python?
No, I do not want to replace the multipart with a set of component parts (like 'explode' ungroup, or multi-to-single). I want to unpack the multigeom into its component parts, apply a user-defined function (in the Q, myFunction() is an identity function as a , but a function like reverse_line(), densify(), simplify()) to the parts to make a set of new parts, then re-pack the new parts into a new, replacement multigeom to use as a replacement. In your answer, how do I append each g.exportToWkt() into a list, and convert the list back into a multipart geometry of the same type as geom?
Jul
25
revised What is the inverse of QgsGeometry.asGeometryCollection() in qgis python?
added 168 characters in body
Jul
25
revised How does one modify selected geometry with a function in qgis python console?
added 46 characters in body
Jul
25
comment What is the inverse of QgsGeometry.asGeometryCollection() in qgis python?
After I iterate, I then want to re-collect the list of paert into a type-conformant replacement for the original geometry. It appears that you can take apart all types of .isMultipart()==True geometries with .asGeometryCollection(), but you can't re-collect them with something like geom.replaceGeometry(geomColl).
Jul
25
revised How does one modify selected geometry with a function in qgis python console?
Fix for WKBMultiLineString / MultiPolyline / multiple-line cases
Jul
25
revised What is the inverse of QgsGeometry.asGeometryCollection() in qgis python?
More trial and error, but found some working code.
Jul
24
revised What is the inverse of QgsGeometry.asGeometryCollection() in qgis python?
Add links to docs.