当前位置: 开发笔记 > 编程语言 > 正文

python与开源gis电子书_python与开源Gis书本知识点测试

作者：兔斯基小兔子_988 | 来源：互联网 | 2023-10-11 07:51

#-*-coding:utf-8-*-print(upython与开源QGis课题研究组)#print(汉字)##创建矢量数据文件#try:fromosgeoimporto

# -*- coding: utf-8 -*-

print(u"python与开源QGis课题研究组")

#print("汉字")

#&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;

#创建矢量数据文件

#&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;

try:

from osgeo import ogr

except:

import ogr

driver &＃61; ogr.GetDriverByName("ESRI Shapefile")

import os

&＃39;&＃39;&＃39;

ds &＃61; driver.CreateDataSource("xx_tesp.shp")

layer &＃61; ds.CreateLayer(&＃39;test&＃39;,geom_type &＃61; ogr.wkbPoint)

fieldDefn &＃61; ogr.FieldDefn(&＃39;id&＃39;,ogr.OFTString)

fieldDefn.SetWidth(4)

layer.CreateField(fieldDefn)

featureDefn &＃61; layer.GetLayerDefn()

print(featureDefn)

feature &＃61; ogr.Feature(featureDefn)

#设定几何形状

point &＃61; ogr.Geometry(ogr.wkbPoint)

point.SetPoint(0,123,123)

#设定字段数值

feature.SetField(&＃39;id&＃39;,23)

#将要素写入到图层

layer.CreateFeature(feature)

ds.Destroy()

import os

out_shp &＃61; "xx_tesp.shp"

dir(out_shp)

if os.path.exists(out_shp):

driver.DeleteDataSource(out_shp)

dir(out_shp)

point &＃61; ogr.Geometry(ogr.wkbPoint)

print(point)

point.AddPoint(10,20)

print(point)

point.AddPoint(30,20)

print(point)

line &＃61; ogr.Geometry(ogr.wkbLineString)

print(line)

line.AddPoint(10,10)

print(line)

line.AddPoint(20,20)

print(line)

line.SetPoint(2,30,30)

print(line)

print(line.GetPointCount())

print(line.GetX(0))

print(line.GetX(1))

print(line.GetX(3))

ring &＃61; ogr.Geometry(ogr.wkbLinearRing)

ring.AddPoint(0,0)

ring.AddPoint(100,0)

ring.AddPoint(100,100)

ring.AddPoint(0,100)

ring.CloseRings()

print(ring)

print(type(ring))

outring &＃61; ogr.Geometry(ogr.wkbLinearRing)

outring.AddPoint(0,0)

outring.AddPoint(100,0)

outring.AddPoint(100,100)

outring.AddPoint(0,100)

outring.AddPoint(0,0)

inring &＃61; ogr.Geometry(ogr.wkbLinearRing)

inring.AddPoint(25,25)

inring.AddPoint(75,25)

inring.AddPoint(75,75)

inring.AddPoint(25,75)

inring.CloseRings()

polygon &＃61; ogr.Geometry(ogr.wkbPolygon)

polygon.AddGeometry(outring)

polygon.AddGeometry(inring)

print(polygon.GetGeometryCount())

#outring2 &＃61; polygon.GetGeometryRef(0)

#inring2 &＃61; polygon.GetGeometryRef(1)

#print(outring2)

#print(inring2)

for ringx in range(polygon.GetGeometryCount()):

print(polygon.GetGeometryRef(ringx))

multipoint &＃61; ogr.Geometry(ogr.wkbMultiPoint)

point &＃61; ogr.Geometry(ogr.wkbPoint)

point.AddPoint(10,10)

multipoint.AddGeometry(point)

print(multipoint)

point.AddPoint(20,20)

multipoint.AddGeometry(point)

print(multipoint)

mp &＃61; multipoint.GetGeometryCount()

print(mp)

for mpx in range(multipoint.GetGeometryCount()):

print(multipoint.GetGeometryRef(mpx))

extfile &＃61; &＃39;xx_data_pt.shp&＃39;

if os.access(extfile,os.F_OK):

driver.DeleteDataSource(extfile)

#1、创建数据源

newds &＃61; driver.CreateDataSource(extfile)

#print(dir(newds))

#print(dir(newds.CreateLayer))

#2、创建数据源图层

lyrn &＃61; newds.CreateLayer(&＃39;point&＃39;,None,ogr.wkbPoint)

#3、定义图层字段&＃xff0c;添加图层字段

fieldcnstr &＃61; ogr.FieldDefn("idx",ogr.OFTInteger)

