Re: [gdal-dev] RE Python ogr GetFieldAsString problem

gt; something like preg_match in php.... > Thanks > steve > Steve Toutant, M. Sc. > Analyste en géomatique > Secteur environnement > Direction de la santé environnementale et de la toxicologie > Institut national de santé publique du Québec > 945, avenue Wolfe > 4e étage, Québec, Qc

Re: [gdal-dev] RE Python ogr GetFieldAsString problem

t is a real, But GetFieldAsDouble also truncate the value > at the comma > > > > Steve Toutant/INSPQ/ > 2015-05-13 15:11 > > A > gdal-dev@lists.osgeo.org > cc > > Objet > Python ogr GetFieldAsString problem > > > > > > > > > I'm try

[gdal-dev] RE Python ogr GetFieldAsString problem

ouble also truncate the value at the comma Steve Toutant/INSPQ/ 2015-05-13 15:11 A gdal-dev@lists.osgeo.org cc Objet Python ogr GetFieldAsString problem I'm trying to parse a GML and read all the values of an attribute. Example -64,2131 Describe

[gdal-dev] Python ogr GetFieldAsString problem

I'm trying to parse a GML and read all the values of an attribute. Example -64,2131 DescribeFeatureType for this layer returns wfs_ds = ogr.Open(tmp_file.name) #tmp_file.name is the response of a WFS getFeature request layer = wfs_ds.GetLayerByName('SMDVDP_ST

[gdal-dev] OGR read WFS collection

Hello I successfully used this code to create a datasource from a URL step 1- wfs_drv = ogr.GetDriverByName('WFS') step 2- wfs_ds = wfs_drv.Open('WFS:' + url) step 3- Then I parsed the result I need to create another script but in input I get the result of a WFS query. So I need to parse this res

Re: [gdal-dev] Parse GML without geometry

Got it! Lots of fun! Good links for python http://pcjericks.github.io/py-gdalogr-cookbook/vector_layers.html#get-wfs-layer-and-iterate-over-features Thousand thanks Even Rouault @lists.osgeo.org Envoyé par : gdal-dev-boun...@lists.osgeo.org 2014-12-12 15:42 A gdal-dev@lists.osgeo.org cc st

Re: [gdal-dev] Parse GML without geometry

thanks Kyle I ended up with geom=osgeo.ogr.CreateGeometryFromGML(wfs_response) Returns ERROR 1: Unrecognised geometry type If I understand well, I think my question now is How do I parse the feature collection with ogr? Kyle Shannon 2014-12-12 14:47 A steve.tout...@inspq.qc.ca cc GDAL

[gdal-dev] Parse GML without geometry

Hi, I need to parse the GML returned by this query, to make a list of unique value of NA2DESC tag using python http://demo.mapserver.org/cgi-bin/wfs?SERVICE=WFS&REQUEST=GetFeature&version=1.1.0&Typename=continents&maxfeatures=10&propertyname=(NA2DESC) I could not make it with a XML parser yet.

[gdal-dev] fatal errr compiling gdal 1.11.0

I already compiled gdal 1.11.0 on centos6 and I need to install on a fedora 20. Using the same method I get this eror libtool: compile: g++ -g -O2 -DHAVE_SSE_AT_COMPILE_TIME -fvisibility=hidden -Wall -I.. -I/opt/geomatique/src/gdal-1.11.0/port -I/opt/geomatique/src/gdal-1.11.0/gcore -I/opt/ge

[gdal-dev] RE [Polluriel potentiel] Re: ogr ODBC problem [SOLVED]

for ogrinfo we were using, as explained here, http://mapserver.org/input/vector/VirtualSpatialData.html#steps-for-display ogrinfo ODBC:User/Pwd@mssql__archives msStatut_station2002_2006 That was working on windows Though, That doesn't work on our linux distro linux sles 11 sp1. It as to be (comma

[gdal-dev] RE [Polluriel potentiel] Re: ogr ODBC problem

Thanks Brent for your answer. But view's primary key? I'm missing something here...What do you mean by that? Do you mean a column with a unique id in the view? There is no such column in my view. Is ogr now needs a unique id in a view? Because it was working before I migrate to linux with a

Re: [gdal-dev] ogr ODBC problem

Thanks Jeff I got now OGR_ODBC: Table ?s???s!.? has no identified FID column. I found that several users had this problem but found no solution. I don't have write access to this MSSQL server. I'm connecting via ODBC to a non spatial table, but it contains latitude and longitude informa

[gdal-dev] RE ogr ODBC problem

Note that it is a non-spatial mssql database Steve Toutant/INSPQ/ 2012-05-30 16:09 A gdal-dev@lists.osgeo.org cc Objet ogr ODBC problem Hi! I use this command to get the tables from an ODBC connection ogrinfo ODBC:User/Pwd@DNS The connection is succesful but I get this error

[gdal-dev] ogr ODBC problem

Hi! I use this command to get the tables from an ODBC connection ogrinfo ODBC:User/Pwd@DNS The connection is succesful but I get this error several times ERROR 1: No column definitions found for table '�s���s!.���', layer not usable. I used OGR ODBC for several months from a Windows serv

[gdal-dev] Linux installation

Hi,Sorry for cross-posting on several lists, but we are in the process of taking some important decisions, on the OS to use...For several years we have mapserver, gdal and Postgis, etc... running on a Windows server.We plan to migrate our GIS tools from the Windows server to SUSE Linux Enterprise S

