316 reputation
29
bio website
location
age
visits member for 2 years, 8 months
seen May 8 at 13:20

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.
May
7
answered Is NAD 83 the same as EPSG:4326?
Apr
5
awarded  Revival
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
answered How to compute area for polygons in an EPSG:4326 shapefile?
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.
Jan
8
answered Viewing/understanding NetCDF files
Jan
7
awarded  Yearling
Jan
7
awarded  Necromancer
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.