你好，我是小张同学！

原创 数据分享：我国ESA_WorldCover_10m土地利用数据分享

该项目的关键成果是于2021年10月发布了2020年可免费获取的10米分辨率全球土地覆盖产品，该产品基于Sentinel-1和Sentinel-2数据，包含11个土地覆盖类别，并独立验证，全球总体准确率为74.4%。该数据源于Sentinel-1/2的10m分辨率图像制作而成，本次分享我国区域土地利用数据。链接：https://pan.baidu.com/s/1ip69EXzJWyop8-8WLU4gAA。共分为11类: 林地、灌木、草地、耕地、建筑、荒漠、雪\冰、水体、湿地、红树林、苔藓\地衣.

原创 2020地理设计组三等奖：基于深度学习的气象干旱时空预测分析

作品简介一、背景与简介从全球范围看，旱灾影响面积最广，造成经济损失最大，被认为是世界上最严重的自然灾害类型之一[1]，而中国干旱地区面积占全国自然灾害影响面积的60%[2]。由于其出现次数多、持续时间长、影响范围大、对农业等经济部门造成的直接损失重，再加上水资源、土地资源及其对社会的潜在影响，成为我国最大的自然灾害之一[3]。因此，量化研究干旱时空变化特征并阐述其形成的机制，对农业生态系统的科学管理和气象灾害预警等方面有重要意义[4]。国际上将干旱分为气象干旱、农业干旱、水文干旱和社会经济干旱四类[5-6]

原创 2020地理设计组三等奖：基于GIS的广州市天河区“两抢一盗”犯罪案件的多尺度时空格局研究

（1）以法院裁判文书网审判文书数据为基础，结合新兴的空间分析技术，探究展开犯罪地理研究的可行性和潜力，通过提取记录、数据清洗、位置定位、整合建库的数据获取过程，并基于Excel对犯罪数据进行时间维度上的分析，基于GIS空间分析方法对犯罪数据进行空间维度上的分析，尝试建立从时空特征分析到挖掘改良区域的整套流程，不仅可用于对犯罪案件发热研究，也可用于其他地理因素或现象的探索分析。犯罪作为人类社会发展的一种社会现象，与人类的生活与活动息息相关，因此打击犯罪，保障人民的财产安全和社区的稳定发展任务仍然紧迫。

原创 ArcMap 制图教程 ： 中国人口分布图制图完整步骤

感觉好像差不多了……因为没有固定的模式，看个人的感觉，所以你可以根据感觉再调整一下，不要受限与我给出的数值）然后在Symbol表头上，点鼠标右键，选择Filp Colors（反色），因为默认是降序的，我们改为升序：？最后的最后，有同学问，说我这个图做出来之后，论文中要求要提交300dpi的图，而截图分辨率不高啊，怎么办。切换回数据视图，然后把比例尺设置到1:800,000,（当然，你也自行调整到合适的大小）其他全部下一步到结束，（大小地图的符号化是一样的色带，所以用同一套图例即可）

原创 ArcGIS小白入门——ArcMap应用基础讲解

在本实验中，通过实际操作，熟悉和掌握了ArcMap中的数据加载、数据显示与查看、简单的符号化以及地图的排版与整饰等应用基础功能。然而，GIS中采集、编辑、存储、分析、表达与制图等信息流内容繁多、ArcMap软件功能丰富，在以后学习中会要更加注重学习方法适当以及学习时间保证。

原创 了解 ArcGIS 编辑器，从标准、高级到拓扑，从此不求人

三种编辑器：标准编辑器工具条的常用工具的使用演示和说明；高级编辑器工具条所有工具的详细演示和使用说明；拓扑编辑工具条中所有涉及拓扑边工具的详细演示和使用说明。编辑器下拉菜单中的通用工具；如何绘制各种垂直和平行的几何、以及流的使用、按指定距离绘制点线以及避免要素被拖动；配置各大类工具条的快捷命令。

原创 数据分享｜MuSyQ 30米/10天叶片叶绿素含量产品 （2019-2020年中国01版））

