zip = Take[EntityValue[Entity["AdministrativeDivision", {"NewYork", "UnitedStates"}],
       EntityProperty["AdministrativeDivision", "ZIPCodes"]], 10];
list = {#, EntityValue[#, {"Polygon", "Population"}]} & /@ zip;

list2 = Complement[list,
    Cases[list, {_, {Missing["NotApplicable"], _}}]] /.
   GeoPosition -> Sequence;
list2[[All, 2, 1]] = First[#] & /@ list2[[All, 2, 1]];
centroids =
 RegionCentroid[DiscretizeGraphics[#]] & /@ list2[[All, 2, 1]]
(*{{41.2703, -71.9886}, {40.7505, -73.9974}, {40.716, -73.9865},
{40.7317, -73.9891}, {40.6944, -74.0158}, {40.7061, -74.0086},
{40.7096, -74.013}, {40.7137, -74.0083}}*)

data = {centroids[[#, 1]], centroids[[#, 2]],
     QuantityMagnitude[list2[[#, 2, 2]]]} & /@ Range@Length@centroids;
clusters = FindClusters[%];
ListPointPlot3D[clusters, PlotStyle -> PointSize[Large]]