问题描述:
当我们计算得到了某一空间范围内的格点值,想可视化或仅对某一指定区域可视化。
输入:大范围的经纬度+相应待可视化的值(excel文件),示例如:
输出:利用克里金(Kriging)插值对指定地理区域的进行填充可视化。如:我想要中国范围的可视化结果:
代码实现:
import pandas as pd
from pykrige.ok import OrdinaryKriging
import copy
import plotnine
from plotnine import *
from matplotlib.colors import Colormap
import geopandas as gpd
import shapefile
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.io.shapereader as shpreader
import cmaps
from matplotlib.path import Path
import numpy as np
from matplotlib.patches import PathPatch