叶片叶绿素含量(leaf chlorophyll content, Chlleaf)直接反映了植被叶片的生理状态，是植被获取光合作用所需太阳辐射信息的重要参数，遥感是获取大范围动态Chlleaf的主要技术手段。本文利用哨兵二号多光谱成像仪(Sentinel-2 MSI)高时空分辨率的特点，基于叶绿素敏感指数（chlorophyll sensitive index, CSI）生产了MuSyQ高分系列中国地区2019 - 2020年30 m/10天分辨率的标准化Chlleaf产品01版。

原创 2020地理设计组三等奖作品：基于新型冠状病毒肺炎(COVID-19)疫情数据的可变面元问题(MAUP)效应分析

进一步测度六个因子对结果的共同作用，发现，随着空间尺度的增大，六因子对结果的拟合效果降低。分异效应探究主要分为省级和市级两种分层方案下，其中，省级分层方案主要是以百度迁徙指数、公路密度、GDP、邻接关系、人口结构、人口密度、铁路密度为自变量，市级分层方案以GDP、百度迁徙指数、城市等级、邻接关系、人口密度、纬度为自变量，均以累计病例数据为因变量，使用地理探测器的方法，研究不同影响因子对疫情的影响，以及在因子两两交互时，对结果的作用，从而找出影响疫情蔓延的存在因子。因此，疫情数据的分异具有时间粒度效应。

原创 2020遥感应用组二等奖：基于LiDAR和资源三号数据的广西高峰林场森林资源三维动态监测

同时，资源三号作为高分辨率的多光谱影像，携带的多光谱信息可以用于辅助反映森林植被信息，且由于其能重复对地进行观测，为林地各类资源的动态变化监测提供了更有力的手段。分别对比2016年和2018年的正视后视和正视前视两种组合所提取的DSM数据，选取误差较小的组合为之后研究所用，借助LiDAR数据提取的DSM检验精度之后，利用ENVI中的Band Math工具将2018年的DSM和2016年的DSM相减，得到2016-2018年的林场树高变化分布图，之后对树高的变化根据树种进行动态变化分析和成因剖析。

原创 2020遥感应用组二等奖：基于多时相主被动遥感的漓江水体监测与叶绿素反演研究

NDWI水体指数对植被与水体有较好的区分能力，MNDWI水体指数适用于增强和提取背景为建成区的水域信息，EWI水体指数则结合了以上两种指数的优点，较好地对植被、土壤、建筑区进行反射抑制，增强了水体亮度，NDPI水体指数不仅可以提取小型池塘和低于0.01公顷的水体，还可以区分水体与周围环境中的植被信息。式中，a、b、c、d为权重系数，BNDWI为研究区NDWI指数数据，BMNDWI为研究区MNDWI指数数据，BEWI为研究区的 EWI指数数据，BS为研究区后向散射系数数据。式中，S为水面周长，A为水体面积；

原创 2020开发竞赛地理设计组二等奖：5G基站选址优化与应急通信调度

总体思路：结合地理建模以及Arcpy二次开发脚本工具，提取DEM候选局部制高点并计算其可视域，根据可视域和预算组合的多目标优化模型进行优化选址，其中最低成本全覆盖模型及最大覆盖有限成本模型均可在考虑区域原有基站的基础上进行优化。根据视域矩阵得到候选山顶点，再删去根据禁止/优先覆盖栅格处理得到的无需覆盖栅格和需要覆盖但有候选点可覆盖栅格得到的候选点，得到实际基站设置点，实际基站设置点乘以平均造价得到成本 object，实现有限成本最大覆盖。根据缓冲区将距离该地区原有基站过近的局部制高点删除，得到候选点数据。

原创 2020竞赛特等奖作品赏析：科技战疫·突发急性传染病多尺度智能监测预警系统

另外，还可以进行CPU的多级分配。疫情监测通过宏观和微观结合进行展示，宏观上，国家、省和直辖市的疫情信息，微观上，市和区的疫情信息。疫情防控之中，最重要的就是对人员轨迹的掌握，传染病传播大多通过人传人，传染源的活动轨迹上遇到的人都有可能被传染，所以病患轨迹是预防、监测的重点，项目采用Hadoop大数据平台发布病患轨迹数据，面对亿级数据量的轨迹信息，以云计算实现海量轨迹的筛查追踪、预警潜在患者、智能分析重点人群轨迹等作用，便于对政府、群众提供预警服务，为城市疫情防控筑起了安全屏障。

