Keep the holes in the polygon while converting Polygon 25D to Polygon 2D using PythonGeopandas/folium map not displaying?Performing Spatial Join / match Points from dataframe to polygons using Python?geopandas geometry.within problem ValueError: Geometry column cannot contain mutiple geometry typesSpatial join in Geopandas is only finding itselfGeopandas not recognizing geometry typeOverlay Union Geopandas improve performanceGetting error while importing GeoPandas in Python consoleShapely : cascaded_union seems to return too many polygonsGeoPandas equivalent of ST_MakePolygonWrite a Spatial Join output to a CSV?

Avoiding the "not like other girls" trope?

How much of data wrangling is a data scientist's job?

How badly should I try to prevent a user from XSSing themselves?

Apex Framework / library for consuming REST services

Extract rows of a table, that include less than x NULLs

What is the most common color to indicate the input-field is disabled?

Assassin's bullet with mercury

How to show a landlord what we have in savings?

How could indestructible materials be used in power generation?

Can I run a new neutral wire to repair a broken circuit?

What about the virus in 12 Monkeys?

Why didn't Miles's spider sense work before?

Is there a hemisphere-neutral way of specifying a season?

What mechanic is there to disable a threat instead of killing it?

What is a romance in Latin?

Can a virus destroy the BIOS of a modern computer?

CAST throwing error when run in stored procedure but not when run as raw query

How can saying a song's name be a copyright violation?

Can compressed videos be decoded back to their uncompresed original format?

How dangerous is XSS?

What's the in-universe reasoning behind sorcerers needing material components?

Should I cover my bicycle overnight while bikepacking?

Can my sorcerer use a spellbook only to collect spells and scribe scrolls, not cast?

Arrow those variables!



Keep the holes in the polygon while converting Polygon 25D to Polygon 2D using Python


Geopandas/folium map not displaying?Performing Spatial Join / match Points from dataframe to polygons using Python?geopandas geometry.within problem ValueError: Geometry column cannot contain mutiple geometry typesSpatial join in Geopandas is only finding itselfGeopandas not recognizing geometry typeOverlay Union Geopandas improve performanceGetting error while importing GeoPandas in Python consoleShapely : cascaded_union seems to return too many polygonsGeoPandas equivalent of ST_MakePolygonWrite a Spatial Join output to a CSV?













0















I converted my input file with geometry Polygon 25D to an output file with geometry Polygon using the below code.
As you can see the below image. Problem is that the holes in the polygon is vanished after applying this code. I would like to keep the holes in the polygon.
To achieve that how can I edit or modify the below code?



def convert_3D_2D(geometry):
import geopandas as gp
new_geo = []
for p in geometry:
if p.has_z:
if p.geom_type == 'Polygon':
lines = [xy[:2] for xy in list(p.exterior.coords)]
new_p = Polygon(lines)
new_geo.append(new_p)
elif p.geom_type == 'MultiPolygon':
new_multi_p = []
for ap in p:
lines = [xy[:2] for xy in list(ap.exterior.coords)]
new_p = Polygon(lines)
new_p
new_multi_p.append(new_p)
new_multi_p
new_geo.append(MultiPolygon(new_multi_p))
new_geo
return new_geo


df = gpd.GeoDataFrame.from_file('input_file_25D.shp')
df.geometry = convert_3D_2D(df.geometry)
df.to_file('output_file_2D.shp')



figure










