Visualize point PMTiles layer

[jinha]

I have a point PMTiles layer uploaded on my server: https://lidar.digitalforestry.org/vt/in/2022/IN2022-Address_Points.pmtiles
When I tried to visualize the point layer using add_pmtiles, it does not show anything. I checked the tutorial and it only showed examples for polygons and lines only. Is there anyway to visualize a point PMTiles layer?

[giswqs]

It should work. Make sure you use the correct the layer name.

https://protomaps.github.io/PMTiles/?url=https%3A%2F%2Flidar.digitalforestry.org%2Fvt%2Fin%2F2022%2FIN2022-Address_Points.pmtiles#map=7.46/39.614/-86.797

[jinha]

I also checked with the PMTiles viewer and confirmed that it is working okay. The problem is that when I tried to load it using leafmap, it is not working. Below is the code I used. I also went to "foliumap" github repository to see if they have any point layer example, but it did not have any point layer example neither.

"""
import leafmap.foliumap as leafmap
ap_url = 'https://lidar.digitalforestry.org/vt/in/2022/IN2022-Address_Points.pmtiles'
ap_metadata = leafmap.pmtiles_metadata(ap_url)
print(f"layer names: {ap_metadata['layer_names']}")
print(f"bounds: {ap_metadata['bounds']}")

m = leafmap.Map(center=[37, -88], zoom=6, height='800px')
m.add_basemap('CartoDB.DarkMatter')
m.add_basemap('Esri.WorldImagery', show=False)

ap_style = {
"version": 8,
"sources": {
"ap_source": {
"type": "vector",
"url": "pmtiles://" + ap_url,
"attribution": 'PMTiles',
}
},
"layers": [
{
"id": "address_point",
"source": "ap_source",
"source-layer": "in2022_address_points",
"type": "point",
"paint": {"fill-color": "#3388ff", "fill-opacity": 0.5},
},
],
}

m.add_pmtiles(
ap_url, style=ap_style, name='Address Points', overlay=True, show=True, zoom_to_layer=False
)

m
"""

[giswqs]

Try this one. The style can be a mapbox style.

import leafmap.foliumap as leafmap

m = leafmap.Map()

style = {
    "version": 8,
    "sources": {
        "example_source": {
            "type": "vector",
            "url": "pmtiles://" + url,
            "attribution": 'PMTiles',
        }
    },
    "layers": [
        {
            "id": "in2022_address_points",
            "source": "example_source",
            "source-layer": "in2022_address_points",
            "type": "circle",
            'paint': {
                'circle-radius': 3,
                'circle-color': '#0000ff'
            }
        },
    ],
}

m.add_pmtiles(
    url, name='PMTiles', style=style, overlay=True, show=True, zoom_to_layer=True
)
m

[giswqs]

More styles available here: https://docs.mapbox.com/style-spec/reference/layers/#circle

[jinha]

Thanks. It is working now. 😄