Re: [gdal-dev] Storing Gdal datasets in Postgis

2010-03-01 Thread Jorge Arevalo
On Fri, Feb 26, 2010 at 2:51 PM, Ricardo Cezar Bonfim Rodrigues rikardoce...@msn.com wrote: Hi, I'm using Gdal to load dted files and get max elevations of a region, but I have constraints concerning speed. Inspite of loading geotifs generated from dted has reduced a lot the processing

RE: [gdal-dev] Storing Gdal datasets in Postgis

2010-03-01 Thread Ricardo Cezar Bonfim Rodrigues
Hi Jorge, I've implemented a customized gdal reader using the method readBlock instead of using RasterIO and now it seems to be much faster and will propably be adequate to my problem. But thanks by the tips, I'll keep it in mind if I need it. Ricardo Rodrigues Brazil From:

[gdal-dev] (no subject)

2010-03-01 Thread Discourse Maps
I am simply trying to open a raster image, manipulate in NumPy and then spit it back out in GDAL. I am having trouble setting null data... I tried: ds.GetRasterBand(1).SetNoDataValue( - ) but nothing seems to work. Any help is appreciated! Below is my code: #! /usr/bin/env python

Re: [gdal-dev] (no subject)

2010-03-01 Thread Frank Warmerdam
Discourse Maps wrote: I am simply trying to open a raster image, manipulate in NumPy and then spit it back out in GDAL. I am having trouble setting null data... I tried: ds.GetRasterBand(1).SetNoDataValue( - ) but nothing seems to work. ... Dear Discourse Maps, Presumably you should be

[gdal-dev] Python Buffer Features

2010-03-01 Thread Spencer Gardner
I'm writing some code that involves buffering a feature in a point layer and then doing a spatial filter with another set of layers. I've got everything working just fine but I'm having a hard time determining how the spatial filter works. Based on the documentation, it looks like all

Re: [gdal-dev] Python Buffer Features

2010-03-01 Thread Frank Warmerdam
Spencer Gardner wrote: I'm writing some code that involves buffering a feature in a point layer and then doing a spatial filter with another set of layers. I've got everything working just fine but I'm having a hard time determining how the spatial filter works. Based on the documentation, it

[gdal-dev] Problems with gdalwarp re-projecting from lat-long to Lambert Conformal

2010-03-01 Thread Bill Hudspeth
Hello, I am having problems re-projecting an ArcAscii Grid from LatLong/WGS84 to Lambert Conformal. The source grid was derived from a GRASS mapset. Using r.out.gdal, a GRASS raster was exported to AAIGrid format using: r.out.gdal input=MODIS_input_raster format=AAIGrid type=UInt16

[gdal-dev] Georeferencing Geostationnary (+proj=geos) projection

2010-03-01 Thread Stéphane
Hello, I'm trying to put georeference information to a raster HDF file which doesn't containt projection nor geo-transformation or any GCP. HDF is not georeferenced, but contain lot of useful metadata, so I manage to get the projection, which is Geostationnary projection ( +proj=geos). However,

[gdal-dev] Transform not working

2010-03-01 Thread Jamie Lahowetz
I'm new to GDAL and am usingthe python bindings. I want to transform a newly created shapefile from WGS84 to a Lambert Conical Conformal so that I can calculate the length of the contained polylines. The script exits with no errors but when I look at the shapefile in arcCatalog I see that the

[gdal-dev] Re: Transform not working

2010-03-01 Thread deadpickle
Not much of a response going on here. In order to change the shapefile projection from geographic to projected do I need to reopen the file, get each feature, and transform the points individually? If so what calls do I need to invoke? On Mar 1, 12:54 pm, Jamie Lahowetz deadpic...@gmail.com

Re: [gdal-dev] Re: Transform not working

2010-03-01 Thread Chris Garrard
Jamie, Try transforming the line geometry before adding it to the feature. That should do it. chris On Mon, Mar 1, 2010 at 4:24 PM, deadpickle deadpic...@gmail.com wrote: Not much of a response going on here. In order to change the shapefile projection from geographic to projected do I need

Re: [gdal-dev] Re: Transform not working

2010-03-01 Thread Jamie Lahowetz
Thanks for the response. I'm sorry I'm really new to this and am not sure where to transform the line geometry, I know it should be before layer.CreateFeature(feature) but after I add the points to the geometry. Is it before or after feature.SetGeometry(line)? And is the right way to call the

Re: [gdal-dev] GDAL for BNG to Lat/Long, WGS84

2010-03-01 Thread scottd777
Please? Anyone? I'm sure the response will be really quick... -- View this message in context: http://n2.nabble.com/gdal-dev-GDAL-for-BNG-to-Lat-Long-WGS84-tp4652000p4659130.html Sent from the GDAL - Dev mailing list archive at Nabble.com. ___

Re: [gdal-dev] GDAL for BNG to Lat/Long, WGS84

2010-03-01 Thread Frank Warmerdam
scottd777 wrote: Hi, and thank heaven for GDAL! I was so happy to find this tool, but I am having some trouble figuring out the exact command I need for an image warping application. If there is someone who doesn't mind a newbie type question, I could really use the help! Essentially, I

Re: [gdal-dev] Re: Transform not working

2010-03-01 Thread Frank Warmerdam
deadpickle wrote: Not much of a response going on here. In order to change the shapefile projection from geographic to projected do I need to reopen the file, get each feature, and transform the points individually? If so what calls do I need to invoke? On Mar 1, 12:54 pm, Jamie Lahowetz

[gdal-dev] Re: Transform not working

2010-03-01 Thread Jamie Lahowetz
I figured it out... with all of your help. I had the coordtrans = osr.CoordinateTransformation switched. I wanted to input lat/lon (WGS84) and transform it to lcc (source) and I also had to set the source prj to lcc. I think you all said this stuff already, just took me a little while. Thanks for

Re: [gdal-dev] Problems with gdalwarp re-projecting from lat-long to Lambert Conformal

2010-03-01 Thread Markus Neteler
On Mon, Mar 1, 2010 at 7:12 PM, Bill Hudspeth w...@unm.edu wrote: ... When I try to re-project the exported AAIGrid to another projection, I use: gdalwarp -s_srs EPSG:4326 -t_srs '+proj=lcc +lat__1=33n +lat_2=45n +lon_0=97w' MODIS_output_latlong.asc MODIS_output_lambert.asc Note that without