GIS hacker

摘要： 都给我用rasterio！快，不止一点点 # Convert bitmask to polygon (zero values only) import rasterio from rasterio.features import shapes from shapely.geometry import 阅读全文

摘要： ArcGIS 自带的 Zonal Statistics As Table 在分区统计栅格面积时很好用，但是这个工具只能针对栅格，不能针对矢量。 最近有这个需求，所以就做了一个，基本都是 ArcPy 的自带函数，输入的参数和 Zonal Statistics As Table 高度类似，在此不作详细解 阅读全文

摘要： import arcpy,os,sys import arcpy.da srcDir = "E:\\xxx\\xxx\\" for file in os.listdir(srcDir): if os.path.splitext(file)[1] == ".shp": print(file) inpu 阅读全文

摘要： 众所周知，GEE JS端，print一下ImageCollection就知道里面是啥 但是，能看到的仅仅是文字信息而已 那么，我们应该怎么把他们批量加载进来呢？ 看代码： var imageList = imageCollection.toList(imageCollection.size()) v 阅读全文

摘要： 今天正式开始整活儿Colab！ Colab自带了OpenJDK，但是如果我为了保证完全的兼容性，可不可以装个Oracle JDK呢？答案是可以！ 话不多说，上代码！ !sudo add-apt-repository ppa:ts.sch.gr/ppa !sudo apt-get update !su 阅读全文

摘要： 前提： 已经在本地配置好GEE Python客户端，方法参见： [GEE Tips-4]本地配置GEE Python API环境（Windows） 以下所有操作都在命令提示符中进行。 查看所有正在进行的任务： earthengine task list 取消所有正在进行的任务 earthengine 阅读全文

摘要： def imageCollectionToDrive(collection,folder,scale): taskList = [] imageList = collection.toList(collection.size()) for n in range(0,collection.size() 阅读全文

摘要： #读csv def readCsvTo2DList(csvFullFileName,skipHeader,ignoreTail): listFile = open(csvFullFileName,'r') context = listFile.read() contextList = context 阅读全文

摘要： #将某一图层某一字段全部值归入一个List的函数 def GetFieldValueList(inTable,inField): value_list=[] rows=arcpy.da.SearchCursor(inTable,inField) for row in rows: value_list 阅读全文

摘要： 主函数最开头加这么一句： 阅读全文