豆芽菜

libpoppler.so.126​文件缺失了，查找网上解决办法知道需要下载正确的poppler​库。的解答办法，下载libpoppler.so.126文件，上传到服务器指定位置/home/XXXX/anaconda3/envs/XXXX/lib/​下，解决问题。原文链接：https://blog.youkuaiyun.com/m0_62443558/article/details/142880775。conda直接安装 gdal后使用报错，下载的gdal版本是3.6.2。在Python环境下导入gdal。

但该方法仅支持GOES-16卫星，笔者以GOES-18为例进行了扩展。GOES卫星产品没有提供经纬度，需要单独计算，这个代码提供了GOES-16卫星获取经纬度的方法。

affine是一个用于描述二维图像仿射变换的包，可以定义变换参数并进行变换，可以几何其他的GIS包使用，比如GDAL, RasterIO等。特别需要注意的是from_gdal中地理参考参数顺序是 c a b f d e。下面给出一个使用affine定义用于GIS包地理参考的例子，

网上流行的多是用gdal进行矢量转栅格，而gdal用起来需要考虑地理坐标信息等，容易出错。通过geopandas和rasterio同样可以实现。需要注意的是对于没有地理信息的矢量和图像，需要特别处理下更新下地理信息，即翻转下y坐标。

需要下些带坐标的高清遥感历史影像，但通过GEE、PIE、地理空间数据云、USGS等途径都找不到符合需求的数据。然后，在网上找了许多教程，发现目前可以通过ersi的World Imagery Wayback和谷歌地图能看到高分辨率的历史影像；但是要下载这两个来源的历史影像还是比较麻烦的，需要其他软件、科学上网等途径作为辅助。接下来我就把这两天学到的下载World Imagery Wayback和谷歌地图带坐标的历史影像分享一下，希望对在为下载高清历史遥感影像烦恼的同学有所帮助。

rpb文件是用来存储用于遥感数据几何校正的RPC（Rational Polynomial Coefficients ）模型的文件。类似的还有RPC文件，rpb与rpc文件只是格式不同，但包含的信息一致。其用于从图像坐标转换到地理坐标，L1级别的遥感数据是没有经过地理编码及几何校正的，因此可以通过数据包内的RPC文件进行数据几何校正（地理编码）。

后续也将基于这个简单的小库，做一些应用的演示。RPC可以将其理解为一个带畸变参数的相机模型，其描述的是从三维的地理坐标到二维的卫星影像坐标之间的转换关系，一般称之为从物方到像方，我们可以理解为3D到2D，这也被称之为正投影，具体为已知lon，lat，h，求得像素坐标s，l。而反向投影的则是从影像坐标系到地理坐标系，即从像方到物方，这里要注意的是，这仍然是一个3D到2D的关系，因为2D是无法升维到3D的，具体为已经s，l，h，求得lon与lat。分别为归一化的行列号，使用归一化值是为了减少计算过程中的误差。

欧空局的卫星任务查询网站，卫星任务十分丰富，介绍也很详尽。

GDAL重采样，可以通过写文件时改变图像尺寸和geo_transformes的分辨率信息实现。GDAL裁切先计算索引，通过索引获取裁切后的图像，同时重新计算geo_transformes里的左上角点坐标。

​美国桑迪亚（sandia）国家实验室提供一系列机载SAR数据，包括MiniSAR、FARAR等。数据分辨率4英寸，约0.1米。原始数据下载地址，数据是复数据，以不同格式存储，数据读取可用网站上提供的matlab代码，例如.gff可以用GFF_READER。

亚米级土耳其地震影像数据下载

SAR是主动式侧视雷达系统，且成像几何属于斜距投影类型。因此SAR图像与光学图像在成像机理、几何特征、辐射特征等方面都有较大的区别。在进行SAR图像处理和应用前，需要了解SAR图像的基本特征。

​目前获取世界港口信息的网站包括searates、marinevesseltraffic。searates提供全球海港数据的地图查询与基本信息查看。​marinevesseltraffic 提供港口查询，航行信息、live照片和video等信息查看。​

​目前比较好的数据源是OurAirports。网站目前（2022.8.25）提供全球73,106个机场的信息，功能主要包括查询、全球机场概览、机场列表、下载、评论信息和地图浏览等。

干旱是一种比较特殊的自然灾害，持续时间长、波及范围广、治理难度大是其主要的特点。干旱作为一种全球普遍存在的自然灾害，有120多个国家和地区每年不同程度地遭受干旱灾害的威胁。我国是饱受干旱灾害之苦的国家，根据南京大学符淙斌院士提供的数据，自上世纪70年代以来每年因干旱所造成的经济损失高达上千亿元人民币。云南省从2009年秋季开始一直延续至今的60年一遇的严重干旱，已经造成超过100亿元的经济损失。世界著名的咨询公司美国麦肯锡公司于2009年11月发布的《从粮仓到旱地？》的报告指出，除非中国加大抗击干旱的投入，