lyrn.CreateField(fieldcnstr)

fieldf &＃61; ogr.FieldDefn("namex",ogr.OFTString)

lyrn.CreateField(fieldf)

point_coors_arr &＃61; [[1,0],[2,0],[3,0],[4,0]]

for idxx,point_coors in enumerate(point_coors_arr):

#print(type(idxx))

#print(type(point_coors))

wkt &＃61; &＃39;POINT (%f %f)&＃39; % (point_coors[0],point_coors[1])

#print(wkt)

geom &＃61; ogr.CreateGeometryFromWkt(wkt)

feat &＃61; ogr.Feature(lyrn.GetLayerDefn())

feat.SetField(&＃39;idx&＃39;,idxx)

feat.SetField(&＃39;namex&＃39;,&＃39;ID{0}&＃39;.format(idxx))

feat.SetGeometry(geom)

lyrn.CreateFeature(feat)

#print(lyrn.GetLayerDefn())

#print(&＃39;idx:%i&＃39; % (idxx))

#print("namex:%s" % &＃39;ID{0}&＃39;.format(idxx))

#x &＃61; &＃39;ID{0}&＃39;.format(idxx)

#print(x)

"""

wkt &＃61; &＃39;POINT (%f %f)&＃39; % (point_coors[0],point_coors[1])

geom &＃61; ogr.CreateGeometryFromWkt(wkt)

feat &＃61; ogr.Feature(lyrn.GetLayerDefn())

feat.SetField(&＃39;idx&＃39;,idxx)

feat.SetField(&＃39;namex&＃39;,&＃39;ID{0}&＃39;.format(idxx))

feat.SetGeometry(geom)

lyrn.CreateFeature(feat)

"""

newds.Destroy()

extfile &＃61; &＃39;xx_data_line.shp&＃39;

driver &＃61; ogr.GetDriverByName("ESRI Shapefile")

if os.access(extfile,os.F_OK):

driver.DeleteDataSource(extfile)

newds &＃61; driver.CreateDataSource(extfile)

lyrn &＃61; newds.CreateLayer(&＃39;line&＃39;,None,ogr.wkbLineString)

#字段定义创建

fieldcnstr &＃61; ogr.FieldDefn("id",ogr.OFTInteger)

fieldf &＃61; ogr.FieldDefn("name",ogr.OFTString)

lyrn.CreateField(fieldcnstr)

lyrn.CreateField(fieldf)

point_coors_arr &＃61; [[0,0,1,2,3,-2,6,0]]

#print(point_coors_arr)

#print(type(point_coors_arr))

#print(len(point_coors_arr))

for idx,point_coors in enumerate(point_coors_arr):

#wkt &＃61; &＃39;LINESTRING(%f %f,%f %f,%f %f,%f %f)&＃39; % (point_coors[len(point_coors_arr) - len(point_coors_arr)],point_coors[len(point_coors_arr) - len(point_coors_arr) &＃43; 1],)

wkt &＃61; &＃39;LINESTRING(%f %f,%f %f,%f %f,%f %f)&＃39; % (point_coors[0],point_coors[1],point_coors[2],point_coors[3],point_coors[4],point_coors[5]

,point_coors[6],point_coors[7])

print(wkt)

geom &＃61; ogr.CreateGeometryFromWkt(wkt)

feat &＃61; ogr.Feature(lyrn.GetLayerDefn())

feat.SetField(&＃39;id&＃39;,idx)

feat.SetField(&＃39;name&＃39;,&＃39;line_one&＃39;)

feat.SetGeometry(geom)

lyrn.CreateFeature(feat)

newds.Destroy()

extfile &＃61; &＃39;xx_data_polygon.shp&＃39;

driver &＃61; ogr.GetDriverByName(&＃39;ESRI Shapefile&＃39;)

if os.access(extfile,os.F_OK):

driver.DeleteDataSource(extfile)

#创建数据源文件

newds &＃61; driver.CreateDataSource(extfile)

#创建几何形状图层

lyrn &＃61; newds.CreateLayer(&＃39;polygon&＃39;,None,ogr.wkbPolygon)

#定义字段

fieldcnstr &＃61; ogr.FieldDefn(&＃39;id&＃39;,ogr.OFTInteger)

fieldf &＃61; ogr.FieldDefn(&＃39;name&＃39;,ogr.OFTString)

#图层添加字段

lyrn.CreateField(fieldcnstr)

lyrn.CreateField(fieldf)

wkt_polygon_1 &＃61; &＃39;POLYGON((2 1,12 1,12 4,2 4,1 2))&＃39;

