按shp统计范围内散点
摘要:import shapely.geometry as geometry from shapely.geometry import Polygon import numpy as np # 数据格式转换 from osgeo import gdal_array, ogr from numpy impo
阅读全文
shp文件按属性删除
摘要:from osgeo import ogr, osr, gdal ,gdalconst import numpy as np import math import sympy as sy import cv2 #注册所有驱动 gdal.AllRegister() #解决中文路径乱码问题 gdal.S
阅读全文
ArcGIS Python 修改面的左上角为第一个点
摘要:面的节点夹角大于30并且小于150,距离外界矩形左上角最近的点为第一个点,对于多部件每一个都修改,修改面的开始点位置,运行界面如图14-10所示。 #coding=utf8 import arcpy import os import sys import math isEdit=False def
阅读全文
使用GDAL/OGR打开矢量并输出每个面外界矩形范围内的point数据
摘要:原文链接:https://blog.csdn.net/weixin_40625478/article/details/106851352 本文主要目的:我们有的时候需要获取矢量数据的外接矩形范围,但是一个图层数据有好几个面要素,如果获取整体的外接矩形进行处理的话,还是会造成数据量较大;那么我们就逐个
阅读全文
基于Python和GDAL提取栅格数据相邻地物的边界
摘要:摘录于 https://blog.csdn.net/weixin_43123242/article/details/93525175 1.下载第三方包 在网址 https://www.lfd.uci.edu/~gohlke/pythonlibs/#lxml下载对应python版本的whl文件。如,G
阅读全文
shp 汉字乱码解决方案
摘要:在ArcGIS 10.2以上版本,汉字经常乱码,解决方法:做一个扩展名为.reg,内容如下: Windows Registry Editor Version 5.00 [HKEY_CURRENT_USER\Software\ESRI\Desktop10.7\Common\CodePage] "dbf
阅读全文
ARGIS栅格计算器常用命令
摘要:1、将特定值(-9999)赋值为0 Con("raster" 9999,0,"raster") 2、将某一范围内(如小于0)赋值为0 Con("raster" <0,0,"raster") 3、将特定的栅格值(如1)赋值为空值 setnull("raster"==1,"raster") 4、给空值赋
阅读全文
线的简化和平滑笔记
摘要:一、简介 在对地理信息数据处理时,常常会遇到对线进行平滑和简化的操作。线的平滑可以使用拟合或插值来完成。多段线(Polyline)简化算法可以帮助我们减少Polyline的点数,从而降低输入规模。对多段线简化算法,通常的做法是在一定的近似精度下,删除一些点或者边。 二、线的平滑 插值和拟合有所不同。
阅读全文
shp平滑处理
摘要:在做图像数据处理时,经常会有栅格数据转矢量数据的操作,转换后的矢量文件会存在锯齿状边缘,不太美观,因此常常需要对矢量(shp)文件做平滑处理。 1 利用arcgis实现shp的平滑和简化 ArcToolbox / Cartography Tool / Generalization / Smooth
阅读全文
矢量裁剪矢量
摘要:近日在博客中看到一篇使用ogr裁剪矢量数据的文章,觉得挺好,就做个笔记来学习。文章链接https://www.cnblogs.com/ljwgis/p/14071604.html 核心代码 from osgeo import ogr, osr import os def createShpByLay
阅读全文
shp的基本操作
摘要:本节将介绍如何利用python完成对shp的基本操作 1.读取shp四至 import shapefile sf = shapefile.Reader(r"E:\shp\1.shp") #读取shp四至 min_x, min_y, max_x, max_y = sf.bbox #读取每个图斑四至 s
阅读全文