python - Bokeh 未从 GeoJson 正确渲染多边形(岛屿)-6ren

python - Bokeh 未从 GeoJson 正确渲染多边形(岛屿)

转载作者：行者123 更新时间：2023-11-28 21:08:18

我正在 Bokeh 上渲染等值线图。我的地理数据是一个包含多边形和多边形的 GeoJSON。

我在渲染多边形时遇到问题:

如果我提取一个要素的所有几何图形(例如，四个岛屿，在一个列表中)，它们的图不会在图形之间“切割”，它们看起来都是一个。它显示了一些“蜘蛛网”的东西，无序地穿过所有点。

如果我为岛屿创建一个列表(我假设这是正确的工作方式)，Bokeh 不会绘制任何东西。甚至没有网格(只有工具栏)......并且没有显示任何错误。

可能是函数“obtCoordMultipoligono”的输出有问题。

我已经在岛上搜索过示例，但没有任何帮助。

提前致谢。

更新:我添加我的片段。它们是原始的片段，但具有功能性。这个想法是像 BokehGallery 上的“德克萨斯失业率”这样的输出，但我的国家上面有岛屿。

我在 argentina.json 中的 GeoJSON(仅提取 1 个多边形；我没有多边形问题):

    {
    "type": "FeatureCollection",
    "features": [

    {
      "geometry": {
        "type": "MultiPolygon",
        "coordinates": [
          [
            [
              [
                -59.68266601562502,
                -52.231640624999976
              ],
              [
                -59.74658203124997,
                -52.25087890624999
              ],
              [
                -59.76445312499996,
                -52.2421875
              ],
              [
                -59.784863281249955,
                -52.2046875
              ],
              [
                -59.78593749999999,
                -52.156152343749966
              ],
              [
                -59.79331054687498,
                -52.134179687500016
              ],
              [
                -59.75322265624999,
                -52.14140624999998
              ],
              [
                -59.681005859375034,
                -52.18007812499995
              ],
              [
                -59.68266601562502,
                -52.231640624999976
              ]
            ]
          ],
          [
            [
              [
                -58.438818359375006,
                -52.011035156249974
              ],
              [
                -58.432714843750006,
                -52.09902343749996
              ],
              [
                -58.512841796874966,
                -52.071093750000045
              ],
              [
                -58.54140625000002,
                -52.02841796874996
              ],
              [
                -58.49707031249997,
                -51.99941406250001
              ],
              [
                -58.46054687499998,
                -52.0015625
              ],
              [
                -58.438818359375006,
                -52.011035156249974
              ]
            ]
          ],
          [
            [
              [
                -61.01875,
                -51.7857421875
              ],
              [
                -60.94726562499997,
                -51.79951171875005
              ],
              [
                -60.87597656250003,
                -51.79423828125004
              ],
              [
                -60.91616210937494,
                -51.89697265625001
              ],
              [
                -60.94755859374996,
                -51.94628906250002
              ],
              [
                -61.031982421875,
                -51.94248046875004
              ],
              [
                -61.11577148437493,
                -51.87529296875003
              ],
              [
                -61.14501953125003,
                -51.83945312500001
              ],
              [
                -61.05166015625002,
                -51.81396484374997
              ],
              [
                -61.01875,
                -51.7857421875
              ]
            ]
          ],
          [
            [
              [
                -60.11171875000002,
                -51.39589843749998
              ],
              [
                -60.24882812499996,
                -51.39599609375
              ],
              [
                -60.27587890624997,
                -51.36318359374997
              ],
              [
                -60.275341796874955,
                -51.28056640625002
              ],
              [
                -60.17138671875,
                -51.273437499999986
              ],
              [
                -60.06982421875,
                -51.307910156249996
              ],
              [
                -60.07646484374993,
                -51.34257812500004
              ],
              [
                -60.11171875000002,
                -51.39589843749998
              ]
            ]
          ]

        ]
      },
      "type": "Feature",
      "properties": {
        "perimeter": 0,
        "vista": "mim",
        "provincia": "Islas Malvinas",
        "objectid": 24,
        "prov": 0,
        "bounds": [
          0,
          0
        ],
        "provif3_": 27.0,
        "ogc_fid": 26,
        "provif3_id": 26.0
      }
    }
    ]
    }

我在 PBIArg.csv 中的数据:

24,AR-V,Islas,13245

我的代码:

<!-- language: lang-py -->

import json,pprint,csv
pp = pprint.PrettyPrinter(indent=4)
from bokeh.io import output_file, show
from bokeh.models import HoverTool
from bokeh.plotting import figure, show, output_file, ColumnDataSource
import pandas as pd
from osgeo import ogr