1 综合报告金融公司报告中金 | 卫星遥感应用：新生的蓝海，崛起中的龙头（2021）太平洋证券 | 遥感智能:来自天空的“海康威视”（2021） 完整版

遥感影像数据集汇总 - 飞桨AI Studio - 人工智能学习实训社区 (baidu.com)

gdal读取大于4GB的bigtiff影像时import gdalgdal.Open("*.tif")出现错误ERROR 4: This is a BigTIFF file. BigTIFF is not supported by this version of GDAL and libtiff.解决方案：使用高版本gdal 2.2.4 +libtiff 4.0.9 +但是某些版本还可能出现问题，使用 -c conda-forge安装gdal而不是其它源，更加稳定..

自然资源部每年7月份发布上一年度的自然资源部卫星遥感应用报告，是对自然资源领域典型应用的总结，对于遥感应用探索有一定参考。自然资源部卫星遥感应用报告-2020版自然资源部卫星遥感应用报告-2019版 摘录2020年版应用简介： 全国应用服务2020 年，面向自然资源调查监测、督察、执法、海洋防灾减灾、海洋生态 保护等需求，持续开展了全国自然资源常态化季度遥感监测、水资源遥感监测、 红树林遥感监测、冰川遥感监测、青藏高原湖泊变化遥...

1 标准规范查询国标、行标、地方规范等可以在全国标准信息公共服务平台查询，例如国标在：国家标准目录查询- 全国标准信息公共服务平台 (samr.gov.cn)2 标准规范下载如果标准规范右侧显示全文链接，则可下载。点击全文点击在线预览打开链接后，使用另存为无法保存pdf，可以使用打印，选择输出pdf可以保存出pdf文件。...