share|improve this question




























    0















    I converted my input file with geometry Polygon 25D to an output file with geometry Polygon using the below code.
    As you can see the below image. Problem is that the holes in the polygon is vanished after applying this code. I would like to keep the holes in the polygon.
    To achieve that how can I edit or modify the below code?



    def convert_3D_2D(geometry):
    import geopandas as gp
    new_geo = []
    for p in geometry:
    if p.has_z:
    if p.geom_type == 'Polygon':
    lines = [xy[:2] for xy in list(p.exterior.coords)]
    new_p = Polygon(lines)
    new_geo.append(new_p)
    elif p.geom_type == 'MultiPolygon':
    new_multi_p = []
    for ap in p:
    lines = [xy[:2] for xy in list(ap.exterior.coords)]
    new_p = Polygon(lines)
    new_p
    new_multi_p.append(new_p)
    new_multi_p
    new_geo.append(MultiPolygon(new_multi_p))
    new_geo
    return new_geo


    df = gpd.GeoDataFrame.from_file('input_file_25D.shp')
    df.geometry = convert_3D_2D(df.geometry)
    df.to_file('output_file_2D.shp')



    figure










    share|improve this question


























      0












      0








      0








      I converted my input file with geometry Polygon 25D to an output file with geometry Polygon using the below code.
      As you can see the below image. Problem is that the holes in the polygon is vanished after applying this code. I would like to keep the holes in the polygon.
      To achieve that how can I edit or modify the below code?



      def convert_3D_2D(geometry):
      import geopandas as gp
      new_geo = []
      for p in geometry:
      if p.has_z:
      if p.geom_type == 'Polygon':
      lines = [xy[:2] for xy in list(p.exterior.coords)]
      new_p = Polygon(lines)
      new_geo.append(new_p)
      elif p.geom_type == 'MultiPolygon':
      new_multi_p = []
      for ap in p:
      lines = [xy[:2] for xy in list(ap.exterior.coords)]
      new_p = Polygon(lines)
      new_p
      new_multi_p.append(new_p)
      new_multi_p
      new_geo.append(MultiPolygon(new_multi_p))
      new_geo
      return new_geo


      df = gpd.GeoDataFrame.from_file('input_file_25D.shp')
      df.geometry = convert_3D_2D(df.geometry)
      df.to_file('output_file_2D.shp')



      figure










      share|improve this question
















      I converted my input file with geometry Polygon 25D to an output file with geometry Polygon using the below code.
      As you can see the below image. Problem is that the holes in the polygon is vanished after applying this code. I would like to keep the holes in the polygon.
      To achieve that how can I edit or modify the below code?



      def convert_3D_2D(geometry):
      import geopandas as gp
      new_geo = []
      for p in geometry:
      if p.has_z:
      if p.geom_type == 'Polygon':
      lines = [xy[:2] for xy in list(p.exterior.coords)]
      new_p = Polygon(lines)
      new_geo.append(new_p)
      elif p.geom_type == 'MultiPolygon':
      new_multi_p = []
      for ap in p:
      lines = [xy[:2] for xy in list(ap.exterior.coords)]
      new_p = Polygon(lines)
      new_p
      new_multi_p.append(new_p)
      new_multi_p
      new_geo.append(MultiPolygon(new_multi_p))
      new_geo
      return new_geo


      df = gpd.GeoDataFrame.from_file('input_file_25D.shp')
      df.geometry = convert_3D_2D(df.geometry)
      df.to_file('output_file_2D.shp')



      figure







      geopandas python-3






      share|improve this question















      share|improve this question













      share|improve this question




      share|improve this question








      edited 2 days ago









      Vince

      14.8k32849




      14.8k32849










      asked 2 days ago









      roshualineroshualine

      297




      297




















          0






          active

          oldest

          votes












          Your Answer








          StackExchange.ready(function()
          var channelOptions =
          tags: "".split(" "),
          id: "79"
          ;
          initTagRenderer("".split(" "), "".split(" "), channelOptions);

          StackExchange.using("externalEditor", function()
          // Have to fire editor after snippets, if snippets enabled
          if (StackExchange.settings.snippets.snippetsEnabled)
          StackExchange.using("snippets", function()
          createEditor();
          );

          else
          createEditor();

          );

          function createEditor()
          StackExchange.prepareEditor(
          heartbeatType: 'answer',
          autoActivateHeartbeat: false,
          convertImagesToLinks: false,
          noModals: true,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: null,
          bindNavPrevention: true,
          postfix: "",
          imageUploader:
          brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
          contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
          allowUrls: true
          ,
          onDemand: true,
          discardSelector: ".discard-answer"
          ,immediatelyShowMarkdownHelp:true
          );



          );













          draft saved

          draft discarded


















          StackExchange.ready(
          function ()
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f317314%2fkeep-the-holes-in-the-polygon-while-converting-polygon-25d-to-polygon-2d-using-p%23new-answer', 'question_page');

          );

          Post as a guest















          Required, but never shown

























          0






          active

          oldest

          votes








          0






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes















          draft saved

          draft discarded
















































          Thanks for contributing an answer to Geographic Information Systems Stack Exchange!


          • Please be sure to answer the question. Provide details and share your research!

          But avoid …


          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.

          To learn more, see our tips on writing great answers.




          draft saved


          draft discarded














          StackExchange.ready(
          function ()
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f317314%2fkeep-the-holes-in-the-polygon-while-converting-polygon-25d-to-polygon-2d-using-p%23new-answer', 'question_page');

          );

          Post as a guest















          Required, but never shown





















































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown

































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown







          Popular posts from this blog

          Magento 1.9 observer on only new product added ( no updated product )Admin Create New Order ObserverObserver checkout_cart_product_add_after get product quantity in cartHow can I differentiate between product creation and saving existing products?How to get the original price Product SaveBeforeMagento Observer for Admin Sales order Address updateHow do i get quantity in a product_save_after event observer magento2?Using observer event to get updated customer infoSave product observer - detect availabilityEvent not triggered when new product added in admin : Magento 2magento 1.9 observer on quantity product change

          Export of reprojected layer from GEE fails Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?Iterate over features and years to export images for each feature-yearEarth engine reproject: why does reprojecting a pixel at 30m scale gives me pixels with an area of ~550?Exporting slope raster for entire country in GEEGEE Landsat SR year composite // mask cloud / shadow out w/ other method than quality pixel"export region contains no valid (un-masked) pixels - Google Earth EgineGoogle Earth Engine, how to distinguish between rivers/streams and ponds/lakes in a water maskPython script tool fails when processing rainfall data from Google Earth EngineGoogle Earth Engine Error: Number of pixels requested from Image.load exceeds the maximum allowedCloudfree images in small area from Sentinel-2Using Image exportToDrive in Google Earth Engine?

          Adding contour labels to SpatialLinesDataFrame or sf object? The Next CEO of Stack OverflowTearing of SpatialLinesDataFrame and problem with gClip functionAutomatically adjusting map labelsExtracting data.frame from simple features object in RR convert sf object back to SpatialPolygonsDataFrameSimplify polygons of sf objectlines of a SpatialLinesDataFrame corresponding to pointsCropping sf object in R?Pattern matching with grepl on SpatialPolygonsDataFrame object in RHow to create and empty SpatialLinesDataFrame in RWarning message in raster::rotate(): “this does not look like an appropriate object for this function”