323 reputation
213
bio website
location
age
visits member for 2 years, 4 months
seen Nov 29 at 12:53

Sep
27
comment Find correct SRID
Right, good catch... I didn't noticed that UNIT isn't closed. Thanks for confirmation.
Sep
27
comment Find correct SRID
I did what was most obvious to me and added new record for SRID 9001 to spatial_ref_sys column, afterwhich all seem to work as expected.
Sep
4
comment How to correctly import GeoJSON in PostGIS?
@user30184 ring errors are because I have many geometries where two separate rings create single feature and they are treated by ogr2ogr as polygons. If I convert these to multipolygon, errors are gone.
Aug
27
comment Crop and reproject shapefile based on GeoTIFF sample
Nice tip. I tried to reproject and clip at the same time, but I put -clipdst before the filenames, so it didn't work, and I didn't read the error message carefully about partial reprojection possiblity with OGR_ENABLE_PARTIAL_REPROJECTION nor I would have known how to define that option ;)
Aug
27
comment Crop and reproject shapefile based on GeoTIFF sample
Thanks, it was countries shapefile, but I managed to get correct output by first cropping to outer rectangle of GeoTIFF file, then reprojecting and then clipping to projected coordinates. I was just hoping that straightforward approach of using georeferenced file as reference exists in gdal/org system.
Aug
27
comment Crop and reproject shapefile based on GeoTIFF sample
Projecting whole globe to destination projection is not possible, so I have to crop the source first, and that is where I stopped in the first place as corner coordinates are not forming rectangle in Natural Earth shapefile source projection. Any suggestions?
Aug
26
comment Query PostGIS database without creating primary key
Sorry, it's my fault. Before issuing PostgreSQL layer has no primary key. there was one more message: Duplicate field Code found. So that was the real reason why the query didn't load as layer - latter message seems like result from the previous fault.
Jun
19
comment Will ogr2ogr correctly merge features from several OSM files to resulting set
osm2pgsql (as least the version I tried on Windows) does not create primary key, and I guess other OS versions does not too, as osm_id is not unique.
Jun
12
comment Geometry to Shapefile
By one of many interfaces to SQL exposed to VB.Net
Jun
4
comment How to migrate PostGIS spatial data to MsSql spatial data?
very cool spartan program, you just copy/paste layers between servers, and it does the rest for you
Jun
4
comment Export GeoJSON to MSSQL with UTF-8 encoding
Well, open ticket you linked and this question show that there is no way to import Unicode with org2org to mssqlserver. UTF-8 data is passed to mssqlserver, but it doesn't recognize it as such because the column type is set to varchar instead nvarchar. If column is set to nvarchar, mssqlserver would apparently import UTF-8 data and convert it on the fly to it's Unicode encoding variant - UTF-16LE.
Jun
3
comment Export GeoJSON to MSSQL with UTF-8 encoding
So this task is impossible currently. Thanks for the information.
Jun
2
comment Export GeoJSON to MSSQL with UTF-8 encoding
It's just the same. Both GeoJSON and now converted GML, render fine in UTF-8 text editor, but not in database table.
May
27
comment How to georeference using GDAL tools?
According documentation and test: -a_ullr ulx uly lrx lry or with lat/lon terms -a_ullr ullon ullat lrlon lrlat is the correct way to use the command.
Jun
8
comment re-project from EPSG:4326 to polar projection
Perhaps input has to be cropped first, like you did. Following your steps (i.e. cropping the input to north hemisphere) I also get expected result. Trying again on whole globe it fails for some reason. Thanks for you answer :)
Jun
8
comment re-project from EPSG:4326 to polar projection
It reaches 90deg north. It's etopo tiff provided here: ngdc.noaa.gov/mgg/global. I'll update question with gdalinfo report. Thanks
Dec
27
comment What's the difference between gdalwarp and gdal_merge for mosaicing
Had no idea about this tuning switch, which indeed speeds up the process. For reference 3000 is 3GB (more info trac.osgeo.org/gdal/wiki/UserDocs/GdalWarp) I have NumPy build with MKL in my system Python, but still gdalwarp must be doing something more, as NumPy can't be faster then C implementation, although can reach it.
Dec
13
comment How to re-project shapefiles?
Thanks I didn't know about that, and Googling it showed me what you meant. It's exactly what I needed. If you could have only posted your last comment yesterday... I did lot of typing to get there ;)
Dec
13
comment How to re-project shapefiles?
I was trying to overlay vector layer from one projection (epsg:31277) to another (epsg:4326). However I did not understood this well, as now I know that I don't need re-projection (changing the shapefile data) but just change to shapefile projection file (.prj) by using -a_srs instead -t_srs. Issue why EPSG:4326 SRS that didn't had projection, was because it's WGS84 geographic coordinate system, so it doesn't need projection, as was shown to me by others.
Dec
13
comment How to re-project shapefiles?
I tried different projections and they all output very different values for coordinates. I'll ask at gdal mailing list, and post results if I get answer. Thanks for your time