Landsat Fractional Snow Covered Area Product是一套30米的雪比例产品（含云掩膜） ​、该产品包括2部分。一是雪比例产品fSCA，产品文件说明LXSS_US_HHHVVV_YYYYMMDD_yyyymmdd_CCC_VVV_PACKAGE.tar(e.g., LE07_CU_003011_200008285_20170918_C01_V01_SNOW.tar)L LandsatX Sensor (“C” = OLI/TIRS, “E” = E

1 简介土地覆盖：地球表面当前所具有的自然和人为影响所形成的覆盖物，是地球表面的自然状态，如森林、草场、农田、土壤、冰川、湖泊、沼泽湿地及道路等。土地利用：是人类在生产活动中为达到一定的经济效益、社会效益和生态效益，对土地资源的开发、经营、使用方式的总称。两者的区别:• 土地利用表示与土地相结合的人类活动而产生的不同利用方式，反映土地的社会和经济属性。• 土地覆盖表示地球表面存在的不同类型的覆盖特征，强调的是土地的表面形状，反映土地的自然属性。土地利用/土地覆盖分类系统LULC分

矢量driver:https://gdal.org/drivers/vector/index.html栅格driver:https://gdal.org/drivers/raster/index.html以矢量driver为例，代码中使用的driver是Short name，例如ESRI Shapefile、GeoJSON

中文核心期刊查询可以在知网的期刊库中查询，地址如下：出版来源导航 (cnki.net)查询时可以按类别导航选择，或直接搜索关键词，通过“核心期刊”筛选出核心期刊

泥石流的解译特征 标准的泥石流由三部分组成：形成区、流通区和堆积区，如图 3-1。其中，堆积区是识别泥石流的重要标志，同时也是评估灾害危害程度的重要指标。 泥石流的形成区作为泥石流固体物源的集聚区，植被覆盖度较差，土壤裸露程度较高，表面岩体风化侵蚀严重，多为三面高地环绕的盆地低谷地带，且多群发小规模崩塌或滑坡等不良地质现象。在遥感影像上，形状多呈漏斗状、瓢形或桃叶形；色调相比周围较浅，且面积较大；由于表面岩体破碎严重，纹理粗糙，且伴有繁杂的雨水冲刷细沟。...

This recipe converts raster pixels with a specified value to vector lines. For example the blue pixels (value = 0) are converted to vector lines.import ogr, gdal, osr, osimport numpy as npimport itertoolsfrom math import sqrt,ceildef pixelOffset2co

pyshp库对shp文件的操作本身不包含投影信息，如果需要增加投影信息的话需要单独写.prj文件1、使用pyshp写shp文件（包含shp, dbf, shx）"""using pyshp create shpfile(1) type:You can reference shape types by the numbers or by constants defined by PyShp:shapefile.NULL = 0 shapefile.POINT = 1 shapefile.POL

欧洲环境局提供了时序10米分辨率的主要城市地表覆盖变化数据集（Urban Atlas LCLU Chang）。The European Urban Atlas provides reliable, inter-comparable, high-resolution land use and land cover data for 785 Functional Urban Area (FUA) for the 2018 reference year in EEA39 countries. The spat

Arcpy批量合并多个gdb中相同的要素类，流程为：1、将所有的gdb放在一个目录下，即代码中的inpath (path_allgdb_except_one)2、剪切一个gdb作为合并后的输出gdb，即代码中的outgdb (one.gdb)3、执行下面脚本 (arcgis10.3成功)# -*- coding: cp936 -*- #import osimport sysimport arcpyinpath = r"D:\path_allgdb_except_one"outgd

Landsat Fractional Snow Covered Area Product是一套30米的雪比例产品（含云掩膜）、该产品包括2部分。一是雪比例产品fSCA，产品文件说明LXSS_US_HHHVVV_YYYYMMDD_yyyymmdd_CCC_VVV_PACKAGE.tar(e.g., LE07_CU_003011_200008285_20170918_C01_V01_SNOW.tar)L LandsatX Sensor (“C” = OLI/TIRS, “E” = ETM+,

Case：需要给一个现有的shp数据创建一个字段，并将属性表中原有的一个文本类型的属性转换为整型后填入新创建的字段。Problem：新字段创建成功，但是赋值操作无效，即无法成功给字段写入值。solution：对字段进行赋值后需要，重新写入Feature，否则赋值无效，即layer.SetFeature(feature)。特别注意：在对数据进行读写操作，一定要以读写的方式打开，即Open(filePath,1),该方法的原型为Open(pszName,int bUpdate = false),并且

获取已有栅格文件的空间信息（geo_transform, projection）import gdal, osrrasterfn = 'rasterfn.tif'newRasterfn = 'newRasterfn.tif'raster = gdal.Open(rasterfn)geotransform = raster.GetGeoTransform() # geotransformprojection= raster.GetProjectionRef() # projection

Create Projectionfrom osgeo import osrspatialRef = osr.SpatialReference()spatialRef.ImportFromEPSG(2927) # from EPSGReproject a Geometryfrom osgeo import ogrfrom osgeo import osrsource = osr.SpatialReference()source.ImportFromEPSG(292

Web GIS系列：1.搭建简易Web GIS网站：使用GeoServer+PostgreSQL+PostGIS+OpenLayers32.使用GeoServer+QGIS发布WMTS服务3.使用GeoServer+OpenLayers发布和调用WMTS、Vector Tile矢量切片服务.4.Leaflet入门：添加点线面并导入GeoJSON数据使用GeoServer发布服务请...

矢量操作可以使用 ogr这类基础方法进行，但同时还有一些高级库支持更便捷的矢量操作，例如Fiona(shapely)、pyshp(shapefile)、geopandas等。教程资料导航ogr教程资料：1.GDAL/OGR例子合集2.GDAL官网API3.网上的学习笔记，教程来自于犹他州立大学4.中国OSGeo教程shapely教程资料：shapely空间分析...

WMTS – OGC Web Map Tile Servicehttps://gdal.org/drivers/raster/wmts.html切瓦片https://www.cnblogs.com/wang985850293/p/5381684.html服务发布QGIShttps://blog.youkuaiyun.com/sf2gis2/article/details/52854059...

EPSG：4490China Geodetic Coordinate System 2000 (CGCS2000)经纬度坐标系Unit:degree (supplier to define representation)Geodetic CRS:China Geodetic Coordinate System 2000Datum:China 2000Ellipsoid:...

语义分割类：SN1_buildings、SN2_buildings、SN3_roads、SN4_buildings、SN5_roads、SN6_buildings通过以下方式查看aws s3 lss3://spacenet-dataset/spacenet/SN1_buildings/aws s3 lss3://spacenet-dataset/spacenet/SN2_build...

本文修改自：【python】geoJson转换Shapefilegdal，ogr实现geojson转shpfileimport jsonimport gdalfrom osgeo import ogrdef create_polygon(coords): ring = ogr.Geometry(ogr.wkbLinearRing) for coord in co...

gdal ReadAsArray支持按块读取影像，有两种方式：（1）针对datasetdataset = gdal.Open(file, gdal.GA_ReadOnly)data = dataset.ReadAsArray(xoff=0, yoff=0, xsize=512, ysize=512)（2）针对banddataset = gdal.Open(file, gdal...