原创 中国各流域分布图

原创 2007~2016 年税调经纬度及其所处的省市区县乡镇数据

之前给大家分享过一份税调企业经纬度及其所处的省市区县数据：2007~2016 年税调企业地理信息数据（含经纬度及其所处的省市区县）：https://rstata.duanshu.com/#/course/76d38022cd004b09b2aa09647936beb0最近有培训班的小伙伴提出是否能根据税调企业经纬度来判断其所属的乡镇信息，于是我就在之前的结果上进行了补充。

原创 基于时空的Ramsar湿地自动淹没映射利用Google Earth Engine

为了保护这8个湿地，可以使用一些机器学习模型来预测湿地的未来淹没情况，并通过未来分析找出湿地淹没面积减少的可能原因，以采取适当的措施保护受威胁的湿地。所有从六月到九月的可用场景都被用来构建合成图像，这个时期对许多区域来说是植被增长和动态变化的季节，监测这个时期的湿地状况可以提供关于植被健康、水源可用性和湿地生态系统整体功能的宝贵洞见。研究的创新性在于提供了一种新的方法和工具，提高了映射结果的准确性和可靠性，并揭示了淹没区域的时空变化对湿地生态系统的影响，为湿地保护和管理提供了科学依据。

原创 ArcGIS点集之间两两连线

比如有4个点,则连3+2+1=6条线;有38个点,则连37+36+35+···+3+2+1=703条线。我看了有一个大佬用python写了代码(不知道为啥我运行一直报错)都快要放弃了,然后翻书在书上找到了。百度发现大佬们的骚操作,

原创 ArcGIS提取路网节点