[gdal-dev] ogr and postgres

Hi,I just download gdalwin32-1.6 and It seems that PostgreSQL driver is not included in it.How can I use ogr with postgis supported on windows?thankssteveSi vous n'êtes pas le destinataire de ce message, veuillez le détruire près avoir informé l'expéditeur de son erreur. Par ailleurs, il est inter

RE: [gdal-dev] Translating between UTM and Geographic

Maybe the problem is in your data. Di you try to reproject zone 11 nad27 with the wkt nad83 below? "Henneke, Amanda M" 23/09/2009 01:44 PM A "steve.tout...@inspq.qc.ca" cc "gdal-dev@lists.osgeo.org" , "gdal-dev-boun...@lists.osgeo.org" Objet RE: [gdal-dev] Translating between UTM and Ge

RE: [gdal-dev] Translating between UTM and Geographic

What is the wkt you are using for zone 12 nad83? "Henneke, Amanda M" @lists.osgeo.org Envoyé par : gdal-dev-boun...@lists.osgeo.org 23/09/2009 01:39 PM A "gdal-dev@lists.osgeo.org" cc Objet RE: [gdal-dev] Translating between UTM and Geographic Here is some more information: So, h

RE [gdal-dev] Translating between UTM and Geographic

The problem might not be in your code but in the WKT projection information (sSpRefParams). Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques biologiques, environnementaux et occupationnels Institut national de santé publique du Québec 945, avenue Wolfe

RE: RE [gdal-dev] Translating between UTM and Geographic

I don't know about the code I'm not a programmer. You said UTM Zone 12 with NAD83 works but this Zone 11 NAD27 doesn't. My guess is you are using Zone 11 NAD83 instead NAD27. I might be wrong though. Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnemen

RE [gdal-dev] Translating between UTM and Geographic

hi, slight shift in the points (to the east) and an extremely slight shift up (to the north). That seems to be a datum error. Make sure you set the good parameters for datum when using nad83 or nad27. Steve "Henneke, Amanda M" @lists.osgeo.org Envoyé par : gdal-dev-boun...@lists.osgeo.org

[gdal-dev] intersection between vector and raster

ue is greater that 205. Is there a gdal utlity I could use to do that without programming? thanks Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques biologiques, environnementaux et occupationnels Institut national de santé publique du Québec 945, a

[gdal-dev] Re: gdalwarp merge problem (SOLVED)

mages. Even if they are in BYTE, gdalinfo returns the nodata value, Min, Max, etc. of the original images in float32. Is that a bug? thanks Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques biologiques, environnementaux et occupationnels Institut nat

[gdal-dev] gdalwarp merge problem

Hi, Still trying to merge images... I need to merge several images. As an example The image EstMtl.tif must be over the image EstGa.tif. Here is the command gdalwarp -srcnodata -3.4028234663852886e+038 -dstnodata 0 -of "GTiff" -te -931841.246 126213 241462 805816 -tr 30.0 30.0 EstGa.tif icu.tif

Re: [gdal-dev] merge with gdalwrap

linfo? - The Pixel Size is not the same for all images. Per example I got (30.099,-30.103) and this for another image Pixel Size = (30.000,-30.000) What will append if I specify -tr 30 -30 as resolution in this case thanks Steve Steve Toutant,

[gdal-dev] merge with gdalwrap

028234663852886e+38 -o icu.tif -of "GTiff" EstRN.img EstSh.img It works but maybe some more parameters are important. Thanks for your help, Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques biologiques, environnementaux et occupationnel

Re: [gdal-dev] Numeric package SOLVEd

Thanks Scott The numpy version installed on the server wasn't for win32. I download the correct version and it is working. Thanks again for your help, Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques biologiques, environnementaux et occupatio

[gdal-dev] Numeric package

numpy and numeric package were installed but it seems not. How can I install the numeric package properly? thanks Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques biologiques, environnementaux et occupationnels Institut national de santé publique du Q

Re: [gdal-dev] Messages: gdal_translate and gdal_rasterize

Thanks Frank, " the command cannot determine it from the .shp file" Why? Is it because the EPSG code of this .shp is unknown to GDAL or GDAL can't get any SRS from a .shp? Should I add the EPSG code somewhere in a configuration file? Best regards, Steve Steve Toutant, M.

[gdal-dev] Messages: gdal_translate and gdal_rasterize

I care about these messages? Is there a parameter missing in my commands? Thanks a lot for your help!! Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques biologiques, environnementaux et occupationnels Institut national de santé publique du Québec 945

Re: [gdal-dev] RE clip intersected

specify a range of value? Example: merge contiguous cells having a value > 16 and < 19 thanks Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques biologiques, environnementaux et occupationnels Institut national de santé publique du Québec 945, avenue

[gdal-dev] RE clip intersected

seems that this script doesn't exist. If the method above is valid, - how can I create a boundary excluding pixels with nodata? - And, of course, where can I get gdal_polygonize.py Thanks! Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risq

[gdal-dev] clip intersected

Hi, I have 2 Erdas (.img) images that overlap. See attached image. I need to clip one of the 2 images by removing the intersected (dark red) part. If it's possible, how should I accomplish this? thanks for your help, Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environn

[gdal-dev] crop an image

polygon is irregular. So the new image contains the region I need with some extra. From this new image Is there a way to assign "nodata" to pixels outside the irregular polygon? Or maybe there is a way to crop an image from an irregular polygon? thanks Steve Steve Toutant, M. Sc. A