fname = r'islas.json' # constante hasta conseguir algo mejor
dname = r'PBIArg.csv' # variable estadística a graficar.



paleta = ["#FFF5F0", "#FEE0D2", "#FCBBA1", "#FC9272", "#FB6A4A", "#EF3B2C", "#CB181D", "#99000D"]


def colorante(rate,max_value,min_value,paleta):
    try:
        intensidad = int(float(len(paleta)-1) * float(rate - min_value) / float(max_value - min_value))
        return paleta[intensidad]
    except:
        intensidad = int(float(len(paleta)-1) * float(rate - min_value) / float(max_value - min_value))
        return paleta[intensidad]


def obtCoordMultipoligono(pcia):
    mpoly = ogr.CreateGeometryFromJson(pcia)

    print('pcia-MPOLY tiene esta cantidad de islas: ', mpoly.GetGeometryCount()) 
    coordX,coordY = [],[]

    # idx poly mpoly
    for idx, poly in enumerate(mpoly): #itero mpoly
        print('POLY tiene esta cantidad de islas: ', poly.GetGeometryCount()) 
        innerX,innerY = [],[]
        # ind 
        for ind in range(0, poly.GetGeometryCount()): #itero poly
            innerPoly = poly.GetGeometryRef(ind)
            print('INNERPOLY tiene esta cantidad de PUNTOS: ', innerPoly.GetPointCount()) 
            lastX,lastY = [],[]
            for i in range(0, innerPoly.GetPointCount()): #itero innerpoly
                # GetPoint returns a tuple not a Geometry
                punto = innerPoly.GetPoint(i)
                print('pto obtenido en X',punto[0])
                # Asigno a lista local
                lastX.append(punto[0])
                lastY.append(punto[1])
            print('LastX:')
            pp.pprint(lastX)
            innerX.append(lastX)
            innerY.append(lastY)
        print('InnerX:')
        pp.pprint(innerX)
        coordX.append(innerX)
        coordY.append(innerY)
        print('CoordX:')
        pp.pprint(coordX)
    dictCoord = dict(coordX=coordX,coordY=coordY)
    print('DictCoord:')
    pp.pprint(dictCoord)
    return dictCoord

############  MAIN  ##################
######## Leo csv estadísticas ########

with open(dname, 'r') as f:
    '''Leo el CSV, creo diccio de pciaID: estadísticaPcial.
    Abajo busco la estadísticaPcial max y min para luego calcular los colores'''
    max_value, min_value = 0,0
    datos = {}
    for row in csv.reader(f):
        estadistica = int(row[3])
        datos[row[0]] = estadistica
        if estadistica > max_value:
            max_value = estadistica
        if estadistica < min_value:
            min_value = estadistica

######## Leo geojson ########            
with open(fname, 'r') as f:
    geojson = f.read()
    geoDict = json.loads(geojson)

######## Parseo geojson ########

dictArg = {}
for pcia in geoDict['features']:
    pciaID = str(pcia['properties']['objectid'])
    nombrePcia = pcia['properties']['provincia']
    coordX = []
    coordY = []

    if pcia['geometry']['type'] == 'Polygon':
        for punto in pcia['geometry']['coordinates'][0]:
            if len(punto) == 2:
                coordX.append(punto[0])
                coordY.append(punto[1])
    elif pcia['geometry']['type'] == 'MultiPolygon': 
        multiJSON = json.dumps(pcia['geometry'])
        dictCoord = obtCoordMultipoligono(multiJSON)
        # print(dictCoord)
        coordX = dictCoord['coordX']
        coordY = dictCoord['coordY']

    # Handling states without data
    try:
        info=int(datos[pciaID])
    except KeyError:
        info = 0    
    color = colorante(info,max_value,min_value,paleta)
    dictPcia = dict(nombre=nombrePcia,coordX=coordX,coordY=coordY, info=info,color=color)
    dictArg[pciaID] = dictPcia
    print('dict',dictArg['19'])

######## saco coord de las pcias ########
provincias = {
    codPcia: nombrePcia for codPcia, nombrePcia in dictArg.items()
}

# print(provincias)
pcia_xs = [provincia['coordX'] for provincia in provincias.values()]
pcia_ys = [provincia['coordY'] for provincia in provincias.values()]

nombres_provincias = [provincia['nombre'] for provincia in provincias.values()]

######## Saco estadísticas de las pcias ########
provincias_datos = [provincia['info'] for provincia in provincias.values()]