将路网通过空间链接到点上。生成一个新的点图层(注意两个图层的顺序,目标要素是{点图层},连接要素是{步骤二中打断线的图层})。这时候就成了带有节点的网络数据集,但这还不是我们需要的点文件,因为这里面有悬挂节点和自相交节点。这是我学校周围的路网,数据就不分享给大家了，看教程知道有这个方法就好了。),删除,剩下的点就是我们要的线的交点,保存文件.先打开编辑器,把他合并成只有一个字段的要素。的文件,打开此文件属性表,发现多了一个字段。这个字段的数值就是点和多少条线相连的数;的点(根据实际情况可以选择小于。

原创 ArcGIS标注的各种用法和示例

上面的例子中出现了【"" + [中类] + ""】这样的语句，这是ArcGIS在标注中的一种格式化标签，用来定义字符的格式，包括文字字体大小，颜色等等。在【方法】栏中选择【定义要素类并且为每个类加不同的标注】，点击【SQL查询】，在弹出的框中输入SQL筛选语句，这里以【面积>10000平方米】为例，只标注面积大于10000平方米的要素，面积小于10000平方米的就不标注。这里表达式中的几个参数可以自己慢慢调整。

原创 arcgis图上添加发光效果！

然后手动画一个大范围的面数据，利用“擦除”工具，将要表达的区域范围擦除掉，调整擦除后留下的区域的透明度，加上标注或注记，得到基础的底图——我们想让边界更加有立体感，这样出图效果更加有质感，可以利用制图表达的“偏移切线”技术，将区域的边界线进行偏移，形成立体阴影效果——二、关键的一步来啦，为每个点创建多环缓冲区，根据小飞飞的经验，最好是创建九级，甚至更多级，因为这样光晕效果会更突出——三、然后，在得到的多环缓冲区属性表内，新建一个“透明度”字段，主要的作用是控制每一级缓冲区的透明度，达到发光的效果——

原创 多维时空数据介绍（4）多维栅格数据构建

上图中包含了经度、纬度和时间三个维度，其变量为经度、纬度、时间以及降雨量，对应的变量单位为东经多少度、北纬多少度、阳历1895-1-1日开始多少天以及每年降雨多少毫升，其中经度数值4个、纬度数值3个、时间数值2个都是一维的数据，而降雨量是一个包含经度、纬度、时间的三维数组，数值有4*3*2=24个。在GeoScene 中NetCDF数据可以用做多种多维时空数据的数据源，进行创建、可视化和分析，它可以是多维栅格数据、时空立方体，还可以是体元图层，区别在于其所要表示的内容和目的。变量又分为坐标变量和属性变量。

原创 多维时空数据介绍（3）局部异常值分析

NUM_OUT、PERC_OUT、N_LOW_CLS、P_LOW_CLS、N_LOW_OUT、P_LOW_OUT、N_HIGH_CLS、P_HIGH_CLS、N_HIGH_OUT、P_HIGH_OUT、NO_SP_NBR字段则分别表示高-高聚类、低-低聚类、高-低异常值、低-高异常值在同一个时间序列中的数量和对应的百分比。CO_TYPE表示结果类别类型，可区分单个条柱具有统计显著性的高值（高-高）聚类、低值（低-低）聚类、高值主要由低值围绕的异常值（高-低）以及低值主要由高值围绕的异常值（低-高）。

原创 多维时空数据介绍（2）新兴时空热点分析

选择热点和冷点趋势主题后，则生成一个新的二维面，其中的属性表包含的字段比新兴时空热点分析结果要少，重要字段为TREND_ZSCORE、TREND_PVALUE、TREND_BIN，仍然表示基于Z得分的趋势分析得到的内容，这与新兴时空热点分析结果中对应的字段内容是完全相同的。（2）在热点分析基础上，按照条柱时间序列执行聚类（值为Z得分）的趋势分析，也使用Mann-Kendall 统计方法，趋势分析的结果为聚类趋势 z 得分、p 值和每个位置的分组类别。在3D视图中，可以查看条柱的冷热点信息，但是不显示趋势。

原创 多维时空数据介绍（1）时空立方体创建及可视化

其生成的立方体的时间间隔与多维栅格图层的时间间隔相同。2D 模式下，用于显示时空立方体工具中的输入参数必须是上述三种创建时空立方体工具中生成的NetCDF文件，或者是上述NetCDF文件做过时空模式挖掘分析后的结果（做完分析的结果以主题的方式保存在原NetCDF文件或者以主题的方式存储在新的NetCDF文件）。通过已定义位置构建的时空立方体，其条柱的形状跟输入要素的形状一致，如果输入是气象点，则生成的立方体为渔网形状的条柱，但与之前聚合点创建的时空立方体紧密相连不同，只在气象点的位置上有立方体，。

原创 明明都是2000坐标的地形图，怎么位置就不一样呢？

前几天有一个同事在工作中遇到一个问题，核心说起来就是，甲方发来两个年份的同一个片区的地形图，并且言之凿凿都是2000坐标的，但是在一个CAD复制后，到另一个CAD中使用“粘贴到原坐标”，就是位置对不上。由于是随手点的“左上角”坐标，所以有误差，这不重要，有没有发现特点：高程点的整体坐标其实是等高线那个CAD的横轴上增加了36000000。在只有CAD的时候不知道带号具体是啥，可以按照我的方法，分别看一下左上角坐标是啥，比较一下，大概就能知道了。所以，其实说是两个2000坐标的CAD，其实。

原创 GIS教程—坝区面积计算

打开编辑器和高级编辑器，选中需要拆解的区域，在高级编辑里面点击拆解，保存编辑，关闭编辑，然后在此筛选出所需要的地区。做坡度分析，打开Slope_create2的图层属性打开，在符号系统里用已分类显示，分成五类，选择合适的色带。在输入要素类时，选择下拉三角选择要素，等高点和等高线的高度字段均要选择BSGC，并得出结果。在符号系统里面将边类型和高程去掉添加具有相同符号的边和具有相同符号的点，放大图层显示。设置好路径原始数据微雨data数据库内，新生成的数据放于shuju这个数据库内。

原创 ArcGIS根据相邻关系提取相邻面&提取面公共线

因此,如果某面的外边界左侧(外部)和孔洞边界左侧(内部)均没有相邻点,则生成线的 LEFT_FID 值为 -1,而 RIGHT_FID 则设为面要素 ID。如果一个面与另一个面重叠,那么将生成两条输出线以便分别表示每个相交边界:第一条线表示其中一个重叠面的外边界,因此该线的 LEFT_FID 为与其相交的面的要素 ID,而。如果面内包含另一个面,则会生成一条顺时针方向的输出线来表示公共边界,该线的 LEFT_FID 设置为外部面要素 ID,而 RIGHT_FID 设置为内部面要素 ID。

原创 山脊图、密度图，最全总结实现方法在这里！

Joy Plot 得名于 1979 年 Joy Division 的专辑封面 Unknown Pleasue，专辑的封面就是现在山脊图的形状。但是直方图有个缺点，设置的分割距离不同，直方图的表现不同。间的分布情况进行比较，是近年来比较新的图绘类型，在一些顶级期刊中经常看到，比如下图就是发表在Science上的一张经典的山脊图，描绘了组间数据的。它能使复杂混乱的数据变得整洁美观，山脊图在 y 轴上有类或组，而在 x 轴上有数字特征。每个类别或类别的组都会产生相互重叠的密度曲线，从而创建一个漂亮的图绘。

原创 警告：范围不一致！一个或多个已添加的图层的范围与关联的空间参考信息不一致。

再看其他数据要素坐标是6位+7位，这就是问题的根源，如果有带号，应该是8位+7位。之所以有警告是因为空间参考是有带号，而实际数据又没有带号。表示6度分带，有带号，带号为13带，中央经度是13X6-3=105。再看其他图层默认参考系统是投影坐标，3度分带，38带，有带号。表示3度分带，有带号，带号为38带，中央经度是38X3=114。尝试用投影工具转换在有带号和无带号的数据都不成功。WKID: 4526 权限: EPSG。表示3度分带，无带号，中央经度是114。表示6度分带，无带号，中央经度是105。

原创 ArcGIS经纬度转平面坐标教程

点击输出坐标系文本框右侧按钮，弹出空间参考属性对话框，依次点击投影坐标系->Gauss Kruger->CGCS2000->CGCS2000_3_Degree_GK_Zone_35，CGCS2000_3_Degree_GK_Zone_35代表3度带、高斯投影、代号为35度带，具体选择多少带号，根据经度值换算，具体换算方式自己网上查询（测绘人都知道）。在输出要素类处，选择输出数据的路径及名称，点击确定，在arcmap右侧内容列表中生成一个新的要素类，此步骤目的是将数据导出为shp文件。

原创 ArcGis中的标注技巧，实现上下标注、换行标注或分式标注的实用技巧

"&[字段]&""&[字段]&""&[字段]&"

原创 Arcgis属性表中如何给某个字段自动编号

我的资源列表里有，小工具成品，直接下载：https://download.csdn.net/download/qq_41127811/21442281?或者导出shape得到FID，字段计算器下计算等于FID即可（此方法较为麻烦）Arcgis属性表中如何给某一字段从“1，2，3，4……因为当你遇上几十万的数据时候，这个方法就显得效率低下，当然用excel表挂接的这个方法就不说了，方法三：通过python制作小工具。属性表中如何给某个字段自动编号。

原创 Python爬取公交线路信息及站点shp数据 文末附数据下载地址

三、现在已经获取到了各个线路名称、链接，以及所途径的上下行站点，那么接下来就应该对开始对第三页url进行信息采集，即本次爬虫的最终目的：爬取所有站点的GPS坐标数据。但坑就坑在，二环一号线属于环线，所以只有左侧有站点数据，而其他非环线车次的站点数据分有上下行，如果用正则匹配，则分不清楚到底哪些是上行站点哪些是下行站点。以前爬过百度地图，当时用的是API平台，加上网站比较规范，所以标签节点什么的都比较清晰，但这次由于特殊原因所选择的网站为公交网，废话不多说直接开始吧。checi=二环一号线。

原创 Python爬取地铁线路及站点数据，文末附上爬取的各地区shp数据下载链接

是字典的格式，其中‘status'键值对代表查询的状态，'total'代表查询的区域内共有多少个POI，如果大于20的话，我们就需要翻页查询了，'result'里就是我们查询的POI的结果，有地址，坐标，甚至还有电话，这里的detail我没有展开，展开的话还有打分评价之类的。按照这个思路，我们就可以查询地铁站周围的各个类型的POI了，具体操作很简单，写个循环就行，唯一需要注意的就是要判别需不需要翻页，此外，百度api限制了并发时间，每次查询完我们需要time.sleep()一下，

原创 全球建筑物轮廓数据

其中包含阿根廷3427787个建筑物轮廓，玻利维亚1015151个建筑物轮廓，巴西18711536个建筑物轮廓，智利2208744个建筑物轮廓，哥伦比亚6083821个建筑物轮廓，厄瓜多尔3674190个建筑物轮廓，圭亚那3339个建筑物轮廓，巴拉圭990756个建筑物轮廓，秘鲁1710431个建筑物轮廓，乌拉圭2656个建筑物轮廓，委内瑞拉6572969个建筑物轮廓。全国范围的50516194个建筑物轮廓数据集，数据集包含尼日利亚35767509个建筑物轮廓，肯尼亚14748685个建筑物轮廓。

原创 中国自然地理分区数据集 （含农业区划、森林工程、生态保护区、九大流域等）

代码名称类型1秦岭山地生态功能保护区水涵养生态功能保护区2若尔盖一玛曲生态功能保护区水涵养生态功能保护区3滇西北生态功能保护区水涵养生态功能保护区4珠江源（云南部分）生态功能保护区水涵养生态功能保护区5雅鲁藏布江源头生态功能保护区水涵养生态功能保护区6长江源生态功能保护区水涵养生态功能保护区7黄河源生态功能保护区水涵养生态功能保护区8阿尔泰山地生态功能保护区水涵养生态功能保护区9天山山地生态功能保护区水涵养生态功能保护区10。

原创 全国土地利用数据10米分辨率

该项目的关键成果是于2021年10月发布了2020年可免费获取的10米分辨率全球土地覆盖产品，该产品基于Sentinel-1和Sentinel-2数据，包含11个土地覆盖类别，并独立验证，全球总体准确率为74.4%。该数据源于Sentinel-1/2的10m分辨率图像制作而成，本次分享我国区域土地利用数据。链接：https://pan.baidu.com/s/1dUix3DG8Nyq6E1YZ-lgSyw。共分为11类: 林地、灌木、草地、耕地、建筑、荒漠、雪\冰、水体、湿地、红树林、苔藓\地衣.

原创 全国分省份水系数据

链接：https://pan.baidu.com/s/12SZNO0QeIic7_0OqBqTO0Q。这次数据太多了，就不一一展示了，大家自行下载去浏览吧。这次给大家分享的数据是全国分省份水系数据。

原创 中国湖泊数据集（1960s-2020）

结合Landsat影像(4215景)、地形图，利用半自动水体提取及人工目视检查编辑，完成了过去60多年来（1960s, 1970s, 1990, 1995, 2000, 2005, 2010, 2015, 2020）详细的中国湖泊(大于1平方公里)数量与面积变化研究。从1960s到2020年，中国湖泊总数量（≥ 1 km^2）从2127个增加到2621个，面积从68537 km^2扩张到82302 km^2。日地空间物理与天文陆地表层。数据文件命名方式和使用方法。日地空间物理与天文陆地表层。

原创 1985-2021年30m全国逐年土地覆被数据（含分省数据）

基于CLCD，揭示了1985年至2019年期间中国土地覆被变化的趋势和模式，如不透水表面（+148.71%）和水（+18.39%）的扩大，耕地（-4.85%）和草地（-3.29%）的减少，以及森林（+4.34%）的增加。由于中国的一系列造林政策，如1980年代启动的三北防护林工程和2000年后启动的绿色增益工程，从1985年到 2019年，森林增加了4.34%，灌木减少了259万公顷，与 ESACCI_LC的减少趋势相似。尽管为了养活不断增长的人口，中国东北和西北地区进行了广泛的开垦，但耕地普遍减少。

原创 中国各城市土地利用类型（城市功能）数据集（shp）

中国各城市土地利用类型（城市功能）数据集时间：2018年全国范围的城市用地类型数据（居住/商业/交通用地等共计11类）分类：居住用地、商业用地、工业用地、医疗设施用地、体育文化设施用地、交通场站用地、绿地等用地类型含城市编码、一级分类5个、二级分类11个1、上海市。