使用 Shapely 从四个经纬度角创建多边形

如何解决使用 Shapely 从四个经纬度角创建多边形

我正在尝试使用来自臭氧监测仪器 (OMI) 的数据文件,并将来自该来源的数据与监测类似数据的地面仪器进行比较。

最终,我想确定是否在特定多边形内。不过,要做到这一点,我需要创建我想象的多边形。

我有以下变量

latmat

latmat[1] = array([-62.2546,-62.371,-62.4871,-62.6032],dtype=float32)
latmat[2] = array([-62.7195,-62.8356,-62.9519,-63.0676],dtype=float32)

lonmat

lonmat[1] = array([135.579,135.606,135.633,135.66 ],dtype=float32)
lonmat[2] = array([135.688,135.717,135.745,135.774],dtype=float32)

latmat 和 lonmat 当然都比两个指数大得多,但我试图保持这些简单。这些值分别代表一个四角卫星像素的纬度和经度点。

最重要的是,我还有表单中的数据

data[1] = 1E+15
data[2] = 3E+15

如何使用 geopandas 创建一个地理数据框,将 4 个纬度/经度点的每个组合识别为多边形?

解决方法

您需要为此使用 zip()

from shapely.geometry import Polygon    
import numpy as np

latmat1 = np.array([-62.2546,-62.371,-62.4871,-62.6032],dtype=np.float32)
latmat2 = np.array([-62.7195,-62.8356,-62.9519,-63.0676],dtype=np.float32)
lonmat1 = np.array([135.579,135.606,135.633,135.66 ],dtype=np.float32)
lonmat2 = np.array([135.688,135.717,135.745,135.774],dtype=np.float32)

pol1 = Polygon([(i,j) for i,j in zip(lonmat1,latmat1 ) ])
pol2 = Polygon([(i,j in zip(lonmat2,latmat2 ) ])
,

这对我有用,我认为您的数据是如何存储的。

from shapely.geometry import Polygon    
import numpy as np

lonmat= [np.array([135.579,135.66 ]),np.array([135.688,135.774])]
latmat = [np.array([-62.2546,-62.6032]),np.array([-62.7195,-63.0676])]


geom_list = [(x,y) for x,y  in zip(lonmat,latmat)]

geom_list_2 = [Polygon(tuple(zip(x,y))) for x,y in geom_list]

enter image description here

,

所以,我最终意识到我可以使用以下代码将我的数组重塑为正确的多边形。

temp = pd.concat([pd.Series(np.ravel(latcorners)),pd.Series(np.ravel(loncorners))],axis = 1)
temp['id'] = pd.Series(np.repeat([x for x in range(int(len(temp)/4))],4))
temp['coord'] = list(zip(temp[1],temp[0]))
temp.columns = ['lat','lon','id','coord']
temp_out = temp.groupby('id').apply(lambda x: Polygon(x['coord'].to_list()))

此代码为数组中的每个变量分配“ID”,这些变量对应于我希望该值成为其中一部分的特定多边形。然后我使用最后一行代码按 ID 分组,并将每组坐标变成一个多边形。

感谢大家的帮助!

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


依赖报错 idea导入项目后依赖报错,解决方案:https://blog.csdn.net/weixin_42420249/article/details/81191861 依赖版本报错:更换其他版本 无法下载依赖可参考:https://blog.csdn.net/weixin_42628809/a
错误1:代码生成器依赖和mybatis依赖冲突 启动项目时报错如下 2021-12-03 13:33:33.927 ERROR 7228 [ main] o.s.b.d.LoggingFailureAnalysisReporter : *************************** APPL
错误1:gradle项目控制台输出为乱码 # 解决方案:https://blog.csdn.net/weixin_43501566/article/details/112482302 # 在gradle-wrapper.properties 添加以下内容 org.gradle.jvmargs=-Df
错误还原:在查询的过程中,传入的workType为0时,该条件不起作用 <select id="xxx"> SELECT di.id, di.name, di.work_type, di.updated... <where> <if test=&qu
报错如下,gcc版本太低 ^ server.c:5346:31: 错误:‘struct redisServer’没有名为‘server_cpulist’的成员 redisSetCpuAffinity(server.server_cpulist); ^ server.c: 在函数‘hasActiveC
解决方案1 1、改项目中.idea/workspace.xml配置文件,增加dynamic.classpath参数 2、搜索PropertiesComponent,添加如下 <property name="dynamic.classpath" value="tru
删除根组件app.vue中的默认代码后报错:Module Error (from ./node_modules/eslint-loader/index.js): 解决方案:关闭ESlint代码检测,在项目根目录创建vue.config.js,在文件中添加 module.exports = { lin
查看spark默认的python版本 [root@master day27]# pyspark /home/software/spark-2.3.4-bin-hadoop2.7/conf/spark-env.sh: line 2: /usr/local/hadoop/bin/hadoop: No s
使用本地python环境可以成功执行 import pandas as pd import matplotlib.pyplot as plt # 设置字体 plt.rcParams['font.sans-serif'] = ['SimHei'] # 能正确显示负号 p
错误1:Request method ‘DELETE‘ not supported 错误还原:controller层有一个接口,访问该接口时报错:Request method ‘DELETE‘ not supported 错误原因:没有接收到前端传入的参数,修改为如下 参考 错误2:cannot r
错误1:启动docker镜像时报错:Error response from daemon: driver failed programming external connectivity on endpoint quirky_allen 解决方法:重启docker -> systemctl r
错误1:private field ‘xxx‘ is never assigned 按Altʾnter快捷键,选择第2项 参考:https://blog.csdn.net/shi_hong_fei_hei/article/details/88814070 错误2:启动时报错,不能找到主启动类 #
报错如下,通过源不能下载,最后警告pip需升级版本 Requirement already satisfied: pip in c:\users\ychen\appdata\local\programs\python\python310\lib\site-packages (22.0.4) Coll
错误1:maven打包报错 错误还原:使用maven打包项目时报错如下 [ERROR] Failed to execute goal org.apache.maven.plugins:maven-resources-plugin:3.2.0:resources (default-resources)
错误1:服务调用时报错 服务消费者模块assess通过openFeign调用服务提供者模块hires 如下为服务提供者模块hires的控制层接口 @RestController @RequestMapping("/hires") public class FeignControl
错误1:运行项目后报如下错误 解决方案 报错2:Failed to execute goal org.apache.maven.plugins:maven-compiler-plugin:3.8.1:compile (default-compile) on project sb 解决方案:在pom.
参考 错误原因 过滤器或拦截器在生效时,redisTemplate还没有注入 解决方案:在注入容器时就生效 @Component //项目运行时就注入Spring容器 public class RedisBean { @Resource private RedisTemplate<String
使用vite构建项目报错 C:\Users\ychen\work>npm init @vitejs/app @vitejs/create-app is deprecated, use npm init vite instead C:\Users\ychen\AppData\Local\npm-