######## Coloreo el mapa a nivel datos ########

provincias_colores = [provincia['color'] for provincia in provincias.values()]

source = ColumnDataSource(data=dict(
    x=pcia_xs,
    y=pcia_ys,
    color=provincias_colores,
    nombre=nombres_provincias,
    dato=provincias_datos,
))

TOOLS="pan,wheel_zoom,box_zoom,reset,hover,save"

p = figure(title="PBI de Argentina por provincia", tools=TOOLS)

p.patches('x', 'y', source=source,
          fill_color='color', fill_alpha=0.9,
          line_color='#767676', line_width=1.5)

hover = p.select_one(HoverTool)
hover.point_policy = "follow_mouse"
hover.tooltips = [
    ("State:", "@nombre"),
    ("Nº:", "@dato"),

]

output_file("PBIar.html", title="Testing islands in bokeh")

show(p)

最佳答案

如果您使用 Bokeh 的 GeoJSONDataSource，而不是遵循原始的“德克萨斯示例”，您可以显着简化您的代码。

使用您的 geojson 的精简示例可能如下所示:

from bokeh.io import show, output_notebook, output_file
from bokeh.models import (
    GeoJSONDataSource,
    HoverTool,
    LinearColorMapper
)
from bokeh.plotting import figure
from bokeh.palettes import Viridis6

with open(r'argentina.geojson', 'r') as f:
    geo_source = GeoJSONDataSource(geojson=f.read())


color_mapper = LinearColorMapper(palette=Viridis6)

TOOLS = "pan,wheel_zoom,box_zoom,reset,hover,save"

p = figure(title="Argentina", tools=TOOLS, x_axis_location=None, y_axis_location=None, width=500, height=300)
p.grid.grid_line_color = None

p.patches('xs', 'ys', fill_alpha=0.7, fill_color={'field': 'objectid', 'transform': color_mapper}, 
          line_color='white', line_width=0.5, source=geo_source)


hover = p.select_one(HoverTool)
hover.point_policy = "follow_mouse"
hover.tooltips = [("Provincia:", "@provincia")]

output_file("PBIar.html", title="Testing islands in bokeh")

show(p)

输出结果如下:

编辑:

这是使用下面评论中提到的整个 geojson 时的输出。

然后放大南部的岛屿:

更新:可变数据的功能

我添加此功能是为了以交互方式编辑 geoJSON，具体取决于传递的 JSON 数据。

geoJSON 现在具有“数据”属性和州国际代码(属性:“ISO_3166-2”)。

JSON数据是这样的:

{
"AR-A": "7",
"AR-B": "53",
"AR-C": "46"
}

该函数读取 geoJSON 并分配数据:

def asignDataToStates(geo,data):
    for pcia in geo['features']:
        codPcia = str(pcia['properties']['ISO_3166-2'])
        
        if codPcia in data.keys():
            if data.values() != 0:
                pcia['properties']['data'] = data[codPcia]
    dataJson = json.dumps(geo,ensure_ascii=True)
    return dataJson

关于python - Bokeh 未从 GeoJson 正确渲染多边形(岛屿)，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/40226189/

文章推荐： javascript - 我在哪里存储 src 方法的文件？

文章推荐： python - 开始使用我自己的数据集启动 google-cloud-ml

