Reputation
340
Top tag
Next privilege 500 Rep.
Access review queues
Badges
2 9
Newest
 Revival
Impact
~22k people reached

Jul
9
comment Does QGIS have WMS-T (WMS with time) support?
It might be rough and awkward, but some of the functionality is available in the python console.
Jul
8
comment Open source web map client with time sliders using WMS source
What does the filter icon look like?
Jul
4
comment Problem with ImageMosaic plugin of Geoserver - invalid shapefile/tif missing in GetCapabilities
Also github.com/geotools/geotools/blob/… details all the time formats it tries.
Jul
4
comment Creating Store in GeoServer using ImageMosaic Plugin gives Error?
Making the directory writable by the tomcat process worked for me. The code attempts to make DIRNAME.(dbf,prj,properties,shp, and shx) filesin the directory, and when it can't it produces this error string.
Jul
4
comment Problem with ImageMosaic plugin of Geoserver - invalid shapefile/tif missing in GetCapabilities
I do not have an answer yet, but I think it is poorly documented. It isn't at all obvious how to mod things. The source code at github.com/geotools/geotools/blob/master/modules/plugin/… and github.com/geotools/geotools/blob/… hint that you could also set a timeregex.properties:format SimpleDateFormat
Jul
2
comment Unable to zoom beyond scale of 1:2,256 in QGIS?
That worked for me. Apparently the plugin ties the zoom level to match pixels with an available zoom from a tile source. Before uninstalling it, I couldn't choose the predefined zooms at the bottom of the UI, now I can set a 1:10,000, etc.. Thanks.
Jun
28
comment How to store and serve the time a GeoTIFF refers to
Geoserver gives some hints near docs.geoserver.org/stable/en/user/tutorials/… Basically one creates a directory of 'temperature_yyyymmddHHMMSS.tif' files and enables the geoserver to parse them into a collection of files, but you could use the same naming convention to help populate and update a 'gdaltindex' shapefile for mapserver.
May
7
comment Proj4 String for NAD83(2011) / Louisiana South (ftUS)
OT, but: If you are using a +k=??? scaling units markedly different than 1.0, then the easting and northing should be similarly scaled. I used a g218gproj4='+proj=lcc +lat_1=25 +lat_2=25 +lat_0=12.19 +lon_0=265 +k_0=8.20463003682238e-05 +x_0=346.906301484053 +y_0=-49.5724920086904 +a=6370000.0,+b=6370000.0 +units=m +no_defs' to translate between geographic coordinates directly into the grid coordinates of a model intended to match the grid from a nco.ncep.noaa.gov/pmb/docs/on388/tableb.html#GRID218 weather model. This scaled to the model's 12km pixels with the SW pixel at 0,0.
Mar
23
comment Convert a KML ImageOverlay to World Files with GDAL
Yes, the WorldFileTool.jar from worldfiletool_0_3_7.zip works.
Mar
23
comment Convert a KML ImageOverlay to World Files with GDAL
The script seems to give an error: '$ java -jar WorldFileTool.jar doc.kml Exception in thread "main" java.lang.NullPointerException at java.io.File.<init>(File.java:277) at ch.hsr.ifs.worldfiletool.ui.WorldFileTool_GUI.<init>(WorldFileTool_GUI.java:98) at ch.hsr.ifs.worldfiletool.StartWorldFileTool.main(StartWorldFileTool.java:40) $ '
Jan
14
comment How to calculate area in km2 Qgis 2.0
There seems to be some room for doubt: 0.004008*cos(latitude)(1.852*60)^2=49.4894km^2 at latitude=0, and there's no latitude where the area will be 1.017* bigger than the equatorial area in decimalDegrees^2.
Nov
19
comment How to easily shift all features in a vector dataset?
Brilliant -- a simple one line CLI solution.
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
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
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
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).