Learner_P-CSDN博客

原创 ArcGIS

ArcGIS10.7/10.8打不开的问题——解决经验&思路一、问题出现1.测试打开其他ArcGIS Desktop2.尝试—重启许可&服务3.删除初始启动文件4. 重新安装5. 最终解决方法6.反思一、问题出现虽然都是软件打不开的问题，但是导致的原因有很多，不建议大家盲目参考，先确认是不是相似操作导致的，解决起来省时省力。之前一直参考麻辣GIS网站安装Arcmap10.7或10.8日常学习，由于工作原因转向Pro学习，这期间做过基于python+QGIS第三方库的平台搭建工作，这个平

2021-07-09 12:58:12 1146 1

原创 Python+OGR库学习（八）：关于面矢量文件的一些小操作

对于OGR库对shp文件的基本读写没有问题，整理一些对geometry简单的小操作代码功能1、面矢量转线：主要理解面由ring构成，提取ring直接写入线矢量文件就OK（或者构建多线几何体，将所有ring都添加进去，完成后可以在arcmap中用拆分多部件要素拆分）2、读取面矢量文件，输出四角多边形：获取图层四角范围，直接写入ring3、输入面矢量，输出它的凸多边形：关键在于创建多面几何体，...

2019-01-05 19:31:57 5507 19

原创 Python+OGR库学习（七）：生成渔网

代码目的输入四角范围，指定单个网格大小，自动生成渔网代码关键点1、整个渔网的生成和数组组成一样，只不过用一网格代替一个数，还是用两个嵌套循环完成2、在1的基础上，从左上角开始，用四角范围中xmin和网格宽度控制列变化，ymax和网格高度控制行变化进行遍历代码思路（）1、导入相关包2、输入渔网范围，格网宽度和高度3、计算行列数4、初始化左上角格网四角范围，创建输出文件5、遍历每一...

2019-01-05 18:36:20 3841

原创 Python+OGR库碎片整理——ExecuteSQL筛选要素

整理一下SQL语句筛选shp要素，关于SQL语句详细使用参考：[OGR：SQL]https://www.gdal.org/ogr_sql.html1、获取cover是grass的要素并id降序返回，然后打印idimport ogr,osos.chdir(r'F:\Python+gdal\7weeks数据\7weeks数据\ospy_data3\ospy3_data')driver = og...

2019-01-01 22:03:16 1175

原创 Python+OGR库学习（六）：批量投影转换（调用自定义函数和模块）

关键思路1、批量转换：定义投影转换函数，参数设置输入输出文件名，输入输出坐标系2、内部功能：针对每一个输入文件，字段属性批量定义（操作对象是layer）、字段值批量写入（操作对象是feature）代码思路1、定义.py文件，含3个函数：字段属性批量定义字段值批量写入投影转换：其中需要调用上面2个函数3、导入以上自定义模块，遍历文件夹，批量投影转换代码——投影转换模块#! usr...

2019-01-01 17:41:49 2075 7

原创 Python+OGR库学习（五）：要素筛选

代码任务多边形是矢量cache_towns.shp,有一个字段name，及点矢量sites.shp，python实现：拿出name为Nibley的要素，并建立1500的缓冲区，用新的要素筛选sites文件，找到在这个缓冲区内的所有点，并打印其id代码思路1、打开两个文件，获取到图层（不管是属性筛选还是空间范围筛选，操作对象那个都是图层）2、layer.SetAttributeFilter(...

2018-12-31 18:32:26 2809 8

原创 Python+OGR库学习（四）：重投影shp文件并另存，属性表保持不变

代码关键点1、首先要定义好转换参数2、主要操作对象是要素，需要提前创建好输出文件，然后遍历所有要素，对每一个几何对象进行坐标转换3、输出文件的字段属性定义需要从输入文件读取代码思路1、导入相关包，切换路径，注册驱动2、定义转换关系3、打开输入文件，读取到图层，创建输出文件，定义字段属性4、遍历要素，循环读取要素、转换投影、创建对应新要素写入几何，并添加要素到图层5、定义空间参考文...

2018-12-31 16:10:18 3155 15

原创 Python+OGR库学习（三）：将含有字段值的TXT文件写入shpfile

代码任务已知有一个TXT文件，每行格式：country_name：X1 Y1,X2 Y2,…Xn,Yn,其中，Xn,Yn指折点坐标，把它导出为shp文件，TXT如下：代码思路1、导入相关库包，切换到当前文件夹2、注册驱动，创建矢量文件，获取图层3、创建输出文件，并获取图层（没有属性定义）4、定义输出图层字段属性：（1）TXT文件中只能得到一个字段，命名为name，设置属性（2）打...

2018-12-29 08:33:27 1918 1

原创 Python+OGR库学习（二）：读取点矢量文件，复制特定属性值点并另存为shp文件

代码思路：1、导入相关库包，切换到当前文件夹2、注册驱动，打开点矢量文件，获取图层3、创建输出文件，并获取图层（没有属性定义）4、定义输出图层字段属性：假设已知文件所有属性字段定义（即ID和cover）（1）读取输入文件中某一要素（2）获取ID、cover字段定义（3）添加定义到输出图层（如果不知道属性表具体字段，可以参考前一篇文章获取：https://blog.csdn.net/...

2018-12-28 08:25:35 4389 2

原创 Python+OGR库学习（一）：读取点矢量文件属性值和坐标，并保存为TXT（一行一个要素值）

代码思路：1、导入相关库包，切换到当前文件夹2、注册驱动，打开点矢量文件，获取图层3、打开待写入TXT文件4、遍历要素：（1）获取当前要素‘ID’和‘cover’字段属性（2）获取当前点要素对应几何对象和其坐标值X,Y（3）将ID、cover、X、Y写入TXT文件（4）清除当前feature缓存并获取下一要素5、清除DataSource并关闭TXT文件关键点：遍历要素，对每一个...

2018-12-26 22:32:12 6145 6

原创 ECMWF ERA数据下载——Python+官网API

利用Python批量下载ERA数据一、注册账号[官网]https://www.ecmwf.int/二、登录后获取下载KEY[获取网址]https://api.ecmwf.int/v1/key/将大括号中内容（包括括号）复制保存为TXT文件，并重命名为.ecmwfapirc.，注意前后都有点，保存在 C:\Users&amp;amp;amp;amp;amp;lt;USERNAME&amp;amp;amp;amp;amp;gt;.ecmwfapirc.，（WINDOW

2018-12-23 19:00:54 10658 6

原创 ECMWF ERA-interim数据下载——手动下载

ECMWF ERA-interim数据下载——手动下载最近研究气溶胶和PM2.5遥感反演，需要用到气象数据，NECP数据分辨率太低，于是尝试ERA-Interim数据，主要用到每天大气边界层高度和相对湿度数据（没有地面相对湿度，可以根据官网公式计算，或用1000pha处相对湿度代替尝试，没有验证精度）进入每日数据下载地址[ERA-Interim,Daily]https://apps.ecmwf...

2018-12-23 18:07:49 14424 40

Winner的博客