wkt_polygon_2 &＃61; &＃39;POLYGON((4 1,8 1,8 3,4 3,3 1))&＃39;

wkt_polygon_3 &＃61; &＃39;POLYGON((8 4,10 4, 10 5,8 5,6 4))&＃39;

#print(type(wkt_polygon_1))

point_coors_arr &＃61; [wkt_polygon_1,wkt_polygon_2,wkt_polygon_3]

#print(type(point_coors_arr))

for idx,point_coors in enumerate(point_coors_arr):

#print(idx)

#print(point_coors)

wkt &＃61; point_coors

#使用wkt创建几何图形

geom &＃61; ogr.CreateGeometryFromWkt(wkt)

#获取图层要素 ogr.GetLayerDefn , ogr.Feature()

feat &＃61; ogr.Feature(lyrn.GetLayerDefn())

feat.SetField(&＃39;id&＃39;,idx)

print(&＃39;poly_{idx}&＃39;.format(idx &＃61; idx))

#feat.SetField(&＃39;name&＃39;,&＃39;poly_{idx}&＃39;.format(idx &＃61; idx))

feat.SetField(&＃39;name&＃39;,&＃39;poly_{idx}&＃39;.format(idx &＃61; idx))

feat.SetGeometry(geom)

lyrn.CreateFeature(feat)

newds.Destroy()

from osgeo import ogr

import os,math

inshp &＃61; "xx_data_polygon.shp"

ds &＃61; ogr.Open(inshp) #打开shp源文件

driver &＃61; ogr.GetDriverByName("ESRI Shapefile")

outputfile &＃61; "xx_data_polygon_copy.shp" #输出shp

if os.access(outputfile,os.F_OK):

driver.DeleteDataSource(outputfile)

pt_cp &＃61; driver.CopyDataSource(ds,outputfile)

pt_cp.Release()

from osgeo import ogr

import os,math

inshp &＃61; &＃39;xx_data_polygon.shp&＃39;

ds &＃61; ogr.Open(inshp)

driver &＃61; ogr.GetDriverByName("ESRI Shapefile")

outputfile &＃61; &＃39;cp_polygon.shp&＃39;

if os.access(outputfile,os.F_OK):

driver.DeleteDataSource(outputfile)

pt_cp &＃61; driver.CopyDataSource(ds,outputfile)

pt_cp.Release()

outputfile &＃61; &＃39;cp2.shp&＃39;

if os.access(outputfile,os.F_OK):

driver.DeleteDataSource(outputfile)

newds &＃61; driver.CreateDataSource(outputfile)

layer &＃61; ds.GetLayer()

pt_layer &＃61; newds.CopyLayer(layer,&＃39;xxx&＃39;)

#newds.Destroy()

outputfile &＃61; &＃39;cp3.shp&＃39;

if os.access(outputfile,os.F_OK):

driver.DeleteDataSource(outputfile)

newds &＃61; driver.CreateDataSource(outputfile)

layernew &＃61; newds.CreateLayer(&＃39;worldcopy&＃39;,None,ogr.wkbLineString)

layer &＃61; ds.GetLayer()

feature &＃61; layer.GetNextFeature()

if feature is not None:

layernew.CreateFeature(feature)

feature &＃61; layer.GetNextFeature()

newds.Destroy()

driver &＃61; ogr.GetDriverByName("ESRI Shapefile")

inshp &＃61; &＃39;xx_data_polygon.shp&＃39;

ds &＃61; ogr.Open(inshp)

outf &＃61; &＃39;cxp.shp&＃39;

if os.access(outf,os.F_OK):

driver.DeleteDataSource(outf)

print("exists")

newds &＃61; driver.CreateDataSource(outf)

layernew &＃61; newds.CreateLayer(&＃39;worldcopy&＃39;,None,ogr.wkbLineString)

layer &＃61; ds.GetLayer()

#print(layer)

feature &＃61; layer.GetNextFeature()

print(feature)

while feature is not None:

layernew.CreateFeature(feature)

feature &＃61; layer.GetNextFeature()

newds.Destroy()

&＃39;&＃39;&＃39;

from osgeo import ogr

ds &＃61; ogr.Open(&＃39;convChk.shp&＃39;)

layer &＃61; ds.GetLayer(0)

spatialRef &＃61; layer.GetSpatialRef()

print(spatialRef)

推荐阅读

text
Web动态服务器Python基本实现

Web动态服务器Python基本实现 ... [详细]

蜡笔小新 2024-11-21 08:01:30
perl
解决iOS应用推送通知错误：未找到有效aps-environment权限

