17,449 reputation
24381
bio website sites.google.com/site/mwtoews
location Wellington, New Zealand
age 35
visits member for 4 years, 5 months
seen 15 hours ago

I'm a hydrogeologist. An expert of many random things.


1d
comment Finding a substring within a table of data between 2 string positions in PostgreSQL
See dba.stackexchange.com/q/105785/1528
Jul
1
comment rasters written with gdal show up with units of meters when they should be degrees
Also the units for EPSG:3410 is metres, not degrees.
Jul
1
comment rasters written with gdal show up with units of meters when they should be degrees
What's the value of projection? My go-to is from osgeo import osr; osr.UseExceptions(); srs = osr.SpatialReference(); srs.ImportFromEPSG(3410); projection = srs.ExportToWkt()
Jun
30
comment geotiff conversion to coordinate system with gdal
Warning: pixelcoord gives the pixel corner coordinate (usually upper left). To get the centre, you need to apply a half-pixel transformation.
Jun
30
comment Downsampling Geotiff using summation - Gdal/Numpy
Yes: avg = sum / count, so sum = avg * count.
Jun
30
answered Downsampling Geotiff using summation - Gdal/Numpy
Jun
29
comment Output of Add XY Coordinates?
@deiniol your data frame possibly has a different coordinate system, such as Spherical Mercator. This is easy to check in ArcGIS.
Jun
29
comment How to build gdal_calc.py from source in MinGW
Were you able to compile without Python (i.e. no --with-python)? What extra library do you need that was not available?
Jun
28
answered Output of Add XY Coordinates?
Jun
28
answered Remove first seven characters from field string using Field Calculator
Jun
28
revised PostGIS Topology Requires Admin
more a PostgreSQL question
Jun
28
answered PostGIS Topology Requires Admin
Jun
28
comment Why can't I import my shp file into a postgis database with SRID 102631 using shp2pgsql?
A really good web resource is epsg.io. For example, see epsg.io/102631 which also has epsg.io/102631.sql for PostGIS
Jun
24
comment H2GIS and POSTGIS
'SRID=25833;POINT(445588.07 5873455.6)' is in Germany (53°0'28.398"N 14°11'20.623"E). Anything else is in the ocean.
Jun
24
comment Using .len() in IF THEN statement in ArcMap field calculator
If there was an integer field called photo, just use str(!photo!).zfill(5)
Jun
23
comment Adding a column of WKT format in the table (postgreSQL)
INSERT works on rows, so you could insert one row with one column of data. But it would make more sense to insert a row of data with other attributes, like magnitude, etc. And if the row already exists, then an UPDATE statement is required. There are too many ways to get data in a table to describe.
Jun
23
comment Adding a column of WKT format in the table (postgreSQL)
You can always add and modify columns to a table. If the WKT is already properly formatted, then you can try to read that. The geometry can be created from many sources, such as WKT, WKB, GeoJSON, etc.
Jun
22
revised Adding a column of WKT format in the table (postgreSQL)
not a pgAdmin question
Jun
22
answered Adding a column of WKT format in the table (postgreSQL)
Jun
22
comment Rotate Polygon to be Parallel with the X Axis
It's missing, but an enhancement ticket is in place