这玩意几年才更新一次，每次更新，每次更新都忘了该怎么算，记一下，以后用。
JMA 震度信息三个层级，区域 -> 城市 -> 观测点，官方的说法叫
“地震情報／細分区域”・”気象・地震・火山情報／市町村等”・”震度観測点”
这个东西可以在気象庁防災情報XMLフォーマット 技術資料下载到
表 24 「AreaForecastLocalE ・ AreaInformationCity ・ PointSeismicIntensity コード表」就是我们需要的数据。
绘制震度图时不需要精确到震度观测点，所以可以写一个脚本对前两个层级聚类。这两个层级也分别对应震度速报和各地区震度信息的电文中的区域代码。
脚本略。
区域这一层就是对都道府县的略微细分，共188个区域，这个几乎不调整，GitHub 上有人算好了坐标信息，在这里。
市町村这一层级隔几年会更新一次，主要都是因为行政区域调整导致的。2024年现在，1894个区域。
这一层级的地理数据可以在 気象データ高度利用ポータルサイト 下载
下载得到的是 shapefile 格式
用 QGIS 打开
图层右键打开属性表，全选要素。
几何图形工具求各要素质心
勾选仅选中的要素，运行。等一会会
质心会计算成一个新的图层，右键导出为 GeoJSON 格式，然后估计就会处理了。
这里要注意，JMA 提供的 GIS 数据里有 1911 个区域，其中没有分配 regionCode 的是无效的，需要过滤。
（自从 JMA 开始提供 GIS 数据之后，相关数据的计算和展示变得简单了许多）