javascript - 如果输入 a 或 b 正确/正确，我如何执行操作？
这个问题已经有答案了: How to do case insensitive string comparison? (23 个回答) 已关闭 3 年前。用户在我的输入栏中写入“足球”，然后执行第 6
javascript - 字符 id= + 是 + 正确= + 正确不正确...我怎样才能使它成为 javascript 中的字符串
啊，不习惯 javascript 中的字符串。 character_id= + id + correct= + correctOrIncorrect 这就是我需要制作成字符串的内容。如果您无法猜测字符
javascript - jQuery计算价格不起作用(正确)
$(function() { var base_price = 0; CalculatePrice(); $(".math1").on('change', function(e) { Calc
kubernetes - 将Spinnaker部署到Spinnaker将管理的同一kubernetes集群是否安全/正确？
我找不到任何文章回答问题:将Spinnaker部署到Spinnaker将管理的同一Kubernetes集群是否安全/正确？我主要是指生产，HA部署。最佳答案我认为Spinnaker和Kuberne
c++ - 正确/快速的方法来更改命令行Qt5源内部版本的配置
我正在使用MSVC在Windows上从源代码(官方源代码发布，而不是从仓库中)构建Qt5(Qt 5.15.0)。我正在设置环境。变量，依赖项等，然后运行具有1600万个选项的configure，最后
java - 计数时数组越界[正确]
我需要打印一个包含重复单词的数组。我的数组已经可以工作，但我不知道如何正确计算单词数。我已经知道，当我的索引计数器 (i) 为 49 时，并且当 (i) 想要计数到 50 时，我会收到错误，但我不知道
javascript - 正确/错误取决于屏幕尺寸动态？
我正在遵循一个指南，该指南允许 Google map 屏幕根据屏幕尺寸禁用滚动。我唯一挣扎的部分是编写一个代码，当我手动调整屏幕大小时动态更改 True/False 值。这是我按照说明操作的网站，但
java - 未调用子类中的方法(正确)
我有一个类“FileButton”。它的目的是将文件链接到 JButton，FileButton 继承自 JButton。子类继承自此以使用链接到按钮的文件做有用的事情。 JingleCardButt
php - 如何仅显示来自好友列表的帖子。 (正确)
我的 friend 数组只返回一个数字而不是所有数字。 ($myfriends = 3) 应该是…… ($myfriends = 3 5 7 8 9 12). 如果我让它进入 while 循环……整个
html - 在这种情况下使用整数作为类名是否可以接受/正确
这个问题在这里已经有了答案: Is there a workaround to make CSS classes with names that start with numbers valid?
javascript - 在窗口更改时自动调整元素大小(正确)
我正在制作一个 JavaScript 函数，当调整窗口大小时，它会自动将 div 的大小调整为与窗口相同的宽度/高度。该功能非常基本，但我注意到在调整窗口大小时出现明显的“绘制”滞后。在 JS fi
javascript - 删除导航栏的类 - 正确
此问题的基本视觉效果可在 http://sevenx.de/demo/bootstrap-carousel/inc.carousel/tabbed-slider.html 获得。 - 如果你想看一看。
c - 从将其内存分配给同一函数的函数返回字符串是否安全/正确？
我明白，如果我想从函数返回一个字符串文字或一个数组，我应该将其声明为静态的，这样当被调用的函数被返回时，内容就不会“消亡”。但我的问题是，当我在函数内部使用 malloc 分配内存时会怎样？在下面
mysql - 正确/错误值的适当数据字段类型？
在 mySQL 数据库中存储 true/false/1/0 值最合适(读取数据消耗最少)的数据字段是什么？我以前使用过一个字符长的 tinyint，但我不确定它是否是最佳解决方案？谢谢! 最佳答案
c++ - 正确，有效地读取文件
我想一次读取并处理CSV文件第一行中的条目(例如打印)。我假设使用Unix风格的\n换行符，没有条目长度超过255个字符，并且(现在)在EOF之前有一个换行符。这意味着它是fgets()后跟strto
c++ - “正确”无符号整数比较
所以，我们都知道 -1 > 2u == true 的 C/C++ 有符号/无符号比较规则，并且我有一种情况，我想有效地实现“正确”比较。我的问题是，考虑到人们熟悉的尽可能多的架构，哪种方法更有效。显
Java异常处理：如何写出“正确”但被编译器认为有语法错误的程序
**摘要：**文章的标题看似自相矛盾。本文分享自华为云社区《Java异常处理：如何写出“正确”但被编译器认为有语法错误的程序》，作者： Jerry Wang 。文章的标题看似自相矛盾，然而我在“正
r - 进行按行替换的“正确”方法
我有一个数据框，看起来像: dataDemo % mutate_each(funs(ifelse(. == '.', REF, as.character(.))), -POS) # POS REF
text - VBScript 正确/重新格式化带分隔符的文本文件？
有人可以帮助我使用 VBScript 重新格式化/正确格式化带分隔符的文本文件吗？我有一个文本文件 ^分界如下: AGREE^NAME^ADD1^ADD2^ADD3^ADD4^PCODE^BAL^A
java - 语言认证以及诸如适当、正确、合法等术语的使用
就目前而言，这个问题不适合我们的问答形式。我们希望答案得到事实、引用或专业知识的支持，但这个问题可能会引起辩论、争论、投票或扩展讨论。如果您觉得这个问题可以改进并可能重新打开，visit the he

行者123

个人简介

我是一名优秀的程序员,十分优秀！

作者热门文章

滴滴打车优惠券免费领取

全站热门文章

首页

博学

6Ren·AI

商城

python - Bokeh 未从 GeoJson 正确渲染多边形(岛屿)

编辑:

更新:可变数据的功能