topic Re: problem with WFS Data, GDAL, Python, ArcGIS, &amp; Defining Projection in Python Questions https://community.esri.com/t5/python-questions/problem-with-wfs-data-gdal-python-arcgis-amp-amp/m-p/239318#M18640 <HTML><HEAD></HEAD><BODY><UL><BR /><LI>You're writing to "E:\\Melo\\USGS\\Roads.gml" but trying to convert "E:\\Roads.gml"</LI><BR /><LI>You don't need the ArcGIS GDAL python bindings, you are not using them in your script.</LI><BR /><LI>You can convert directly from the WFS (including defining the projection) using:</LI><BR /></UL><PRE __default_attr="plain" __jive_macro_name="code" class="jive_macro_code jive_text_macro">ogr2ogr -a_srs EPSG:4267 -f "ESRI Shapefile" E:\\Roads.shp WFS:http://path/to/WFS/service?required_or_optional_parameters_etc...</PRE></BODY></HTML> Wed, 07 Aug 2013 05:35:10 GMT Luke_Pinner 2013-08-07T05:35:10Z problem with WFS Data, GDAL, Python, ArcGIS, &amp;amp; Defining Projection https://community.esri.com/t5/python-questions/problem-with-wfs-data-gdal-python-arcgis-amp-amp/m-p/239317#M18639 <HTML><HEAD></HEAD><BODY><SPAN>I'm pulling in data in gml2 format with a EPSG#4267 from a wfs.</SPAN><BR /><BR /><SPAN>I create a datastream in python and then run the following from the python window in ArcMap:</SPAN><BR /><BR /><PRE class="lia-code-sample line-numbers language-none">data = datastream.read() #write the gml data to a gml file rasterdata = open("E:\\Melo\\USGS" + "\\Roads.gml", 'w') rasterdata.write(data) rasterdata.close() #Run the ogr2ogr to convert to shp cmd = ["ogr2ogr", "-f", "ESRI Shapefile", "E:\\Roads.shp", "E:\\Roads.gml"] subprocess.call(cmd) #Define the projection infc ="E:\\Roads.shp" prjfile = "E:\\4267.prj" arcpy.DefineProjection_management(infc, prjfile) </PRE><BR /><BR /><SPAN>Unfortunately, when i try to add the layer in ArcMap 10.1 it "appears" blank.&nbsp; Typically, I've only seen this happen when their was something wrong with the projection.&nbsp; I know there is data in the file because I've converted it successfully using alternate methods.</SPAN><BR /><BR /><SPAN>I installd GDAL according to </SPAN><A href="http://stackoverflow.com/questions/6009506/unable-to-install-python-and-gdal-dll-load-failed" rel="nofollow noopener noreferrer" target="_blank">these</A><SPAN> instructions exactly (1500/32 bit/2.7).</SPAN><BR /><BR /><SPAN>I did NOT download and install the </SPAN><A href="http://www.arcgis.com/home/item.html?id=1eec30bf5fa042a5a2227b094db89441" rel="nofollow noopener noreferrer" target="_blank">ArcGIS bindings </A><SPAN>.&nbsp; Do I need them?</SPAN><BR /><BR /><SPAN>Any ideas on what could be causing the layer to appear blank in ArcMap?</SPAN></BODY></HTML> Sat, 11 Dec 2021 12:02:57 GMT https://community.esri.com/t5/python-questions/problem-with-wfs-data-gdal-python-arcgis-amp-amp/m-p/239317#M18639 Town_ofSnowflake 2021-12-11T12:02:57Z Re: problem with WFS Data, GDAL, Python, ArcGIS, & Defining Projection https://community.esri.com/t5/python-questions/problem-with-wfs-data-gdal-python-arcgis-amp-amp/m-p/239318#M18640 <HTML><HEAD></HEAD><BODY><UL><BR /><LI>You're writing to "E:\\Melo\\USGS\\Roads.gml" but trying to convert "E:\\Roads.gml"</LI><BR /><LI>You don't need the ArcGIS GDAL python bindings, you are not using them in your script.</LI><BR /><LI>You can convert directly from the WFS (including defining the projection) using:</LI><BR /></UL><PRE __default_attr="plain" __jive_macro_name="code" class="jive_macro_code jive_text_macro">ogr2ogr -a_srs EPSG:4267 -f "ESRI Shapefile" E:\\Roads.shp WFS:http://path/to/WFS/service?required_or_optional_parameters_etc...</PRE></BODY></HTML> Wed, 07 Aug 2013 05:35:10 GMT https://community.esri.com/t5/python-questions/problem-with-wfs-data-gdal-python-arcgis-amp-amp/m-p/239318#M18640 Luke_Pinner 2013-08-07T05:35:10Z