[gdal-dev] Missing geometry type when using ogr2ogr to write FileGDB

[Stefan Keller]

I got the real problem: It occurs with this WFS source:
http://maps.zh.ch/wfs/FnsNSWFS?SERVICE=WFS&REQUEST=GetCapabilities

=> ogrinfo reports no geometry with this WFS:

$ ogrinfo -ro WFS:http://maps.zh.ch/wfs/FnsNSWFS?TYPENAME=svo_zonen
INFO: Open of `WFS:http://maps.zh.ch/wfs/FnsNSWFS?TYPENAME=svo_zonen'
      using driver `WFS' successful.
1: svo_zonen

Should be
1: svo_zonen (Polygon)

No geometry reported here in parantheses. So it's a configuration or
installation problem of the WFS server (ArcGIS)

In QGIS this works (QGIS reader seems to guess somehow the geometry type).

Yours, Stefan

P.S. Still unsure why the OGR WFS reader stumbles over additional WFS URL
parameters appended. Could be because of the key parameter?

$ ogrinfo -ro WFS:
http://openpoimap.ch/featureserver/workspace.cgi?key=HG4kFX2kEEqBSqvcswuh2E&SERVICE=WFS
INFO: Open of `WFS:
http://openpoimap.ch/featureserver/workspace.cgi?key=HG4kFX2kEEqBSqvcswuh2E'
      using driver `WFS' successful.
1: pois (Point)
'SERVICE' is not recognized as an internal or external command,
operable program or batch file.



2013/4/21 Stefan Keller <sfkeller at gmail.com>

> Weird enough, this ogrinfo command works
>
> $ ogrinfo -ro WFS:
> http://openpoimap.ch/featureserver/workspace.cgi?key=HG4kFX2kEEqBSqvcswuh2E
> INFO: Open of `WFS:
> http://openpoimap.ch/featureserver/workspace.cgi?key=HG4kFX2kEEqBSqvcswuh2E
> '
>       using driver `WFS' successful.
> 1: pois (Point)
>
> This works too but stumbles over any additional URL parameters:
> $ ogrinfo -ro WFS:
> http://openpoimap.ch/featureserver/workspace.cgi?key=HG4kFX2kEEqBSqvcswuh2E&SERVICE=WFS&TYPENAME=pois
> INFO: Open of `WFS:
> http://openpoimap.ch/featureserver/workspace.cgi?key=HG4kFX2kEEqBSqvcswuh2E
> '
>       using driver `WFS' successful.
> 1: pois (Point)
> 'SERVICE' is not recognized as an internal or external command,
> operable program or batch file.
> 'TYPENAME' is not recognized as an internal or external command,
> operable program or batch file.
>
> Yours, Stefan
>
>
> 2013/4/21 Stefan Keller <sfkeller at gmail.com>
>
>> Many thanks for the quick reply, Even!
>>
>> You're right, it seems to work with Shapefile - but not for the WFS
>> reader (and own server) we are testing...
>>
>> Below you find the GML output of our WFS source.
>>
>> So is the "fault" in the WFS reader - or even in our WFS Response XML?
>>
>> Yours, Stefan
>>
>>
>> Testing OpenPOIMap WFS Service
>> --------------------------------------------------
>>
>> Get all hostels of Switzerland from OpenStreetMap database
>> ("tourism=hostel"). Response edited:
>>
>> GetCapabilities:
>> http://openpoimap.ch/featureserver/workspace.cgi?key=HG4kFX2kEEqBSqvcswuh2E&SERVICE=WFS&REQUEST=GetCapabilities
>> Response => (...)
>>
>> GetFeature:
>> http://openpoimap.ch/featureserver/workspace.cgi?key=HG4kFX2kEEqBSqvcswuh2E&SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=pois&SRSNAME=EPSG:4326
>>
>> Response =>
>> <wfs:FeatureCollection
>>   xmlns:fs="http://example.com/featureserver"
>>   xmlns:wfs="http://www.opengis.net/wfs"
>>   xmlns:gml="http://www.opengis.net/gml"
>>   xmlns:ogc="http://www.opengis.net/ogc"
>>   xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
>>   xsi:schemaLocation="http://www.opengis.net/wfs
>> http://schemas.opengeospatial.net//wfs/1.0.0/WFS-basic.xsd">
>>
>>   <gml:featureMember gml:id="-324390">
>>     <fs:pois fid="-324390">
>>       <fs:way>
>>         <gml:Point srsName="EPSG:4326">
>>           <gml:coordinates decimal="." cs="," ts="
>> ">6.59728931318,46.5195096642</gml:coordinates>
>>         </gml:Point>
>>       </fs:way>
>>        <fs:shop>None</fs:shop>
>>       <fs:amenity>None</fs:amenity>
>>       <fs:surface>None</fs:surface>
>>       <fs:street>None</fs:street>
>>       ...
>>     </fs:pois>
>>   </gml:featureMember>
>>   ...
>> </wfs:FeatureCollection>
>>
>>
>> 2013/4/21 Even Rouault <even.rouault at mines-paris.org>
>>
>>> Le dimanche 21 avril 2013 17:49:28, Stefan Keller a écrit :
>>> > When trying to use ogr2ogr (ver. 1.9) to convert a Shapefile (geometry
>>> type
>>> > Polygon) to file geodatabase (FileGDB) I'm getting an error like "FGDB
>>> > layers cannot be created with a wkbUnknown layer geometry type.":
>>> >
>>> > $ ogr2ogr -f filegdb out.gdb polygon.shp
>>> >
>>> > When adding -nlt it seems to work:
>>> >
>>> > $ ogr2ogr -f filegdb out.gdb polygon.shp -nlt POLYGON
>>> >
>>> > This is to a) an usual OGR behavior, and b) a show stopper for my
>>> script
>>> > because I don't know before hand (and won't check with ogrinfo) what
>>> > geometry type the input source is. Seems like there is a problem in the
>>> > FileGDB driver.
>>> >
>>> > Can anybody confirm and explain this bug?
>>>
>>> Stefan,
>>>
>>> I don't confirm that behaviour with either 1.9.2 or 1.10.0. Conversion
>>> works
>>> fine for me without specifying -nlt. What is the geometry type reported
>>> by
>>> ogrinfo on polygon.shp ?
>>>
>>> Even
>>>
>>
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20130421/9547383f/attachment.html>