- 博客(29)
- 收藏
- 关注
RSS订阅原创 GDAL Python 过滤Shape Polygon中的面积小于某个阈值的小图斑
【代码】GDAL Python 过滤Shape Polygon中的面积小于某个阈值的小图斑。
2023-09-06 11:10:31
240
原创 GDAL循环读取矢量文件中的geometry
GetArea和Area在geom尺度上是同一个函数。按照feature数量进行循环。3就代表polygon。
2023-09-06 10:12:32
224
原创 关于gdal numpy计算NDVI等指数的问题
在利用python gdal numpy计算NDVI等运算栅格数组的指数时出现nan和inf的问题。会出现nan和inf，在QGIS中加载出现问题。去除nan和inf，并进行二值化。以上得到np多维数组。
2023-09-05 10:24:33
178
原创 pycharm添加conda虚拟环境
使用prefix在指定文件夹创建Anaconda虚拟环境。添加本地解释器，找到虚拟环境目录后添加解释器。使用conda安装所需要的包。使用pycharm创建项目。
2022-11-30 17:42:16
1288
1
原创 12_OSR模块和空间参考系
定义在WGS84基准面下的北半球UTM17带。设置一个定义用户名字的投影坐标系统并确定系统被投影过，应该是参数是用户随意输入，只表示用于定义坐标系的名称；设定一个坐标系统，EPSG中已有的各种参数；设置投影转换参数细节，包含两个参数，zone表示投影带号，north表示是否为北半球，默认是True。
2022-11-30 15:41:17
1172
原创 10_空间过滤器Spatial filters
使用来实现根据空间位置的筛选。先设置两个shp文件，其中的要素需要相交。首先需要获取重叠两个shp文件的图层。获得叠加图层的空间坐标使用函数对原有图层进行筛选，此时的是已经被筛选过后的。直接CopyLayer这个图层形成新的shp文件。这个函数可以设定周围四个边界点形成矩形框。效果和shp筛选一致。注意：经过后的layer并不是完整的。若需要对layer执行其他操作，那么需要将其进行筛选清空。
2022-11-30 15:38:11
351
原创 9_OGR创建Shapefile
使用ogr创建一个shp文件，并且生成一个点坐标，添加两个字段并填入内容。整个逻辑过程如下：首先，生成SHP的Driver；判断当前路径是否已经有同名的SHP文件存在，若存在则通过os模块删除；通过Driver创建数据集，即创建矢量文件；为矢量文件新建图层，同时在参数中设定图层名称、空间参考系、几何类型；设定字段，并在矢量文件的图层中添加字段；为图层创建feature；创建几何图形，并设定几何图形的坐标，将其添加到feature中；为feature设定字段属性；将feature写入。
2022-11-30 15:36:57
1024
原创 7_使用OGR模块操作矢量数据
其实参数都是一样的，第一个参数表示传入的路径，第二个表示读取方式。在实际的编程过程中，最好是先声明打开文件类型的。文件而言，这个是缺省的打开方式。函数之间有点小小的不同。一样的操作，但是对于。
2022-11-09 14:56:00
239
原创 6_GDAL和Pillow互操作
GDAL和PIL很相像，它们处理和操作的对象都是栅格图像。区别在于GDAL主要重点放在地理或遥感数据的读写和数据建模以及地理定位和转换；PIL的重点是放在图像本身处理上的。在底层逻辑上，两个库都是对numpy的数据进行处理，因此理论上可以相互共享和交换数据。GDAL的核心在波段band， 一切操作的基础和核心都在波段。波段可以单独拿出来操作，至于波段在数据集中的顺序无关紧要。遥感图像大多比RGB图像的波段要多，而每个波段单独都是一个完整的整体， 每个波段单独拿出来都是一个数据集。Pillow。
2022-11-06 20:54:30
253
原创 5_GDAL创建影像
该函数用于创建影像。指保存的路径和文件名称；指生成影像的长和宽；int bands指生成影像的波段数目，默认是1；指影像的DataType，一般与原始影像的DT相同。生成影像的基本流程是：先声明一个指定图像格式的Driver；通过函数生成初始化的影像，同时指明路径、大小、波段和DT；为Create后的影像设置Projection、；写入具有值的矩阵到新影像中。gdal。
2022-11-06 20:53:28
4028
原创 4_获取影像像元数据
本文中的方法可以对第三章进行补充，能够访问像元灰度值，并对其进行处理。GDAL提供了下面两个函数来访问影像的数值。读取图像数据(以二进制的形式)读取图像数据(以数组的形式)xoff, yoff指定想要读取的部分原点位置在整张图像中距离全图原点的位置；指定要读取部分图像的矩形的长和宽；~~放缩比例，最后读取影像的长宽会根据这个比例因子进行缩放；~~测试发现这两个参数并不是缩放比例因子，并不是缩放比例因子，而是新的指定区域，比如指定为(20, 20)，那么。
2022-11-05 21:05:43
856
原创 2_开始使用GDAL
GDAL和OGR都归于osgeo除了GDAL包，还有一个gdalconst包最好也导入。gdalconst也是osgeo的一个包，它只是在代码中对GDAL中用到的一些常量进行了绑定。其中gdalconst中的常量都加了前缀，力图与其他模块冲突最小。
2022-11-05 21:03:47
636
原创 1_GDAL简介
包含GDAL的许多基本概念，感觉是机翻，但都是GIS和RS的部分基础衍生。GDAL——处理栅格空间数据的类库OGR——处理矢量数据的类库。
2022-11-05 21:02:11
429
原创 cesium中调整图层顺序
cesium中图层由_layers存储viewer.imageryLayers._layers其中默认图层存在时数组中只有一个元素，默认为下图2的ion添加其他图层可以追加在该数组之后当前有个项目要求sample始终在basemap之上，这里需要一个简单的交换位置。// 全局变量 用于存储底图源 和 底图层let basemapSource = '';let basemap = '';function imageryBasemapLayers(layers) { if(base
2022-04-29 15:02:43
1976
原创 cesium加载图层
// 全局变量 用于存储底图源 和 底图层let basemapSource = '';let basemap = '';function imageryBasemapLayers(layers) { if(basemapSource == ''){ basemapSource = new Cesium.WebMapServiceImageryProvider({ // 这里是你的 geoserver服务点击查看图层的 url u
2022-04-29 10:37:11
658
原创 jQuery数组函数和工具函数
jQuery数组函数$.each(arr, function(index, item)) parameter：Array 回调函数function(index, item) 可以遍历纯数组对象（array），类数组对象（arguments对象，DOM对象）;$.grep(arr, function(item, index)) 可以过滤数组中的元素，只保留符合特定条件的元素并返回一个新数组。 let result = $.grep(list, function (ite
2020-11-08 20:36:59
241
1
原创 Geoserver+Postgis+Mapbox实现路径规划
Geoserver+Postgis+Mapbox实现路径规划 2020/10/14，利用Postgis数据库对LineString的Shp文件进行路径规划，利用Geoserver对规划的结果进行发布，在Mapbox中使用WMS服务调用并呈现在地图上。 本文是以经纬度坐标为依据进行计算，参考其他网上已有教程可以选择以节点进行计算，其他链接放置本文底部。软件版本Geoserver：geoserver-2.17.2-bin；PostGIS：2.1.7PostgreSQL：9.3QGIS：3
2020-10-14 21:03:57
3490
2
原创 Geoserver调用Arcgis切片
Geoserver调用Arcgis切片 2020\09\18 对Arcgis Server发布江苏海洋大学矢量化地图切片进行Geoserver（简称GS）本地发布。本质是调用第三方切片，在GS中发布，从而绕开Arcgis Server。软件版本：java 1.8.0_261geoserver 2.17.2geowebcache 1.10.1（原本采用1.15.2，但会在进入http://localhost:8080/geowebcache/demo报错拒绝访问，更换低版本可以进入）操
2020-10-14 19:39:06
1079
1
原创 c#属性的定义
1、什么是属性属性是一种用于访问对象或类的特性的成员。属性有访问器，这些访问器指定在他们的值被读取或者写入时所需要执行的语句。我们知道的访问器有get访问器和set访问器两种2、属性的声明public class Person{ private string name;//定义个私有的字段，外部不能对私有变量name直接进行读写，但可以通过其属性进行。 定义字段用小写 ...
2019-07-25 21:45:54
5272
原创 输入一段中英文混合的字符串 自动统计其中的英文字母出现频率 忽略大小写
import java.util.Scanner;public class Main { public static void main(String args[]) { String str1,str2; //定义两个String的变量 str1和str2 int ary[]=new int[26]; //定义一个新的数组ary[] 其中有26个字母 Scann...
2019-03-19 14:30:53
4227
原创 简单的java华氏转摄氏温度
import java.util.Scanner;public class Transform { public static void main(String args[]) { int x; Scanner i=new Scanner(System.in); x=i.nextInt(); System.out.println((x-32)*5/9); }}
2019-03-18 15:29:32
873
原创 Java计算100到1000以内两个数之间的prime之和
在中定义两个数n和m 而且符合100<n<m<1000import java.util.Scanner;public class Prime { public static void main(String args[]) { int n,m; int sum=0; Scanner i=new Scanner(System.in); System.out.p...
2019-03-18 15:28:36
452
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人