When reading a PointZ shapefile (<a href="https://www.arcgis.com/home/item.html?id=1f8

Hi <a class="user-mention notranslate" data-hovercard-type="user" data-hovercard-url="

PointZ shapefiles are not properly read about nettopologysuite.io.esri HOT 2 CLOSED

nettopologysuite commented on September 25, 2024

PointZ shapefiles are not properly read

from nettopologysuite.io.esri.

Comments (2)

KubaSzostak commented on September 25, 2024 1

Hi @ShinNoNoir, I checked the layer you provided. It has a geometry of type Point, not PointZM. In that case if you try to read Z coordinate from geometry then you get NaN. There is a valid Z value, but it's not a coordinate, but an attribute.

That's how it looks in QGIS:

You can read the same data using C# code:

var features = Shapefile.ReadAllFeatures("2021_Paspunten_BM5.shp");
foreach (var feature in features)
{
    Console.WriteLine($"Geometry: {feature.Geometry}");
    Console.WriteLine($"XYMZ: {feature.Geometry.Coordinate.X} {feature.Geometry.Coordinate.Y} {feature.Geometry.Coordinate.M} {feature.Geometry.Coordinate.Z}");
                Console.WriteLine("Attributes:");
    foreach( var name in feature.Attributes.GetNames())
    {
        var value = feature.Attributes[name];
        Console.WriteLine($"  {name}: {value}");
    }
    Console.WriteLine();
}

Output:

Geometry: POINT (172354.84 534578.751)
XYMZ: 172354.84 534578.751 NaN NaN
Attributes:
  Puntnummer: 580_1010
  X: 172354.84
  Y: 534578.751
  Z: -3.055
  Bron: 2021_paspuntcoordinaten_05_80_hrl
  POINT_X: 
  POINT_Y: 

...

I hope this explanation solves your issue. If you need more support on this issue provide further details, please.

from nettopologysuite.io.esri.