- 博客(50)
- 收藏
- 关注
RSS订阅原创 如何使用ArcGIS将自己的研究区放在三维地球上展示
以中国为例，右键点击Vertical Perspective (world)，选择复制并修改，对Longtitude_of_center和Latitude_of_center的值修改，这里分别修改为112和15，可以根据自的需要进行设置，修改后，点击应用，即可的得到中国朝向的伪三维地球展示。2. 接着在数据框或数据的坐标系属性里，搜索【54049】坐标系，可以看到World_Vertical_Perspective投影坐标系，然后应用该坐标系，即可得到“三维”视角的地球。
2024-08-28 22:42:16
225
原创 主成分分析学习
8 主成分分析 | 多元统计分析讲义 (pku.edu.cn)清风数学建模学习笔记——主成分分析(PCA)原理详解及案例分析_x10为生均教育经费对以上指标数据做主成分分析，并提取主成分-CSDN博客PCA（主成分分析）-------原理，推导，步骤、实例、代码_pca一般步骤-CSDN博客主成分分析例题_主成分分析例题手算-CSDN博客主成分分析(PCA)算法介绍及matlab实现案例 - XieHS - 博客园 (cnblogs.com)
2024-06-09 20:19:11
340
原创 Arcgis pro报错：000472: 单波段格网的名称不能超过 13 个字符
人家arcgis 10.8都没这个限制，你一个pro版本还搞这个？
2024-06-03 17:10:26
467
原创 MATLAB中自定义栅格数据地理坐标R，利用geotifwrite写入tif
以全球数据为例，今天发现在matlab中对矩阵进行显示后，该数据是正常的，但是添加了空间参考输出后，数据的南北维度颠倒了。如图1是正常显示的结果，图2是输出后的南北纬被颠倒的结果。图下方是一段示例代码。有时候将nc格式的数据转成tiff，或者是将一个矩阵输出成带有地理坐标信息tiff数据时，常常涉及到空间参考的定义和geotiffwrite()函数。在完成了矩阵A的赋值后，应该在后面添加一个rot90()函数将原来的数据倒置一下，所以纠正后的代码应该如下。
2024-05-03 14:44:33
1088
原创 《托斯卡纳艳阳》
有爱情不断流淌的地方，当然充满了诱惑。蛇的设定自然源于亚当夏娃，而蛇的就此消失，当然告诉我们浪漫不代表淫乱，人需要有足够坚定的内心来对抗战胜诱惑。我觉得这就是托斯卡纳艳阳下用细节用镜头用语言所告诉我们的，艳阳下的真实的，有美好的爱情的，人的生活。这爱情可能短暂，可能只是荷尔蒙作祟，可能艰难曲折，可能突如其来，可能柏拉图，但是一直是爱情。此时佛兰西斯说，她是为了生活买了房子，但是，房子里却没有生活，她渴望在这所房子里，能够举行一场婚礼。我们知道，她渴望的是丰盈的内心，充满爱意的人生，但是此时，她什么也没有。
2024-04-26 21:56:52
249
原创 好用的在线Matlab编辑器，风格很炫
Octave Online · Cloud IDE compatible with MATLAB (octave-online.net)
2024-04-24 16:20:43
1179
原创 Matlab 对nc文件进行处理
主要是需要将 nc 数据转为带有地理坐标系的 Geotif 文件，可以很方便的用 matlab 或者 python 完成单个或者批处理的转换，随后，利用感兴趣区的shapefile文件对 tif 文件进行裁剪从而达到提取数据的目的，裁剪我还未在matlab中实现过，批量裁剪有多种方法，我使用过Arcgis pro的模型构建器、python的库有GDAL、Arcpy等。•维(dimension)：一个维对应着函数中的某个自变量，或者说函数图象中的一个坐标轴，在线性代数中就是一个N维向量的一个分量。
2024-04-24 01:52:04
1474
原创 关于审图号的说明
为加强地图审核管理，保证地图质量，根据《中华人民共和国测绘法》等有关法律、法规规定：在中华人民共和国境内公开出版地图、引进地图、展示、登载地图以及在生产加工的产品上附加的地图图形都需要经审核，审核通过之后编发审图号。省级测绘行政主管部门审图号为：省（自治区、直辖市）简称S（××××年）××××号〔如：京S（2004）006号〕国务院测绘行政主管部门审图号为：GS（××××年）××××号〔如：GS（2004）001号〕发文章切记使用带有审图号的底图，不然基本都会被退回修改。（在标准地图官网可查）
2024-04-15 19:34:57
1214
原创 李连杰版《霍元甲》
决定茶品的上下高低，是品茶的人，不同的人当然会有不同的选择，同样是说：“功夫没有高低”，只有习武的人才有强弱之别，武术也不是为了战胜对手，而是发现自己，我们真正的对手是自己，即最可怕的敌人就在自己内心。田中认为比武就像品茶，同样都是茶，就绝对有高低和不同品性之分，意思是世上无数门派众多，自然功夫也有高低强弱，而霍元甲却说：什么是高，什么是低？这些茶都是生在自然中，所以没有高低之分。田中回应：“看来阁下真的不懂”。
2024-04-07 13:10:43
149
原创 混淆矩阵计算遥感分类精度（总体精度、Kappa系数、用户精度、生产者精度等）
假设分类结果的某类别的生产者精度为 80%，那么如果该类别有10个像素的话，则有2个没有分到该类别（丢失了）。假设分类结果的某类别的用户精度为 80%，则用户拿到分类结果后，10个像元中有8个是对的，2个是错的；
2024-04-06 14:07:54
2237
原创 Geopandas的介绍（1）
GeoDataFrame：一个包含GeoSeries的列表数据结构。重要属性：geometry。可以通过geometry属性(gdf.geometry)访问，并且可以通过输入gdf.geometry.name找到geometry列的名称。GeoSeries：本质上是一个矢量，矢量中的每个条目都是一组对应于一个观测结果的形状。一个条目可以只包含一个形状(比如一个单一的多边形)，也可以包含多个形状。Geopandas有三种基本的几何对象：点、线、面。
2024-03-29 20:47:35
254
原创 GDAL的数据类型（9）
GDT_Int32 = 5(C++中对应int 或 long 或 long int)GDT_Int16 = 3(C++中对应 short 或 short int)GDT_UInt16 = 2(C++中对应 unsigned short)GDT_UInt32 = 4(C++中对应unsigned long)GDT_Float64 = 7 (C++中对应double)GDT_Byte = 1(C++中对应unsigned char)GDT_Float32 = 6(C++中对应float)
2024-03-26 22:10:20
413
原创 GDAL 镶嵌至新图像（7）
其中的SetProjection()会为影像添加上 地图的投影、地理坐标系和投影坐标系。out_ds.SetGeoTransform(gt)方法重新计算地理转换。
2024-03-26 21:08:50
190
原创 联想拯救者Y9000P 2023款 内存占比巨高解决办法
Win+R打开运行输入“MdSched”即可使用内存诊断工具，用了直接降很多！电脑刚开机，内存就占用巨高，啥都还没开呢，先给我个下马威。
2024-03-21 14:12:02
832
原创 学习工作日记
批量读取nc并转成tif，计算青藏高原的平均蒸散发。2、LAI的反演（利用植被指数经验法）1、下午需要完成两笔单子。3、还有个万恶的英语作业。
2024-03-19 13:40:49
99
原创 使用Arcpy对gee上得到的数据进行拼接（6）
原因：数据太大，在gee中自动把一张tif给拆成了两张输出到云盘中。拼接闫凯教授的LAI数据时的代码（Arcpy）
2024-03-17 18:54:44
258
原创 GEE 解析 ee.Reducer（四）
2、ee.Reducer.toList() 用于将像素值列表归并到一个列表中的 reducer，会主动排除Nodata。1、ee.Reducer.count() 获取非Nodata值的个数，以计算指定区域内的非Nodata像素的数量。
2024-03-16 20:23:21
525
原创 GEE 导出至Google drive（三）
evaluate：主要用来异步导出和显示，加快速度，一般没有返回值，结合 for 和 map 一起用，这里的 map 里可以有 print,Map,toDrive 方法。之前只知道怎么Copy代码，但是不会去深究每行代码的含义，今天是周六，正好闲来无事打算瞅瞅这么个事。
2024-03-16 20:23:15
595
原创 GDAL(python) 之GeoTransform（3）
使用GDAL库读出的dataset带有两个重要的地理参数，分别是Projection（暂未搞清这个是干啥的）和GeoTransform（用作仿射变换的）。有了这两个参数，就确定了影像的地理位置。注意：本代码仅测试输出的top和botton可以上下接起来，并不讨论top.tif和bottopm.tif的真实投影位置。左上角x坐标， 水平分辨率，旋转参数， 左上角y坐标，旋转参数，竖直分辨率（负值，原点坐标在左上角）。再GDAL for Python中，GeoTransform是一个六个元素的元组。
2024-03-08 10:49:22
510
原创 GEE批量下载ImageCollection且无需呆点run！（一）
其中，collection是你需要导出的影像集合名称, folder是你导出到drive的文件夹名称, options是你需要设置的参数（scale、region、单张影像name、maxPixels以及数据类型等）。批量执行GEE导出任务：首先，您需要生成“”导出任务”（全部在task中显示run才行）。然后按F12进入控制台，然后将这些脚本粘贴到其中，然后按enter。需要下载某位大佬生产的LAI/FPAR数据，且数据都在GEE上，数据量达到约4500张tif，要是单点run就太呆了。
2024-03-07 13:41:47
585
原创 教你如何安装GDAL（1）
打开python的IDLE查看python的版本，需要安装与python版本一致的GDAL 下载与刚才确认的两个参数一致的GDAL文件包。下载链接如下：Archived: Python Extension Packages for Windows - Christoph Gohlke (uci.edu) 这里看两个参数，一个是cp屁股后面的数字。比如，35就是python版本3.5的意思，至于是3.5.1还是3.5.9都不重要，下载3.5就行了。第二个是看后面amd64还是win32，这个容易理解（和p
2024-03-04 14:30:31
2326
6
原创 IGBP classification (with code)
Derived from IGBP classificationsBSV Barren Sparse Vegetation: Lands exposed soil, sand, or rocks and has less than 10% vegetative cover during any time of the year. CRO Croplands: Lands covered with temporary crops followed by harvest and a
2023-12-15 15:02:17
988
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人