如果知道中心点，如何计算网格的角点

问题描述

我正在处理一个坐标为 wsg84 格式的文件。从这个文件中，我使用以下代码计算了笛卡尔坐标：

import numpy as np
import pyproj
from netCDF4 import Dataset
_projections = {}


def unproject(z,l,x,y):
    if z not in _projections:
        _projections[z] = pyproj.Proj(proj='utm',zone=z,ellps='wgs84')
    if l < 'N':
        y -= 10000000
    lng,lat = _projections[z](x,y,inverse=True)
    return (lng,lat)

def calculate_grid_deFinition(in_path):
    """
    Calculate the latitude and longitude coordinates of the cell.

    :param in_path: Path to the file that contains all the information.
    :type in_path: str

    :return: Latitudes array,Longitudes array,Latitude interval,Longitude interval.
    :rtype: numpy.array,numpy.array,float,float
    """
    with Dataset(in_path) as f:
        utm_lat = f.variables['LAT'][:]
        utm_lon = f.variables['LON'][:]
        geolat = np.zeros((len(utm_lon),len(utm_lat)))
        geolon = np.zeros((len(utm_lon),len(utm_lat)))
        for x,u_lon in enumerate(utm_lon):
            for y,u_lat in enumerate(utm_lat):
                lon_Now,lat_Now = unproject(32,'N',u_lon,u_lat)
                geolat[x,y] = lat_Now
                geolon[x,y] = lon_Now

        return geolat,geolon

inpfile = 'path/to/filename.nc'
geolat,geolon = calculate_grid_deFinition(inpfile)

现在我还想计算网格单元格的角点。假设一个规则的矩形网格，使用以下代码将非常简单：

def create_bounds(coordinates):
    interval_x = coordinates[1,:] - coordinates[0,:]
    interval_y = coordinates[:,1] - coordinates[:,0]
    interval_x = np.tile(interval_x,(len(interval_y),1))
    interval_y = np.tile(interval_y,(interval_x.shape[-1],1)).T

    coords_left = coordinates - interval_x / 2
    coords_right = coordinates + interval_x / 2
    coords_bottom = coordinates - interval_y / 2
    coords_top = coordinates + interval_y / 2

    bound_coords = np.dstack((coords_left,coords_top,coords_right,coords_bottom))
    return bound_coords


lat_bnds = create_bounds(geolat)
lon_bnds = create_bounds(geolon)

但是，这是不对的，因为网格倾斜并导致重叠。我将不胜感激任何线索。 为了您的信息，您可以从这个 file

下载原始网格坐标

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）