在尝试加载支持推送通知的iOS应用程序的Ad Hoc构建时，遇到了‘no valid aps-environment entitlement found for application’的错误提示。本文将探讨此错误的原因及多种可能的解决方案。 ... [详细]

蜡笔小新 2024-11-21 19:26:31
input
UVALive 8201 - BBP 公式计算圆周率

在1995年，Simon Plouffe 发现了一种特殊的求和方法来表示某些常数。两年后，Bailey 和 Borwein 在他们的论文中发表了这一发现，这种方法被命名为 Bailey-Borwein-Plouffe (BBP) 公式。该问题要求计算圆周率 π 的第 n 个十六进制数字。 ... [详细]

蜡笔小新 2024-11-21 18:32:57
text
优化 DOM 以提升 JavaScript 性能

本文探讨了如何通过优化 DOM 操作来提升 JavaScript 的性能，包括使用 `createElement` 函数、动画元素、理解重绘事件及处理鼠标滚动事件等关键主题。 ... [详细]

蜡笔小新 2024-11-21 18:16:19
bit
二维码的实现与应用

本文介绍了二维码的基本概念、分类及其优缺点，并详细描述了如何使用Java编程语言结合第三方库（如ZXing和qrcode.jar）来实现二维码的生成与解析。 ... [详细]

蜡笔小新 2024-11-21 17:10:15
main
Java 中的十进制样式 getZeroDigit()方法，示例

Java 中的十进制样式 getZeroDigit()方法，示例 ... [详细]

蜡笔小新 2024-11-21 16:53:03
text
Struts2 + json+ jquery 实现三级联动action和jsp代码竟然有小红叉，提示缺双引号，检查了转义符号也没缺啊，求解

publicclassBindActionextendsActionSupport{privateStringproString;privateStringcitString; ... [详细]

蜡笔小新 2024-11-21 16:25:41
join
Spring AOP学习笔记Advice执行顺序

一、Advice执行顺序二、Advice在同一个Aspect中三、Advice在不同的Aspect中一、Advice执行顺序如果多个Advice和同一个JointPoint连接& ... [详细]

蜡笔小新 2024-11-21 15:28:36
main
Singleton单例模式和DoubleChecked Locking双重检查锁定模式

问题描述现在，不管开发一个多大的系统（至少我现在的部门是这样的），都会带一个日志功能；在实际开发过程中 ... [详细]

蜡笔小新 2024-11-21 15:14:45
text
IC卡操作功能实现

本文介绍了如何通过C#语言调用动态链接库（DLL）中的函数来实现IC卡的基本操作，包括初始化设备、设置密码模式、获取设备状态等，并详细展示了将TextBox中的数据写入IC卡的具体实现方法。 ... [详细]

蜡笔小新 2024-11-21 11:02:19
text
SIP基础概览

本文介绍了SIP（Session Initiation Protocol，会话发起协议）的基本概念、功能、消息格式及其实现机制。SIP是一种在IP网络上用于建立、管理和终止多媒体通信会话的应用层协议。 ... [详细]

蜡笔小新 2024-11-21 17:42:08
list
Java 实现生成指定数量的不重复随机数

本文将详细介绍如何使用Java编程语言生成指定数量的不重复随机数，包括具体的实现方法和代码示例。适合初学者和有一定基础的开发者参考。 ... [详细]

蜡笔小新 2024-11-21 12:08:06
buffer
深入解析 org.apache.jena.atlas.lib.ByteBufferLib 的 acopyArray 方法

本文详细介绍了 Java 中 org.apache.jena.atlas.lib.ByteBufferLib 类下的 acopyArray 方法，并提供了多个实际应用中的代码示例，帮助开发者更好地理解和使用该方法。 ... [详细]

蜡笔小新 2024-11-21 10:26:11
range
解决JavaScript中法语字符排序问题

在开发一个使用JavaScript、HTML和CSS的Web应用时，遇到从SQLite数据库中提取的法语词汇排序不正确的问题，特别是带重音符号的字母未按预期排序。 ... [详细]

蜡笔小新 2024-11-21 09:08:57
text
设置Shadowsocks公共代理的关键步骤

本文详细介绍了如何正确设置Shadowsocks公共代理，包括调整超时设置、检查系统限制、防止滥用及遵守DMCA法规等关键步骤。 ... [详细]

蜡笔小新 2024-11-20 20:41:33

兔斯基小兔子_988

这个家伙很懒，什么也没留下！

Tags | 热门标签

RankList | 热门文章